solve_linear_system Interface

public interface solve_linear_system

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

Type IntentOptional Attributes Name
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

Type IntentOptional Attributes Name
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)