MathUtil.F90 File Reference

Modules
module	mathutil
	Math functions.

Functions/Subroutines
subroutine	mathutil::computequadraturepoints (x1, x2, points, weights, nbq, method)
subroutine	mathutil::computematrixtransform_mapping21 (mat, Tl, Tr, ni, nj, side)
subroutine	mathutil::zspecialmatrixtransform (mat, coeffs, nx, ny, leftside)
subroutine	mathutil::computevectortransform_mapping21 (vec, T, ni, nj, side)
subroutine	mathutil::zcomputevectortransform_mapping21 (vec, T, ni, nj, side)
subroutine	mathutil::computevelocityspacediscretization (dx, nqw, nqa, lw, nw, na)
real function	mathutil::uniform ()
real function	mathutil::gaussian_number ()
real function	mathutil::plusminus ()

Variables
double complex, parameter	mathutil::cpone =cmplx(1d0, 0d0)
double complex, parameter	mathutil::cpzero =cmplx(0d0, 0d0)
integer, parameter	mathutil::quad_gauss =1
integer, parameter	mathutil::quad_equid =2
double precision, dimension(2), save	mathutil::coeffsq = (/1d0,-1d0/)
double precision, dimension(2, 2), save	mathutil::tlq = reshape( (/ 2d0, -1d0, 0d0, 1d0 /), (/2, 2/) )
double precision, dimension(2, 2), save	mathutil::tlq_inv = reshape( (/ 0.5d0, 0.5d0, 0d0, 1d0 /), (/2, 2/) )
double precision, dimension(2, 2), save	mathutil::tlq_tr = reshape( (/ 2d0, 0d0,-1d0,1d0 /), (/2, 2/) )
double precision, dimension(2, 2), save	mathutil::trq = reshape( (/ 1d0, 0d0, -1d0, 2d0 /), (/2, 2/) )
double precision, dimension(2, 2), save	mathutil::trq_inv = reshape( (/ 1d0, 0d0, 0.5d0,0.5d0 /), (/2, 2/) )
double precision, dimension(2, 2), save	mathutil::trq_tr = reshape( (/ 1d0, -1d0,0d0, 2d0 /), (/2, 2/) )
double precision, dimension(3), save	mathutil::coeffsc = (/5d0,-4d0,-1d0/)
double precision, dimension(3, 3), save	mathutil::tlc = reshape( (/ 6d0, -4d0, -1d0, 0d0, 1d0, 0d0, 0d0, 0d0, 1d0 /), (/3, 3/) )
double precision, dimension(3, 3), save	mathutil::tlc_inv = reshape( (/ 1d0/6d0, 2d0/3d0, 1d0/6d0, 0d0, 1d0, 0d0, 0d0, 0d0, 1d0 /), (/3, 3/) )
double precision, dimension(3, 3), save	mathutil::tlc_tr = reshape( (/ 6d0, 0d0, 0d0, -4d0, 1d0, 0d0, -1d0, 0d0, 1d0 /), (/3, 3/) )
double precision, dimension(3, 3), save	mathutil::trc = reshape( (/ 1d0, 0d0, 0d0, 0d0, 1d0, 0d0, -1d0, -4d0, 6d0 /), (/3, 3/) )
double precision, dimension(3, 3), save	mathutil::trc_inv = reshape( (/ 1d0, 0d0, 0d0, 0d0, 1d0, 0d0, 1d0/6d0, 2d0/3d0, 1d0/6d0 /), (/3, 3/) )
double precision, dimension(3, 3), save	mathutil::trc_tr = reshape( (/ 1d0, 0d0, -1d0, 0d0, 1d0, -4d0, 0d0, 0d0, 6d0 /), (/3, 3/) )
double precision, dimension(3, 3), save	mathutil::trc_deboor = reshape( (/ 1.0, 0.0, 0.0, 1.0/3.0, 2.0/3.0, 0.0, 1.0/6.0, 2.0/3.0, 1.0/6.0 /), (/3, 3/) )
double precision, dimension(3, 3), save	mathutil::tlc_deboor = reshape( (/ 1.0/6.0, 2.0/3.0, 1.0/6.0, 0.0, 2.0/3.0, 1.0/3.0, 0.0, 0.0, 1.0 /), (/3, 3/) )