CG  Version 25
Macros | Functions | Variables
insp.temp.C File Reference
#include "OB_CompositeGridSolver.h"
#include "OB_MappedGridSolver.h"
#include "Parameters.h"
#include "MappedGridOperators.h"
#include "Ogmg.h"
#include "ParallelUtility.h"
#include "SparseRep.h"
#include <float.h>
#include "turbulenceParameters.h"
#include "viscoPlasticMacrosCpp.h"
#include "EquationDomain.h"
#include "SurfaceEquation.h"
#include "Cgins.h"
Include dependency graph for insp.temp.C:

Macros

#define POW2(x)   pow((x),2)
 
#define ForBoundary(side, axis)
 
#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 PN1(I1, I2, I3)
 
#define PXB2(I1, I2, I3)
 
#define PYB2(I1, I2, I3)
 
#define PN2(I1, I2, I3)
 
#define DELTAU(I1, I2, I3, dir)   (uxx(I1,I2,I3,dir)+uyy(I1,I2,I3,dir)+uzz(I1,I2,I3,dir))
 
#define P3B(I1, I2, I3, dir)
 
#define PN3(I1, I2, I3)
 
#define realSmall   (REAL_MIN*1.e5)
 
#define DAI1(cd, I1, I2, I3)
 
#define DAI2(cd, I1, I2, I3)
 
#define DAI3(cd, I1, I2, I3)
 
#define D2UW2(I1, I2, I3)
 
#define AD2(kd)
 
#define AD23(kd)
 
#define AD4(kd)
 
#define AD43(kd)
 
#define assignPressureRHSOpt   assignpressurerhsopt_
 
#define P0(c, I1, I2, I3, t)   e(c,I1,I2,I3,pc,t)
 
#define P02N(c, I1, I2, I3, t)
 
#define P03N(c, I1, I2, I3, t)
 
#define PF2(c, I1, I2, I3, t)
 
#define PF3(c, I1, I2, I3, t)
 
#define FB21(c, I1, I2, I3, t)
 
#define FB22(c, I1, I2, I3, t)
 
#define FB2N(c, I1, I2, I3, t)
 
#define FB31(c, I1, I2, I3, t)
 
#define FB32(c, I1, I2, I3, t)
 
#define FB33(c, I1, I2, I3, t)
 
#define FB3N(c, I1, I2, I3, t)
 

Functions

void getLocalBoundsAndBoundaryConditions (const realMappedGridFunction &a, IntegerArray &gidLocal, IntegerArray &dimensionLocal, IntegerArray &bcLocal)
 
void assignPressureRHSOpt (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 &rx, const real &u, const real &uu, real &f, const real &gv, const real &divDamping, const int &bc, const int &indexRange, const int &ndp, const real &pressureValue, const int &nr1a, const int &nr1b, const int &nr2a, const int &nr2b, const int &nr3a, const int &nr3b, const real &normal00, const real &normal10, const real &normal01, const real &normal11, const real &normal02, const real &normal12, const int &ipar, const real &rpar, const int &ierr)
 

Variables

ListOfEquationDomains equationDomainList
 
SurfaceEquation surfaceEquation
 
DomainSolver * pDomainSolver
 

Macro Definition Documentation

#define AD2 (   kd)
Value:
( \
(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)
Value:
( \
(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)
Value:
( \
(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)
Value:
( \
(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 assignPressureRHSOpt   assignpressurerhsopt_

Referenced by Cgins::assignPressureRHS().

#define D2UW2 (   I1,
  I2,
  I3 
)
Value:
(UW(I1+1,I2,I3)+UW(I1-1,I2,I3) \
+UW(I1,I2+1,I3)+UW(I1,I2-1,I3)-4.*UW(I1,I2,I3))
#define DAI1 (   cd,
  I1,
  I2,
  I3 
)
Value:
( \
cd/ max(realSmall,\
( POW2(xy(I1+1,I2 ,I3 ,0)-xy(I1-1,I2 ,I3 ,0)) \
+POW2(xy(I1+1,I2 ,I3 ,1)-xy(I1-1,I2 ,I3 ,1)) ) ) )
#define DAI2 (   cd,
  I1,
  I2,
  I3 
)
Value:
( \
cd/max(realSmall,\
( POW2(xy(I1+1,I2 ,I3 ,0)-xy(I1-1,I2 ,I3 ,0)) \
+POW2(xy(I1+1,I2 ,I3 ,1)-xy(I1-1,I2 ,I3 ,1)) ) ) \
+ cd/max(realSmall, \
( POW2(xy(I1 ,I2+1,I3 ,0)-xy(I1 ,I2-1,I3 ,0)) \
+POW2(xy(I1 ,I2+1,I3 ,1)-xy(I1 ,I2-1,I3 ,1)) ) ) )
#define DAI3 (   cd,
  I1,
  I2,
  I3 
)
Value:
( \
cd/max(realSmall, \
( POW2(xy(I1+1,I2 ,I3 ,0)-xy(I1-1,I2 ,I3 ,0)) \
+POW2(xy(I1+1,I2 ,I3 ,1)-xy(I1-1,I2 ,I3 ,1)) \
+POW2(xy(I1+1,I2 ,I3 ,2)-xy(I1-1,I2 ,I3 ,2)) ) ) \
+ cd/max(realSmall, \
( POW2(xy(I1 ,I2+1,I3 ,0)-xy(I1 ,I2-1,I3 ,0)) \
+POW2(xy(I1 ,I2+1,I3 ,1)-xy(I1 ,I2-1,I3 ,1)) \
+POW2(xy(I1 ,I2+1,I3 ,2)-xy(I1 ,I2-1,I3 ,2)) ) ) \
+ cd/max(realSmall, \
( POW2(xy(I1 ,I2 ,I3+1,0)-xy(I1 ,I2 ,I3-1,0)) \
+POW2(xy(I1 ,I2 ,I3+1,1)-xy(I1 ,I2 ,I3-1,1)) \
+POW2(xy(I1 ,I2 ,I3+1,2)-xy(I1 ,I2 ,I3-1,2)) ) ) )
#define DELTAU (   I1,
  I2,
  I3,
  dir 
)    (uxx(I1,I2,I3,dir)+uyy(I1,I2,I3,dir)+uzz(I1,I2,I3,dir))
#define FB21 (   c,
  I1,
  I2,
  I3,
  t 
)
Value:
( \
u0t(I1,I2,I3) \
+uuLocal(I1,I2,I3,uc)*u0x(I1,I2,I3) \
+uuLocal(I1,I2,I3,vc)*u0y(I1,I2,I3) \
+p0x(I1,I2,I3) \
-nu*(u0xx+u0yy) \
)
#define FB22 (   c,
  I1,
  I2,
  I3,
  t 
)
Value:
( \
v0t(I1,I2,I3) \
+uuLocal(I1,I2,I3,uc)*v0x(I1,I2,I3) \
+uuLocal(I1,I2,I3,vc)*v0y(I1,I2,I3) \
+p0y(I1,I2,I3) \
-nu*(v0xx+v0yy) \
)
#define FB2N (   c,
  I1,
  I2,
  I3,
  t 
)
Value:
( normal(I1,I2,I3,0)*FB21(c,I1,I2,I3,t) \
+normal(I1,I2,I3,1)*FB22(c,I1,I2,I3,t) )
#define FB31 (   c,
  I1,
  I2,
  I3,
  t 
)
Value:
( \
u0t(I1,I2,I3) \
+uuLocal(I1,I2,I3,uc)*u0x(I1,I2,I3) \
+uuLocal(I1,I2,I3,vc)*u0y(I1,I2,I3) \
+uuLocal(I1,I2,I3,wc)*u0z(I1,I2,I3) \
+p0x(I1,I2,I3) \
-nu*(u0xx+u0yy+u0zz) \
)
#define FB32 (   c,
  I1,
  I2,
  I3,
  t 
)
Value:
( \
v0t(I1,I2,I3) \
+uuLocal(I1,I2,I3,uc)*v0x(I1,I2,I3) \
+uuLocal(I1,I2,I3,vc)*v0y(I1,I2,I3) \
+uuLocal(I1,I2,I3,wc)*v0z(I1,I2,I3) \
+p0y(I1,I2,I3) \
-nu*(v0xx+v0yy+v0zz) \
)
#define FB33 (   c,
  I1,
  I2,
  I3,
  t 
)
Value:
( \
w0t(I1,I2,I3) \
+uuLocal(I1,I2,I3,uc)*w0x(I1,I2,I3) \
+uuLocal(I1,I2,I3,vc)*w0y(I1,I2,I3) \
+uuLocal(I1,I2,I3,wc)*w0z(I1,I2,I3) \
+p0z(I1,I2,I3) \
-nu*(w0xx+w0yy+w0zz) \
)
#define FB3N (   c,
  I1,
  I2,
  I3,
  t 
)
Value:
( normal(I1,I2,I3,0)*FB31(c,I1,I2,I3,t) \
+normal(I1,I2,I3,1)*FB32(c,I1,I2,I3,t) \
+normal(I1,I2,I3,2)*FB33(c,I1,I2,I3,t) )
#define ForBoundary (   side,
  axis 
)
Value:
for( int axis=0; axis<c.numberOfDimensions(); axis++ ) \
for( int side=0; side<=1; side++ )
#define P0 (   c,
  I1,
  I2,
  I3,
  t 
)    e(c,I1,I2,I3,pc,t)

Referenced by Cgsm::assignHempInitialConditions(), F90_ID(), and Cgasf::getUt().

#define P02N (   c,
  I1,
  I2,
  I3,
  t 
)
Value:
( normal(I1,I2,I3,0)*p0x(I1,I2,I3) \
+normal(I1,I2,I3,1)*p0y(I1,I2,I3) )
#define P03N (   c,
  I1,
  I2,
  I3,
  t 
)
Value:
( normal(I1,I2,I3,0)*p0x(I1,I2,I3) \
+normal(I1,I2,I3,1)*p0y(I1,I2,I3) \
+normal(I1,I2,I3,2)*p0z(I1,I2,I3) )
#define P3B (   I1,
  I2,
  I3,
  dir 
)
Value:
( nu*DELTAU(I1,I2,I3,dir) \
+advectionCoefficient*( uu(I1,I2,I3,uc)*ux(I1,I2,I3,dir) \
+ uu(I1,I2,I3,vc)*uy(I1,I2,I3,dir) \
+ uu(I1,I2,I3,wc)*uz(I1,I2,I3,dir)) )
#define PF2 (   c,
  I1,
  I2,
  I3,
  t 
)
Value:
( \
advectionCoefficient*( u0x*u0x+2.*u0y*v0x+v0y*v0y ) + p0xx + p0yy )
#define PF3 (   c,
  I1,
  I2,
  I3,
  t 
)
Value:
( \
advectionCoefficient*( u0x*u0x+v0y*v0y+w0z*w0z+2.*(u0y*v0x+u0z*w0x+v0z*w0y) )+p0xx+p0yy+p0zz \
)
#define PN1 (   I1,
  I2,
  I3 
)
Value:
( (2*side-1)*( nu*uxx(I1,I2,I3,uc) \
+ advectionCoefficient*uu(I1,I2,I3,uc)*ux(I1,I2,I3,uc) ) )

Referenced by Cgasf::solveForTimeIndependentVariables().

#define PN2 (   I1,
  I2,
  I3 
)
Value:
( normal(I1,I2,I3,0)*PXB2(I1,I2,I3) \
+normal(I1,I2,I3,1)*PYB2(I1,I2,I3) )

Referenced by Cgasf::addForcingToPressureEquation(), and Cgasf::solveForTimeIndependentVariables().

#define PN3 (   I1,
  I2,
  I3 
)
Value:
( normal(I1,I2,I3,0)*P3B(I1,I2,I3,uc) \
+normal(I1,I2,I3,1)*P3B(I1,I2,I3,vc) \
+normal(I1,I2,I3,2)*P3B(I1,I2,I3,wc) )

Referenced by Cgasf::solveForTimeIndependentVariables().

#define POW2 (   x)    pow((x),2)
#define PXB2 (   I1,
  I2,
  I3 
)
Value:
( nu*(uxx(I1,I2,I3,uc)+uyy(I1,I2,I3,uc)) \
+advectionCoefficient*( uu(I1,I2,I3,uc)*ux(I1,I2,I3,uc) + uu(I1,I2,I3,vc)*uy(I1,I2,I3,uc)) )
#define PYB2 (   I1,
  I2,
  I3 
)
Value:
( nu*(uxx(I1,I2,I3,vc)+uyy(I1,I2,I3,vc)) \
+advectionCoefficient*( uu(I1,I2,I3,uc)*ux(I1,I2,I3,vc) + uu(I1,I2,I3,vc)*uy(I1,I2,I3,vc)) )
#define realSmall   (REAL_MIN*1.e5)

Referenced by Cgins::updateDivergenceDamping().

#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 UXX (   c)    uxx(I1,I2,I3,c)
#define UXY (   c)    uxy(I1,I2,I3,c)
#define UXZ (   c)    uxz(I1,I2,I3,c)
#define UY (   c)    uy(I1,I2,I3,c)
#define UYY (   c)    uyy(I1,I2,I3,c)
#define UYZ (   c)    uyz(I1,I2,I3,c)
#define UZ (   c)    uz(I1,I2,I3,c)
#define UZZ (   c)    uzz(I1,I2,I3,c)

Function Documentation

void assignPressureRHSOpt ( 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 &  rx,
const real &  u,
const real &  uu,
real &  f,
const real &  gv,
const real &  divDamping,
const int &  bc,
const int &  indexRange,
const int &  ndp,
const real &  pressureValue,
const int &  nr1a,
const int &  nr1b,
const int &  nr2a,
const int &  nr2b,
const int &  nr3a,
const int &  nr3b,
const real &  normal00,
const real &  normal10,
const real &  normal01,
const real &  normal11,
const real &  normal02,
const real &  normal12,
const int &  ipar,
const real &  rpar,
const int &  ierr 
)
void getLocalBoundsAndBoundaryConditions ( const realMappedGridFunction &  a,
IntegerArray &  gidLocal,
IntegerArray &  dimensionLocal,
IntegerArray &  bcLocal 
)

Variable Documentation

ListOfEquationDomains equationDomainList
DomainSolver* pDomainSolver
SurfaceEquation surfaceEquation
Generated on Fri Jan 4 2013 10:19:32 for CG by   doxygen 1.8.3