SICOPOLIS V5-dev  Revision 1368
Data Types | Functions/Subroutines
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]a0a0(j) is element A_(j,j-1) of Matrix A
[in]a1a1(j) is element A_(j,j) of Matrix A
[in]a2a2(j) is element A_(j,j+1) of Matrix A
[in]binhomogeneity vector
[in]nrowssize of matrix A (indices run from 0 (!!!) to nrows)
[out]xSolution vector.

Definition at line 124 of file sico_maths_m_stub.F90.