SICOPOLIS V5-dev
Revision 1420
|
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... | |
Solvers for systems of linear equations used by SICOPOLIS.
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.
[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.