LIBRARY ROUTINES FOR SPECTRAL METHODS. More...

Functions/Subroutines
subroutine	zwgl (Z, W, NP)
	Generate `NP` Gauss Legendre points `Z` and weights `W` associated with Jacobi polynomial \( P(N)(\alpha=0, \beta=0) \). The polynomial degree `N = NP-1`. More...

subroutine	zwgll (Z, W, NP)

subroutine	zwgj (Z, W, NP, ALPHA, BETA)

subroutine	zwgjd (Z, W, NP, ALPHA, BETA)

subroutine	zwglj (Z, W, NP, ALPHA, BETA)

subroutine	zwgljd (Z, W, NP, ALPHA, BETA)

real(kind=rp) function	endw1 (N, ALPHA, BETA)

real(kind=rp) function	endw2 (N, ALPHA, BETA)

real(kind=rp) function	gammaf (X)

real(kind=rp) function	pnormj (N, ALPHA, BETA)

subroutine	jacg (XJAC, NP, ALPHA, BETA)

subroutine	jacobf (POLY, PDER, POLYM1, PDERM1, POLYM2, PDERM2, N, ALP, BET, X)

real(kind=rp) function	hgj (II, Z, ZGJ, NP, ALPHA, BETA)

real(kind=rp) function	hgjd (II, Z, ZGJ, NP, ALPHA, BETA)

real(kind=rp) function	hglj (II, Z, ZGLJ, NP, ALPHA, BETA)

real(kind=rp) function	hgljd (I, Z, ZGLJ, NP, ALPHA, BETA)

subroutine	dgj (D, DT, Z, NZ, NZD, ALPHA, BETA)

subroutine	dgjd (D, DT, Z, NZ, NZD, ALPHA, BETA)

subroutine	dglj (D, DT, Z, NZ, NZD, ALPHA, BETA)

subroutine	dgljd (D, DT, Z, NZ, NZD, ALPHA, BETA)

subroutine	dgll (D, DT, Z, NZ, NZD)

real(kind=rp) function	hgll (I, Z, ZGLL, NZ)

real(kind=rp) function	hgl (I, Z, ZGL, NZ)

real(kind=rp) function	pnleg (Z, N)

subroutine	legendre_poly (L, x, N)

real(kind=rp) function	pndleg (Z, N)

subroutine	dgllgl (D, DT, ZM1, ZM2, IM12, NZM1, NZM2, ND1, ND2)

subroutine	dgljgj (D, DT, ZGL, ZG, IGLG, NPGL, NPG, ND1, ND2, ALPHA, BETA)

subroutine	dgljgjd (D, DT, ZGL, ZG, IGLG, NPGL, NPG, ND1, ND2, ALPHA, BETA)

subroutine	iglm (I12, IT12, Z1, Z2, NZ1, NZ2, ND1, ND2)

subroutine	igllm (I12, IT12, Z1, Z2, NZ1, NZ2, ND1, ND2)

subroutine	igjm (I12, IT12, Z1, Z2, NZ1, NZ2, ND1, ND2, ALPHA, BETA)

subroutine	igljm (I12, IT12, Z1, Z2, NZ1, NZ2, ND1, ND2, ALPHA, BETA)

Detailed Description

March 1989

For questions, comments or suggestions, please contact:

Einar Malvin Ronquist Room 3-243 Department of Mechanical Engineering Massachusetts Institute of Technology 77 Massachusetts Avenue Cambridge, MA 0299 U.S.A.

Function/Subroutine Documentation

◆ dgj()

subroutine speclib::dgj	(	real(kind=rp), dimension(nzd,nzd)	D,
		real(kind=rp), dimension(nzd,nzd)	DT,
		real(kind=rp), dimension(1)	Z,
			NZ,
			NZD,
		real(kind=rp)	ALPHA,
		real(kind=rp)	BETA
	)

Definition at line 690 of file speclib.f90.

Here is the call graph for this function:

◆ dgjd()

subroutine speclib::dgjd	(	real(kind=rp), dimension(nzd,nzd)	D,
		real(kind=rp), dimension(nzd,nzd)	DT,
		real(kind=rp), dimension(1)	Z,
			NZ,
			NZD,
		real(kind=rp)	ALPHA,
		real(kind=rp)	BETA
	)

Definition at line 734 of file speclib.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

◆ dglj()

subroutine speclib::dglj	(	real(kind=rp), dimension(nzd,nzd)	D,
		real(kind=rp), dimension(nzd,nzd)	DT,
		real(kind=rp), dimension(1)	Z,
			NZ,
			NZD,
		real(kind=rp)	ALPHA,
		real(kind=rp)	BETA
	)

Definition at line 773 of file speclib.f90.

Here is the call graph for this function:

◆ dgljd()

subroutine speclib::dgljd	(	real(kind=rp), dimension(nzd,nzd)	D,
		real(kind=rp), dimension(nzd,nzd)	DT,
		real(kind=rp), dimension(1)	Z,
			NZ,
			NZD,
		real(kind=rp)	ALPHA,
		real(kind=rp)	BETA
	)

Definition at line 817 of file speclib.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

◆ dgljgj()

subroutine speclib::dgljgj	(	real(kind=rp), dimension(nd2,nd1)	D,
		real(kind=rp), dimension(nd1,nd2)	DT,
		real(kind=rp), dimension(nd1)	ZGL,
		real(kind=rp), dimension(nd2)	ZG,
		real(kind=rp), dimension(nd2,nd1)	IGLG,
			NPGL,
			NPG,
			ND1,
			ND2,
			ALPHA,
			BETA
	)

Definition at line 1057 of file speclib.f90.

Here is the call graph for this function:

◆ dgljgjd()

subroutine speclib::dgljgjd	(	real(kind=rp), dimension(nd2,nd1)	D,
		real(kind=rp), dimension(nd1,nd2)	DT,
		real(kind=rp), dimension(nd1)	ZGL,
		real(kind=rp), dimension(nd2)	ZG,
		real(kind=rp), dimension(nd2,nd1)	IGLG,
			NPGL,
			NPG,
			ND1,
			ND2,
		real(kind=rp)	ALPHA,
		real(kind=rp)	BETA
	)

Definition at line 1113 of file speclib.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

◆ dgll()

subroutine speclib::dgll	(	real(kind=rp), dimension(nzd,nzd)	D,
		real(kind=rp), dimension(nzd,nzd)	DT,
		real(kind=rp), dimension(1)	Z,
			NZ,
			NZD
	)

Definition at line 864 of file speclib.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

◆ dgllgl()

subroutine speclib::dgllgl	(	real(kind=rp), dimension(nd2,nd1)	D,
		real(kind=rp), dimension(nd1,nd2)	DT,
		real(kind=rp), dimension(nd1)	ZM1,
		real(kind=rp), dimension(nd2)	ZM2,
		real(kind=rp), dimension(nd2,nd1)	IM12,
			NZM1,
			NZM2,
			ND1,
			ND2
	)

Definition at line 1020 of file speclib.f90.

Here is the call graph for this function:

◆ endw1()

real(kind=rp) function speclib::endw1	(		N,
		real(kind=rp)	ALPHA,
		real(kind=rp)	BETA
	)

Definition at line 344 of file speclib.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

◆ endw2()

real(kind=rp) function speclib::endw2	(		N,
		real(kind=rp)	ALPHA,
		real(kind=rp)	BETA
	)

Definition at line 388 of file speclib.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

◆ gammaf()

real(kind=rp) function speclib::gammaf ( real(kind=rp) X )

Definition at line 432 of file speclib.f90.

Here is the caller graph for this function:

◆ hgj()

real(kind=rp) function speclib::hgj	(		II,
		real(kind=rp)	Z,
		real(kind=rp), dimension(1)	ZGJ,
			NP,
		real(kind=rp)	ALPHA,
		real(kind=rp)	BETA
	)

Definition at line 580 of file speclib.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

◆ hgjd()

real(kind=rp) function speclib::hgjd	(		II,
		real(kind=rp)	Z,
		real(kind=rp), dimension(1)	ZGJ,
			NP,
		real(kind=rp)	ALPHA,
		real(kind=rp)	BETA
	)

Definition at line 609 of file speclib.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

◆ hgl()

real(kind=rp) function speclib::hgl	(		I,
		real(kind=rp)	Z,
		real(kind=rp), dimension(1)	ZGL,
			NZ
	)

Definition at line 923 of file speclib.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

◆ hglj()

real(kind=rp) function speclib::hglj	(		II,
		real(kind=rp)	Z,
		real(kind=rp), dimension(1)	ZGLJ,
			NP,
		real(kind=rp)	ALPHA,
		real(kind=rp)	BETA
	)

Definition at line 633 of file speclib.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

◆ hgljd()

real(kind=rp) function speclib::hgljd	(		I,
		real(kind=rp)	Z,
		real(kind=rp), dimension(1)	ZGLJ,
			NP,
		real(kind=rp)	ALPHA,
		real(kind=rp)	BETA
	)

Definition at line 662 of file speclib.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

◆ hgll()

real(kind=rp) function speclib::hgll	(		I,
		real(kind=rp)	Z,
		real(kind=rp), dimension(1)	ZGLL,
			NZ
	)

Definition at line 901 of file speclib.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

◆ igjm()

subroutine speclib::igjm	(	real(kind=rp), dimension(nd2,nd1)	I12,
		real(kind=rp), dimension(nd1,nd2)	IT12,
		real(kind=rp), dimension(nd1)	Z1,
		real(kind=rp), dimension(nd2)	Z2,
			NZ1,
			NZ2,
			ND1,
			ND2,
		real(kind=rp)	ALPHA,
		real(kind=rp)	BETA
	)

Definition at line 1219 of file speclib.f90.

Here is the call graph for this function:

◆ igljm()

subroutine speclib::igljm	(	real(kind=rp), dimension(nd2,nd1)	I12,
		real(kind=rp), dimension(nd1,nd2)	IT12,
		real(kind=rp), dimension(nd1)	Z1,
		real(kind=rp), dimension(nd2)	Z2,
			NZ1,
			NZ2,
			ND1,
			ND2,
		real(kind=rp)	ALPHA,
		real(kind=rp)	BETA
	)

Definition at line 1246 of file speclib.f90.

Here is the call graph for this function:

◆ igllm()

subroutine speclib::igllm	(	real(kind=rp), dimension(nd2,nd1)	I12,
		real(kind=rp), dimension(nd1,nd2)	IT12,
		real(kind=rp), dimension(nd1)	Z1,
		real(kind=rp), dimension(nd2)	Z2,
			NZ1,
			NZ2,
			ND1,
			ND2
	)

Definition at line 1193 of file speclib.f90.

Here is the call graph for this function:

◆ iglm()

subroutine speclib::iglm	(	real(kind=rp), dimension(nd2,nd1)	I12,
		real(kind=rp), dimension(nd1,nd2)	IT12,
		real(kind=rp), dimension(nd1)	Z1,
		real(kind=rp), dimension(nd2)	Z2,
			NZ1,
			NZ2,
			ND1,
			ND2
	)