LinearSolversModule Module

Solve the matrix equation AX = B



Contents


Subroutines

public subroutine CGSolver(AMat, XMat, BMat, solver_parameters_in)

Solve the matrix equation AX = B using the conjugate gradient method.

Arguments

Type IntentOptional AttributesName
type(Matrix_ps), intent(in) :: AMat

The matrix A, must be hermitian, positive definite.

type(Matrix_ps), intent(inout) :: XMat

The solved for matrix X.

type(Matrix_ps), intent(in) :: BMat

The right hand side.

type(SolverParameters_t), intent(in), optional :: solver_parameters_in

Parameters for the solver

public subroutine CholeskyDecomposition(AMat, LMat, solver_parameters_in)

Compute The Cholesky Decomposition of a Hermitian Positive Definite matrix. This is a really naive implementation, that might be worth revisiting.

Arguments

Type IntentOptional AttributesName
type(Matrix_ps), intent(in) :: AMat

The matrix A, must be hermitian, positive definite.

type(Matrix_ps), intent(inout) :: LMat

The lower diagonal matrix computed.

type(SolverParameters_t), intent(in), optional :: solver_parameters_in

Parameters for the solver