CG  Version 25
Macros | Functions
ins.temp.C File Reference
#include "OB_CompositeGridSolver.h"
#include "OB_MappedGridSolver.h"
#include "Parameters.h"
#include "MappedGridOperators.h"
#include "SparseRep.h"
#include "Ogmg.h"
#include "display.h"
#include "HDF_DataBase.h"
#include "ParallelUtility.h"
#include "turbulenceModels.h"
#include "turbulenceParameters.h"
#include "viscoPlasticMacrosCpp.h"
Include dependency graph for ins.temp.C:

Macros

#define insdt   insdt_
 
#define insdts   insdts_
 
#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 AD2(kd)
 
#define AD2N(kd)
 
#define AD23(kd)
 
#define AD23N(kd)
 
#define AD4(kd)
 
#define AD43(kd)
 
#define RX(i1, i2, i3, m, n)   inverseVertexDerivative(i1,i2,i3,m+numberOfDimensions*(n))
 
#define uxmzzR(i1, i2, i3, c)   (u(i1,i2,i3,c)-u(i1-1,i2,i3,c))*dxi
 
#define uymzzR(i1, i2, i3, c)   (u(i1,i2+1,i3,c)-u(i1,i2-1,i3,c)+u(i1-1,i2+1,i3,c)-u(i1-1,i2-1,i3,c))*dyi*.25
 
#define uzmzzR(i1, i2, i3, c)   (u(i1,i2,i3+1,c)-u(i1,i2,i3-1,c)+u(i1-1,i2,i3+1,c)-u(i1-1,i2,i3-1,c))*dzi*.25
 
#define uxzmzR(i1, i2, i3, c)   (u(i1+1,i2,i3,c)-u(i1-1,i2,i3,c)+u(i1+1,i2-1,i3,c)-u(i1-1,i2-1,i3,c))*dxi*.25
 
#define uyzmzR(i1, i2, i3, c)   (u(i1,i2,i3,c)-u(i1,i2-1,i3,c))*dyi
 
#define uzzmzR(i1, i2, i3, c)   (u(i1,i2,i3+1,c)-u(i1,i2,i3-1,c)+u(i1,i2-1,i3+1,c)-u(i1,i2-1,i3-1,c))*dzi*.25
 
#define uxzzmR(i1, i2, i3, c)   (u(i1+1,i2,i3,c)-u(i1-1,i2,i3,c)+u(i1+1,i2,i3-1,c)-u(i1-1,i2,i3-1,c))*dxi*.25
 
#define uyzzmR(i1, i2, i3, c)   (u(i1,i2+1,i3,c)-u(i1,i2-1,i3,c)+u(i1,i2+1,i3-1,c)-u(i1,i2-1,i3-1,c))*dyi*.25
 
#define uzzzmR(i1, i2, i3, c)   (u(i1,i2,i3,c)-u(i1,i2,i3-1,c))*dzi
 
#define udmzC(i1, i2, i3, m, c)
 
#define udzmC(i1, i2, i3, m, c)
 
#define udmzzC(i1, i2, i3, m, c)
 
#define udzmzC(i1, i2, i3, m, c)
 
#define udzzmC(i1, i2, i3, m, c)
 
#define admzR(i1, i2, i3)
 
#define adzmR(i1, i2, i3)
 
#define admzzR(i1, i2, i3)
 
#define adzmzR(i1, i2, i3)
 
#define adzzmR(i1, i2, i3)
 
#define admzC(i1, i2, i3)
 
#define adzmC(i1, i2, i3)
 
#define admzzC(i1, i2, i3)
 
#define adzmzC(i1, i2, i3)
 
#define adzzmC(i1, i2, i3)
 
#define admzRSA(i1, i2, i3)   ad21+cd22*( fabs(uxmzzR(i1,i2,i3,nc))+fabs(uymzzR(i1,i2,i3,nc)) )
 
#define adzmRSA(i1, i2, i3)   ad21+cd22*( fabs(uxzmzR(i1,i2,i3,nc))+fabs(uyzmzR(i1,i2,i3,nc)) )
 
#define admzzRSA(i1, i2, i3)   ad21+cd22*( fabs(uxmzzR(i1,i2,i3,nc))+fabs(uymzzR(i1,i2,i3,nc))+fabs(uzmzzR(i1,i2,i3,nc)) )
 
#define adzmzRSA(i1, i2, i3)   ad21+cd22*( fabs(uxzmzR(i1,i2,i3,nc))+fabs(uyzmzR(i1,i2,i3,nc))+fabs(uzzmzR(i1,i2,i3,nc)) )
 
#define adzzmRSA(i1, i2, i3)   ad21+cd22*( fabs(uxzzmR(i1,i2,i3,nc))+fabs(uyzzmR(i1,i2,i3,nc))+fabs(uzzzmR(i1,i2,i3,nc)) )
 
#define admzCSA(i1, i2, i3)   ad21+cd22*( fabs(udmzC(i1,i2,i3,0,nc))+fabs(udmzC(i1,i2,i3,1,nc)) )
 
#define adzmCSA(i1, i2, i3)   ad21+cd22*( fabs(udzmC(i1,i2,i3,0,nc))+fabs(udzmC(i1,i2,i3,1,nc)) )
 
#define admzzCSA(i1, i2, i3)   ad21+cd22*( fabs(udmzzC(i1,i2,i3,0,nc))+fabs(udmzzC(i1,i2,i3,1,nc))+fabs(udmzzC(i1,i2,i3,2,nc)) )
 
#define adzmzCSA(i1, i2, i3)   ad21+cd22*( fabs(udzmzC(i1,i2,i3,0,nc))+fabs(udzmzC(i1,i2,i3,1,nc))+fabs(udzmzC(i1,i2,i3,2,nc)) )
 
#define adzzmCSA(i1, i2, i3)   ad21+cd22*( fabs(udzzmC(i1,i2,i3,0,nc))+fabs(udzzmC(i1,i2,i3,1,nc))+fabs(udzzmC(i1,i2,i3,2,nc)) )
 
#define AD2F(m)
 
#define AD3F(m)
 
#define uv(m)   e (mg,I1,I2,I3,m,t)
 
#define uvt(m)   e.t (mg,I1,I2,I3,m,t)
 
#define uvx(m)   e.x (mg,I1,I2,I3,m,t)
 
#define uvy(m)   e.y (mg,I1,I2,I3,m,t)
 
#define uvz(m)   e.z (mg,I1,I2,I3,m,t)
 
#define uvxx(m)   e.xx(mg,I1,I2,I3,m,t)
 
#define uvxy(m)   e.xy(mg,I1,I2,I3,m,t)
 
#define uvxz(m)   e.xz(mg,I1,I2,I3,m,t)
 
#define uvyy(m)   e.yy(mg,I1,I2,I3,m,t)
 
#define uvyz(m)   e.yz(mg,I1,I2,I3,m,t)
 
#define uvzz(m)   e.zz(mg,I1,I2,I3,m,t)
 
#define ADTZ2(cc, q0)
 
#define ADTZ4(cc, q0)
 
#define ADTZ2(cc, q0)
 
#define ADTZ4(cc, q0)
 
#define M123(m1, m2, m3)   (m1+halfWidth1+width*(m2+halfWidth2+width*(m3+halfWidth3)))
 
#define M123N(m1, m2, m3, n)   (M123(m1,m2,m3)+stencilSize*(n))
 
#define CE(c, e)   (stencilSize*((c)+numberOfComponentsForCoefficients*(e)))
 
#define M123CE(m1, m2, m3, c, e)   (M123(m1,m2,m3)+CE(c,e))
 
#define MCE(m1, m2, m3)   (M123(m1,m2,m3)+CE(c0,e0))
 

Functions

void insdt (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 &gv, const real &dw, const int &bc, const int &ipar, const real &rpar, const int &ierr)
 
void insdts (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, const real &gv, const real &dw, const real &divDamp, const real &dtVar, const int &bc, const int &ipar, const real &rpar, const int &ierr)
 
int outputCompositeGrid (CompositeGrid &cg, const aString &gridFileName)
 

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 AD23N (   kd)
Value:
( \
(ad21 + cd22*( fabs(UX(nc))+fabs(UY(nc))+fabs(UZ(nc)) ) )\
*(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 AD2F (   m)
Value:
adf0(I1+1,I2,I3)*(u(I1+1,I2,I3,m)-u(I1,I2,I3,m))-adf0(I1,I2,I3)*(u(I1,I2,I3,m)-u(I1-1,I2,I3,m))+\
adf1(I1,I2+1,I3)*(u(I1,I2+1,I3,m)-u(I1,I2,I3,m))-adf1(I1,I2,I3)*(u(I1,I2,I3,m)-u(I1,I2-1,I3,m))
#define AD2N (   kd)
Value:
( \
(ad21 + cd22*( fabs(UX(nc))+fabs(UY(nc)) ) )\
*(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 AD3F (   m)
Value:
adf0(I1+1,I2,I3)*(u(I1+1,I2,I3,m)-u(I1,I2,I3,m))-adf0(I1,I2,I3)*(u(I1,I2,I3,m)-u(I1-1,I2,I3,m))+\
adf1(I1,I2+1,I3)*(u(I1,I2+1,I3,m)-u(I1,I2,I3,m))-adf1(I1,I2,I3)*(u(I1,I2,I3,m)-u(I1,I2-1,I3,m))+\
adf2(I1,I2,I3+1)*(u(I1,I2,I3+1,m)-u(I1,I2,I3,m))-adf2(I1,I2,I3)*(u(I1,I2,I3,m)-u(I1,I2,I3-1,m))
#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 admzC (   i1,
  i2,
  i3 
)
Value:
ad21+cd22*( fabs(udmzC(i1,i2,i3,0,uc))+fabs(udmzC(i1,i2,i3,0,vc))+ \
fabs(udmzC(i1,i2,i3,1,uc))+fabs(udmzC(i1,i2,i3,1,vc)) )
#define admzCSA (   i1,
  i2,
  i3 
)    ad21+cd22*( fabs(udmzC(i1,i2,i3,0,nc))+fabs(udmzC(i1,i2,i3,1,nc)) )
#define admzR (   i1,
  i2,
  i3 
)
Value:
ad21+cd22*( fabs(uxmzzR(i1,i2,i3,uc))+fabs(uxmzzR(i1,i2,i3,vc))+\
fabs(uymzzR(i1,i2,i3,uc))+fabs(uymzzR(i1,i2,i3,vc)) )
#define admzRSA (   i1,
  i2,
  i3 
)    ad21+cd22*( fabs(uxmzzR(i1,i2,i3,nc))+fabs(uymzzR(i1,i2,i3,nc)) )
#define admzzC (   i1,
  i2,
  i3 
)
Value:
ad21+cd22*( fabs(udmzzC(i1,i2,i3,0,uc))+fabs(udmzzC(i1,i2,i3,0,vc))+fabs(udmzzC(i1,i2,i3,0,wc))+\
fabs(udmzzC(i1,i2,i3,1,uc))+fabs(udmzzC(i1,i2,i3,1,vc))+fabs(udmzzC(i1,i2,i3,1,wc))+ \
fabs(udmzzC(i1,i2,i3,2,uc))+fabs(udmzzC(i1,i2,i3,2,vc))+fabs(udmzzC(i1,i2,i3,2,wc)) )
#define admzzCSA (   i1,
  i2,
  i3 
)    ad21+cd22*( fabs(udmzzC(i1,i2,i3,0,nc))+fabs(udmzzC(i1,i2,i3,1,nc))+fabs(udmzzC(i1,i2,i3,2,nc)) )
#define admzzR (   i1,
  i2,
  i3 
)
Value:
ad21+cd22*( fabs(uxmzzR(i1,i2,i3,uc))+fabs(uxmzzR(i1,i2,i3,vc))+fabs(uxmzzR(i1,i2,i3,wc))+\
fabs(uymzzR(i1,i2,i3,uc))+fabs(uymzzR(i1,i2,i3,vc))+fabs(uymzzR(i1,i2,i3,wc))+\
fabs(uzmzzR(i1,i2,i3,uc))+fabs(uzmzzR(i1,i2,i3,vc))+fabs(uzmzzR(i1,i2,i3,wc)) )
#define admzzRSA (   i1,
  i2,
  i3 
)    ad21+cd22*( fabs(uxmzzR(i1,i2,i3,nc))+fabs(uymzzR(i1,i2,i3,nc))+fabs(uzmzzR(i1,i2,i3,nc)) )
#define ADTZ2 (   cc,
  q0 
)
Value:
adc2*(e(mg,I1-1,I2,I3,cc,t)+e(mg,I1+1,I2,I3,cc,t)+\
e(mg,I1,I2-1,I3,cc,t)+e(mg,I1,I2+1,I3,cc,t)+\
-4.*q0 )
#define ADTZ2 (   cc,
  q0 
)
Value:
adc2*(e(mg,I1-1,I2,I3,cc,t)+e(mg,I1+1,I2,I3,cc,t)+\
e(mg,I1,I2-1,I3,cc,t)+e(mg,I1,I2+1,I3,cc,t)+\
e(mg,I1,I2,I3-1,cc,t)+e(mg,I1,I2,I3+1,cc,t)\
-6.*q0 )
#define ADTZ4 (   cc,
  q0 
)
Value:
adc4*(-(e(mg,I1-2,I2,I3,cc,t)+e(mg,I1+2,I2,I3,cc,t)+\
e(mg,I1,I2-2,I3,cc,t)+e(mg,I1,I2+2,I3,cc,t))\
+4.*(e(mg,I1-1,I2,I3,cc,t)+e(mg,I1+1,I2,I3,cc,t)+\
e(mg,I1,I2-1,I3,cc,t)+e(mg,I1,I2+1,I3,cc,t))\
-12.*q0 )
#define ADTZ4 (   cc,
  q0 
)
Value:
adc4*(-(e(mg,I1-2,I2,I3,cc,t)+e(mg,I1+2,I2,I3,cc,t)+\
e(mg,I1,I2-2,I3,cc,t)+e(mg,I1,I2+2,I3,cc,t)+\
e(mg,I1,I2,I3-2,cc,t)+e(mg,I1,I2,I3+2,cc,t))\
+4.*(e(mg,I1-1,I2,I3,cc,t)+e(mg,I1+1,I2,I3,cc,t)+\
e(mg,I1,I2-1,I3,cc,t)+e(mg,I1,I2+1,I3,cc,t)+\
e(mg,I1,I2,I3-1,cc,t)+e(mg,I1,I2,I3+1,cc,t))\
-18.*q0 )
#define adzmC (   i1,
  i2,
  i3 
)
Value:
ad21+cd22*( fabs(udzmC(i1,i2,i3,0,uc))+fabs(udzmC(i1,i2,i3,0,vc))+ \
fabs(udzmC(i1,i2,i3,1,uc))+fabs(udzmC(i1,i2,i3,1,vc)) )
#define adzmCSA (   i1,
  i2,
  i3 
)    ad21+cd22*( fabs(udzmC(i1,i2,i3,0,nc))+fabs(udzmC(i1,i2,i3,1,nc)) )
#define adzmR (   i1,
  i2,
  i3 
)
Value:
ad21+cd22*( fabs(uxzmzR(i1,i2,i3,uc))+fabs(uxzmzR(i1,i2,i3,vc))+\
fabs(uyzmzR(i1,i2,i3,uc))+fabs(uyzmzR(i1,i2,i3,vc)) )
#define adzmRSA (   i1,
  i2,
  i3 
)    ad21+cd22*( fabs(uxzmzR(i1,i2,i3,nc))+fabs(uyzmzR(i1,i2,i3,nc)) )
#define adzmzC (   i1,
  i2,
  i3 
)
Value:
ad21+cd22*( fabs(udzmzC(i1,i2,i3,0,uc))+fabs(udzmzC(i1,i2,i3,0,vc))+fabs(udzmzC(i1,i2,i3,0,wc))+\
fabs(udzmzC(i1,i2,i3,1,uc))+fabs(udzmzC(i1,i2,i3,1,vc))+fabs(udzmzC(i1,i2,i3,1,wc))+ \
fabs(udzmzC(i1,i2,i3,2,uc))+fabs(udzmzC(i1,i2,i3,2,vc))+fabs(udzmzC(i1,i2,i3,2,wc)) )
#define adzmzCSA (   i1,
  i2,
  i3 
)    ad21+cd22*( fabs(udzmzC(i1,i2,i3,0,nc))+fabs(udzmzC(i1,i2,i3,1,nc))+fabs(udzmzC(i1,i2,i3,2,nc)) )
#define adzmzR (   i1,
  i2,
  i3 
)
Value:
ad21+cd22*( fabs(uxzmzR(i1,i2,i3,uc))+fabs(uxzmzR(i1,i2,i3,vc))+fabs(uxzmzR(i1,i2,i3,wc))+\
fabs(uyzmzR(i1,i2,i3,uc))+fabs(uyzmzR(i1,i2,i3,vc))+fabs(uyzmzR(i1,i2,i3,wc))+\
fabs(uzzmzR(i1,i2,i3,uc))+fabs(uzzmzR(i1,i2,i3,vc))+fabs(uzzmzR(i1,i2,i3,wc)) )
#define adzmzRSA (   i1,
  i2,
  i3 
)    ad21+cd22*( fabs(uxzmzR(i1,i2,i3,nc))+fabs(uyzmzR(i1,i2,i3,nc))+fabs(uzzmzR(i1,i2,i3,nc)) )
#define adzzmC (   i1,
  i2,
  i3 
)
Value:
ad21+cd22*( fabs(udzzmC(i1,i2,i3,0,uc))+fabs(udzzmC(i1,i2,i3,0,vc))+fabs(udzzmC(i1,i2,i3,0,wc))+\
fabs(udzzmC(i1,i2,i3,1,uc))+fabs(udzzmC(i1,i2,i3,1,vc))+fabs(udzzmC(i1,i2,i3,1,wc))+ \
fabs(udzzmC(i1,i2,i3,2,uc))+fabs(udzzmC(i1,i2,i3,2,vc))+fabs(udzzmC(i1,i2,i3,2,wc)) )
#define adzzmCSA (   i1,
  i2,
  i3 
)    ad21+cd22*( fabs(udzzmC(i1,i2,i3,0,nc))+fabs(udzzmC(i1,i2,i3,1,nc))+fabs(udzzmC(i1,i2,i3,2,nc)) )
#define adzzmR (   i1,
  i2,
  i3 
)
Value:
ad21+cd22*( fabs(uxzzmR(i1,i2,i3,uc))+fabs(uxzzmR(i1,i2,i3,vc))+fabs(uxzzmR(i1,i2,i3,wc))+\
fabs(uyzzmR(i1,i2,i3,uc))+fabs(uyzzmR(i1,i2,i3,vc))+fabs(uyzzmR(i1,i2,i3,wc))+\
fabs(uzzzmR(i1,i2,i3,uc))+fabs(uzzzmR(i1,i2,i3,vc))+fabs(uzzzmR(i1,i2,i3,wc)) )
#define adzzmRSA (   i1,
  i2,
  i3 
)    ad21+cd22*( fabs(uxzzmR(i1,i2,i3,nc))+fabs(uyzzmR(i1,i2,i3,nc))+fabs(uzzzmR(i1,i2,i3,nc)) )
#define CE (   c,
  e 
)    (stencilSize*((c)+numberOfComponentsForCoefficients*(e)))

Referenced by Maxwell::assignInitialConditions().

#define ForBoundary (   side,
  axis 
)
Value:
for( axis=0; axis<mg.numberOfDimensions(); axis++ ) \
for( side=0; side<=1; side++ )
#define insdt   insdt_
#define insdts   insdts_
#define M123 (   m1,
  m2,
  m3 
)    (m1+halfWidth1+width*(m2+halfWidth2+width*(m3+halfWidth3)))
#define M123CE (   m1,
  m2,
  m3,
  c,
  e 
)    (M123(m1,m2,m3)+CE(c,e))
#define M123N (   m1,
  m2,
  m3,
  n 
)    (M123(m1,m2,m3)+stencilSize*(n))
#define MCE (   m1,
  m2,
  m3 
)    (M123(m1,m2,m3)+CE(c0,e0))
#define RX (   i1,
  i2,
  i3,
  m,
  n 
)    inverseVertexDerivative(i1,i2,i3,m+numberOfDimensions*(n))
#define U (   c)    u(I1,I2,I3,c)
#define udmzC (   i1,
  i2,
  i3,
  m,
  c 
)
Value:
(RX(i1,i2,i3,0,m)+RX(i1-1,i2,i3,0,m))*(u(i1,i2,i3,c)-u(i1-1,i2 ,i3,c))*dr2i +\
(RX(i1,i2,i3,1,m)+RX(i1-1,i2,i3,1,m))*(\
u(i1,i2+1,i3,c)-u(i1,i2-1,i3,c)+ u(i1-1,i2+1,i3,c)-u(i1-1,i2-1,i3,c))*dsi*.125
#define udmzzC (   i1,
  i2,
  i3,
  m,
  c 
)
Value:
(RX(i1,i2,i3,0,m)+RX(i1-1,i2,i3,0,m))*(u(i1,i2,i3,c)-u(i1-1,i2 ,i3,c))*dr2i +\
(RX(i1,i2,i3,1,m)+RX(i1-1,i2,i3,1,m))*(\
u(i1,i2+1,i3,c)-u(i1,i2-1,i3,c)+ u(i1-1,i2+1,i3,c)-u(i1-1,i2-1,i3,c))*dsi*.125+\
(RX(i1,i2,i3,2,m)+RX(i1-1,i2,i3,2,m))*(\
u(i1,i2,i3+1,c)-u(i1,i2,i3-1,c)+ u(i1-1,i2,i3+1,c)-u(i1-1,i2,i3-1,c))*dti*.125
#define udzmC (   i1,
  i2,
  i3,
  m,
  c 
)
Value:
(RX(i1,i2,i3,1,m)+RX(i1,i2-1,i3,1,m))*(u(i1,i2,i3,c)-u(i1,i2-1,i3,c))*ds2i +\
(RX(i1,i2,i3,0,m)+RX(i1,i2-1,i3,0,m))*(\
u(i1+1,i2,i3,c)-u(i1-1,i2,i3,c)+ u(i1+1,i2-1,i3,c)-u(i1-1,i2-1,i3,c))*dri*.125
#define udzmzC (   i1,
  i2,
  i3,
  m,
  c 
)
Value:
(RX(i1,i2,i3,1,m)+RX(i1,i2-1,i3,1,m))*(u(i1,i2,i3,c)-u(i1,i2-1,i3,c))*ds2i +\
(RX(i1,i2,i3,0,m)+RX(i1,i2-1,i3,0,m))*(\
u(i1+1,i2,i3,c)-u(i1-1,i2,i3,c)+ u(i1+1,i2-1,i3,c)-u(i1-1,i2-1,i3,c))*dri*.125+\
(RX(i1,i2,i3,2,m)+RX(i1,i2-1,i3,2,m))*(\
u(i1,i2,i3+1,c)-u(i1,i2,i3-1,c)+ u(i1,i2-1,i3+1,c)-u(i1,i2-1,i3-1,c))*dti*.125
#define udzzmC (   i1,
  i2,
  i3,
  m,
  c 
)
Value:
(RX(i1,i2,i3,2,m)+RX(i1,i2,i3-1,2,m))*(u(i1,i2,i3,c)-u(i1,i2,i3-1,c))*dt2i +\
(RX(i1,i2,i3,0,m)+RX(i1,i2,i3-1,0,m))*(\
u(i1+1,i2,i3,c)-u(i1-1,i2,i3,c)+ u(i1+1,i2,i3-1,c)-u(i1-1,i2,i3-1,c))*dri*.125+\
(RX(i1,i2,i3,1,m)+RX(i1,i2,i3-1,1,m))*(\
u(i1,i2+1,i3,c)-u(i1,i2-1,i3,c)+ u(i1,i2+1,i3-1,c)-u(i1,i2-1,i3-1,c))*dsi*.125
#define UU (   c)    uu(I1,I2,I3,c)
#define uv (   m)    e (mg,I1,I2,I3,m,t)
#define uvt (   m)    e.t (mg,I1,I2,I3,m,t)
#define uvx (   m)    e.x (mg,I1,I2,I3,m,t)
#define uvxx (   m)    e.xx(mg,I1,I2,I3,m,t)
#define uvxy (   m)    e.xy(mg,I1,I2,I3,m,t)
#define uvxz (   m)    e.xz(mg,I1,I2,I3,m,t)
#define uvy (   m)    e.y (mg,I1,I2,I3,m,t)
#define uvyy (   m)    e.yy(mg,I1,I2,I3,m,t)
#define uvyz (   m)    e.yz(mg,I1,I2,I3,m,t)
#define uvz (   m)    e.z (mg,I1,I2,I3,m,t)
#define uvzz (   m)    e.zz(mg,I1,I2,I3,m,t)
#define UX (   c)    ux(I1,I2,I3,c)
#define uxmzzR (   i1,
  i2,
  i3,
  c 
)    (u(i1,i2,i3,c)-u(i1-1,i2,i3,c))*dxi
#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 uxzmzR (   i1,
  i2,
  i3,
  c 
)    (u(i1+1,i2,i3,c)-u(i1-1,i2,i3,c)+u(i1+1,i2-1,i3,c)-u(i1-1,i2-1,i3,c))*dxi*.25
#define uxzzmR (   i1,
  i2,
  i3,
  c 
)    (u(i1+1,i2,i3,c)-u(i1-1,i2,i3,c)+u(i1+1,i2,i3-1,c)-u(i1-1,i2,i3-1,c))*dxi*.25
#define UY (   c)    uy(I1,I2,I3,c)
#define uymzzR (   i1,
  i2,
  i3,
  c 
)    (u(i1,i2+1,i3,c)-u(i1,i2-1,i3,c)+u(i1-1,i2+1,i3,c)-u(i1-1,i2-1,i3,c))*dyi*.25
#define UYY (   c)    uyy(I1,I2,I3,c)
#define UYZ (   c)    uyz(I1,I2,I3,c)
#define uyzmzR (   i1,
  i2,
  i3,
  c 
)    (u(i1,i2,i3,c)-u(i1,i2-1,i3,c))*dyi
#define uyzzmR (   i1,
  i2,
  i3,
  c 
)    (u(i1,i2+1,i3,c)-u(i1,i2-1,i3,c)+u(i1,i2+1,i3-1,c)-u(i1,i2-1,i3-1,c))*dyi*.25
#define UZ (   c)    uz(I1,I2,I3,c)
#define uzmzzR (   i1,
  i2,
  i3,
  c 
)    (u(i1,i2,i3+1,c)-u(i1,i2,i3-1,c)+u(i1-1,i2,i3+1,c)-u(i1-1,i2,i3-1,c))*dzi*.25
#define UZZ (   c)    uzz(I1,I2,I3,c)
#define uzzmzR (   i1,
  i2,
  i3,
  c 
)    (u(i1,i2,i3+1,c)-u(i1,i2,i3-1,c)+u(i1,i2-1,i3+1,c)-u(i1,i2-1,i3-1,c))*dzi*.25
#define uzzzmR (   i1,
  i2,
  i3,
  c 
)    (u(i1,i2,i3,c)-u(i1,i2,i3-1,c))*dzi

Function Documentation

void insdt ( 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 &  gv,
const real &  dw,
const int &  bc,
const int &  ipar,
const real &  rpar,
const int &  ierr 
)
void insdts ( 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,
const real &  gv,
const real &  dw,
const real &  divDamp,
const real &  dtVar,
const int &  bc,
const int &  ipar,
const real &  rpar,
const int &  ierr 
)
int outputCompositeGrid ( CompositeGrid &  cg,
const aString &  gridFileName 
)

Referenced by Cgins::solveForTimeIndependentVariables().

Generated on Fri Jan 4 2013 10:19:32 for CG by   doxygen 1.8.3