/external/eigen/doc/snippets/ |
MatrixBase_adjoint.cpp | 3 cout << "Here is the adjoint of m:" << endl << m.adjoint() << endl;
|
Jacobi_makeJacobi.cpp | 2 m = (m + m.adjoint()).eval(); 6 m.applyOnTheLeft(0, 1, J.adjoint());
|
LLT_solve.cpp | 7 = (samples.adjoint() * samples).llt().solve((samples.adjoint()*elevations));
|
tut_arithmetic_transpose_conjugate.cpp | 10 cout << "Here is the matrix a^*\n" << a.adjoint() << endl;
|
Jacobi_makeGivens.cpp | 5 v.applyOnTheLeft(0, 1, G.adjoint());
|
RealQZ_compute.cpp | 15 << "\n|QQ* - I|: " << (qz.matrixQ()*qz.matrixQ().adjoint() - MatrixXf::Identity(4,4)).norm() 16 << ", |ZZ* - I|: " << (qz.matrixZ()*qz.matrixZ().adjoint() - MatrixXf::Identity(4,4)).norm()
|
HouseholderSequence_HouseholderSequence.cpp | 15 Matrix3d H0 = Matrix3d::Identity() - h(0) * v0 * v0.adjoint(); 18 Matrix3d H1 = Matrix3d::Identity() - h(1) * v1 * v1.adjoint(); 21 Matrix3d H2 = Matrix3d::Identity() - h(2) * v2 * v2.adjoint();
|
Tridiagonalization_diagonal.cpp | 2 MatrixXcd A = X + X.adjoint(); 3 cout << "Here is a random self-adjoint 4x4 matrix:" << endl << A << endl << endl;
|
/external/eigen/test/ |
product_syrk.cpp | 39 ((s1 * rhs2 * rhs2.adjoint()).eval().template triangularView<Lower>().toDenseMatrix())); 41 VERIFY_IS_APPROX(((m2.template triangularView<Lower>() += s1 * rhs2 * rhs22.adjoint()).nestedExpression()), 42 ((s1 * rhs2 * rhs22.adjoint()).eval().template triangularView<Lower>().toDenseMatrix())); 47 (s1 * rhs2 * rhs2.adjoint()).eval().template triangularView<Upper>().toDenseMatrix()); 49 VERIFY_IS_APPROX((m2.template triangularView<Upper>() += s1 * rhs22 * rhs2.adjoint()).nestedExpression(), 50 (s1 * rhs22 * rhs2.adjoint()).eval().template triangularView<Upper>().toDenseMatrix()); 54 VERIFY_IS_APPROX(m2.template selfadjointView<Lower>().rankUpdate(rhs1.adjoint(),s1)._expression(), 55 (s1 * rhs1.adjoint() * rhs1).eval().template triangularView<Lower>().toDenseMatrix()); 57 VERIFY_IS_APPROX((m2.template triangularView<Lower>() += s1 * rhs11.adjoint() * rhs1).nestedExpression(), 58 (s1 * rhs11.adjoint() * rhs1).eval().template triangularView<Lower>().toDenseMatrix()) [all...] |
product_selfadjoint.cpp | 38 m1 = (m1.adjoint() + m1).eval(); 43 VERIFY_IS_APPROX(m2, (m1 + v1 * v2.adjoint()+ v2 * v1.adjoint()).template triangularView<Lower>().toDenseMatrix()); 47 VERIFY_IS_APPROX(m2, (m1 + (s3*(-v1)*(s2*v2).adjoint()+numext::conj(s3)*(s2*v2)*(-v1).adjoint())).template triangularView<Upper>().toDenseMatrix()); 50 m2.template selfadjointView<Upper>().rankUpdate(-s2*r1.adjoint(),r2.adjoint()*s3,s1); 51 VERIFY_IS_APPROX(m2, (m1 + s1*(-s2*r1.adjoint())*(r2.adjoint()*s3).adjoint() + numext::conj(s1)*(r2.adjoint()*s3) * (-s2*r1.adjoint()).adjoint()).template triangularView<Upp (…) [all...] |
eigensolver_generalized_real.cpp | 33 MatrixType spdA = a.adjoint() * a + a1.adjoint() * a1; 34 MatrixType spdB = b.adjoint() * b + b1.adjoint() * b1; 75 GeneralizedSelfAdjointEigenSolver<MatrixType> eig1(a.adjoint() * a,b.adjoint() * b); 76 eig1.compute(a.adjoint() * a,b.adjoint() * b); 77 GeneralizedEigenSolver<MatrixType> eig2(a.adjoint() * a,b.adjoint() * b) [all...] |
product_notemporary.cpp | 45 VERIFY_EVALUATION_COUNT( m3 = (m1 * m2.adjoint()), 1); 46 VERIFY_EVALUATION_COUNT( m3 = (m1 * m2.adjoint()).transpose(), 1); 47 VERIFY_EVALUATION_COUNT( m3.noalias() = m1 * m2.adjoint(), 0); 53 VERIFY_EVALUATION_COUNT( m3 = m3 + (m1 * m2.adjoint()), 1); 54 VERIFY_EVALUATION_COUNT( m3 = m3 - (m1 * m2.adjoint()), 1); 56 VERIFY_EVALUATION_COUNT( m3 = m3 + (m1 * m2.adjoint()).transpose(), 1); 64 VERIFY_EVALUATION_COUNT( m3.noalias() = s1 * m1 * s2 * m2.adjoint(), 0); 65 VERIFY_EVALUATION_COUNT( m3.noalias() = s1 * m1 * s2 * (m1*s3+m2*s2).adjoint(), 1); 66 VERIFY_EVALUATION_COUNT( m3.noalias() = (s1 * m1).adjoint() * s2 * m2, 0); 67 VERIFY_EVALUATION_COUNT( m3.noalias() += s1 * (-m1*s3).adjoint() * (s2 * m2 * s3), 0) [all...] |
adjoint.cpp | 22 // check compatibility of dot and adjoint 23 VERIFY(test_isApproxWithRef(v1.dot(square * v2), (square.adjoint() * v1).dot(v2), 0)); 57 // check compatibility of dot and adjoint 58 ref = NumTraits<Scalar>::IsInteger ? 0 : (std::max)((std::max)(v1.norm(),v2.norm()),(std::max)((square * v2).norm(),(square.adjoint() * v1).norm())); 59 VERIFY(internal::isMuchSmallerThan(abs(v1.dot(square * v2) - (square.adjoint() * v1).dot(v2)), ref, test_precision<Scalar>())); 67 template<typename MatrixType> void adjoint(const MatrixType& m) function 95 // check basic compatibility of adjoint, transpose, conjugate 96 VERIFY_IS_APPROX(m1.transpose().conjugate().adjoint(), m1); 97 VERIFY_IS_APPROX(m1.adjoint().conjugate().transpose(), m1); 100 VERIFY_IS_APPROX((m1.adjoint() * m2).adjoint(), m2.adjoint() * m1) [all...] |
product_symm.cpp | 27 m1 = (m1+m1.adjoint()).eval(); 54 VERIFY_IS_APPROX(rhs12 = (s1*m2).adjoint().template selfadjointView<Upper>() * (s2*rhs1), 55 rhs13 = (s1*m1).adjoint() * (s2*rhs1)); 57 VERIFY_IS_APPROX(rhs12 = (s1*m2).template selfadjointView<Lower>().adjoint() * (s2*rhs1), 58 rhs13 = (s1*m1).adjoint() * (s2*rhs1)); 67 VERIFY_IS_APPROX(rhs12 = (s1*m2).template selfadjointView<Lower>() * (s2*rhs2.adjoint()), 68 rhs13 = (s1*m1) * (s2*rhs2.adjoint())); 71 VERIFY_IS_APPROX(rhs12 = (s1*m2).template selfadjointView<Upper>() * (s2*rhs2.adjoint()), 72 rhs13 = (s1*m1) * (s2*rhs2.adjoint())); 75 VERIFY_IS_APPROX(rhs12 = (s1*m2.adjoint()).template selfadjointView<Lower>() * (s2*rhs2.adjoint()) [all...] |
product_trmm.cpp | 54 VERIFY_IS_APPROX( ge_xs.noalias() = (s1*mat.adjoint()).template triangularView<Mode>() * (s2*ge_left.transpose()), s1*triTr.conjugate() * (s2*ge_left.transpose())); 55 VERIFY_IS_APPROX( ge_sx.noalias() = ge_right.transpose() * mat.adjoint().template triangularView<Mode>(), ge_right.transpose() * triTr.conjugate()); 57 VERIFY_IS_APPROX( ge_xs.noalias() = (s1*mat.adjoint()).template triangularView<Mode>() * (s2*ge_left.adjoint()), s1*triTr.conjugate() * (s2*ge_left.adjoint())); 58 VERIFY_IS_APPROX( ge_sx.noalias() = ge_right.adjoint() * mat.adjoint().template triangularView<Mode>(), ge_right.adjoint() * triTr.conjugate()); 61 VERIFY_IS_APPROX( (ge_xs_save + s1*triTr.conjugate() * (s2*ge_left.adjoint())).eval(), ge_xs.noalias() += (s1*mat.adjoint()).template triangularView<Mode>() * (s2*ge_left.adjoint()) ) [all...] |
mixingtypes.cpp | 160 VERIFY_IS_APPROX(sd*md.adjoint()*mcd, (sd*md).template cast<CD>().eval().adjoint()*mcd); 161 VERIFY_IS_APPROX(sd*mcd.adjoint()*md, sd*mcd.adjoint()*md.template cast<CD>()); 162 VERIFY_IS_APPROX(sd*md.adjoint()*mcd.adjoint(), (sd*md).template cast<CD>().eval().adjoint()*mcd.adjoint()); 163 VERIFY_IS_APPROX(sd*mcd.adjoint()*md.adjoint(), sd*mcd.adjoint()*md.template cast<CD>().adjoint()) [all...] |
upperbidiagonalization.cpp | 26 MatrixType c = ubd.householderU() * b * ubd.householderV().adjoint(); 28 TransposeMatrixType d = ubd.householderV() * b.adjoint() * ubd.householderU().adjoint(); 29 VERIFY_IS_APPROX(a.adjoint(),d);
|
nomalloc.cpp | 48 m2.col(0).noalias() -= m1.adjoint() * m1.col(0); 49 m2.col(0).noalias() -= m1 * m1.row(0).adjoint(); 50 m2.col(0).noalias() -= m1.adjoint() * m1.row(0).adjoint(); 53 m2.row(0).noalias() -= m1.row(0) * m1.adjoint(); 54 m2.row(0).noalias() -= m1.col(0).adjoint() * m1; 55 m2.row(0).noalias() -= m1.col(0).adjoint() * m1.adjoint(); 59 m2.col(0).noalias() -= m1.adjoint().template triangularView<Upper>() * m1.col(0); 60 m2.col(0).noalias() -= m1.template triangularView<Upper>() * m1.row(0).adjoint(); [all...] |
product_mmtr.cpp | 48 CHECK_MMTR(matc, Lower, = s*soc*sor.adjoint()); 49 CHECK_MMTR(matc, Upper, = s*(soc*soc.adjoint())); 50 CHECK_MMTR(matr, Lower, = s*soc*soc.adjoint()); 51 CHECK_MMTR(matr, Upper, = soc*(s*sor.adjoint())); 53 CHECK_MMTR(matc, Lower, += s*soc*soc.adjoint()); 55 CHECK_MMTR(matr, Lower, += s*sor*soc.adjoint()); 56 CHECK_MMTR(matr, Upper, += soc*(s*soc.adjoint())); 58 CHECK_MMTR(matc, Lower, -= s*soc*soc.adjoint()); 60 CHECK_MMTR(matr, Lower, -= s*soc*soc.adjoint()); 61 CHECK_MMTR(matr, Upper, -= soc*(s*soc.adjoint())); [all...] |
product_extra.cpp | 41 VERIFY_IS_APPROX(m3.noalias() = m1 * m2.adjoint(), m1 * m2.adjoint().eval()); 42 VERIFY_IS_APPROX(m3.noalias() = m1.adjoint() * square.adjoint(), m1.adjoint().eval() * square.adjoint().eval()); 43 VERIFY_IS_APPROX(m3.noalias() = m1.adjoint() * m2, m1.adjoint().eval() * m2); 44 VERIFY_IS_APPROX(m3.noalias() = (s1 * m1.adjoint()) * m2, (s1 * m1.adjoint()).eval() * m2) [all...] |
product_trmv.cpp | 60 VERIFY((m3.adjoint() * v1).isApprox(m1.adjoint().template triangularView<Eigen::Lower>() * v1, largerEps)); 62 VERIFY((m3.adjoint() * (s1*v1.conjugate())).isApprox(m1.adjoint().template triangularView<Eigen::Upper>() * (s1*v1.conjugate()), largerEps)); 68 VERIFY((v1.adjoint() * m3).isApprox(v1.adjoint() * m1.template triangularView<Eigen::Lower>(), largerEps)); 69 VERIFY((v1.adjoint() * m3.adjoint()).isApprox(v1.adjoint() * m1.template triangularView<Eigen::Lower>().adjoint(), largerEps)) [all...] |
lu.cpp | 92 VERIFY_IS_APPROX(m1 * m1.adjoint() * m1image, m1image); 113 m2 = m1.adjoint()*m3; 116 VERIFY_IS_APPROX(m2, m1.adjoint()*m3); 118 m3 = lu.adjoint().solve(m2); 119 VERIFY_IS_APPROX(m2, m1.adjoint()*m3); 169 VERIFY_IS_APPROX(m3, m1.adjoint()*m2); 171 m3 = lu.adjoint().solve(m2); 172 VERIFY_IS_APPROX(m2, m1.adjoint()*m3); 214 VERIFY_IS_APPROX(m3, m1.adjoint()*m2); 216 m3 = plu.adjoint().solve(m2) [all...] |
/external/eigen/doc/examples/ |
tut_arithmetic_dot_cross.cpp | 12 double dp = v.adjoint()*w; // automatic conversion of the inner product to a scalar
|
/external/eigen/lapack/ |
svd.cpp | 63 matrix(vt,*n,*n,*ldvt) = svd.matrixV().adjoint(); 68 matrix(vt,diag_size,*n,*ldvt) = svd.matrixV().adjoint(); 73 matrix(vt,*n,*n,*ldvt) = svd.matrixV().adjoint(); 78 matrix(a,diag_size,*n,*lda) = svd.matrixV().adjoint(); 133 if(*jobv=='A') matrix(vt,*n,*n,*ldvt) = svd.matrixV().adjoint(); 134 else if(*jobv=='S') matrix(vt,diag_size,*n,*ldvt) = svd.matrixV().adjoint(); 135 else if(*jobv=='O') matrix(a,diag_size,*n,*lda) = svd.matrixV().adjoint();
|
/external/eigen/Eigen/src/Householder/ |
BlockHouseholder.h | 38 // triFactor.col(i).head(i).noalias() = -h * vectors.block(i, 0, rs, i).adjoint() 63 triFactor.row(i).tail(rt).noalias() = -hCoeffs(i) * vectors.col(i).tail(rs).adjoint() 92 VectorsType::MaxColsAtCompileTime,MatrixType::MaxColsAtCompileTime> tmp = V.adjoint() * mat; 95 else tmp = T.template triangularView<Upper>().adjoint() * tmp;
|