95    constexpr static float abserr_cholesky = 3e-4;
96    constexpr static float relerr_cholesky = 1e-4;
112    constexpr static float abserr_cholesky = 3e-4;
113    constexpr static float relerr_cholesky = 1e-4;
120    constexpr static float abserr_inv = 4.0e-6;
121    constexpr static float relerr_inv = 5.0e-6;
130    constexpr static float abserr_cholesky = 2e-1;
131    constexpr static float relerr_cholesky = 2e-1;
133    constexpr static float abserr_householder = 2e-4;
134    constexpr static float relerr_householder = 2e-3;
138    constexpr static float relerr_qr = 1e-2;
140    constexpr static float abserr_inv = 3e-2;
141    constexpr static float relerr_inv = 3e-2;
362          b.row(r,c) -= b.row(c,c)*tmp;  in _matinv()
363          res.row(r) -= res.row(c)*tmp;  in _matinv()
369          b.row(r,c) -= b.row(c,c)*tmp;  in _matinv()
370          res.row(r) -= res.row(c)*tmp;  in _matinv()
413    std::cout << "----\r\n";  in testinv()
463    std::cout << "----\r\n";  in testadd()
528    std::cout << "----\r\n";  in testdiag()
594    std::cout << "----\r\n";  in testouter()
640    std::cout << "----\r\n";  in testview()
712    std::cout << "----\r\n";  in testmatvec()
778    std::cout << "----\r\n";  in testcomplexmatvec()
857    std::cout << "----\r\n";  in testmatmult()
952    std::cout << "----\r\n";  in testsubmatmult()
1050    std::cout << "----\r\n";  in testmattranspose()
1157          beta = -_gen_sqrt(alpha*alpha+xnorm2);  in _householder()
1159       T r = number_traits<T>::one() / (alpha - beta);  in _householder()
1161       tau = (beta - alpha)/beta;  in _householder()
1231    std::cout << "----\r\n" << "N = " << NB << "\r\n";  in testHouseholder()
1309    for(index_t c=0;c<NBC-1;c++)  in QR()
1314       MatrixView<T> vt(tmpvec,1,NBR-c);  in QR()
1315       dot(tmpmat.sub(0,1,0,NBC-c),vt,RM.sub(c,c));  in QR()
1318         RM.sub(c,c) - beta * outer(tmpvec.sub(0,NBR-c),tmpmat.row(0,0,NBC-c));  in QR()
1321       // Not valid when c == C-1  in QR()
1322       // We don't want to use a  test since CMSIS-DSP is not using  in QR()
1324       RM.col(c,c+1) = copy(tmpvec.sub(1,NBR-c));  in QR()
1329    auto beta = householder(RM.col(NBC-1,NBC-1),eps,tmpvec);  in QR()
1330    tau[NBC-1] = beta;  in QR()
1332    MatrixView<T> vt(tmpvec,1,NBR-(NBC-1));  in QR()
1333    dot(tmpmat.sub(0,1,0,NBC-(NBC-1)),vt,RM.sub(NBC-1,NBC-1));  in QR()
1335    RM.sub(NBC-1,NBC-1) =   in QR()
1336         RM.sub(NBC-1,NBC-1) - beta * outer(tmpvec.sub(0,NBR-(NBC-1)),tmpmat.row(0,0,NBC-(NBC-1)));  in QR()
1343       for(index_t c=NBC-1;c>=0;c--)  in QR()
1348          MatrixView<T> vt(tmpvec,1,NBR-c);  in QR()
1349          dot(tmpmat.sub(0,1,0,NBR-c),vt,Q.sub(c,c));  in QR()
1352            Q.sub(c,c) - tau[c] * outer(tmpvec.sub(0,NBR-c),tmpmat.row(0,0,NBR-c));  in QR()
1364    std::cout << "----\r\n";  in testQR()
1427    //std::cout << "-------\r\n";  in testQR()
1474    Matrix<T,R,R,TMP_ALLOC> g = a;  in cholesky()  local
1477    Matrix<T> g = a;  in cholesky()
1483    g.col(0,0) = g.col(0,0) * (T)(number_traits<T>::one() / _gen_sqrt(g(0,0)));  in cholesky()
1487       dot(tmp.sub(j),g.sub(j,NBR,0,j) , g.row(j,0,j));  in cholesky()
1489 … g.col(j,j) = (g.col(j,j) - tmp.sub(j)) * (T)(number_traits<T>::one() / _gen_sqrt(g(j,j)- tmp[j]));  in cholesky()
1492    return(g);  in cholesky()
1499    std::cout << "----\r\n";  in testCholesky()
1516    // temporary buffer like CMSIS-DSP  in testCholesky()
1537    //std::cout << "-------\r\n";  in testCholesky()