solve_linear_system Interface

  1. solver_module.f90
  2. solver
  3. solve_linear_system

public interface solve_linear_system

Contents

Module Procedures


Module Procedures

private function solve_linear_system_dp(A, b, n, lerr) result(x)

Solves the linear equation of the form A*x = b for x. A is an (n,n) arrays x and b are (n) arrays Uses Gaussian elimination, so will have issues if nbody_system is ill-conditioned. Uses quad precision intermidiate values, so works best on small arrays.

Arguments

TypeIntentOptionalAttributesName
real(kind=DP), intent(in), dimension(:,:):: A
real(kind=DP), intent(in), dimension(:):: b
integer(kind=I4B), intent(in) :: n
logical, intent(out) :: lerr

Return Value real(kind=DP),dimension(n)

private function solve_linear_system_qp(A, b, n, lerr) result(x)

Solves the linear equation of the form A*x = b for x. A is an (n,n) arrays x and b are (n) arrays Uses Gaussian elimination, so will have issues if nbody_system is ill-conditioned. Uses quad precision intermidiate values, so works best on small arrays.

Arguments

TypeIntentOptionalAttributesName
real(kind=QP), intent(in), dimension(:,:):: A
real(kind=QP), intent(in), dimension(:):: b
integer(kind=I4B), intent(in) :: n
logical, intent(out) :: lerr

Return Value real(kind=QP),dimension(n)