mathutil Module Reference

Math functions. More...

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

Variables
double complex, parameter	cpone =cmplx(1d0, 0d0)
double complex, parameter	cpzero =cmplx(0d0, 0d0)
integer, parameter	quad_gauss =1
integer, parameter	quad_equid =2
double precision, dimension(2), save	coeffsq = (/1d0,-1d0/)
double precision, dimension(2, 2), save	tlq = reshape( (/ 2d0, -1d0, 0d0, 1d0 /), (/2, 2/) )
double precision, dimension(2, 2), save	tlq_inv = reshape( (/ 0.5d0, 0.5d0, 0d0, 1d0 /), (/2, 2/) )
double precision, dimension(2, 2), save	tlq_tr = reshape( (/ 2d0, 0d0,-1d0,1d0 /), (/2, 2/) )
double precision, dimension(2, 2), save	trq = reshape( (/ 1d0, 0d0, -1d0, 2d0 /), (/2, 2/) )
double precision, dimension(2, 2), save	trq_inv = reshape( (/ 1d0, 0d0, 0.5d0,0.5d0 /), (/2, 2/) )
double precision, dimension(2, 2), save	trq_tr = reshape( (/ 1d0, -1d0,0d0, 2d0 /), (/2, 2/) )
double precision, dimension(3), save	coeffsc = (/5d0,-4d0,-1d0/)
double precision, dimension(3, 3), save	tlc = reshape( (/ 6d0, -4d0, -1d0, 0d0, 1d0, 0d0, 0d0, 0d0, 1d0 /), (/3, 3/) )
double precision, dimension(3, 3), save	tlc_inv = reshape( (/ 1d0/6d0, 2d0/3d0, 1d0/6d0, 0d0, 1d0, 0d0, 0d0, 0d0, 1d0 /), (/3, 3/) )
double precision, dimension(3, 3), save	tlc_tr = reshape( (/ 6d0, 0d0, 0d0, -4d0, 1d0, 0d0, -1d0, 0d0, 1d0 /), (/3, 3/) )
double precision, dimension(3, 3), save	trc = reshape( (/ 1d0, 0d0, 0d0, 0d0, 1d0, 0d0, -1d0, -4d0, 6d0 /), (/3, 3/) )
double precision, dimension(3, 3), save	trc_inv = reshape( (/ 1d0, 0d0, 0d0, 0d0, 1d0, 0d0, 1d0/6d0, 2d0/3d0, 1d0/6d0 /), (/3, 3/) )
double precision, dimension(3, 3), save	trc_tr = reshape( (/ 1d0, 0d0, -1d0, 0d0, 1d0, -4d0, 0d0, 0d0, 6d0 /), (/3, 3/) )
double precision, dimension(3, 3), save	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	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/) )

Detailed Description

Math functions.

Author

Julien Dominski

Date

January 2013

Function/Subroutine Documentation

computematrixtransform_mapping21()

subroutine mathutil::computematrixtransform_mapping21	(	type(gemat), intent(inout)	mat,
		double precision, dimension(:,:), intent(in)	Tl,
		double precision, dimension(:,:), intent(in)	Tr,
		integer, intent(in)	ni,
		integer, intent(in)	nj,
		integer, intent(in)	side
	)

computequadraturepoints()

subroutine mathutil::computequadraturepoints	(	double precision, intent(in)	x1,
		double precision, intent(in)	x2,
		double precision, dimension(:), intent(inout)	points,
		double precision, dimension(:), intent(inout)	weights,
		integer, intent(in)	nbq,
		integer, intent(in)	method
	)

Here is the caller graph for this function:

computevectortransform_mapping21()

subroutine mathutil::computevectortransform_mapping21	(	double precision, dimension(:), intent(inout)	vec,
		double precision, dimension(:,:), intent(in)	T,
		integer, intent(in)	ni,
		integer, intent(in)	nj,
		integer, intent(in)	side
	)

computevelocityspacediscretization()

subroutine mathutil::computevelocityspacediscretization	(	double precision, intent(in)	dx,
		integer, intent(in)	nqw,
		integer, intent(in)	nqa,
		double precision, intent(in)	lw,
		integer, intent(out)	nw,
		integer, intent(out)	na
	)

Here is the caller graph for this function:

gaussian_number()

real function mathutil::gaussian_number

Here is the call graph for this function:

Here is the caller graph for this function:

plusminus()

real function mathutil::plusminus

Here is the call graph for this function:

Here is the caller graph for this function:

uniform()

real function mathutil::uniform

Here is the call graph for this function:

Here is the caller graph for this function:

zcomputevectortransform_mapping21()

subroutine mathutil::zcomputevectortransform_mapping21	(	double complex, dimension(:), intent(inout)	vec,
		double precision, dimension(:,:), intent(in)	T,
		integer, intent(in)	ni,
		integer, intent(in)	nj,
		integer, intent(in)	side
	)

zspecialmatrixtransform()

subroutine mathutil::zspecialmatrixtransform	(	type(zpbmat), intent(inout)	mat,
		double precision, dimension(:), intent(in)	coeffs,
		integer, intent(in)	nx,
		integer, intent(in)	ny,
		logical, intent(in)	leftside
	)

Variable Documentation

coeffsc

double precision, dimension(3), save mathutil::coeffsc = (/5d0,-4d0,-1d0/)

coeffsq

double precision, dimension(2), save mathutil::coeffsq = (/1d0,-1d0/)

cpone

double complex, parameter mathutil::cpone =cmplx(1d0, 0d0)

cpzero

double complex, parameter mathutil::cpzero =cmplx(0d0, 0d0)

quad_equid

integer, parameter mathutil::quad_equid =2

quad_gauss

integer, parameter mathutil::quad_gauss =1

tlc

double precision, dimension(3,3), save mathutil::tlc = reshape( (/ 6d0, -4d0, -1d0, 0d0, 1d0, 0d0, 0d0, 0d0, 1d0 /), (/3, 3/) )

tlc_deboor

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/) )

tlc_inv

double precision, dimension(3,3), save mathutil::tlc_inv = reshape( (/ 1d0/6d0, 2d0/3d0, 1d0/6d0, 0d0, 1d0, 0d0, 0d0, 0d0, 1d0 /), (/3, 3/) )

tlc_tr

double precision, dimension(3,3), save mathutil::tlc_tr = reshape( (/ 6d0, 0d0, 0d0, -4d0, 1d0, 0d0, -1d0, 0d0, 1d0 /), (/3, 3/) )

tlq

double precision, dimension(2,2), save mathutil::tlq = reshape( (/ 2d0, -1d0, 0d0, 1d0 /), (/2, 2/) )

tlq_inv

double precision, dimension(2,2), save mathutil::tlq_inv = reshape( (/ 0.5d0, 0.5d0, 0d0, 1d0 /), (/2, 2/) )

tlq_tr

double precision, dimension(2,2), save mathutil::tlq_tr = reshape( (/ 2d0, 0d0,-1d0,1d0 /), (/2, 2/) )

trc

double precision, dimension(3,3), save mathutil::trc = reshape( (/ 1d0, 0d0, 0d0, 0d0, 1d0, 0d0, -1d0, -4d0, 6d0 /), (/3, 3/) )

trc_deboor

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/) )

trc_inv

double precision, dimension(3,3), save mathutil::trc_inv = reshape( (/ 1d0, 0d0, 0d0, 0d0, 1d0, 0d0, 1d0/6d0, 2d0/3d0, 1d0/6d0 /), (/3, 3/) )

trc_tr

double precision, dimension(3,3), save mathutil::trc_tr = reshape( (/ 1d0, 0d0, -1d0, 0d0, 1d0, -4d0, 0d0, 0d0, 6d0 /), (/3, 3/) )

trq

double precision, dimension(2,2), save mathutil::trq = reshape( (/ 1d0, 0d0, -1d0, 2d0 /), (/2, 2/) )

trq_inv

double precision, dimension(2,2), save mathutil::trq_inv = reshape( (/ 1d0, 0d0, 0.5d0,0.5d0 /), (/2, 2/) )

trq_tr

double precision, dimension(2,2), save mathutil::trq_tr = reshape( (/ 1d0, -1d0,0d0, 2d0 /), (/2, 2/) )