7template <
typename Scalar>
8void r1mpyq(DenseIndex m, DenseIndex n, Scalar *a,
const std::vector<JacobiRotation<Scalar> > &v_givens,
const std::vector<JacobiRotation<Scalar> > &w_givens)
10 typedef DenseIndex
Index;
13 for (
Index j = n-2; j>=0; --j)
14 for (
Index i = 0; i<m; ++i) {
15 Scalar temp = v_givens[j].c() * a[i+m*j] - v_givens[j].s() * a[i+m*(n-1)];
16 a[i+m*(n-1)] = v_givens[j].s() * a[i+m*j] + v_givens[j].c() * a[i+m*(n-1)];
20 for (
Index j = 0; j<n-1; ++j)
21 for (
Index i = 0; i<m; ++i) {
22 Scalar temp = w_givens[j].c() * a[i+m*j] + w_givens[j].s() * a[i+m*(n-1)];
23 a[i+m*(n-1)] = -w_givens[j].s() * a[i+m*j] + w_givens[j].c() * a[i+m*(n-1)];
Namespace containing all symbols from the Eigen library.
EIGEN_DEFAULT_DENSE_INDEX_TYPE Index