CG  Version 25
Macros | Functions
pressureEquation.C File Reference
#include "Cgins.h"
#include "Parameters.h"
#include "MappedGridOperators.h"
#include "Ogmg.h"
#include "ParallelUtility.h"
#include "ParallelGridUtility.h"
#include "SparseRep.h"
#include "HDF_DataBase.h"
#include "App.h"
#include "turbulenceParameters.h"
#include "EquationDomain.h"
Include dependency graph for pressureEquation.C:

Macros

#define POW2(x)   pow((x),2)
 
#define ForBoundary(side, axis)
 
#define assignPressureRHSOpt   EXTERN_C_NAME(assignpressurerhsopt)
 
#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 AD2(kd)
 
#define mixedRHS(component, side, axis, grid)   bcData(component+numberOfComponents*(0),side,axis,grid)
 
#define mixedCoeff(component, side, axis, grid)   bcData(component+numberOfComponents*(1),side,axis,grid)
 
#define mixedNormalCoeff(component, side, axis, grid)   bcData(component+numberOfComponents*(2),side,axis,grid)
 
#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 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 &xy, const real &rx, const real &radiusInverse, const real &u, const real &uu, real &f, const real &gv, const real &divDamping, const real &userDefinedForce, 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)
 

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 assignPressureRHSOpt   EXTERN_C_NAME(assignpressurerhsopt)
#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) )

Referenced by Cgins::addForcingToPressureEquation().

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

Referenced by Cgins::addForcingToPressureEquation().

#define ForBoundary (   side,
  axis 
)
Value:
for( int axis=0; axis<c.numberOfDimensions(); axis++ ) \
for( int side=0; side<=1; side++ )

Referenced by Cgins::addForcingToPressureEquation(), and Cgins::assignPressureRHS().

#define mixedCoeff (   component,
  side,
  axis,
  grid 
)    bcData(component+numberOfComponents*(1),side,axis,grid)
#define mixedNormalCoeff (   component,
  side,
  axis,
  grid 
)    bcData(component+numberOfComponents*(2),side,axis,grid)
#define mixedRHS (   component,
  side,
  axis,
  grid 
)    bcData(component+numberOfComponents*(0),side,axis,grid)
#define P0 (   c,
  I1,
  I2,
  I3,
  t 
)    e(c,I1,I2,I3,pc,t)
#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) )

Referenced by Cgins::addForcingToPressureEquation().

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

Referenced by Cgins::addForcingToPressureEquation().

#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) ) )
#define PN2 (   I1,
  I2,
  I3 
)
Value:
( normal(I1,I2,I3,0)*PXB2(I1,I2,I3) \
+normal(I1,I2,I3,1)*PYB2(I1,I2,I3) )
#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) )
#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 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 &  xy,
const real &  rx,
const real &  radiusInverse,
const real &  u,
const real &  uu,
real &  f,
const real &  gv,
const real &  divDamping,
const real &  userDefinedForce,
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 
)
Generated on Fri Jan 4 2013 10:19:32 for CG by   doxygen 1.8.3