tlapack/potri_8hpp_source.html

//

// Copyright (c) 2025, University of Colorado Denver. All rights reserved.

//

// This file is part of <T>LAPACK.

// <T>LAPACK is free software: you can redistribute it and/or modify it under

// the terms of the BSD 3-Clause license. See the accompanying LICENSE file.


#ifndef TLAPACK_POTRI_HH

#define TLAPACK_POTRI_HH


#include "tlapack/base/utils.hpp"

#include "tlapack/lapack/lauum_recursive.hpp"

#include "tlapack/lapack/potrf2.hpp"

#include "tlapack/lapack/trtri_recursive.hpp"


namespace tlapack {


template <TLAPACK_SMATRIX matrix_t>


void potri(Uplo uplo, matrix_t& A)

{

    tlapack_check(uplo == LOWER_TRIANGLE || uplo == UPPER_TRIANGLE);

    tlapack_check(nrows(A) == ncols(A));


    potrf2(uplo, A);


    trtri_recursive(uplo, NON_UNIT_DIAG, A);


    lauum_recursive(uplo, A);

}


}  // namespace tlapack


#endif  // TLAPACK_POTRI_HH

tlapack::LOWER_TRIANGLE
constexpr internal::LowerTriangle LOWER_TRIANGLE
Lower Triangle access.
Definition types.hpp:188

tlapack::UPPER_TRIANGLE
constexpr internal::UpperTriangle UPPER_TRIANGLE
Upper Triangle access.
Definition types.hpp:186

tlapack::NON_UNIT_DIAG
constexpr internal::NonUnitDiagonal NON_UNIT_DIAG
The main diagonal is not assumed to consist of 1's.
Definition types.hpp:220

tlapack::Uplo
Uplo
Definition types.hpp:50

utils.hpp

tlapack::potri
void potri(Uplo uplo, matrix_t &A)
Computes the Inverse of a Hermitian positive definite matrix A using recursive algorithms.
Definition potri.hpp:47

tlapack::potrf2
int potrf2(uplo_t uplo, matrix_t &A, const EcOpts &opts={})
Computes the Cholesky factorization of a Hermitian positive definite matrix A using the recursive alg...
Definition potrf2.hpp:73

tlapack::trtri_recursive
int trtri_recursive(uplo_t uplo, Diag diag, matrix_t &C, const EcOpts &opts={})
TRTRI computes the inverse of a triangular matrix in-place Input is a triangular matrix,...
Definition trtri_recursive.hpp:48

tlapack_check
#define tlapack_check(cond)
Throw an error if cond is false.
Definition exceptionHandling.hpp:98

lauum_recursive.hpp

tlapack::lauum_recursive
int lauum_recursive(const Uplo &uplo, matrix_t &C)
LAUUM is a specific type of inplace HERK.
Definition lauum_recursive.hpp:45

potrf2.hpp
Computes the Cholesky factorization of a Hermitian positive definite matrix A using the recursive alg...

tlapack::real_type
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

trtri_recursive.hpp