Lines Matching refs:shifts

177 …const ArrayRef& diag, const IndicesRef& perm, VectorType& singVals, ArrayRef shifts, ArrayRef mus);
178 …iag, const IndicesRef& perm, const VectorType& singVals, const ArrayRef& shifts, const ArrayRef& m…
179 …iag, const IndicesRef& perm, const VectorType& singVals, const ArrayRef& shifts, const ArrayRef& m…
601   Map<ArrayXr> shifts(m_workspace.data()+1*n, n);
612   computeSingVals(col0, diag, perm, singVals, shifts, mus);
618   std::cout << "  shift:    " << shifts.transpose() << "\n";
624 …std::cout << "    check1 (expect0) : " << ((singVals.array()-(shifts+mus)) / singVals.array()).hea…
634   assert(shifts.allFinite());
638   perturbCol0(col0, diag, perm, singVals, shifts, mus, zhat);
647   computeSingVecs(zhat, diag, perm, singVals, shifts, mus, U, V);
707                                          VectorType& singVals, ArrayRef shifts, ArrayRef mus)
724       shifts(k) = k==0 ? col0(0) : diag(k);
866     shifts[k] = shift;
882     const ArrayRef& shifts, const ArrayRef& mus, ArrayRef zhat)
902       RealScalar prod = (singVals(last) + dk) * (mus(last) + (shifts(last) - dk));
910           prod *= ((singVals(j)+dk) / ((diag(i)+dk))) * ((mus(j)+(shifts(j)-dk)) / ((diag(i)-dk)));
912 …if(i!=k && std::abs(((singVals(j)+dk)*(mus(j)+(shifts(j)-dk)))/((diag(i)+dk)*(diag(i)-dk)) - 1) > …
913 …singVals(j)+dk)*(mus(j)+(shifts(j)-dk)))/((diag(i)+dk)*(diag(i)-dk)) << " == (" << (singVals(j)+dk…
919 …( " << prod << ")  ;  " << (singVals(last) + dk) << " * " << mus(last) + shifts(last) << " - " << …
931     const ArrayRef& shifts, const ArrayRef& mus, MatrixXr& U, MatrixXr& V)
949         U(i,k) = zhat(i)/(((diag(i) - shifts(k)) - mus(k)) )/( (diag(i) + singVals[k]));
960 …       V(i,k) = diag(i) * zhat(i) / (((diag(i) - shifts(k)) - mus(k)) )/( (diag(i) + singVals[k]));