CG  Version 25
Macros | Functions
insBC4.C File Reference
#include "Cgins.h"
#include "display.h"
#include "Parameters.h"
#include "ParallelUtility.h"
#include "ParallelGridUtility.h"
#include "InsParameters.h"
#include "Insbc4WorkSpace.h"
#include "App.h"
#include "cgux4a.h"
Include dependency graph for insBC4.C:

Macros

#define inscr   EXTERN_C_NAME(inscr)
 
#define insbc4   EXTERN_C_NAME(insbc4)
 
#define bcOptTemperature   EXTERN_C_NAME(bcopttemperature)
 
#define U   u
 
#define UR4   UR4A
 
#define US4   US4A
 
#define UT4   UT4A
 
#define fu1(m, ib)   pfu1[m+ndf*(ib)]
 
#define fu2(m, ib)   pfu2[m+ndf*(ib)]
 
#define ibe(m, ib)   pibe[m+3*(ib)]
 
#define u(i1, i2, i3, n)   pu[i1-ndra+ndr*(i2-ndsa+nds*(i3-ndta+ndt*(n)))]
 
#define u0(i1, i2, i3, n)   pu0[i1-ndra0+ndr0*(i2-ndsa0+nds0*(i3-ndta0+ndt0*(n)))]
 
#define cx(i, ord)   cxp[ord-2][i]
 
#define cxn(i, ord)   cxnp[ord-2][i]
 
#define XX(n1, n2, i1, i2, i3, kd, iord)   (cx(0,iord)*u(i1 ,i2 ,i3,kd)+cx(1,iord)*u(i1+ n1,i2+ n2,i3,kd)+cx(2,iord)*u(i1+2*n1,i2+2*n2,i3,kd)+cx(3,iord)*u(i1+3*n1,i2+3*n2,i3,kd)+cx(4,iord)*u(i1+4*n1,i2+4*n2,i3,kd)+cx(5,iord)*u(i1+5*n1,i2+5*n2,i3,kd)+cx(6,iord)*u(i1+6*n1,i2+6*n2,i3,kd) )
 
#define XXN(n1, n2, i1, i2, i3, kd, iord)   (cxn(0,iord)*u(i1 ,i2 ,i3,kd)+cxn(1,iord)*u(i1+ n1,i2+ n2,i3,kd)+cxn(2,iord)*u(i1+2*n1,i2+2*n2,i3,kd)+cxn(3,iord)*u(i1+3*n1,i2+3*n2,i3,kd)+cxn(4,iord)*u(i1+4*n1,i2+4*n2,i3,kd))
 
#define X2(n1, n2, i1, i2, i3, kd)   (-u(i1 ,i2 ,i3,kd)+ 2.*u(i1+ n1,i2+ n2,i3,kd)- u(i1+2*n1,i2+2*n2,i3,kd) )
 
#define X3(n1, n2, i1, i2, i3, kd)   (-u(i1 ,i2 ,i3,kd)+ 3.*u(i1+ n1,i2+ n2,i3,kd)- 3.*u(i1+2*n1,i2+2*n2,i3,kd)+ u(i1+3*n1,i2+3*n2,i3,kd))
 
#define X4(n1, n2, i1, i2, i3, kd)   (- (u(i1 ,i2 ,i3,kd)+u(i1+4*n1,i2+4*n2,i3,kd))+ 4.*(u(i1+ n1,i2+ n2,i3,kd)+u(i1+3*n1,i2+3*n2,i3,kd))- 6.* u(i1+2*n1,i2+2*n2,i3,kd) )
 
#define XT3(n1, n2, n3, i1, i2, i3, kd)   ( - u(i1 ,i2 ,i3 ,kd) + 3.*u(i1+ n1,i2+ n2,i3+ n3,kd) - 3.*u(i1+2*n1,i2+2*n2,i3+2*n3,kd) + u(i1+3*n1,i2+3*n2,i3+3*n3,kd))
 
#define XT6(n1, n2, n3, i1, i2, i3, kd)   ( - u(i1 ,i2 ,i3 ,kd) + 6.*u(i1+ n1,i2+ n2,i3+ n3,kd) -15.*u(i1+2*n1,i2+2*n2,i3+2*n3,kd) +20.*u(i1+3*n1,i2+3*n2,i3+3*n3,kd) -15.*u(i1+4*n1,i2+4*n2,i3+4*n3,kd) + 6.*u(i1+5*n1,i2+5*n2,i3+5*n3,kd) - u(i1+6*n1,i2+6*n2,i3+6*n3,kd))
 
#define XT6N(n1, n2, n3, i1, i2, i3, kd)   ( - u(i1 ,i2 ,i3 ,kd) + 6.*u(i1+ n1,i2+ n2,i3+ n3,kd) -15.*u(i1+2*n1,i2+2*n2,i3+2*n3,kd) +20.*u(i1+3*n1,i2+3*n2,i3+3*n3,kd) -15.*u(i1+4*n1,i2+4*n2,i3+4*n3,kd))
 
#define u0bt(i1, i2, i3, kd, t)   e.t(xyz(i1,i2,i3,0),xyz(i1,i2,i3,1),z0,kd,t)
 
#define u0bt3d(i1, i2, i3, kd, t)   e.t(xyz(i1,i2,i3,0),xyz(i1,i2,i3,1),xyz(i1,i2,i3,2),kd,t)
 
#define DU2(i1, i2, i3, kd)   ( nu*(LAPLACIAN42(i1,i2,i3,kd))- anl*(u0(i1,i2,i3,uc)*UX42(i1,i2,i3,kd) +u0(i1,i2,i3,vc)*UY42(i1,i2,i3,kd)) )
 
#define DU2F1(i1, i2, i3, kd)   ( UX42(i1,i2,i3,pc) +u0bt(i1,i2,i3,kd,t) )
 
#define DU2F2(i1, i2, i3, kd)   ( UY42(i1,i2,i3,pc) +u0bt(i1,i2,i3,kd,t) )
 
#define DU2A(i1, i2, i3, kd)   ( nu*LAPLACIAN42(i1,i2,i3,kd) )
 
#define DU2F1A(i1, i2, i3, kd)   ( UX42(i1,i2,i3,pc) +anl*(u0(i1,i2,i3,uc)*UX42(i1,i2,i3,kd) +u0(i1,i2,i3,vc)*UY42(i1,i2,i3,kd)) )
 
#define DU2F2A(i1, i2, i3, kd)   ( UY42(i1,i2,i3,pc) +anl*(u0(i1,i2,i3,uc)*UX42(i1,i2,i3,kd) +u0(i1,i2,i3,vc)*UY42(i1,i2,i3,kd)) )
 
#define X0(i1, i2, i3)   xyz(i1,i2,i3,0)
 
#define Y0(i1, i2, i3)   xyz(i1,i2,i3,1)
 
#define Z0(i1, i2, i3)   xyz(i1,i2,i3,2)
 
#define DU2B(i1, i2, i3, kd)   ( nu*LAPLACIAN42(i1,i2,i3,kd) )
 
#define DU2F1B(i1, i2, i3, kd)   ( UX42(i1,i2,i3,pc) +anl*(e(X0(i1,i2,i3),Y0(i1,i2,i3),z0,uc,t)*e.x(X0(i1,i2,i3),Y0(i1,i2,i3),z0,kd,t) +e(X0(i1,i2,i3),Y0(i1,i2,i3),z0,vc,t)*e.y(X0(i1,i2,i3),Y0(i1,i2,i3),z0,kd,t) ) +u0bt(i1,i2,i3,kd,t) )
 
#define DU2F2B(i1, i2, i3, kd)   ( UY42(i1,i2,i3,pc) +anl*(e(X0(i1,i2,i3),Y0(i1,i2,i3),z0,uc,t)*e.x(X0(i1,i2,i3),Y0(i1,i2,i3),z0,kd,t) +e(X0(i1,i2,i3),Y0(i1,i2,i3),z0,vc,t)*e.y(X0(i1,i2,i3),Y0(i1,i2,i3),z0,kd,t) ) +u0bt(i1,i2,i3,kd,t) )
 
#define FB1(x0, y0, z0, t0)   ( -nu*(e.xx(x0,y0,z0,uc,t0)+e.yy(x0,y0,z0,uc,t0)) +anl*(e(x0,y0,z0,uc,t0)*e.x(x0,y0,z0,uc,t0)+e(x0,y0,z0,vc,t0)*e.y(x0,y0,z0,uc,t0)) +e.x(x0,y0,z0,pc,t0)+e.t(x0,y0,z0,uc,t0) )
 
#define FB2(x0, y0, z0, t0)   ( -nu*(e.xx(x0,y0,z0,vc,t0)+e.yy(x0,y0,z0,vc,t0)) +anl*(e(x0,y0,z0,uc,t0)*e.x(x0,y0,z0,vc,t0)+e(x0,y0,z0,vc,t0)*e.y(x0,y0,z0,vc,t0)) +e.y(x0,y0,z0,pc,t0)+e.t(x0,y0,z0,vc,t0) )
 
#define DU3(i1, i2, i3, kd)   ( nu*LAPLACIAN43(i1,i2,i3,kd) )
 
#define DU3F1(i1, i2, i3, kd)   ( UX43(i1,i2,i3,pc)+anl*( u0(i1,i2,i3,uc)*UX43(i1,i2,i3,kd)+ u0(i1,i2,i3,vc)*UY43(i1,i2,i3,kd)+ u0(i1,i2,i3,wc)*UZ43(i1,i2,i3,kd) ) )
 
#define DU3F2(i1, i2, i3, kd)   ( UY43(i1,i2,i3,pc) +anl*( u0(i1,i2,i3,uc)*UX43(i1,i2,i3,kd)+ u0(i1,i2,i3,vc)*UY43(i1,i2,i3,kd)+ u0(i1,i2,i3,wc)*UZ43(i1,i2,i3,kd) ) )
 
#define DU3F3(i1, i2, i3, kd)   ( UZ43(i1,i2,i3,pc) +anl*( u0(i1,i2,i3,uc)*UX43(i1,i2,i3,kd)+ u0(i1,i2,i3,vc)*UY43(i1,i2,i3,kd)+ u0(i1,i2,i3,wc)*UZ43(i1,i2,i3,kd) ) )
 
#define DU3B(i1, i2, i3, kd)   ( nu*LAPLACIAN43(i1,i2,i3,kd) )
 
#define DU3F1B(i1, i2, i3, kd)   ( UX43(i1,i2,i3,pc)+anl*( e(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),uc,t)*e.x(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),kd,t)+ e(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),vc,t)*e.y(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),kd,t)+ e(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),wc,t)*e.z(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),kd,t) ) +u0bt3d(i1,i2,i3,kd,t) )
 
#define DU3F2B(i1, i2, i3, kd)   ( UY43(i1,i2,i3,pc) +anl*( e(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),uc,t)*e.x(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),kd,t)+ e(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),vc,t)*e.y(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),kd,t)+ e(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),wc,t)*e.z(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),kd,t) ) +u0bt3d(i1,i2,i3,kd,t) )
 
#define DU3F3B(i1, i2, i3, kd)   ( UZ43(i1,i2,i3,pc) +anl*( e(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),uc,t)*e.x(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),kd,t)+ e(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),vc,t)*e.y(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),kd,t)+ e(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),wc,t)*e.z(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),kd,t) ) +u0bt3d(i1,i2,i3,kd,t) )
 
#define F31(x0, y0, z0, t0)   ( e.t(x0,y0,z0,uc,t0) +anl*(e(x0,y0,z0,uc,t0)*e.x(x0,y0,z0,uc,t0) +e(x0,y0,z0,vc,t0)*e.y(x0,y0,z0,uc,t0) +e(x0,y0,z0,wc,t0)*e.z(x0,y0,z0,uc,t0)) + e.x(x0,y0,z0,pc,t0) - nu*( e.xx(x0,y0,z0,uc,t0)+e.yy(x0,y0,z0,uc,t0)+e.zz(x0,y0,z0,uc,t0) ) )
 
#define F32(x0, y0, z0, t0)   ( e.t(x0,y0,z0,vc,t0) +anl*(e(x0,y0,z0,uc,t0)*e.x(x0,y0,z0,vc,t0) +e(x0,y0,z0,vc,t0)*e.y(x0,y0,z0,vc,t0) +e(x0,y0,z0,wc,t0)*e.z(x0,y0,z0,vc,t0)) + e.y(x0,y0,z0,pc,t0) - nu*( e.xx(x0,y0,z0,vc,t0)+e.yy(x0,y0,z0,vc,t0)+e.zz(x0,y0,z0,vc,t0) ) )
 
#define F33(x0, y0, z0, t0)   ( e.t(x0,y0,z0,wc,t0) +anl*(e(x0,y0,z0,uc,t0)*e.x(x0,y0,z0,wc,t0) +e(x0,y0,z0,vc,t0)*e.y(x0,y0,z0,wc,t0) +e(x0,y0,z0,wc,t0)*e.z(x0,y0,z0,wc,t0)) + e.z(x0,y0,z0,pc,t0) - nu*( e.xx(x0,y0,z0,wc,t0)+e.yy(x0,y0,z0,wc,t0)+e.zz(x0,y0,z0,wc,t0) ) )
 
#define d12(m)   d12p[m]
 
#define d22(m)   d22p[m]
 
#define d14(m)   d14p[m]
 
#define d24(m)   d24p[m]
 
#define xv(m)   xvp[m]
 
#define df(m, ib)   pdf[m+ndf*(ib)]
 
#define DIV42(i1, i2, i3)   (UX42(i1,i2,i3,uc)+UY42(i1,i2,i3,vc))
 
#define DIV43(i1, i2, i3)   (UX43(i1,i2,i3,uc)+UY43(i1,i2,i3,vc)+UZ43(i1,i2,i3,wc))
 
#define DIV42_R(i1, i2, i3)   ( RX (i1,i2,i3)*URR4(i1,i2,i3,uc) +SX (i1,i2,i3)*URS4(i1,i2,i3,uc) +RXR4(i1,i2,i3)*UR4 (i1,i2,i3,uc) +SXR4(i1,i2,i3)*US4 (i1,i2,i3,uc) +RY (i1,i2,i3)*URR4(i1,i2,i3,vc) +SY (i1,i2,i3)*URS4(i1,i2,i3,vc) +RYR4(i1,i2,i3)*UR4 (i1,i2,i3,vc) +SYR4(i1,i2,i3)*US4 (i1,i2,i3,vc) )
 
#define DIV42_S(i1, i2, i3)   ( RX (i1,i2,i3)*URS4(i1,i2,i3,uc) +SX (i1,i2,i3)*USS4(i1,i2,i3,uc) +RXS4(i1,i2,i3)*UR4 (i1,i2,i3,uc) +SXS4(i1,i2,i3)*US4 (i1,i2,i3,uc) +RY (i1,i2,i3)*URS4(i1,i2,i3,vc) +SY (i1,i2,i3)*USS4(i1,i2,i3,vc) +RYS4(i1,i2,i3)*UR4 (i1,i2,i3,vc) +SYS4(i1,i2,i3)*US4 (i1,i2,i3,vc) )
 
#define DIV43_R(i1, i2, i3)   ( RX (i1,i2,i3)*URR4(i1,i2,i3,uc) +SX (i1,i2,i3)*URS4(i1,i2,i3,uc) +TX (i1,i2,i3)*URT4(i1,i2,i3,uc) +RXR4(i1,i2,i3)*UR4 (i1,i2,i3,uc) +SXR4(i1,i2,i3)*US4 (i1,i2,i3,uc) +TXR4(i1,i2,i3)*UT4 (i1,i2,i3,uc) +RY (i1,i2,i3)*URR4(i1,i2,i3,vc) +SY (i1,i2,i3)*URS4(i1,i2,i3,vc) +TY (i1,i2,i3)*URT4(i1,i2,i3,vc) +RYR4(i1,i2,i3)*UR4 (i1,i2,i3,vc) +SYR4(i1,i2,i3)*US4 (i1,i2,i3,vc) +TYR4(i1,i2,i3)*UT4 (i1,i2,i3,vc) +RZ (i1,i2,i3)*URR4(i1,i2,i3,wc) +SZ (i1,i2,i3)*URS4(i1,i2,i3,wc) +TZ (i1,i2,i3)*URT4(i1,i2,i3,wc) +RZR4(i1,i2,i3)*UR4 (i1,i2,i3,wc) +SZR4(i1,i2,i3)*US4 (i1,i2,i3,wc) +TZR4(i1,i2,i3)*UT4 (i1,i2,i3,wc) )
 
#define DIV43_S(i1, i2, i3)   ( RX (i1,i2,i3)*URS4(i1,i2,i3,uc) +SX (i1,i2,i3)*USS4(i1,i2,i3,uc) +TX (i1,i2,i3)*UST4(i1,i2,i3,uc) +RXS4(i1,i2,i3)*UR4 (i1,i2,i3,uc) +SXS4(i1,i2,i3)*US4 (i1,i2,i3,uc) +TXS4(i1,i2,i3)*UT4 (i1,i2,i3,uc) +RY (i1,i2,i3)*URS4(i1,i2,i3,vc) +SY (i1,i2,i3)*USS4(i1,i2,i3,vc) +TY (i1,i2,i3)*UST4(i1,i2,i3,vc) +RYS4(i1,i2,i3)*UR4 (i1,i2,i3,vc) +SYS4(i1,i2,i3)*US4 (i1,i2,i3,vc) +TYS4(i1,i2,i3)*UT4 (i1,i2,i3,vc) +RZ (i1,i2,i3)*URS4(i1,i2,i3,wc) +SZ (i1,i2,i3)*USS4(i1,i2,i3,wc) +TZ (i1,i2,i3)*UST4(i1,i2,i3,wc) +RZS4(i1,i2,i3)*UR4 (i1,i2,i3,wc) +SZS4(i1,i2,i3)*US4 (i1,i2,i3,wc) +TZS4(i1,i2,i3)*UT4 (i1,i2,i3,wc) )
 
#define DIV43_T(i1, i2, i3)   ( RX (i1,i2,i3)*URT4(i1,i2,i3,uc) +SX (i1,i2,i3)*UST4(i1,i2,i3,uc) +TX (i1,i2,i3)*UTT4(i1,i2,i3,uc) +RXT4(i1,i2,i3)*UR4 (i1,i2,i3,uc) +SXT4(i1,i2,i3)*US4 (i1,i2,i3,uc) +TXT4(i1,i2,i3)*UT4 (i1,i2,i3,uc) +RY (i1,i2,i3)*URT4(i1,i2,i3,vc) +SY (i1,i2,i3)*UST4(i1,i2,i3,vc) +TY (i1,i2,i3)*UTT4(i1,i2,i3,vc) +RYT4(i1,i2,i3)*UR4 (i1,i2,i3,vc) +SYT4(i1,i2,i3)*US4 (i1,i2,i3,vc) +TYT4(i1,i2,i3)*UT4 (i1,i2,i3,vc) +RZ (i1,i2,i3)*URT4(i1,i2,i3,wc) +SZ (i1,i2,i3)*UST4(i1,i2,i3,wc) +TZ (i1,i2,i3)*UTT4(i1,i2,i3,wc) +RZT4(i1,i2,i3)*UR4 (i1,i2,i3,wc) +SZT4(i1,i2,i3)*US4 (i1,i2,i3,wc) +TZT4(i1,i2,i3)*UT4 (i1,i2,i3,wc) )
 
#define u(i1, i2, i3, n)   pu[i1-ndra+ndr*(i2-ndsa+nds*(i3-ndta+ndt*(n)))]
 
#define d12(m)   d12p[m]
 
#define d22(m)   d22p[m]
 
#define d14(m)   d14p[m]
 
#define d24(m)   d24p[m]
 
#define d14(m)   d14p[m]
 
#define d24(m)   d24p[m]
 
#define dbc(s, a, side, axis)   (pdbc[(s)+2*((a)+3*((side)+2*(axis)))])
 
#define addBoundaryForcing(side, axis)   (pAddBoundaryForcing[(side)+2*(axis)])
 
#define bcfOffset(side, axis)   pbcfOffset[(side)+2*(axis)]
 
#define u(i1, i2, i3, n)   pu[i1-ndra+ndr*(i2-ndsa+nds*(i3-ndta+ndt*(n)))]
 
#define d14(m)   d14p[m]
 
#define d24(m)   d24p[m]
 
#define DIVBC(j1)   u0(j1+2*is1,i2,i3,uc)=8.*(u0(j1+is1,i2,i3,uc)-u0(j1-is1,i2-2*is2,i3,uc))+u0(j1-2*is1,i2,i3,uc) +is1*dx[0]/dx[1]*( 8.*(u0(j1,i2+1,i3,vc)-u0(j1,i2-1,i3,vc))+ u0(j1,i2-2,i3,vc)-u0(j1,i2+2,i3,vc));
 
#define mr(m, kd, ks)   mrp[m+3*(kd)+9*(ks)]
 
#define nrm(kd, ks, kd1, ks1)   nrmp[kd+3*(ks)+6*(kd1)+18*(ks1)]
 
#define nib(kd, ks, m, k)   pnib[kd+3*(ks)+6*(m)+12*(k)]
 
#define rx3(kd1, kd2)   (prx[i1-ndra+ndr*(i2-ndsa+nds*(i3-ndta+ndt*(kd1+nd*(kd2))))])
 
#define cf1(n1, n2, n3)   pcf1[n1+ndcf*(n2+ndcf2*(n3))]
 
#define cf2(n1, n2, n3)   pcf2[n1+ndcf*(n2+ndcf2*(n3))]
 
#define a0(n1, n2)   pa0[n1+2*(n2)]
 
#define u4(i1, i2, i3, m)   pu4[i1+2+5*(i2+2+5*(i3+2+5*(m)))]
 

Functions

void inscr (const int &kd1, const int &ks1, const int &kd2, const int &ks2, const int &nd, const int &gridIndexRange, const int &bc, const int &ndra, const int &ndrb, const int &ndsa, const int &ndsb, const int &ndta, const int &ndtb, real &u, const real &t, const real &d14, const real &d24, const real &ajs, const real &xy, const real &rsxy)
 
void insbc4 (const int &bcOption, const int &nd, 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 &ipar, const real &rpar, real &u, const int &mask, const real &x, const real &rx, const real &gridVelocity, const real &gtt, const int &bc, const int &indexRange, const int &ierr)
 
void bcOptTemperature (const int &nd, const int &nd1a, const int &nd1b, const int &nd2a, const int &nd2b, const int &nd3a, const int &nd3b, const int &gridIndexRange, const int &dimRange, const int &isPeriodic, real &u, const int &mask, const real &rsxy, const real &xy, const int &boundaryCondition, const int &addBoundaryForcing, const int &interfaceType, const int &ndb, const real &bcData, const int &dim, const real &bcf0, const int64_t &bcfOffset, const int &ipar, const real &rpar, int &ierr)
 
void insbe1 (const int iopt, int *pibe, const int nib, const int kd, const int ks, const int ndf, real *pfu1, real *pfu2, real t, realMappedGridFunction &ua, const realArray &inverseVertexDerivative, const int ndra, const int ndrb, const int ndsa, const int ndsb, const int ndta, const int ndtb, real *pu, const IntegerArray &ipar, const RealArray &rpar, OGFunction &e)
 
void insbe2 (const int iopt, const int ks, const int kd, const int ndf, real *pdf, const int nib, int *pibe, const int ndra, const int ndrb, const int ndsa, const int ndsb, const int ndta, const int ndtb, real *pu, realMappedGridFunction &ua, const realArray &inverseVertexDerivative, const IntegerArray &ipar, const RealArray &rpar)
 
void applyFourthOrderBoundaryConditionsOld (realMappedGridFunction &u0, real t, int grid, Insbc4WorkSpace &ws, IntegerArray &ipar, RealArray &rpar, OGFunction &exact, Parameters &parameters)
 

Macro Definition Documentation

#define a0 (   n1,
  n2 
)    pa0[n1+2*(n2)]

Referenced by applyFourthOrderBoundaryConditionsOld().

#define addBoundaryForcing (   side,
  axis 
)    (pAddBoundaryForcing[(side)+2*(axis)])
#define bcfOffset (   side,
  axis 
)    pbcfOffset[(side)+2*(axis)]
#define bcOptTemperature   EXTERN_C_NAME(bcopttemperature)

Referenced by Cgins::applyFourthOrderBoundaryConditions().

#define cf1 (   n1,
  n2,
  n3 
)    pcf1[n1+ndcf*(n2+ndcf2*(n3))]

Referenced by applyFourthOrderBoundaryConditionsOld().

#define cf2 (   n1,
  n2,
  n3 
)    pcf2[n1+ndcf*(n2+ndcf2*(n3))]

Referenced by applyFourthOrderBoundaryConditionsOld().

#define cx (   i,
  ord 
)    cxp[ord-2][i]

Referenced by Parameters::setTwilightZoneParameters().

#define cxn (   i,
  ord 
)    cxnp[ord-2][i]
#define d12 (   m)    d12p[m]
#define d12 (   m)    d12p[m]
#define d14 (   m)    d14p[m]

Referenced by Cgins::applyFourthOrderBoundaryConditions(), applyFourthOrderBoundaryConditionsOld(), insbe1(), and insbe2().

#define d14 (   m)    d14p[m]
#define d14 (   m)    d14p[m]
#define d14 (   m)    d14p[m]
#define d22 (   m)    d22p[m]
#define d22 (   m)    d22p[m]
#define d24 (   m)    d24p[m]

Referenced by Cgins::applyFourthOrderBoundaryConditions(), applyFourthOrderBoundaryConditionsOld(), insbe1(), and insbe2().

#define d24 (   m)    d24p[m]
#define d24 (   m)    d24p[m]
#define d24 (   m)    d24p[m]
#define dbc (   s,
  a,
  side,
  axis 
)    (pdbc[(s)+2*((a)+3*((side)+2*(axis)))])
#define df (   m,
  ib 
)    pdf[m+ndf*(ib)]

Referenced by RigidBodyMotion::getTimeStepEstimate(), and insbe2().

#define DIV42 (   i1,
  i2,
  i3 
)    (UX42(i1,i2,i3,uc)+UY42(i1,i2,i3,vc))

Referenced by insbe2().

#define DIV42_R (   i1,
  i2,
  i3 
)    ( RX (i1,i2,i3)*URR4(i1,i2,i3,uc) +SX (i1,i2,i3)*URS4(i1,i2,i3,uc) +RXR4(i1,i2,i3)*UR4 (i1,i2,i3,uc) +SXR4(i1,i2,i3)*US4 (i1,i2,i3,uc) +RY (i1,i2,i3)*URR4(i1,i2,i3,vc) +SY (i1,i2,i3)*URS4(i1,i2,i3,vc) +RYR4(i1,i2,i3)*UR4 (i1,i2,i3,vc) +SYR4(i1,i2,i3)*US4 (i1,i2,i3,vc) )

Referenced by insbe2().

#define DIV42_S (   i1,
  i2,
  i3 
)    ( RX (i1,i2,i3)*URS4(i1,i2,i3,uc) +SX (i1,i2,i3)*USS4(i1,i2,i3,uc) +RXS4(i1,i2,i3)*UR4 (i1,i2,i3,uc) +SXS4(i1,i2,i3)*US4 (i1,i2,i3,uc) +RY (i1,i2,i3)*URS4(i1,i2,i3,vc) +SY (i1,i2,i3)*USS4(i1,i2,i3,vc) +RYS4(i1,i2,i3)*UR4 (i1,i2,i3,vc) +SYS4(i1,i2,i3)*US4 (i1,i2,i3,vc) )

Referenced by insbe2().

#define DIV43 (   i1,
  i2,
  i3 
)    (UX43(i1,i2,i3,uc)+UY43(i1,i2,i3,vc)+UZ43(i1,i2,i3,wc))

Referenced by insbe2().

#define DIV43_R (   i1,
  i2,
  i3 
)    ( RX (i1,i2,i3)*URR4(i1,i2,i3,uc) +SX (i1,i2,i3)*URS4(i1,i2,i3,uc) +TX (i1,i2,i3)*URT4(i1,i2,i3,uc) +RXR4(i1,i2,i3)*UR4 (i1,i2,i3,uc) +SXR4(i1,i2,i3)*US4 (i1,i2,i3,uc) +TXR4(i1,i2,i3)*UT4 (i1,i2,i3,uc) +RY (i1,i2,i3)*URR4(i1,i2,i3,vc) +SY (i1,i2,i3)*URS4(i1,i2,i3,vc) +TY (i1,i2,i3)*URT4(i1,i2,i3,vc) +RYR4(i1,i2,i3)*UR4 (i1,i2,i3,vc) +SYR4(i1,i2,i3)*US4 (i1,i2,i3,vc) +TYR4(i1,i2,i3)*UT4 (i1,i2,i3,vc) +RZ (i1,i2,i3)*URR4(i1,i2,i3,wc) +SZ (i1,i2,i3)*URS4(i1,i2,i3,wc) +TZ (i1,i2,i3)*URT4(i1,i2,i3,wc) +RZR4(i1,i2,i3)*UR4 (i1,i2,i3,wc) +SZR4(i1,i2,i3)*US4 (i1,i2,i3,wc) +TZR4(i1,i2,i3)*UT4 (i1,i2,i3,wc) )

Referenced by insbe2().

#define DIV43_S (   i1,
  i2,
  i3 
)    ( RX (i1,i2,i3)*URS4(i1,i2,i3,uc) +SX (i1,i2,i3)*USS4(i1,i2,i3,uc) +TX (i1,i2,i3)*UST4(i1,i2,i3,uc) +RXS4(i1,i2,i3)*UR4 (i1,i2,i3,uc) +SXS4(i1,i2,i3)*US4 (i1,i2,i3,uc) +TXS4(i1,i2,i3)*UT4 (i1,i2,i3,uc) +RY (i1,i2,i3)*URS4(i1,i2,i3,vc) +SY (i1,i2,i3)*USS4(i1,i2,i3,vc) +TY (i1,i2,i3)*UST4(i1,i2,i3,vc) +RYS4(i1,i2,i3)*UR4 (i1,i2,i3,vc) +SYS4(i1,i2,i3)*US4 (i1,i2,i3,vc) +TYS4(i1,i2,i3)*UT4 (i1,i2,i3,vc) +RZ (i1,i2,i3)*URS4(i1,i2,i3,wc) +SZ (i1,i2,i3)*USS4(i1,i2,i3,wc) +TZ (i1,i2,i3)*UST4(i1,i2,i3,wc) +RZS4(i1,i2,i3)*UR4 (i1,i2,i3,wc) +SZS4(i1,i2,i3)*US4 (i1,i2,i3,wc) +TZS4(i1,i2,i3)*UT4 (i1,i2,i3,wc) )

Referenced by insbe2().

#define DIV43_T (   i1,
  i2,
  i3 
)    ( RX (i1,i2,i3)*URT4(i1,i2,i3,uc) +SX (i1,i2,i3)*UST4(i1,i2,i3,uc) +TX (i1,i2,i3)*UTT4(i1,i2,i3,uc) +RXT4(i1,i2,i3)*UR4 (i1,i2,i3,uc) +SXT4(i1,i2,i3)*US4 (i1,i2,i3,uc) +TXT4(i1,i2,i3)*UT4 (i1,i2,i3,uc) +RY (i1,i2,i3)*URT4(i1,i2,i3,vc) +SY (i1,i2,i3)*UST4(i1,i2,i3,vc) +TY (i1,i2,i3)*UTT4(i1,i2,i3,vc) +RYT4(i1,i2,i3)*UR4 (i1,i2,i3,vc) +SYT4(i1,i2,i3)*US4 (i1,i2,i3,vc) +TYT4(i1,i2,i3)*UT4 (i1,i2,i3,vc) +RZ (i1,i2,i3)*URT4(i1,i2,i3,wc) +SZ (i1,i2,i3)*UST4(i1,i2,i3,wc) +TZ (i1,i2,i3)*UTT4(i1,i2,i3,wc) +RZT4(i1,i2,i3)*UR4 (i1,i2,i3,wc) +SZT4(i1,i2,i3)*US4 (i1,i2,i3,wc) +TZT4(i1,i2,i3)*UT4 (i1,i2,i3,wc) )

Referenced by insbe2().

#define DIVBC (   j1)    u0(j1+2*is1,i2,i3,uc)=8.*(u0(j1+is1,i2,i3,uc)-u0(j1-is1,i2-2*is2,i3,uc))+u0(j1-2*is1,i2,i3,uc) +is1*dx[0]/dx[1]*( 8.*(u0(j1,i2+1,i3,vc)-u0(j1,i2-1,i3,vc))+ u0(j1,i2-2,i3,vc)-u0(j1,i2+2,i3,vc));

Referenced by applyFourthOrderBoundaryConditionsOld().

#define DU2 (   i1,
  i2,
  i3,
  kd 
)    ( nu*(LAPLACIAN42(i1,i2,i3,kd))- anl*(u0(i1,i2,i3,uc)*UX42(i1,i2,i3,kd) +u0(i1,i2,i3,vc)*UY42(i1,i2,i3,kd)) )
#define DU2A (   i1,
  i2,
  i3,
  kd 
)    ( nu*LAPLACIAN42(i1,i2,i3,kd) )

Referenced by insbe1().

#define DU2B (   i1,
  i2,
  i3,
  kd 
)    ( nu*LAPLACIAN42(i1,i2,i3,kd) )

Referenced by insbe1().

#define DU2F1 (   i1,
  i2,
  i3,
  kd 
)    ( UX42(i1,i2,i3,pc) +u0bt(i1,i2,i3,kd,t) )
#define DU2F1A (   i1,
  i2,
  i3,
  kd 
)    ( UX42(i1,i2,i3,pc) +anl*(u0(i1,i2,i3,uc)*UX42(i1,i2,i3,kd) +u0(i1,i2,i3,vc)*UY42(i1,i2,i3,kd)) )

Referenced by insbe1().

#define DU2F1B (   i1,
  i2,
  i3,
  kd 
)    ( UX42(i1,i2,i3,pc) +anl*(e(X0(i1,i2,i3),Y0(i1,i2,i3),z0,uc,t)*e.x(X0(i1,i2,i3),Y0(i1,i2,i3),z0,kd,t) +e(X0(i1,i2,i3),Y0(i1,i2,i3),z0,vc,t)*e.y(X0(i1,i2,i3),Y0(i1,i2,i3),z0,kd,t) ) +u0bt(i1,i2,i3,kd,t) )

Referenced by insbe1().

#define DU2F2 (   i1,
  i2,
  i3,
  kd 
)    ( UY42(i1,i2,i3,pc) +u0bt(i1,i2,i3,kd,t) )
#define DU2F2A (   i1,
  i2,
  i3,
  kd 
)    ( UY42(i1,i2,i3,pc) +anl*(u0(i1,i2,i3,uc)*UX42(i1,i2,i3,kd) +u0(i1,i2,i3,vc)*UY42(i1,i2,i3,kd)) )

Referenced by insbe1().

#define DU2F2B (   i1,
  i2,
  i3,
  kd 
)    ( UY42(i1,i2,i3,pc) +anl*(e(X0(i1,i2,i3),Y0(i1,i2,i3),z0,uc,t)*e.x(X0(i1,i2,i3),Y0(i1,i2,i3),z0,kd,t) +e(X0(i1,i2,i3),Y0(i1,i2,i3),z0,vc,t)*e.y(X0(i1,i2,i3),Y0(i1,i2,i3),z0,kd,t) ) +u0bt(i1,i2,i3,kd,t) )

Referenced by insbe1().

#define DU3 (   i1,
  i2,
  i3,
  kd 
)    ( nu*LAPLACIAN43(i1,i2,i3,kd) )

Referenced by insbe1().

#define DU3B (   i1,
  i2,
  i3,
  kd 
)    ( nu*LAPLACIAN43(i1,i2,i3,kd) )

Referenced by insbe1().

#define DU3F1 (   i1,
  i2,
  i3,
  kd 
)    ( UX43(i1,i2,i3,pc)+anl*( u0(i1,i2,i3,uc)*UX43(i1,i2,i3,kd)+ u0(i1,i2,i3,vc)*UY43(i1,i2,i3,kd)+ u0(i1,i2,i3,wc)*UZ43(i1,i2,i3,kd) ) )

Referenced by insbe1().

#define DU3F1B (   i1,
  i2,
  i3,
  kd 
)    ( UX43(i1,i2,i3,pc)+anl*( e(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),uc,t)*e.x(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),kd,t)+ e(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),vc,t)*e.y(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),kd,t)+ e(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),wc,t)*e.z(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),kd,t) ) +u0bt3d(i1,i2,i3,kd,t) )

Referenced by insbe1().

#define DU3F2 (   i1,
  i2,
  i3,
  kd 
)    ( UY43(i1,i2,i3,pc) +anl*( u0(i1,i2,i3,uc)*UX43(i1,i2,i3,kd)+ u0(i1,i2,i3,vc)*UY43(i1,i2,i3,kd)+ u0(i1,i2,i3,wc)*UZ43(i1,i2,i3,kd) ) )

Referenced by insbe1().

#define DU3F2B (   i1,
  i2,
  i3,
  kd 
)    ( UY43(i1,i2,i3,pc) +anl*( e(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),uc,t)*e.x(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),kd,t)+ e(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),vc,t)*e.y(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),kd,t)+ e(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),wc,t)*e.z(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),kd,t) ) +u0bt3d(i1,i2,i3,kd,t) )

Referenced by insbe1().

#define DU3F3 (   i1,
  i2,
  i3,
  kd 
)    ( UZ43(i1,i2,i3,pc) +anl*( u0(i1,i2,i3,uc)*UX43(i1,i2,i3,kd)+ u0(i1,i2,i3,vc)*UY43(i1,i2,i3,kd)+ u0(i1,i2,i3,wc)*UZ43(i1,i2,i3,kd) ) )

Referenced by insbe1().

#define DU3F3B (   i1,
  i2,
  i3,
  kd 
)    ( UZ43(i1,i2,i3,pc) +anl*( e(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),uc,t)*e.x(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),kd,t)+ e(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),vc,t)*e.y(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),kd,t)+ e(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),wc,t)*e.z(X0(i1,i2,i3),Y0(i1,i2,i3),Z0(i1,i2,i3),kd,t) ) +u0bt3d(i1,i2,i3,kd,t) )

Referenced by insbe1().

#define F31 (   x0,
  y0,
  z0,
  t0 
)    ( e.t(x0,y0,z0,uc,t0) +anl*(e(x0,y0,z0,uc,t0)*e.x(x0,y0,z0,uc,t0) +e(x0,y0,z0,vc,t0)*e.y(x0,y0,z0,uc,t0) +e(x0,y0,z0,wc,t0)*e.z(x0,y0,z0,uc,t0)) + e.x(x0,y0,z0,pc,t0) - nu*( e.xx(x0,y0,z0,uc,t0)+e.yy(x0,y0,z0,uc,t0)+e.zz(x0,y0,z0,uc,t0) ) )

Referenced by insbe1().

#define F32 (   x0,
  y0,
  z0,
  t0 
)    ( e.t(x0,y0,z0,vc,t0) +anl*(e(x0,y0,z0,uc,t0)*e.x(x0,y0,z0,vc,t0) +e(x0,y0,z0,vc,t0)*e.y(x0,y0,z0,vc,t0) +e(x0,y0,z0,wc,t0)*e.z(x0,y0,z0,vc,t0)) + e.y(x0,y0,z0,pc,t0) - nu*( e.xx(x0,y0,z0,vc,t0)+e.yy(x0,y0,z0,vc,t0)+e.zz(x0,y0,z0,vc,t0) ) )

Referenced by insbe1().

#define F33 (   x0,
  y0,
  z0,
  t0 
)    ( e.t(x0,y0,z0,wc,t0) +anl*(e(x0,y0,z0,uc,t0)*e.x(x0,y0,z0,wc,t0) +e(x0,y0,z0,vc,t0)*e.y(x0,y0,z0,wc,t0) +e(x0,y0,z0,wc,t0)*e.z(x0,y0,z0,wc,t0)) + e.z(x0,y0,z0,pc,t0) - nu*( e.xx(x0,y0,z0,wc,t0)+e.yy(x0,y0,z0,wc,t0)+e.zz(x0,y0,z0,wc,t0) ) )

Referenced by insbe1().

#define FB1 (   x0,
  y0,
  z0,
  t0 
)    ( -nu*(e.xx(x0,y0,z0,uc,t0)+e.yy(x0,y0,z0,uc,t0)) +anl*(e(x0,y0,z0,uc,t0)*e.x(x0,y0,z0,uc,t0)+e(x0,y0,z0,vc,t0)*e.y(x0,y0,z0,uc,t0)) +e.x(x0,y0,z0,pc,t0)+e.t(x0,y0,z0,uc,t0) )

Referenced by insbe1().

#define FB2 (   x0,
  y0,
  z0,
  t0 
)    ( -nu*(e.xx(x0,y0,z0,vc,t0)+e.yy(x0,y0,z0,vc,t0)) +anl*(e(x0,y0,z0,uc,t0)*e.x(x0,y0,z0,vc,t0)+e(x0,y0,z0,vc,t0)*e.y(x0,y0,z0,vc,t0)) +e.y(x0,y0,z0,pc,t0)+e.t(x0,y0,z0,vc,t0) )

Referenced by insbe1().

#define fu1 (   m,
  ib 
)    pfu1[m+ndf*(ib)]

Referenced by insbe1().

#define fu2 (   m,
  ib 
)    pfu2[m+ndf*(ib)]

Referenced by insbe1().

#define ibe (   m,
  ib 
)    pibe[m+3*(ib)]

Referenced by applyFourthOrderBoundaryConditionsOld(), insbe1(), and insbe2().

#define insbc4   EXTERN_C_NAME(insbc4)

Referenced by Cgins::applyFourthOrderBoundaryConditions().

#define inscr   EXTERN_C_NAME(inscr)

Referenced by applyFourthOrderBoundaryConditionsOld().

#define mr (   m,
  kd,
  ks 
)    mrp[m+3*(kd)+9*(ks)]

Referenced by applyFourthOrderBoundaryConditionsOld().

#define nib (   kd,
  ks,
  m,
  k 
)    pnib[kd+3*(ks)+6*(m)+12*(k)]

Referenced by applyFourthOrderBoundaryConditionsOld().

#define nrm (   kd,
  ks,
  kd1,
  ks1 
)    nrmp[kd+3*(ks)+6*(kd1)+18*(ks1)]

Referenced by PenaltyWallFunctionBC::applyBC(), and applyFourthOrderBoundaryConditionsOld().

#define rx3 (   kd1,
  kd2 
)    (prx[i1-ndra+ndr*(i2-ndsa+nds*(i3-ndta+ndt*(kd1+nd*(kd2))))])

Referenced by applyFourthOrderBoundaryConditionsOld().

#define U   u
#define u (   i1,
  i2,
  i3,
  n 
)    pu[i1-ndra+ndr*(i2-ndsa+nds*(i3-ndta+ndt*(n)))]

Referenced by DomainSolver::adaptGrids(), Cgsm::addDissipation(), Maxwell::addDissipation(), Maxwell::addFilter(), DomainSolver::addGrids(), DomainSolver::advance(), Cgsm::advance(), DomainSolver::advanceAdamsPredictorCorrector(), Cgad::advanceADI(), Maxwell::advanceFDTD(), Cgsm::advanceFOS(), DomainSolver::advanceImplicitMultiStep(), Cgins::advanceLineSolveNew(), Cgins::advanceLineSolveOld(), Cgsm::advanceMethodOfLines(), Maxwell::advanceNew(), Maxwell::advanceNFDTD(), DomainSolver::advanceSecondOrderSystem(), Cgsm::advanceSOS(), Maxwell::advanceSOSUP(), DomainSolver::advanceSteadyStateRungeKutta(), DomainSolver::advanceTrapezoidal(), DomainSolver::advanceVariableTimeStepAdamsPredictorCorrector(), Cgasf::allSpeedImplicitTimeStep(), PenaltyWallFunctionBC::applyBC(), Cgcns::applyBoundaryConditions(), Cgins::applyBoundaryConditions(), Cgsm::applyBoundaryConditions(), Cgcns::applyBoundaryConditionsForImplicitTimeStepping(), Cgad::applyBoundaryConditionsForImplicitTimeStepping(), Cgins::applyBoundaryConditionsForImplicitTimeStepping(), GridFunctionFilter::applyFilter(), DomainSolver::applyFilter(), applyFourthOrderBoundaryConditionsOld(), Cgsm::assignAnnulusEigenfunction(), RadiationBoundaryCondition::assignBoundaryConditions(), Maxwell::assignBoundaryConditions(), Cgsm::assignBoundaryConditionsFOS(), Cgsm::assignBoundaryConditionsSOS(), Cgsm::assignGaussianPulseInitialConditions(), Cgsm::assignHempInitialConditions(), Cgsm::assignInitialConditions(), DomainSolver::assignInitialConditions(), Cgmp::assignInterfaceRightHandSide(), Cgins::assignLineSolverBoundaryConditions(), Cgsm::assignParabolicInitialConditions(), Cgins::assignPressureRHS(), Cgsm::assignSpecialInitialConditions(), Cgcns::assignTestProblem(), Cgasf::assignTestProblem(), Cgsm::assignTwilightZoneInitialConditions(), DomainSolver::boundaryConditionPredictor(), DomainSolver::buildAmrGridsForInitialConditions(), DomainSolver::buildGrid(), Cgmp::buildRunTimeDialog(), Cgsm::buildRunTimeDialog(), Maxwell::buildRunTimeDialog(), DomainSolver::buildRunTimeDialog(), Cgsm::checkDisplacementAndStress(), checkGridFunction(), checkNormalBC(), DomainSolver::checkSolution(), DomainSolver::computeBodyForcing(), Cgsm::computeDissipation(), Maxwell::computeDissipation(), Maxwell::computeIntensity(), CnsParameters::conservativeToPrimitive(), convertToVertexCentered(), Cgins::determineErrors(), DomainSolver::determineErrors(), DomainSolver::endTimeStepFE(), DomainSolver::endTimeStepIM(), DomainSolver::endTimeStepPC(), DomainSolver::extrapolateInterpolationNeighbours(), for(), Cgcns::formImplicitTimeSteppingMatrix(), Cgcns::getAugmentedSolution(), Cgsm::getAugmentedSolution(), Maxwell::getAugmentedSolution(), Cgins::getAugmentedSolution(), DomainSolver::getAugmentedSolution(), Maxwell::getChargeDensity(), FlowSolutions::getDeformingDiffuser(), AsfParameters::getDerivedFunction(), CnsParameters::getDerivedFunction(), InsParameters::getDerivedFunction(), Cgsm::getEnergy(), Maxwell::getEnergy(), Cgsm::getErrors(), FlowSolutions::getFlowSolution(), Cgsm::getForcing(), Maxwell::getForcing(), getGridMotion(), Cgsm::getInitialConditions(), DomainSolver::getInitialConditions(), Cgmp::getInterfaceResiduals(), Cgins::getLineSolverBoundaryConditions(), Cgsm::getMaxDivAndCurl(), Maxwell::getMaxDivergence(), FlowSolutions::getObliqueShockFlow(), FlowSolutions::getPistonFlow(), FlowSolutions::getShockElasticPistonFlow(), DomainSolver::getSolutionBounds(), FlowSolutions::getSupersonicExpandingFlow(), Cgcns::getTimeSteppingEigenvalue(), Cgasf::getTimeSteppingEigenvalue(), InsParameters::getTurbulenceModelVariables(), InsParameters::getTwoPhaseFlowVariables(), Parameters::getUserDefinedKnownSolution(), Cgcns::getUt(), Cgasf::getUt(), Cgins::getUt(), Maxwell::getValuesFDTD(), Cgsm::getVelocityAndStress(), InsParameters::getViscoPlasticVariables(), MovingGrids::gridAccelerationBC(), Cgasf::gridAccelerationBC(), Cgins::gridAccelerationBC(), if(), DomainSolver::initializeSolution(), DomainSolver::initializeTimeSteppingIM(), DomainSolver::initializeTimeSteppingPC(), Cgins::initializeTurbulenceModels(), DeformingBodyMotion::integrate(), DomainSolver::interpolate(), Cgins::lineSolverBoundaryConditions(), main(), Cgasf::maxMachNumber(), ogfTaylor(), DomainSolver::output(), DomainSolver::outputProbes(), Maxwell::outputResultsAfterEachTimeStep(), Cgsm::outputResultsAfterEachTimeStep(), DomainSolver::parabolicInflow(), Cgmp::plot(), Maxwell::plot(), Cgsm::plot(), DomainSolver::plot(), CnsParameters::primitiveToConservative(), Cgins::printTimeStepInfo(), Maxwell::project(), Cgins::project(), Maxwell::projectInterpolationPoints(), GridFunction::put(), MovingGrids::rigidBodyMotion(), Maxwell::saveShow(), Cgsm::saveShow(), DomainSolver::saveShow(), DomainSolver::setParametersInteractively(), Cgmp::setup(), DomainSolver::setup(), Cgsm::setupGridFunctions(), DomainSolver::setupGridFunctions(), Cgad::setupPde(), Cgcns::setupPde(), Cgasf::setupPde(), Cgins::setupPde(), Cgsm::setupPde(), DomainSolver::setVariableBoundaryValues(), DomainSolver::smoothVelocity(), Cgsm::solve(), Cgasf::solveForAllSpeedPressure(), Cgasf::solveForTimeIndependentVariables(), Cgins::solveForTimeIndependentVariables(), Reactions::solveImplicitForPTYGivenR(), Reactions::solveImplicitForRTYGivenP(), Reactions::solveImplicitForYGivenRTP(), DomainSolver::startTimeStepAF(), DomainSolver::takeTimeStepAF(), DomainSolver::takeTimeStepIM(), DomainSolver::takeTimeStepPC(), Cgins::turbulenceModels(), MovingGrids::update(), Cgins::updatePressureEquation(), Cgcns::updateStateVariables(), DomainSolver::updateToMatchNewGrid(), DomainSolver::updateVariableTimeInterpolation(), Cgsm::userDefinedBoundaryValues(), DomainSolver::userDefinedBoundaryValues(), DeformingBodyMotion::userDefinedDeformingSurface(), Maxwell::userDefinedForcing(), Cgsm::userDefinedForcing(), DomainSolver::userDefinedForcing(), DomainSolver::userDefinedGrid(), Maxwell::userDefinedInitialConditions(), DomainSolver::userDefinedOutput(), and DomainSolver::variableTimeStepBoundaryInterpolation().

#define u (   i1,
  i2,
  i3,
  n 
)    pu[i1-ndra+ndr*(i2-ndsa+nds*(i3-ndta+ndt*(n)))]
#define u (   i1,
  i2,
  i3,
  n 
)    pu[i1-ndra+ndr*(i2-ndsa+nds*(i3-ndta+ndt*(n)))]
#define u0 (   i1,
  i2,
  i3,
  n 
)    pu0[i1-ndra0+ndr0*(i2-ndsa0+nds0*(i3-ndta0+ndt0*(n)))]

Referenced by Cgasf::addAllSpeedImplicitForcing(), Cgins::addForcing(), Cgasf::addForcing(), Cgcns::addForcing(), Cgins::addForcingToPressureEquation(), Cgasf::addForcingToPressureEquation(), DomainSolver::advanceAdamsPredictorCorrector(), DomainSolver::advanceImplicitMultiStep(), DomainSolver::advanceSecondOrderSystem(), DomainSolver::advanceSteadyStateRungeKutta(), Cgasf::allSpeedImplicitTimeStep(), Cgins::applyFourthOrderBoundaryConditions(), applyFourthOrderBoundaryConditionsOld(), Cgins::assignPressureRHS(), Cgcns::assignTestProblem(), Cgasf::assignTestProblem(), DomainSolver::boundaryConditionPredictor(), Parameters::chooseUserDefinedBoundaryValues(), Cgins::computeAxisymmetricDivergence(), Cgcns::getAugmentedSolution(), Cgins::getAugmentedSolution(), DomainSolver::getAugmentedSolution(), FlowSolutions::getPistonFlow(), Cgins::getResidual(), FlowSolutions::getSupersonicExpandingFlow(), DomainSolver::getTimeDependentBoundaryConditions(), Cgcns::getTimeSteppingEigenvalue(), Cgasf::getTimeSteppingEigenvalue(), Cgins::getTimeSteppingEigenvalue(), Parameters::getUserDefinedKnownSolution(), Cgcns::getUt(), Cgasf::getUt(), Cgasf::initializeSolution(), Cgins::insImplicitMatrix(), DeformingBodyMotion::integrate(), DomainSolver::parabolicInflow(), Cgasf::solveForAllSpeedPressure(), Cgasf::solveForTimeIndependentVariables(), DomainSolver::takeTimeStepIM(), DomainSolver::takeTimeStepPC(), Parameters::updateUserDefinedKnownSolution(), DomainSolver::userDefinedBoundaryValues(), DeformingBodyMotion::userDefinedDeformingSurface(), and DomainSolver::userDefinedInitialConditions().

#define u0bt (   i1,
  i2,
  i3,
  kd,
  t 
)    e.t(xyz(i1,i2,i3,0),xyz(i1,i2,i3,1),z0,kd,t)
#define u0bt3d (   i1,
  i2,
  i3,
  kd,
  t 
)    e.t(xyz(i1,i2,i3,0),xyz(i1,i2,i3,1),xyz(i1,i2,i3,2),kd,t)
#define u4 (   i1,
  i2,
  i3,
  m 
)    pu4[i1+2+5*(i2+2+5*(i3+2+5*(m)))]

Referenced by applyFourthOrderBoundaryConditionsOld().

#define UR4   UR4A
#define US4   US4A
#define UT4   UT4A
#define X0 (   i1,
  i2,
  i3 
)    xyz(i1,i2,i3,0)
#define X2 (   n1,
  n2,
  i1,
  i2,
  i3,
  kd 
)    (-u(i1 ,i2 ,i3,kd)+ 2.*u(i1+ n1,i2+ n2,i3,kd)- u(i1+2*n1,i2+2*n2,i3,kd) )
#define X3 (   n1,
  n2,
  i1,
  i2,
  i3,
  kd 
)    (-u(i1 ,i2 ,i3,kd)+ 3.*u(i1+ n1,i2+ n2,i3,kd)- 3.*u(i1+2*n1,i2+2*n2,i3,kd)+ u(i1+3*n1,i2+3*n2,i3,kd))

Referenced by insbe1(), and insbe2().

#define X4 (   n1,
  n2,
  i1,
  i2,
  i3,
  kd 
)    (- (u(i1 ,i2 ,i3,kd)+u(i1+4*n1,i2+4*n2,i3,kd))+ 4.*(u(i1+ n1,i2+ n2,i3,kd)+u(i1+3*n1,i2+3*n2,i3,kd))- 6.* u(i1+2*n1,i2+2*n2,i3,kd) )

Referenced by insbe2().

#define XT3 (   n1,
  n2,
  n3,
  i1,
  i2,
  i3,
  kd 
)    ( - u(i1 ,i2 ,i3 ,kd) + 3.*u(i1+ n1,i2+ n2,i3+ n3,kd) - 3.*u(i1+2*n1,i2+2*n2,i3+2*n3,kd) + u(i1+3*n1,i2+3*n2,i3+3*n3,kd))

Referenced by insbe1().

#define XT6 (   n1,
  n2,
  n3,
  i1,
  i2,
  i3,
  kd 
)    ( - u(i1 ,i2 ,i3 ,kd) + 6.*u(i1+ n1,i2+ n2,i3+ n3,kd) -15.*u(i1+2*n1,i2+2*n2,i3+2*n3,kd) +20.*u(i1+3*n1,i2+3*n2,i3+3*n3,kd) -15.*u(i1+4*n1,i2+4*n2,i3+4*n3,kd) + 6.*u(i1+5*n1,i2+5*n2,i3+5*n3,kd) - u(i1+6*n1,i2+6*n2,i3+6*n3,kd))

Referenced by insbe1().

#define XT6N (   n1,
  n2,
  n3,
  i1,
  i2,
  i3,
  kd 
)    ( - u(i1 ,i2 ,i3 ,kd) + 6.*u(i1+ n1,i2+ n2,i3+ n3,kd) -15.*u(i1+2*n1,i2+2*n2,i3+2*n3,kd) +20.*u(i1+3*n1,i2+3*n2,i3+3*n3,kd) -15.*u(i1+4*n1,i2+4*n2,i3+4*n3,kd))

Referenced by insbe1().

#define xv (   m)    xvp[m]

Referenced by boundaryDistance(), DomainSolver::computeBodyForcing(), Parameters::defineVariableBoundaryValues(), Cgins::determineErrors(), DomainSolver::determineErrors(), MovingGrids::getUserDefinedBoundaryAcceleration(), MovingGrids::getUserDefinedGridVelocity(), if(), insbe1(), main(), DomainSolver::outputProbes(), DomainSolver::setVariableBoundaryValues(), DomainSolver::setVariableMaterialProperties(), TimeFunction::update(), DomainSolver::userDefinedForcing(), MovingGrids::userDefinedGridAccelerationBC(), and MovingGrids::userDefinedTransformMotion().

#define XX (   n1,
  n2,
  i1,
  i2,
  i3,
  kd,
  iord 
)    (cx(0,iord)*u(i1 ,i2 ,i3,kd)+cx(1,iord)*u(i1+ n1,i2+ n2,i3,kd)+cx(2,iord)*u(i1+2*n1,i2+2*n2,i3,kd)+cx(3,iord)*u(i1+3*n1,i2+3*n2,i3,kd)+cx(4,iord)*u(i1+4*n1,i2+4*n2,i3,kd)+cx(5,iord)*u(i1+5*n1,i2+5*n2,i3,kd)+cx(6,iord)*u(i1+6*n1,i2+6*n2,i3,kd) )
#define XXN (   n1,
  n2,
  i1,
  i2,
  i3,
  kd,
  iord 
)    (cxn(0,iord)*u(i1 ,i2 ,i3,kd)+cxn(1,iord)*u(i1+ n1,i2+ n2,i3,kd)+cxn(2,iord)*u(i1+2*n1,i2+2*n2,i3,kd)+cxn(3,iord)*u(i1+3*n1,i2+3*n2,i3,kd)+cxn(4,iord)*u(i1+4*n1,i2+4*n2,i3,kd))

Referenced by insbe1().

#define Y0 (   i1,
  i2,
  i3 
)    xyz(i1,i2,i3,1)
#define Z0 (   i1,
  i2,
  i3 
)    xyz(i1,i2,i3,2)

Function Documentation

void applyFourthOrderBoundaryConditionsOld ( realMappedGridFunction &  u0,
real  t,
int  grid,
Insbc4WorkSpace &  ws,
IntegerArray &  ipar,
RealArray &  rpar,
OGFunction &  exact,
Parameters &  parameters 
)

References a0, assert(), axis, boundaryCondition(), c, cf1, cf2, d, d14, d24, Parameters::dbase, Parameters::dirichletBoundaryCondition, display(), DIVBC, Insbc4WorkSpace::first, getIndex(), grid, i, i1, I1, i2, I2, i3, I3, ibe, InsParameters::inflowWithVelocityGiven, insbe1(), insbe2(), inscr, is2, is3, m, m1, m2, m3, mask, mr, N(), nib, normal, Parameters::noSlipWall, nrm, InsParameters::outflow, Insbc4WorkSpace::pcn, Insbc4WorkSpace::pct, Insbc4WorkSpace::pibe, Insbc4WorkSpace::pnib, Insbc4WorkSpace::pnibt, Insbc4WorkSpace::ptg, rx, rx3, side, Parameters::slipWall, u, u0, u4, v, V, and x2.

Referenced by Cgins::applyFourthOrderBoundaryConditions().

void bcOptTemperature ( const int &  nd,
const int &  nd1a,
const int &  nd1b,
const int &  nd2a,
const int &  nd2b,
const int &  nd3a,
const int &  nd3b,
const int &  gridIndexRange,
const int &  dimRange,
const int &  isPeriodic,
real &  u,
const int &  mask,
const real &  rsxy,
const real &  xy,
const int &  boundaryCondition,
const int &  addBoundaryForcing,
const int &  interfaceType,
const int &  ndb,
const real &  bcData,
const int &  dim,
const real &  bcf0,
const int64_t &  bcfOffset,
const int &  ipar,
const real &  rpar,
int &  ierr 
)
void insbc4 ( const int &  bcOption,
const int &  nd,
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 &  ipar,
const real &  rpar,
real &  u,
const int &  mask,
const real &  x,
const real &  rx,
const real &  gridVelocity,
const real &  gtt,
const int &  bc,
const int &  indexRange,
const int &  ierr 
)
void insbe1 ( const int  iopt,
int *  pibe,
const int  nib,
const int  kd,
const int  ks,
const int  ndf,
real *  pfu1,
real *  pfu2,
real  t,
realMappedGridFunction &  ua,
const realArray &  inverseVertexDerivative,
const int  ndra,
const int  ndrb,
const int  ndsa,
const int  ndsb,
const int  ndta,
const int  ndtb,
real *  pu,
const IntegerArray &  ipar,
const RealArray &  rpar,
OGFunction &  e 
)

References assert(), axis, d12, d14, d22, d24, DU2A, DU2B, DU2F1A, DU2F1B, DU2F2A, DU2F2B, DU3, DU3B, DU3F1, DU3F1B, DU3F2, DU3F2B, DU3F3, DU3F3B, e, F31, F32, F33, FB1, FB2, fu1, fu2, i1, i2, i3, ibe, InsParameters::inflowWithVelocityGiven, mg, normal, InsParameters::outflow, pc, Parameters::slipWall, uc, vc, wc, X3, XT3, XT6, XT6N, xv, XX, and XXN.

Referenced by applyFourthOrderBoundaryConditionsOld().

void insbe2 ( const int  iopt,
const int  ks,
const int  kd,
const int  ndf,
real *  pdf,
const int  nib,
int *  pibe,
const int  ndra,
const int  ndrb,
const int  ndsa,
const int  ndsb,
const int  ndta,
const int  ndtb,
real *  pu,
realMappedGridFunction &  ua,
const realArray &  inverseVertexDerivative,
const IntegerArray &  ipar,
const RealArray &  rpar 
)

References axis, d12, d14, d22, d24, df, DIV42, DIV42_R, DIV42_S, DIV43, DIV43_R, DIV43_S, DIV43_T, i1, i2, i3, ibe, InsParameters::inflowWithVelocityGiven, is2, mg, normal, InsParameters::outflow, pc, uc, vc, wc, X3, and X4.

Referenced by applyFourthOrderBoundaryConditionsOld().

void inscr ( const int &  kd1,
const int &  ks1,
const int &  kd2,
const int &  ks2,
const int &  nd,
const int &  gridIndexRange,
const int &  bc,
const int &  ndra,
const int &  ndrb,
const int &  ndsa,
const int &  ndsb,
const int &  ndta,
const int &  ndtb,
real &  u,
const real &  t,
const real &  d14,
const real &  d24,
const real &  ajs,
const real &  xy,
const real &  rsxy 
)
Generated on Fri Jan 4 2013 10:19:32 for CG by   doxygen 1.8.3