sico_maths_m_stub Module Reference

Solvers for systems of linear equations used by SICOPOLIS. More...

Data Types
interface	my_erfc
interface	sor_sprs
interface	tri_sle

Functions/Subroutines
subroutine	sor_sprs_stub (lgs_a_value, lgs_a_index, lgs_a_diag_index, lgs_a_ptr, lgs_b_value, nnz, nmax,
	SOR solver for a system of linear equations lgs_a*lgs_x=lgs_b [matrix storage: compressed sparse row CSR, represented by arrays lgs_a_value(values), lgs_a_index (indices) and lgs_a_ptr (pointers)]. More...
subroutine	tri_sle_stub (a0, a1, a2, x, b, nrows)
	Solution of a system of linear equations Ax=b with tridiagonal matrix A. More...
subroutine	my_erfc_stub (x, retval)
	Bilinear interpolation. More...

Detailed Description

Solvers for systems of linear equations used by SICOPOLIS.

Function/Subroutine Documentation

subroutine sico_maths_m_stub::my_erfc_stub	(	real(dp), intent(in)	x,
		real(dp), intent(out)	retval
	)

Bilinear interpolation.

Computation of the complementary error function erfc(x) = 1-erf(x) with a fractional error everywhere less than 1.2 x 10^(-7) (formula by Press et al., 'Numerical Recipes in Fortran 77').

Definition at line 201 of file sico_maths_m_stub.F90.

subroutine sico_maths_m_stub::sor_sprs_stub	(	real(dp), dimension(nnz), intent(inout)	lgs_a_value,
		integer(i4b), dimension(nnz), intent(inout)	lgs_a_index,
		integer(i4b), dimension(nmax), intent(inout)	lgs_a_diag_index,
		integer(i4b), dimension(nmax+1), intent(inout)	lgs_a_ptr,
		real(dp), dimension(nmax), intent(inout)	lgs_b_value,
		integer(i4b), intent(in)	nnz,
		integer(i4b), intent(in)	nmax
	)

SOR solver for a system of linear equations lgs_a*lgs_x=lgs_b [matrix storage: compressed sparse row CSR, represented by arrays lgs_a_value(values), lgs_a_index (indices) and lgs_a_ptr (pointers)].

Definition at line 84 of file sico_maths_m_stub.F90.

subroutine sico_maths_m_stub::tri_sle_stub	(	real(dp), dimension(0:kcmax+ktmax+krmax+imax+jmax), intent(in)	a0,
		real(dp), dimension(0:kcmax+ktmax+krmax+imax+jmax), intent(inout)	a1,
		real(dp), dimension(0:kcmax+ktmax+krmax+imax+jmax), intent(in)	a2,
		real(dp), dimension(0:kcmax+ktmax+krmax+imax+jmax), intent(out)	x,
		real(dp), dimension(0:kcmax+ktmax+krmax+imax+jmax), intent(inout)	b,
		integer(i4b), intent(in)	nrows
	)

Solution of a system of linear equations Ax=b with tridiagonal matrix A.

Parameters

[in]	a0	a0(j) is element A_(j,j-1) of Matrix A
[in]	a1	a1(j) is element A_(j,j) of Matrix A
[in]	a2	a2(j) is element A_(j,j+1) of Matrix A
[in]	b	inhomogeneity vector
[in]	nrows	size of matrix A (indices run from 0 (!!!) to nrows)
[out]	x	Solution vector.

Definition at line 124 of file sico_maths_m_stub.F90.