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(kind=r8) function	mathutil::uniform ()
real(kind=r8) function	mathutil::gaussian_number ()
real(kind=r8) function	mathutil::plusminus ()

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