CG  Version 25
Macros | Functions
asf.C File Reference
#include "Cgasf.h"
#include "AsfParameters.h"
#include "MappedGridOperators.h"
#include "interpPoints.h"
#include "Chemkin.h"
#include "ParallelUtility.h"
Include dependency graph for asf.C:

Macros

#define ForBoundary(side, axis)   for( axis=0; axis<c.numberOfDimensions(); axis++ ) for( side=0; side<=1; side++ )
 
#define asfdt   EXTERN_C_NAME(asfdt)
 
#define U(c)   u(I1,I2,I3,c)
 
#define UU(c)   uu(I1,I2,I3,c)
 
#define UX(c)   ux(I1,I2,I3,c)
 
#define UY(c)   uy(I1,I2,I3,c)
 
#define UZ(c)   uz(I1,I2,I3,c)
 
#define UXX(c)   uxx(I1,I2,I3,c)
 
#define UXY(c)   uxy(I1,I2,I3,c)
 
#define UXZ(c)   uxz(I1,I2,I3,c)
 
#define UYY(c)   uyy(I1,I2,I3,c)
 
#define UYZ(c)   uyz(I1,I2,I3,c)
 
#define UZZ(c)   uzz(I1,I2,I3,c)
 
#define AD2(kd)   ( (ad21 + cd22* ( fabs(UX(uc))+fabs(UY(uc)) +fabs(UX(vc))+fabs(UY(vc)) ) ) *(u(I1+1,I2,I3,kd)-4.*u(I1,I2,I3,kd)+u(I1-1,I2,I3,kd) +u(I1,I2+1,I3,kd) +u(I1,I2-1,I3,kd)) )
 
#define AD23(kd)   ( (ad21 + cd22* ( fabs(UX(uc))+fabs(UY(uc))+fabs(UZ(uc)) +fabs(UX(vc))+fabs(UY(vc))+fabs(UZ(vc)) +fabs(UX(wc))+fabs(UY(wc))+fabs(UZ(wc)) ) ) *(u(I1+1,I2,I3,kd)-6.*u(I1,I2,I3,kd)+u(I1-1,I2,I3,kd) +u(I1,I2+1,I3,kd) +u(I1,I2-1,I3,kd) +u(I1,I2,I3+1,kd) +u(I1,I2,I3-1,kd)) )
 
#define AD4(kd)   ( (ad41 + cd42* ( fabs(UX(uc))+fabs(UY(uc)) +fabs(UX(vc))+fabs(UY(vc)) ) ) *( -u(I1+2,I2,I3,kd)-u(I1-2,I2,I3,kd) -u(I1,I2+2,I3,kd)-u(I1,I2-2,I3,kd) +4.*(u(I1+1,I2,I3,kd)+u(I1-1,I2,I3,kd) +u(I1,I2+1,I3,kd)+u(I1,I2-1,I3,kd)) -12.*u(I1,I2,I3,kd) ) )
 
#define AD43(kd)   ( (ad41 + cd42* ( fabs(UX(uc))+fabs(UY(uc))+fabs(UZ(uc)) +fabs(UX(vc))+fabs(UY(vc))+fabs(UZ(vc)) +fabs(UX(wc))+fabs(UY(wc))+fabs(UZ(wc)) ) ) *( -u(I1+2,I2,I3,kd)-u(I1-2,I2,I3,kd) -u(I1,I2+2,I3,kd)-u(I1,I2-2,I3,kd) -u(I1,I2,I3+2,kd)-u(I1,I2,I3-2,kd) +4.*(u(I1+1,I2,I3,kd)+u(I1-1,I2,I3,kd) +u(I1,I2+1,I3,kd)+u(I1,I2-1,I3,kd) +u(I1,I2,I3+1,kd)+u(I1,I2,I3-1,kd)) -18.*u(I1,I2,I3,kd) ) )
 
#define R(t)   e (mg,I1,I2,I3,rc,t)
 
#define RT(t)   e.t (mg,I1,I2,I3,rc,t)
 
#define RX(t)   e.x (mg,I1,I2,I3,rc,t)
 
#define RY(t)   e.y (mg,I1,I2,I3,rc,t)
 
#define RZ(t)   e.z (mg,I1,I2,I3,rc,t)
 
#define RXX(t)   e.xx(mg,I1,I2,I3,rc,t)
 
#define RYY(t)   e.yy(mg,I1,I2,I3,rc,t)
 
#define RZZ(t)   e.zz(mg,I1,I2,I3,rc,t)
 
#define U(t)   e (mg,I1,I2,I3,uc,t)
 
#define UT(t)   e.t (mg,I1,I2,I3,uc,t)
 
#define UX(t)   e.x (mg,I1,I2,I3,uc,t)
 
#define UY(t)   e.y (mg,I1,I2,I3,uc,t)
 
#define UZ(t)   e.z (mg,I1,I2,I3,uc,t)
 
#define UXX(t)   e.xx(mg,I1,I2,I3,uc,t)
 
#define UXY(t)   e.xy(mg,I1,I2,I3,uc,t)
 
#define UXZ(t)   e.xz(mg,I1,I2,I3,uc,t)
 
#define UYY(t)   e.yy(mg,I1,I2,I3,uc,t)
 
#define UYZ(t)   e.yz(mg,I1,I2,I3,uc,t)
 
#define UZZ(t)   e.zz(mg,I1,I2,I3,uc,t)
 
#define V(t)   e (mg,I1,I2,I3,vc,t)
 
#define VT(t)   e.t (mg,I1,I2,I3,vc,t)
 
#define VX(t)   e.x (mg,I1,I2,I3,vc,t)
 
#define VY(t)   e.y (mg,I1,I2,I3,vc,t)
 
#define VZ(t)   e.z (mg,I1,I2,I3,vc,t)
 
#define VXX(t)   e.xx(mg,I1,I2,I3,vc,t)
 
#define VXY(t)   e.xy(mg,I1,I2,I3,vc,t)
 
#define VXZ(t)   e.xz(mg,I1,I2,I3,vc,t)
 
#define VYY(t)   e.yy(mg,I1,I2,I3,vc,t)
 
#define VYZ(t)   e.yz(mg,I1,I2,I3,vc,t)
 
#define VZZ(t)   e.zz(mg,I1,I2,I3,vc,t)
 
#define W(t)   e (mg,I1,I2,I3,wc,t)
 
#define WT(t)   e.t (mg,I1,I2,I3,wc,t)
 
#define WX(t)   e.x (mg,I1,I2,I3,wc,t)
 
#define WY(t)   e.y (mg,I1,I2,I3,wc,t)
 
#define WZ(t)   e.z (mg,I1,I2,I3,wc,t)
 
#define WXX(t)   e.xx(mg,I1,I2,I3,wc,t)
 
#define WXY(t)   e.xy(mg,I1,I2,I3,wc,t)
 
#define WXZ(t)   e.xz(mg,I1,I2,I3,wc,t)
 
#define WYY(t)   e.yy(mg,I1,I2,I3,wc,t)
 
#define WYZ(t)   e.yz(mg,I1,I2,I3,wc,t)
 
#define WZZ(t)   e.zz(mg,I1,I2,I3,wc,t)
 
#define P(t)   (e (mg,I1,I2,I3,pc,t)+pressureLevel)
 
#define PT(t)   e.t (mg,I1,I2,I3,pc,t)
 
#define PX(t)   e.x (mg,I1,I2,I3,pc,t)
 
#define PY(t)   e.y (mg,I1,I2,I3,pc,t)
 
#define PZ(t)   e.z (mg,I1,I2,I3,pc,t)
 
#define PXX(t)   e.xx(mg,I1,I2,I3,pc,t)
 
#define PXY(t)   e.xy(mg,I1,I2,I3,pc,t)
 
#define PXZ(t)   e.xz(mg,I1,I2,I3,pc,t)
 
#define PYY(t)   e.yy(mg,I1,I2,I3,pc,t)
 
#define PYZ(t)   e.yz(mg,I1,I2,I3,pc,t)
 
#define PZZ(t)   e.zz(mg,I1,I2,I3,pc,t)
 
#define T(t)   e (mg,I1,I2,I3,tc,t)
 
#define TT(t)   e.t (mg,I1,I2,I3,tc,t)
 
#define TX(t)   e.x (mg,I1,I2,I3,tc,t)
 
#define TY(t)   e.y (mg,I1,I2,I3,tc,t)
 
#define TZ(t)   e.z (mg,I1,I2,I3,tc,t)
 
#define TXX(t)   e.xx(mg,I1,I2,I3,tc,t)
 
#define TXY(t)   e.xy(mg,I1,I2,I3,tc,t)
 
#define TXZ(t)   e.xz(mg,I1,I2,I3,tc,t)
 
#define TYY(t)   e.yy(mg,I1,I2,I3,tc,t)
 
#define TYZ(t)   e.yz(mg,I1,I2,I3,tc,t)
 
#define TZZ(t)   e.zz(mg,I1,I2,I3,tc,t)
 
#define S(t, s)   e (mg,I1,I2,I3,s,t)
 
#define ST(t, s)   e.t (mg,I1,I2,I3,s,t)
 
#define SX(t, s)   e.x (mg,I1,I2,I3,s,t)
 
#define SY(t, s)   e.y (mg,I1,I2,I3,s,t)
 
#define SZ(t, s)   e.z (mg,I1,I2,I3,s,t)
 
#define SXX(t, s)   e.xx(mg,I1,I2,I3,s,t)
 
#define SXY(t, s)   e.xy(mg,I1,I2,I3,s,t)
 
#define SXZ(t, s)   e.xz(mg,I1,I2,I3,s,t)
 
#define SYY(t, s)   e.yy(mg,I1,I2,I3,s,t)
 
#define SYZ(t, s)   e.yz(mg,I1,I2,I3,s,t)
 
#define SZZ(t, s)   e.zz(mg,I1,I2,I3,s,t)
 
#define E(t)   e (mg,I1,I2,I3,ec,t)
 
#define ET(t)   e.t (mg,I1,I2,I3,ec,t)
 
#define EX(t)   e.x (mg,I1,I2,I3,ec,t)
 
#define EXX(t)   e.xx(mg,I1,I2,I3,ec,t)
 

Functions

void asfdt (const int &nd, const int &n1a, const int &n1b, const int &n2a, const int &n2b, const int &n3a, const int &n3b, const int &nd1a, const int &nd1b, const int &nd2a, const int &nd2b, const int &nd3a, const int &nd3b, const int &nd4a, const int &nd4b, const int &mask, const real &xy, const real &rx, const real &u, const real &uu, real &ut, real &uti, const real &rhoL, const real &pL, const real &gv, const real &dw, const int &bc, const int &ipar, const real &rpar, const int &ierr)
 

Macro Definition Documentation

#define AD2 (   kd)    ( (ad21 + cd22* ( fabs(UX(uc))+fabs(UY(uc)) +fabs(UX(vc))+fabs(UY(vc)) ) ) *(u(I1+1,I2,I3,kd)-4.*u(I1,I2,I3,kd)+u(I1-1,I2,I3,kd) +u(I1,I2+1,I3,kd) +u(I1,I2-1,I3,kd)) )

Referenced by Cgasf::getUt().

#define AD23 (   kd)    ( (ad21 + cd22* ( fabs(UX(uc))+fabs(UY(uc))+fabs(UZ(uc)) +fabs(UX(vc))+fabs(UY(vc))+fabs(UZ(vc)) +fabs(UX(wc))+fabs(UY(wc))+fabs(UZ(wc)) ) ) *(u(I1+1,I2,I3,kd)-6.*u(I1,I2,I3,kd)+u(I1-1,I2,I3,kd) +u(I1,I2+1,I3,kd) +u(I1,I2-1,I3,kd) +u(I1,I2,I3+1,kd) +u(I1,I2,I3-1,kd)) )
#define AD4 (   kd)    ( (ad41 + cd42* ( fabs(UX(uc))+fabs(UY(uc)) +fabs(UX(vc))+fabs(UY(vc)) ) ) *( -u(I1+2,I2,I3,kd)-u(I1-2,I2,I3,kd) -u(I1,I2+2,I3,kd)-u(I1,I2-2,I3,kd) +4.*(u(I1+1,I2,I3,kd)+u(I1-1,I2,I3,kd) +u(I1,I2+1,I3,kd)+u(I1,I2-1,I3,kd)) -12.*u(I1,I2,I3,kd) ) )
#define AD43 (   kd)    ( (ad41 + cd42* ( fabs(UX(uc))+fabs(UY(uc))+fabs(UZ(uc)) +fabs(UX(vc))+fabs(UY(vc))+fabs(UZ(vc)) +fabs(UX(wc))+fabs(UY(wc))+fabs(UZ(wc)) ) ) *( -u(I1+2,I2,I3,kd)-u(I1-2,I2,I3,kd) -u(I1,I2+2,I3,kd)-u(I1,I2-2,I3,kd) -u(I1,I2,I3+2,kd)-u(I1,I2,I3-2,kd) +4.*(u(I1+1,I2,I3,kd)+u(I1-1,I2,I3,kd) +u(I1,I2+1,I3,kd)+u(I1,I2-1,I3,kd) +u(I1,I2,I3+1,kd)+u(I1,I2,I3-1,kd)) -18.*u(I1,I2,I3,kd) ) )
#define asfdt   EXTERN_C_NAME(asfdt)

Referenced by Cgasf::getUt().

#define E (   t)    e (mg,I1,I2,I3,ec,t)
#define ET (   t)    e.t (mg,I1,I2,I3,ec,t)
#define EX (   t)    e.x (mg,I1,I2,I3,ec,t)
#define EXX (   t)    e.xx(mg,I1,I2,I3,ec,t)
#define ForBoundary (   side,
  axis 
)    for( axis=0; axis<c.numberOfDimensions(); axis++ ) for( side=0; side<=1; side++ )
#define P (   t)    (e (mg,I1,I2,I3,pc,t)+pressureLevel)

Referenced by Cgasf::addAllSpeedImplicitForcing(), Cgasf::addForcing(), F90_ID(), and RigidBodyMotion::getBodyForceOption().

#define PT (   t)    e.t (mg,I1,I2,I3,pc,t)

Referenced by Cgasf::addForcing().

#define PX (   t)    e.x (mg,I1,I2,I3,pc,t)

Referenced by Cgasf::addAllSpeedImplicitForcing(), and Cgasf::addForcing().

#define PXX (   t)    e.xx(mg,I1,I2,I3,pc,t)
#define PXY (   t)    e.xy(mg,I1,I2,I3,pc,t)
#define PXZ (   t)    e.xz(mg,I1,I2,I3,pc,t)
#define PY (   t)    e.y (mg,I1,I2,I3,pc,t)

Referenced by Cgasf::addAllSpeedImplicitForcing(), and Cgasf::addForcing().

#define PYY (   t)    e.yy(mg,I1,I2,I3,pc,t)
#define PYZ (   t)    e.yz(mg,I1,I2,I3,pc,t)
#define PZ (   t)    e.z (mg,I1,I2,I3,pc,t)

Referenced by Cgasf::addAllSpeedImplicitForcing(), and Cgasf::addForcing().

#define PZZ (   t)    e.zz(mg,I1,I2,I3,pc,t)
#define R (   t)    e (mg,I1,I2,I3,rc,t)

Referenced by Cgasf::addAllSpeedImplicitForcing(), and Cgasf::addForcing().

#define RT (   t)    e.t (mg,I1,I2,I3,rc,t)

Referenced by Cgasf::addForcing().

#define RX (   t)    e.x (mg,I1,I2,I3,rc,t)

Referenced by Cgasf::addForcing().

#define RXX (   t)    e.xx(mg,I1,I2,I3,rc,t)

Referenced by Cgasf::addForcing().

#define RY (   t)    e.y (mg,I1,I2,I3,rc,t)

Referenced by Cgasf::addForcing().

#define RYY (   t)    e.yy(mg,I1,I2,I3,rc,t)

Referenced by Cgasf::addForcing().

#define RZ (   t)    e.z (mg,I1,I2,I3,rc,t)
#define RZZ (   t)    e.zz(mg,I1,I2,I3,rc,t)
#define S (   t,
  s 
)    e (mg,I1,I2,I3,s,t)

Referenced by Cgasf::allSpeedImplicitTimeStep(), Cgcns::applyBoundaryConditions(), Reactions::checkChemicalSourceJacobian(), Cgsm::checkDisplacementAndStress(), Cgasf::computeSource(), Cgins::getLineSolverBoundaryConditions(), Parameters::getUserDefinedKnownSolution(), Cgasf::getUt(), Cgsm::getVelocityAndStress(), Cgins::lineSolverBoundaryConditions(), Parameters::readFromAShowFile(), DomainSolver::saveSequencesToShowFile(), Maxwell::setBoundaryCondition(), Cgsm::setBoundaryCondition(), DomainSolver::setBoundaryConditionsInteractively(), Parameters::setBoundaryConditionValues(), Reactions::solveImplicitForPTYGivenR(), Reactions::solveImplicitForRTYGivenP(), Reactions::solveImplicitForYGivenRTP(), Parameters::thereAreTimeDependentUserBoundaryConditions(), SurfaceEquation::update(), and DomainSolver::userDefinedInitialConditions().

#define ST (   t,
  s 
)    e.t (mg,I1,I2,I3,s,t)
#define SX (   t,
  s 
)    e.x (mg,I1,I2,I3,s,t)
#define SXX (   t,
  s 
)    e.xx(mg,I1,I2,I3,s,t)
#define SXY (   t,
  s 
)    e.xy(mg,I1,I2,I3,s,t)
#define SXZ (   t,
  s 
)    e.xz(mg,I1,I2,I3,s,t)
#define SY (   t,
  s 
)    e.y (mg,I1,I2,I3,s,t)
#define SYY (   t,
  s 
)    e.yy(mg,I1,I2,I3,s,t)
#define SYZ (   t,
  s 
)    e.yz(mg,I1,I2,I3,s,t)
#define SZ (   t,
  s 
)    e.z (mg,I1,I2,I3,s,t)
#define SZZ (   t,
  s 
)    e.zz(mg,I1,I2,I3,s,t)
#define T (   t)    e (mg,I1,I2,I3,tc,t)

Referenced by Cgasf::addForcing(), and MovingGrids::rigidBodyMotion().

#define TT (   t)    e.t (mg,I1,I2,I3,tc,t)

Referenced by Cgasf::addForcing().

#define TX (   t)    e.x (mg,I1,I2,I3,tc,t)

Referenced by Cgasf::addForcing().

#define TXX (   t)    e.xx(mg,I1,I2,I3,tc,t)

Referenced by Cgasf::addForcing().

#define TXY (   t)    e.xy(mg,I1,I2,I3,tc,t)
#define TXZ (   t)    e.xz(mg,I1,I2,I3,tc,t)
#define TY (   t)    e.y (mg,I1,I2,I3,tc,t)

Referenced by Cgasf::addForcing().

#define TYY (   t)    e.yy(mg,I1,I2,I3,tc,t)

Referenced by Cgasf::addForcing().

#define TYZ (   t)    e.yz(mg,I1,I2,I3,tc,t)
#define TZ (   t)    e.z (mg,I1,I2,I3,tc,t)

Referenced by Cgasf::addForcing().

#define TZZ (   t)    e.zz(mg,I1,I2,I3,tc,t)

Referenced by Cgasf::addForcing().

#define U (   c)    u(I1,I2,I3,c)
#define U (   t)    e (mg,I1,I2,I3,uc,t)
#define UT (   t)    e.t (mg,I1,I2,I3,uc,t)

Referenced by Cgasf::addForcing().

#define UU (   c)    uu(I1,I2,I3,c)
#define UX (   c)    ux(I1,I2,I3,c)
#define UX (   t)    e.x (mg,I1,I2,I3,uc,t)
#define UXX (   c)    uxx(I1,I2,I3,c)
#define UXX (   t)    e.xx(mg,I1,I2,I3,uc,t)
#define UXY (   c)    uxy(I1,I2,I3,c)
#define UXY (   t)    e.xy(mg,I1,I2,I3,uc,t)
#define UXZ (   c)    uxz(I1,I2,I3,c)
#define UXZ (   t)    e.xz(mg,I1,I2,I3,uc,t)
#define UY (   c)    uy(I1,I2,I3,c)
#define UY (   t)    e.y (mg,I1,I2,I3,uc,t)
#define UYY (   c)    uyy(I1,I2,I3,c)
#define UYY (   t)    e.yy(mg,I1,I2,I3,uc,t)
#define UYZ (   c)    uyz(I1,I2,I3,c)
#define UYZ (   t)    e.yz(mg,I1,I2,I3,uc,t)
#define UZ (   c)    uz(I1,I2,I3,c)
#define UZ (   t)    e.z (mg,I1,I2,I3,uc,t)
#define UZZ (   c)    uzz(I1,I2,I3,c)
#define UZZ (   t)    e.zz(mg,I1,I2,I3,uc,t)
#define V (   t)    e (mg,I1,I2,I3,vc,t)

Referenced by Cgasf::addForcing().

#define VT (   t)    e.t (mg,I1,I2,I3,vc,t)

Referenced by Cgasf::addForcing().

#define VX (   t)    e.x (mg,I1,I2,I3,vc,t)

Referenced by Cgasf::addForcing().

#define VXX (   t)    e.xx(mg,I1,I2,I3,vc,t)

Referenced by Cgasf::addForcing().

#define VXY (   t)    e.xy(mg,I1,I2,I3,vc,t)

Referenced by Cgasf::addForcing().

#define VXZ (   t)    e.xz(mg,I1,I2,I3,vc,t)
#define VY (   t)    e.y (mg,I1,I2,I3,vc,t)

Referenced by Cgasf::addAllSpeedImplicitForcing(), and Cgasf::addForcing().

#define VYY (   t)    e.yy(mg,I1,I2,I3,vc,t)

Referenced by Cgasf::addForcing().

#define VYZ (   t)    e.yz(mg,I1,I2,I3,vc,t)

Referenced by Cgasf::addForcing().

#define VZ (   t)    e.z (mg,I1,I2,I3,vc,t)

Referenced by Cgasf::addForcing().

#define VZZ (   t)    e.zz(mg,I1,I2,I3,vc,t)

Referenced by Cgasf::addForcing().

#define W (   t)    e (mg,I1,I2,I3,wc,t)

Referenced by Cgasf::addForcing(), and getCrossProductMatrix().

#define WT (   t)    e.t (mg,I1,I2,I3,wc,t)

Referenced by Cgasf::addForcing().

#define WX (   t)    e.x (mg,I1,I2,I3,wc,t)

Referenced by Cgasf::addForcing().

#define WXX (   t)    e.xx(mg,I1,I2,I3,wc,t)

Referenced by Cgasf::addForcing().

#define WXY (   t)    e.xy(mg,I1,I2,I3,wc,t)
#define WXZ (   t)    e.xz(mg,I1,I2,I3,wc,t)

Referenced by Cgasf::addForcing().

#define WY (   t)    e.y (mg,I1,I2,I3,wc,t)

Referenced by Cgasf::addForcing().

#define WYY (   t)    e.yy(mg,I1,I2,I3,wc,t)

Referenced by Cgasf::addForcing().

#define WYZ (   t)    e.yz(mg,I1,I2,I3,wc,t)

Referenced by Cgasf::addForcing().

#define WZ (   t)    e.z (mg,I1,I2,I3,wc,t)

Referenced by Cgasf::addAllSpeedImplicitForcing(), and Cgasf::addForcing().

#define WZZ (   t)    e.zz(mg,I1,I2,I3,wc,t)

Referenced by Cgasf::addForcing().

Function Documentation

void asfdt ( const int &  nd,
const int &  n1a,
const int &  n1b,
const int &  n2a,
const int &  n2b,
const int &  n3a,
const int &  n3b,
const int &  nd1a,
const int &  nd1b,
const int &  nd2a,
const int &  nd2b,
const int &  nd3a,
const int &  nd3b,
const int &  nd4a,
const int &  nd4b,
const int &  mask,
const real &  xy,
const real &  rx,
const real &  u,
const real &  uu,
real &  ut,
real &  uti,
const real &  rhoL,
const real &  pL,
const real &  gv,
const real &  dw,
const int &  bc,
const int &  ipar,
const real &  rpar,
const int &  ierr 
)
Generated on Fri Jan 4 2013 10:19:30 for CG by   doxygen 1.8.3