annulusEigenFunction.h File Reference

#include "besselPrimeZeros.h"
#include "besselZeros.h"

Include dependency graph for annulusEigenFunction.h:

Go to the source code of this file.

Functions
	if (numberOfDimensions==2)
	assert (m< mdbz &&n< ndbz)
	printF ("***Cylinder: Bessel function soln: n=%i, m=%i, k=%i, lambda=%e, omega=%e (c=%8.2e) [za,zb]=[%4.2f,%4.2f]\n", n, m, k, lambda, omega, c, cylinderAxisStart, cylinderAxisEnd)
	for (int k=2;k<=n+1;k++) np1Factorial *
	FOR_3D (i1, i2, i3, J1, J2, J3)

Variables
Index	J1 = Range(max(I1.getBase(),uLocal.getBase(0)),min(I1.getBound(),uLocal.getBound(0)))
Index	J2 = Range(max(I2.getBase(),uLocal.getBase(1)),min(I2.getBound(),uLocal.getBound(1)))
Index	J3 = Range(max(I3.getBase(),uLocal.getBase(2)),min(I3.getBound(),uLocal.getBound(2)))
else const int	ndbz =10
else const int	mdbz =10
real	besselZeros [ndbz][mdbz]
const real	cylinderLength =cylinderAxisEnd-cylinderAxisStart
const int	n = int(initialConditionParameters[0]+.5)
const int	m = int(initialConditionParameters[1]+.5)
const int	k = int(initialConditionParameters[2]+.5)
real	lambda = besselZeros[n][m]
real	omega = sqrt( SQR(k*Pi/cylinderLength) + lambda*lambda )
const real	eps =sqrt(REAL_EPSILON)
real	np1Factorial =1.
int	i1
int	i2
int	i3
real	r
real	gr
real	xd
real	yd
real	zd
real	bj
real	bjp
real	rx
real	ry
real	theta
real	thetax
real	thetay
real	cosTheta
real	sinTheta
real	bjThetax
real	bjThetay
real	uex
real	uey
real	cosn
real	sinn
real	sinkz
real	coskz
real	cost
real	sint

Function Documentation

assert

(

)

for

(

)

References a, addBoundaryForcing, axis, bcfOffset, cg, dbc, fn, grid, i, mg, n, ndp, pbcf, side, timeCoefficientsForTZ, and u.

FOR_3D	(	i1	,
		i2	,
		i3	,
		J1	,
		J2	,
		J3
	)

References cost, cosTheta, ERREX, ERREY, ERREZ, gr, lambda, r, sinn, sinTheta, UEX, UEY, UEZ, UMEX, UMEY, UMEZ, and X.

if

(

numberOfDimensions

= =2

)

References assert(), besselPrimeZeros, bj, bjp, bjThetax, bjThetay, cosn, cost, cosTheta, eps, ERREX, ERREY, ERRHZ, FOR_3D, gr, i1, i2, i3, J1, k, m, mdbpz, n, ndbpz, np1Factorial, omega, r, rx, ry, sinn, sint, sinTheta, theta, thetax, thetay, uex, UEX, uey, UEY, UHZ, uLocal, UMEX, UMEY, UMHZ, X, xd, yd, and zd.

printF	(	"***Cylinder: Bessel function soln:	n = %i,
		m	= %i,
		k	= %i,
		lambda	= %e,
		omega	= %e (c=%8.2e) [za,
		zb]	= [%4.2f,
		%4.2f]\n"	,
		n	,
		m	,
		k	,
		lambda	,
		omega	,
		c	,
		cylinderAxisStart	,
		cylinderAxisEnd
	)

Variable Documentation

real besselZeros[ndbz][mdbz]

Referenced by Maxwell::assignInitialConditions(), and Maxwell::getErrors().

real bj

Referenced by Maxwell::assignInitialConditions(), Maxwell::getErrors(), and if().

real bjp

Referenced by Maxwell::assignInitialConditions(), Maxwell::getErrors(), and if().

real bjThetax

Referenced by Maxwell::assignInitialConditions(), Maxwell::getErrors(), and if().

real bjThetay

Referenced by Maxwell::assignInitialConditions(), Maxwell::getErrors(), and if().

real coskz

Referenced by Maxwell::assignInitialConditions(), and Maxwell::getErrors().

real cosn

Referenced by Maxwell::assignInitialConditions(), Maxwell::getErrors(), and if().

real cost

Referenced by Cgsm::assignAnnulusEigenfunction(), Maxwell::assignBoundaryConditions(), Maxwell::assignInitialConditions(), Cgsm::assignSpecialInitialConditions(), FOR_3D(), Maxwell::getErrors(), MovingGrids::getGridVelocity(), if(), Maxwell::userDefinedForcing(), Cgsm::userDefinedForcing(), and Maxwell::userDefinedInitialConditions().

real cosTheta

Referenced by Cgsm::assignAnnulusEigenfunction(), Maxwell::assignInitialConditions(), FOR_3D(), Maxwell::getErrors(), if(), and DomainSolver::userDefinedInitialConditions().

const real cylinderLength =cylinderAxisEnd-cylinderAxisStart

Referenced by Maxwell::assignInitialConditions(), and Maxwell::getErrors().

const real eps =sqrt(REAL_EPSILON)

Referenced by RadiationBoundaryCondition::assignBoundaryConditions(), DomainSolver::buildGrid(), Reactions::chemistrySourceAndJacobian(), Cgins::computeTurbulenceQuantities(), getExactSolution(), DomainSolver::getTimeStep(), RigidBodyMotion::getTimeStepEstimate(), if(), InterfaceTransfer::initialize(), DomainSolver::saveShow(), and DomainSolver::userDefinedInitialConditions().

real gr

Referenced by Maxwell::assignInitialConditions(), FOR_3D(), Maxwell::getErrors(), and if().

int i1

int i2

int i3

Index J1 = Range(max(I1.getBase(),uLocal.getBase(0)),min(I1.getBound(),uLocal.getBound(0)))

Referenced by DomainSolver::addArtificialDissipation(), Maxwell::advanceC(), Cgins::advanceLineSolveNew(), Cgins::advanceLineSolveOld(), Maxwell::advanceNew(), Maxwell::advanceNFDTD(), DomainSolver::advanceSecondOrderSystem(), Cgcns::applyBoundaryConditions(), Maxwell::assignInitialConditions(), Cgmp::assignInterfaceBoundaryConditions(), Maxwell::assignInterfaceBoundaryConditions(), Cgmp::assignInterfaceRightHandSideOld(), Cgsm::getEnergy(), Maxwell::getEnergy(), Cgsm::getErrors(), Maxwell::getErrors(), Maxwell::getForcing(), Cgmp::getInterfaceResidualsOld(), Cgcns::getTimeSteppingEigenvalue(), Cgcns::getUt(), DeformingBodyMotion::getVelocity(), MovingGrids::gridAccelerationBC(), if(), InterfaceTransfer::initialize(), Cgmp::initializeInterfaces(), Maxwell::initializeInterfaces(), Cgsm::initializeInterfaces(), Cgmp::interfaceProjection(), Cgins::lineSolverBoundaryConditions(), DomainSolver::outputProbes(), DomainSolver::takeTimeStepPC(), and MovingGrids::userDefinedGridAccelerationBC().

Index J2 = Range(max(I2.getBase(),uLocal.getBase(1)),min(I2.getBound(),uLocal.getBound(1)))

Referenced by DomainSolver::addArtificialDissipation(), Maxwell::advanceC(), Cgins::advanceLineSolveNew(), Cgins::advanceLineSolveOld(), Maxwell::advanceNew(), Maxwell::advanceNFDTD(), DomainSolver::advanceSecondOrderSystem(), Cgcns::applyBoundaryConditions(), Maxwell::assignInitialConditions(), Cgmp::assignInterfaceBoundaryConditions(), Maxwell::assignInterfaceBoundaryConditions(), Cgmp::assignInterfaceRightHandSideOld(), Cgsm::getEnergy(), Maxwell::getEnergy(), Cgsm::getErrors(), Maxwell::getErrors(), Maxwell::getForcing(), Cgmp::getInterfaceResidualsOld(), Cgcns::getTimeSteppingEigenvalue(), Cgcns::getUt(), DeformingBodyMotion::getVelocity(), MovingGrids::gridAccelerationBC(), if(), InterfaceTransfer::initialize(), Cgmp::initializeInterfaces(), Maxwell::initializeInterfaces(), Cgsm::initializeInterfaces(), Cgmp::interfaceProjection(), Cgins::lineSolverBoundaryConditions(), DomainSolver::outputProbes(), DomainSolver::takeTimeStepPC(), and MovingGrids::userDefinedGridAccelerationBC().

Index J3 = Range(max(I3.getBase(),uLocal.getBase(2)),min(I3.getBound(),uLocal.getBound(2)))

Referenced by DomainSolver::addArtificialDissipation(), Maxwell::advanceC(), Cgins::advanceLineSolveNew(), Cgins::advanceLineSolveOld(), Maxwell::advanceNew(), Maxwell::advanceNFDTD(), DomainSolver::advanceSecondOrderSystem(), Cgcns::applyBoundaryConditions(), Maxwell::assignInitialConditions(), Cgmp::assignInterfaceBoundaryConditions(), Maxwell::assignInterfaceBoundaryConditions(), Cgmp::assignInterfaceRightHandSideOld(), Cgsm::getEnergy(), Maxwell::getEnergy(), Cgsm::getErrors(), Maxwell::getErrors(), Maxwell::getForcing(), Cgmp::getInterfaceResidualsOld(), Cgcns::getTimeSteppingEigenvalue(), Cgcns::getUt(), DeformingBodyMotion::getVelocity(), MovingGrids::gridAccelerationBC(), if(), InterfaceTransfer::initialize(), Cgmp::initializeInterfaces(), Maxwell::initializeInterfaces(), Cgsm::initializeInterfaces(), Cgmp::interfaceProjection(), Cgins::lineSolverBoundaryConditions(), DomainSolver::outputProbes(), DomainSolver::takeTimeStepPC(), and MovingGrids::userDefinedGridAccelerationBC().

const int k = int(initialConditionParameters[2]+.5)

real lambda = besselZeros[n][m]

const int m = int(initialConditionParameters[1]+.5)

else const int mdbz =10

Referenced by Maxwell::assignInitialConditions(), and Maxwell::getErrors().

const int n = int(initialConditionParameters[0]+.5)

else const int ndbz =10

Referenced by Maxwell::assignInitialConditions(), and Maxwell::getErrors().

real np1Factorial =1.

Referenced by Maxwell::assignInitialConditions(), Maxwell::getErrors(), and if().

real omega = sqrt( SQR(k*Pi/cylinderLength) + lambda*lambda )

real r

Referenced by Cgins::addForcing(), Cgad::advanceADI(), Cgsm::assignAnnulusEigenfunction(), Maxwell::assignInitialConditions(), Cgsm::assignSpecialInitialConditions(), DomainSolver::buildGrid(), Cgmp::checkIfInterfacesMatch(), Reactions::chemistrySource(), Cgasf::computeSource(), Cgins::determineErrors(), DomainSolver::determineErrors(), RigidBodyMotion::dirkImplicitSolve(), TimeFunction::evalDerivative(), evaluateBoundaryCurve(), F90_ID(), FOR_3D(), getClosestGridPoint(), RigidBodyMotion::getCoordinates(), Maxwell::getErrors(), Cgsm::getInputOutputOption(), MatrixMotion::getMotion(), RigidBodyMotion::getPointTransformationMatrix(), RigidBodyMotion::getRotationMatrix(), Parameters::getUserDefinedKnownSolution(), Cgasf::getUt(), if(), initialConditionsFromAProfile(), RadiationBoundaryCondition::initialize(), DeformingBodyMotion::initialize(), Maxwell::initializePlaneMaterialInterface(), Maxwell::interactiveUpdate(), main(), MovingGrids::moveGrids(), ogfTaylor(), Maxwell::outputHeader(), DomainSolver::parabolicInflow(), Reactions::pFromRTY(), DeformingBodyMotion::regenerateComponentGrids(), Reactions::rFromPTY(), MovingGrids::rigidBodyMotion(), Maxwell::setupDSICoefficients(), Maxwell::setupGrids(), TestRigidBody::solve(), Reactions::solveImplicitForPTYGivenR(), Reactions::solveImplicitForRTYGivenP(), Reactions::solveImplicitForYGivenRTP(), Cgins::turbulenceModels(), Cgsm::userDefinedBoundaryValues(), and DomainSolver::userDefinedGrid().

real rx

real ry

real sinkz

Referenced by Maxwell::assignInitialConditions(), and Maxwell::getErrors().

real sinn

Referenced by Maxwell::assignInitialConditions(), FOR_3D(), Maxwell::getErrors(), and if().

real sint

Referenced by Maxwell::assignBoundaryConditions(), Maxwell::assignInitialConditions(), Cgsm::assignSpecialInitialConditions(), Maxwell::getErrors(), MovingGrids::getGridVelocity(), if(), Maxwell::userDefinedForcing(), and Cgsm::userDefinedForcing().

real sinTheta

Referenced by Cgsm::assignAnnulusEigenfunction(), Maxwell::assignInitialConditions(), FOR_3D(), Maxwell::getErrors(), if(), and DomainSolver::userDefinedInitialConditions().

real theta

Referenced by Cgcns::applyBoundaryConditionsForImplicitTimeStepping(), Maxwell::assignInitialConditions(), Cgcns::formImplicitTimeSteppingMatrix(), Maxwell::getErrors(), MovingGrids::getGridVelocity(), MatrixMotion::getMotion(), FlowSolutions::getSupersonicExpandingFlow(), Parameters::getUserDefinedKnownSolution(), Cgcns::getUt(), if(), DeformingBodyMotion::initialize(), DeformingBodyMotion::regenerateComponentGrids(), Cgsm::userDefinedBoundaryValues(), and DomainSolver::userDefinedInitialConditions().

real thetax

Referenced by Maxwell::assignInitialConditions(), Maxwell::getErrors(), and if().

real thetay

Referenced by Maxwell::assignInitialConditions(), Maxwell::getErrors(), and if().

real uex

Referenced by Cgad::advanceADI(), Cgcns::applyBoundaryConditions(), Cgins::applyBoundaryConditions(), Cgad::applyBoundaryConditionsForImplicitTimeStepping(), Cgins::applyBoundaryConditionsForImplicitTimeStepping(), Maxwell::assignInitialConditions(), Cgins::assignLineSolverBoundaryConditions(), Maxwell::getErrors(), if(), Cgins::lineSolverBoundaryConditions(), and Maxwell::outputResultsAfterEachTimeStep().

real uey

Referenced by Cgad::advanceADI(), Cgcns::applyBoundaryConditions(), Cgins::applyBoundaryConditions(), Cgad::applyBoundaryConditionsForImplicitTimeStepping(), Cgins::applyBoundaryConditionsForImplicitTimeStepping(), Maxwell::assignInitialConditions(), Cgins::assignLineSolverBoundaryConditions(), Maxwell::getErrors(), if(), Cgins::lineSolverBoundaryConditions(), and Maxwell::outputResultsAfterEachTimeStep().

real xd

Referenced by Maxwell::assignInitialConditions(), Maxwell::getErrors(), Maxwell::getForcing(), and if().

real yd

Referenced by Maxwell::assignInitialConditions(), Maxwell::getErrors(), Maxwell::getForcing(), and if().

real zd

Referenced by Maxwell::assignInitialConditions(), Maxwell::getErrors(), Maxwell::getForcing(), and if().