Searched refs:A_jj (Results 1 – 5 of 5) sorted by relevance
| /dports/comms/liquid-dsp/liquid-dsp-1.3.2/src/matrix/src/ |
| H A D | matrix.chol.c | 48 T A_jj; in MATRIX() local
56 A_jj = matrix_access(_A,_n,_n,j,j); in MATRIX()
57 if ( creal(A_jj) < 0.0 ) { in MATRIX()
58 …: matrix_chol(), matrix is not positive definite (real{A[%u,%u]} = %12.4e < 0)\n",j,j,creal(A_jj)); in MATRIX()
62 if ( fabs(cimag(A_jj)) > 0.0 ) { in MATRIX()
63 …_chol(), matrix is not positive definite (|imag{A[%u,%u]}| = %12.4e > 0)\n",j,j,fabs(cimag(A_jj))); in MATRIX()
79 if ( creal(A_jj) < t0 ) { in MATRIX()
80 …_chol(), matrix is not positive definite (real{A[%u,%u]} = %12.4e < %12.4e)\n",j,j,creal(A_jj),t0); in MATRIX()
83 L_jj = sqrt( A_jj - t0 ); in MATRIX()
|
| /dports/math/armadillo/armadillo-10.7.1/include/armadillo_bits/ |
| H A D | sympd_helper.hpp | 57 const eT A_jj = A_col[j]; in guess_sympd_worker() local
59 if(A_jj <= eT(0)) { return false; } in guess_sympd_worker()
61 max_diag = (A_jj > max_diag) ? A_jj : max_diag; in guess_sympd_worker()
73 const eT A_jj = A_col[j]; in guess_sympd_worker() local
97 if( (A_ij_abs + A_ij_abs) >= (A_ii + A_jj) ) { return false; } in guess_sympd_worker()
133 const eT& A_jj = A_col[j]; in guess_sympd_worker() local
134 const T A_jj_real = std::real(A_jj); in guess_sympd_worker()
135 const T A_jj_imag = std::imag(A_jj); in guess_sympd_worker()
|
| /dports/math/R-cran-RcppArmadillo/RcppArmadillo/inst/include/armadillo_bits/ |
| H A D | sympd_helper.hpp | 57 const eT A_jj = A_col[j]; in guess_sympd_worker() local
59 if(A_jj <= eT(0)) { return false; } in guess_sympd_worker()
61 max_diag = (A_jj > max_diag) ? A_jj : max_diag; in guess_sympd_worker()
73 const eT A_jj = A_col[j]; in guess_sympd_worker() local
97 if( (A_ij_abs + A_ij_abs) >= (A_ii + A_jj) ) { return false; } in guess_sympd_worker()
133 const eT& A_jj = A_col[j]; in guess_sympd_worker() local
134 const T A_jj_real = std::real(A_jj); in guess_sympd_worker()
135 const T A_jj_imag = std::imag(A_jj); in guess_sympd_worker()
|
| /dports/math/viennacl/ViennaCL-1.7.1/viennacl/linalg/ |
| H A D | qr.hpp | 66 ScalarType A_jj = A(j,j); in setup_householder_vector_ublas() local
78 ScalarType mu = std::sqrt(sigma + A_jj*A_jj); in setup_householder_vector_ublas()
80 ScalarType v1 = (A_jj <= 0) ? (A_jj - mu) : (-sigma / (A_jj + mu)); in setup_householder_vector_ublas()
106 ScalarType A_jj = A(j,j); in setup_householder_vector_viennacl() local
118 ScalarType mu = std::sqrt(sigma + A_jj*A_jj); in setup_householder_vector_viennacl()
120 ScalarType v1 = (A_jj <= 0) ? (A_jj - mu) : (-sigma / (A_jj + mu)); in setup_householder_vector_viennacl()
|
| /dports/science/dakota/dakota-6.13.0-release-public.src-UI/packages/external/trilinos/packages/teuchos/numerics/src/ |
| H A D | Teuchos_Details_Lapack128.cpp | 93 __float128* const A_jj = A + (j-1)*LDA + (j-1); in GETRF() local
96 const int jp = (j - 1) + blas.IAMAX (M - j + 1, A_jj, 1); in GETRF()
113 if (fabsq (*A_jj) >= sfmin) { in GETRF()
114 blas.SCAL (M-j, one / *A_jj, A_j1_j, 1); in GETRF()
118 *A_jpi_j /= *A_jj; in GETRF()
|