12#ifndef TLAPACK_POTRF_BLOCKED_RL_HH
13#define TLAPACK_POTRF_BLOCKED_RL_HH
60template <TLAPACK_UPLO uplo_t, TLAPACK_SMATRIX matrix_t>
70 const idx_t n = nrows(
A);
71 const idx_t nb =
opts.nb;
87 for (idx_t j = 0; j < n; j += nb) {
88 idx_t
jb = min(nb, n - j);
97 "The leading minor of the reported order is not "
99 " and the factorization could not be completed.");
114 for (idx_t j = 0; j < n; j += nb) {
115 idx_t
jb = min(nb, n - j);
124 "The leading minor of the reported order is not "
126 " and the factorization could not be completed.");
void herk(Uplo uplo, Op trans, const alpha_t &alpha, const matrixA_t &A, const beta_t &beta, matrixC_t &C)
Hermitian rank-k update:
Definition herk.hpp:68
void trsm(Side side, Uplo uplo, Op trans, Diag diag, const alpha_t &alpha, const matrixA_t &A, matrixB_t &B)
Solve the triangular matrix-vector equation.
Definition trsm.hpp:76
int potrf_rl(uplo_t uplo, matrix_t &A, const BlockedCholeskyOpts &opts={})
Computes the Cholesky factorization of a Hermitian positive definite matrix A using a blocked algorit...
Definition potrf_blocked_right_looking.hpp:61
int potf2(uplo_t uplo, matrix_t &A)
Computes the Cholesky factorization of a Hermitian positive definite matrix A using a level-2 algorit...
Definition potf2.hpp:58
#define tlapack_check(cond)
Throw an error if cond is false.
Definition exceptionHandling.hpp:98
#define tlapack_error(info, detailedInfo)
Error handler.
Definition exceptionHandling.hpp:142
Computes the Cholesky factorization of a Hermitian positive definite matrix A using a level-2 algorit...
Sort the numbers in D in increasing order (if ID = 'I') or in decreasing order (if ID = 'D' ).
Definition arrayTraits.hpp:15
typename traits::real_type_traits< Types..., int >::type real_type
The common real type of the list of types.
Definition scalar_type_traits.hpp:113
constexpr internal::LowerTriangle LOWER_TRIANGLE
Lower Triangle access.
Definition types.hpp:188
constexpr internal::UpperTriangle UPPER_TRIANGLE
Upper Triangle access.
Definition types.hpp:186
constexpr internal::RightSide RIGHT_SIDE
right side
Definition types.hpp:296
constexpr internal::NonUnitDiagonal NON_UNIT_DIAG
The main diagonal is not assumed to consist of 1's.
Definition types.hpp:220
constexpr internal::ConjTranspose CONJ_TRANS
conjugate transpose
Definition types.hpp:264
@ Upper
0 <= i <= j, 0 <= j <= n.
@ Lower
0 <= i <= m, 0 <= j <= i.
constexpr internal::NoTranspose NO_TRANS
no transpose
Definition types.hpp:260
constexpr internal::LeftSide LEFT_SIDE
left side
Definition types.hpp:294
Computes the Cholesky factorization of a Hermitian positive definite matrix A using a blocked algorit...
Definition potrf_blocked.hpp:22