Blender V2.61 - r43446
Defines

solver_relax.h File Reference

Go to the source code of this file.

Defines

#define CAUSE_PANIC   { this->mPanic=1; }
#define CSRC_C   RAC(ccel , dC )
#define CSRC_E   RAC(ccel + (-1) *(dTotalNum), dE )
#define CSRC_W   RAC(ccel + (+1) *(dTotalNum), dW )
#define CSRC_N   RAC(ccel + (-mLevel[lev].lOffsx) *(dTotalNum), dN )
#define CSRC_S   RAC(ccel + (+mLevel[lev].lOffsx) *(dTotalNum), dS )
#define CSRC_NE   RAC(ccel + (-mLevel[lev].lOffsx-1) *(dTotalNum), dNE)
#define CSRC_NW   RAC(ccel + (-mLevel[lev].lOffsx+1) *(dTotalNum), dNW)
#define CSRC_SE   RAC(ccel + (+mLevel[lev].lOffsx-1) *(dTotalNum), dSE)
#define CSRC_SW   RAC(ccel + (+mLevel[lev].lOffsx+1) *(dTotalNum), dSW)
#define CSRC_T   RAC(ccel + (-mLevel[lev].lOffsy) *(dTotalNum), dT )
#define CSRC_B   RAC(ccel + (+mLevel[lev].lOffsy) *(dTotalNum), dB )
#define CSRC_ET   RAC(ccel + (-mLevel[lev].lOffsy-1) *(dTotalNum), dET)
#define CSRC_EB   RAC(ccel + (+mLevel[lev].lOffsy-1) *(dTotalNum), dEB)
#define CSRC_WT   RAC(ccel + (-mLevel[lev].lOffsy+1) *(dTotalNum), dWT)
#define CSRC_WB   RAC(ccel + (+mLevel[lev].lOffsy+1) *(dTotalNum), dWB)
#define CSRC_NT   RAC(ccel + (-mLevel[lev].lOffsy-mLevel[lev].lOffsx) *(dTotalNum), dNT)
#define CSRC_NB   RAC(ccel + (+mLevel[lev].lOffsy-mLevel[lev].lOffsx) *(dTotalNum), dNB)
#define CSRC_ST   RAC(ccel + (-mLevel[lev].lOffsy+mLevel[lev].lOffsx) *(dTotalNum), dST)
#define CSRC_SB   RAC(ccel + (+mLevel[lev].lOffsy+mLevel[lev].lOffsx) *(dTotalNum), dSB)
#define XSRC_C(x)   RAC(ccel + (x) *dTotalNum, dC )
#define XSRC_E(x)   RAC(ccel + ((x)-1) *dTotalNum, dE )
#define XSRC_W(x)   RAC(ccel + ((x)+1) *dTotalNum, dW )
#define XSRC_N(x)   RAC(ccel + ((x)-mLevel[lev].lOffsx) *dTotalNum, dN )
#define XSRC_S(x)   RAC(ccel + ((x)+mLevel[lev].lOffsx) *dTotalNum, dS )
#define XSRC_NE(x)   RAC(ccel + ((x)-mLevel[lev].lOffsx-1) *dTotalNum, dNE)
#define XSRC_NW(x)   RAC(ccel + ((x)-mLevel[lev].lOffsx+1) *dTotalNum, dNW)
#define XSRC_SE(x)   RAC(ccel + ((x)+mLevel[lev].lOffsx-1) *dTotalNum, dSE)
#define XSRC_SW(x)   RAC(ccel + ((x)+mLevel[lev].lOffsx+1) *dTotalNum, dSW)
#define XSRC_T(x)   RAC(ccel + ((x)-mLevel[lev].lOffsy) *dTotalNum, dT )
#define XSRC_B(x)   RAC(ccel + ((x)+mLevel[lev].lOffsy) *dTotalNum, dB )
#define XSRC_ET(x)   RAC(ccel + ((x)-mLevel[lev].lOffsy-1) *dTotalNum, dET)
#define XSRC_EB(x)   RAC(ccel + ((x)+mLevel[lev].lOffsy-1) *dTotalNum, dEB)
#define XSRC_WT(x)   RAC(ccel + ((x)-mLevel[lev].lOffsy+1) *dTotalNum, dWT)
#define XSRC_WB(x)   RAC(ccel + ((x)+mLevel[lev].lOffsy+1) *dTotalNum, dWB)
#define XSRC_NT(x)   RAC(ccel + ((x)-mLevel[lev].lOffsy-mLevel[lev].lOffsx) *dTotalNum, dNT)
#define XSRC_NB(x)   RAC(ccel + ((x)+mLevel[lev].lOffsy-mLevel[lev].lOffsx) *dTotalNum, dNB)
#define XSRC_ST(x)   RAC(ccel + ((x)-mLevel[lev].lOffsy+mLevel[lev].lOffsx) *dTotalNum, dST)
#define XSRC_SB(x)   RAC(ccel + ((x)+mLevel[lev].lOffsy+mLevel[lev].lOffsx) *dTotalNum, dSB)
#define OMEGA(l)   mLevel[(l)].omega
#define EQC   ( DFL1*(rho - usqr))
#define EQN   ( DFL2*(rho + uy*(4.5*uy + 3.0) - usqr))
#define EQS   ( DFL2*(rho + uy*(4.5*uy - 3.0) - usqr))
#define EQE   ( DFL2*(rho + ux*(4.5*ux + 3.0) - usqr))
#define EQW   ( DFL2*(rho + ux*(4.5*ux - 3.0) - usqr))
#define EQT   ( DFL2*(rho + uz*(4.5*uz + 3.0) - usqr))
#define EQB   ( DFL2*(rho + uz*(4.5*uz - 3.0) - usqr))
#define EQNE   ( DFL3*(rho + (+ux+uy)*(4.5*(+ux+uy) + 3.0) - usqr))
#define EQNW   ( DFL3*(rho + (-ux+uy)*(4.5*(-ux+uy) + 3.0) - usqr))
#define EQSE   ( DFL3*(rho + (+ux-uy)*(4.5*(+ux-uy) + 3.0) - usqr))
#define EQSW   ( DFL3*(rho + (-ux-uy)*(4.5*(-ux-uy) + 3.0) - usqr))
#define EQNT   ( DFL3*(rho + (+uy+uz)*(4.5*(+uy+uz) + 3.0) - usqr))
#define EQNB   ( DFL3*(rho + (+uy-uz)*(4.5*(+uy-uz) + 3.0) - usqr))
#define EQST   ( DFL3*(rho + (-uy+uz)*(4.5*(-uy+uz) + 3.0) - usqr))
#define EQSB   ( DFL3*(rho + (-uy-uz)*(4.5*(-uy-uz) + 3.0) - usqr))
#define EQET   ( DFL3*(rho + (+ux+uz)*(4.5*(+ux+uz) + 3.0) - usqr))
#define EQEB   ( DFL3*(rho + (+ux-uz)*(4.5*(+ux-uz) + 3.0) - usqr))
#define EQWT   ( DFL3*(rho + (-ux+uz)*(4.5*(-ux+uz) + 3.0) - usqr))
#define EQWB   ( DFL3*(rho + (-ux-uz)*(4.5*(-ux-uz) + 3.0) - usqr))
#define MSRC_C   m[dC ]
#define MSRC_N   m[dN ]
#define MSRC_S   m[dS ]
#define MSRC_E   m[dE ]
#define MSRC_W   m[dW ]
#define MSRC_T   m[dT ]
#define MSRC_B   m[dB ]
#define MSRC_NE   m[dNE]
#define MSRC_NW   m[dNW]
#define MSRC_SE   m[dSE]
#define MSRC_SW   m[dSW]
#define MSRC_NT   m[dNT]
#define MSRC_NB   m[dNB]
#define MSRC_ST   m[dST]
#define MSRC_SB   m[dSB]
#define MSRC_ET   m[dET]
#define MSRC_EB   m[dEB]
#define MSRC_WT   m[dWT]
#define MSRC_WB   m[dWB]
#define CCEL_C   RAC(ccel, dC )
#define CCEL_N   RAC(ccel, dN )
#define CCEL_S   RAC(ccel, dS )
#define CCEL_E   RAC(ccel, dE )
#define CCEL_W   RAC(ccel, dW )
#define CCEL_T   RAC(ccel, dT )
#define CCEL_B   RAC(ccel, dB )
#define CCEL_NE   RAC(ccel, dNE)
#define CCEL_NW   RAC(ccel, dNW)
#define CCEL_SE   RAC(ccel, dSE)
#define CCEL_SW   RAC(ccel, dSW)
#define CCEL_NT   RAC(ccel, dNT)
#define CCEL_NB   RAC(ccel, dNB)
#define CCEL_ST   RAC(ccel, dST)
#define CCEL_SB   RAC(ccel, dSB)
#define CCEL_ET   RAC(ccel, dET)
#define CCEL_EB   RAC(ccel, dEB)
#define CCEL_WT   RAC(ccel, dWT)
#define CCEL_WB   RAC(ccel, dWB)
#define CCELG_C(f)   (RAC(ccel, dC )*mGaussw[(f)])
#define CCELG_N(f)   (RAC(ccel, dN )*mGaussw[(f)])
#define CCELG_S(f)   (RAC(ccel, dS )*mGaussw[(f)])
#define CCELG_E(f)   (RAC(ccel, dE )*mGaussw[(f)])
#define CCELG_W(f)   (RAC(ccel, dW )*mGaussw[(f)])
#define CCELG_T(f)   (RAC(ccel, dT )*mGaussw[(f)])
#define CCELG_B(f)   (RAC(ccel, dB )*mGaussw[(f)])
#define CCELG_NE(f)   (RAC(ccel, dNE)*mGaussw[(f)])
#define CCELG_NW(f)   (RAC(ccel, dNW)*mGaussw[(f)])
#define CCELG_SE(f)   (RAC(ccel, dSE)*mGaussw[(f)])
#define CCELG_SW(f)   (RAC(ccel, dSW)*mGaussw[(f)])
#define CCELG_NT(f)   (RAC(ccel, dNT)*mGaussw[(f)])
#define CCELG_NB(f)   (RAC(ccel, dNB)*mGaussw[(f)])
#define CCELG_ST(f)   (RAC(ccel, dST)*mGaussw[(f)])
#define CCELG_SB(f)   (RAC(ccel, dSB)*mGaussw[(f)])
#define CCELG_ET(f)   (RAC(ccel, dET)*mGaussw[(f)])
#define CCELG_EB(f)   (RAC(ccel, dEB)*mGaussw[(f)])
#define CCELG_WT(f)   (RAC(ccel, dWT)*mGaussw[(f)])
#define CCELG_WB(f)   (RAC(ccel, dWB)*mGaussw[(f)])
#define CSMOMEGA_STATS(dlev, domega)
#define SRCS(l)   mLevel[(l)].setCurr
#define TSET(l)   mLevel[(l)].setOther
#define LBMDS_ADDMOV(linv, l)
#define DEFAULT_STREAM_FREESLIP(l, invl, mnbf)
#define STREAMCHECK(id, ni, nj, nk, nl)
#define COLLCHECK
#define DEFAULT_STREAM
#define DEFAULT_COLLIDEG(grav)
#define OPTIMIZED_STREAMCOLLIDE
#define OPTIMIZED_STREAMCOLLIDE___UNUSED
#define USQRMAXCHECK(Cusqr, Cux, Cuy, Cuz,CmMaxVlen, CmMxvx, CmMxvy, CmMxvz)
#define WO1D1   ( 1.0/ 2.0)
#define WO1D2   ( 1.0/ 4.0)
#define WO1D3   ( 1.0/ 8.0)
#define WO2D1_1   (-1.0/16.0)
#define WO2D1_9   ( 9.0/16.0)
#define WO2D2_11   (WO2D1_1 * WO2D1_1)
#define WO2D2_19   (WO2D1_9 * WO2D1_1)
#define WO2D2_91   (WO2D1_9 * WO2D1_1)
#define WO2D2_99   (WO2D1_9 * WO2D1_9)
#define WO2D3_111   (WO2D1_1 * WO2D1_1 * WO2D1_1)
#define WO2D3_191   (WO2D1_9 * WO2D1_1 * WO2D1_1)
#define WO2D3_911   (WO2D1_9 * WO2D1_1 * WO2D1_1)
#define WO2D3_991   (WO2D1_9 * WO2D1_9 * WO2D1_1)
#define WO2D3_119   (WO2D1_1 * WO2D1_1 * WO2D1_9)
#define WO2D3_199   (WO2D1_9 * WO2D1_1 * WO2D1_9)
#define WO2D3_919   (WO2D1_9 * WO2D1_1 * WO2D1_9)
#define WO2D3_999   (WO2D1_9 * WO2D1_9 * WO2D1_9)
#define ADD_INT_FLAGCHECK(alev, ai, aj, ak, at, afac)
#define ADD_INT_DFSCHECK(alev, ai, aj, ak, at, afac, l)
#define INTSTRICTCHECK(x, y, z, caseId)
#define INTUNUTCHECK(ix, iy, iz)
#define INTDEBOUT
#define ADD_INT_DFS(alev, ai, aj, ak, at, afac)
#define IDF_WRITEBACK_
#define IDF_WRITEBACK
#define DEBUG_CALCPRINTCELL(str, df)

Detailed Description

Definition in file solver_relax.h.


Define Documentation

#define ADD_INT_DFS (   alev,
  ai,
  aj,
  ak,
  at,
  afac 
)
Value:
ADD_INT_FLAGCHECK(alev, ai,aj,ak, at, afac); \
                        FORDF0{ \
                            LbmFloat df = ( \
                                    QCELL((alev), (ai),(aj),(ak),mLevel[(alev)].setCurr , l)*(1.0-(at)) + \
                                    QCELL((alev), (ai),(aj),(ak),mLevel[(alev)].setOther, l)*(    (at)) \
                                    ) ; \
                            ADD_INT_DFSCHECK(alev, ai,aj,ak, at, afac, l); \
                            df *= (afac); \
                            rho += df;  \
                            ux  += (this->dfDvecX[l]*df);  \
                            uy  += (this->dfDvecY[l]*df);   \
                            uz  += (this->dfDvecZ[l]*df);   \
                            intDf[l] += df; \
                        }

Definition at line 876 of file solver_relax.h.

Referenced by LbmFsgrSolver::interpolateCellFromCoarse().

#define ADD_INT_DFSCHECK (   alev,
  ai,
  aj,
  ak,
  at,
  afac,
 
)

Definition at line 839 of file solver_relax.h.

#define ADD_INT_FLAGCHECK (   alev,
  ai,
  aj,
  ak,
  at,
  afac 
)

Definition at line 838 of file solver_relax.h.

#define CAUSE_PANIC   { this->mPanic=1; }
#define CCEL_B   RAC(ccel, dB )

Definition at line 120 of file solver_relax.h.

#define CCEL_C   RAC(ccel, dC )

Definition at line 114 of file solver_relax.h.

#define CCEL_E   RAC(ccel, dE )

Definition at line 117 of file solver_relax.h.

#define CCEL_EB   RAC(ccel, dEB)

Definition at line 130 of file solver_relax.h.

#define CCEL_ET   RAC(ccel, dET)

Definition at line 129 of file solver_relax.h.

#define CCEL_N   RAC(ccel, dN )

Definition at line 115 of file solver_relax.h.

#define CCEL_NB   RAC(ccel, dNB)

Definition at line 126 of file solver_relax.h.

#define CCEL_NE   RAC(ccel, dNE)

Definition at line 121 of file solver_relax.h.

#define CCEL_NT   RAC(ccel, dNT)

Definition at line 125 of file solver_relax.h.

#define CCEL_NW   RAC(ccel, dNW)

Definition at line 122 of file solver_relax.h.

#define CCEL_S   RAC(ccel, dS )

Definition at line 116 of file solver_relax.h.

#define CCEL_SB   RAC(ccel, dSB)

Definition at line 128 of file solver_relax.h.

#define CCEL_SE   RAC(ccel, dSE)

Definition at line 123 of file solver_relax.h.

#define CCEL_ST   RAC(ccel, dST)

Definition at line 127 of file solver_relax.h.

#define CCEL_SW   RAC(ccel, dSW)

Definition at line 124 of file solver_relax.h.

#define CCEL_T   RAC(ccel, dT )

Definition at line 119 of file solver_relax.h.

#define CCEL_W   RAC(ccel, dW )

Definition at line 118 of file solver_relax.h.

#define CCEL_WB   RAC(ccel, dWB)

Definition at line 132 of file solver_relax.h.

#define CCEL_WT   RAC(ccel, dWT)

Definition at line 131 of file solver_relax.h.

#define CCELG_B (   f)    (RAC(ccel, dB )*mGaussw[(f)])

Definition at line 140 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define CCELG_C (   f)    (RAC(ccel, dC )*mGaussw[(f)])

Definition at line 134 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define CCELG_E (   f)    (RAC(ccel, dE )*mGaussw[(f)])

Definition at line 137 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define CCELG_EB (   f)    (RAC(ccel, dEB)*mGaussw[(f)])

Definition at line 150 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define CCELG_ET (   f)    (RAC(ccel, dET)*mGaussw[(f)])

Definition at line 149 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define CCELG_N (   f)    (RAC(ccel, dN )*mGaussw[(f)])

Definition at line 135 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define CCELG_NB (   f)    (RAC(ccel, dNB)*mGaussw[(f)])

Definition at line 146 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define CCELG_NE (   f)    (RAC(ccel, dNE)*mGaussw[(f)])

Definition at line 141 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define CCELG_NT (   f)    (RAC(ccel, dNT)*mGaussw[(f)])

Definition at line 145 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define CCELG_NW (   f)    (RAC(ccel, dNW)*mGaussw[(f)])

Definition at line 142 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define CCELG_S (   f)    (RAC(ccel, dS )*mGaussw[(f)])

Definition at line 136 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define CCELG_SB (   f)    (RAC(ccel, dSB)*mGaussw[(f)])

Definition at line 148 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define CCELG_SE (   f)    (RAC(ccel, dSE)*mGaussw[(f)])

Definition at line 143 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define CCELG_ST (   f)    (RAC(ccel, dST)*mGaussw[(f)])

Definition at line 147 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define CCELG_SW (   f)    (RAC(ccel, dSW)*mGaussw[(f)])

Definition at line 144 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define CCELG_T (   f)    (RAC(ccel, dT )*mGaussw[(f)])

Definition at line 139 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define CCELG_W (   f)    (RAC(ccel, dW )*mGaussw[(f)])

Definition at line 138 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define CCELG_WB (   f)    (RAC(ccel, dWB)*mGaussw[(f)])

Definition at line 152 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define CCELG_WT (   f)    (RAC(ccel, dWT)*mGaussw[(f)])

Definition at line 151 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define COLLCHECK

Definition at line 306 of file solver_relax.h.

#define CSMOMEGA_STATS (   dlev,
  domega 
)

Definition at line 162 of file solver_relax.h.

#define CSRC_B   RAC(ccel + (+mLevel[lev].lOffsy) *(dTotalNum), dB )

Definition at line 36 of file solver_relax.h.

#define CSRC_C   RAC(ccel , dC )

Definition at line 26 of file solver_relax.h.

#define CSRC_E   RAC(ccel + (-1) *(dTotalNum), dE )

Definition at line 27 of file solver_relax.h.

#define CSRC_EB   RAC(ccel + (+mLevel[lev].lOffsy-1) *(dTotalNum), dEB)

Definition at line 38 of file solver_relax.h.

#define CSRC_ET   RAC(ccel + (-mLevel[lev].lOffsy-1) *(dTotalNum), dET)

Definition at line 37 of file solver_relax.h.

#define CSRC_N   RAC(ccel + (-mLevel[lev].lOffsx) *(dTotalNum), dN )

Definition at line 29 of file solver_relax.h.

#define CSRC_NB   RAC(ccel + (+mLevel[lev].lOffsy-mLevel[lev].lOffsx) *(dTotalNum), dNB)

Definition at line 42 of file solver_relax.h.

#define CSRC_NE   RAC(ccel + (-mLevel[lev].lOffsx-1) *(dTotalNum), dNE)

Definition at line 31 of file solver_relax.h.

#define CSRC_NT   RAC(ccel + (-mLevel[lev].lOffsy-mLevel[lev].lOffsx) *(dTotalNum), dNT)

Definition at line 41 of file solver_relax.h.

#define CSRC_NW   RAC(ccel + (-mLevel[lev].lOffsx+1) *(dTotalNum), dNW)

Definition at line 32 of file solver_relax.h.

#define CSRC_S   RAC(ccel + (+mLevel[lev].lOffsx) *(dTotalNum), dS )

Definition at line 30 of file solver_relax.h.

#define CSRC_SB   RAC(ccel + (+mLevel[lev].lOffsy+mLevel[lev].lOffsx) *(dTotalNum), dSB)

Definition at line 44 of file solver_relax.h.

#define CSRC_SE   RAC(ccel + (+mLevel[lev].lOffsx-1) *(dTotalNum), dSE)

Definition at line 33 of file solver_relax.h.

#define CSRC_ST   RAC(ccel + (-mLevel[lev].lOffsy+mLevel[lev].lOffsx) *(dTotalNum), dST)

Definition at line 43 of file solver_relax.h.

#define CSRC_SW   RAC(ccel + (+mLevel[lev].lOffsx+1) *(dTotalNum), dSW)

Definition at line 34 of file solver_relax.h.

#define CSRC_T   RAC(ccel + (-mLevel[lev].lOffsy) *(dTotalNum), dT )

Definition at line 35 of file solver_relax.h.

#define CSRC_W   RAC(ccel + (+1) *(dTotalNum), dW )

Definition at line 28 of file solver_relax.h.

#define CSRC_WB   RAC(ccel + (+mLevel[lev].lOffsy+1) *(dTotalNum), dWB)

Definition at line 40 of file solver_relax.h.

#define CSRC_WT   RAC(ccel + (-mLevel[lev].lOffsy+1) *(dTotalNum), dWT)

Definition at line 39 of file solver_relax.h.

#define DEBUG_CALCPRINTCELL (   str,
  df 
)
Value:
{\
        LbmFloat prho=df[0], pux=0., puy=0., puz=0.; \
        for(int dfl=1; dfl<this->cDfNum; dfl++) { \
            prho += df[dfl];  \
            pux  += (this->dfDvecX[dfl]*df[dfl]);  \
            puy  += (this->dfDvecY[dfl]*df[dfl]);  \
            puz  += (this->dfDvecZ[dfl]*df[dfl]);  \
        } \
        errMsg("DEBUG_CALCPRINTCELL",">"<<str<<" rho="<<prho<<" vel="<<ntlVec3Gfx(pux,puy,puz) ); \
    }

Definition at line 1103 of file solver_relax.h.

#define DEFAULT_COLLIDEG (   grav)
Value:
this->collideArrays(lev, i,j,k, m, rho,ux,uy,uz, OMEGA(lev), grav, mLevel[lev].lcsmago, &mDebugOmegaRet, &lcsmqo ); \
    CSMOMEGA_STATS(lev,mDebugOmegaRet); \
    FORDF0 { RAC(tcel,l) = m[l]; } \
    usqr = 1.5 * (ux*ux + uy*uy + uz*uz); \
    COLLCHECK; \

Definition at line 354 of file solver_relax.h.

Referenced by LbmFsgrSolver::mainLoop(), and LbmFsgrSolver::standingFluidPreinit().

#define DEFAULT_STREAM

Definition at line 310 of file solver_relax.h.

Referenced by LbmFsgrSolver::mainLoop(), and LbmFsgrSolver::standingFluidPreinit().

#define DEFAULT_STREAM_FREESLIP (   l,
  invl,
  mnbf 
)

Definition at line 209 of file solver_relax.h.

#define EQB   ( DFL2*(rho + uz*(4.5*uz - 3.0) - usqr))

Definition at line 76 of file solver_relax.h.

Referenced by LbmFsgrSolver::mainLoop().

#define EQC   ( DFL1*(rho - usqr))

Definition at line 70 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define EQE   ( DFL2*(rho + ux*(4.5*ux + 3.0) - usqr))

Definition at line 73 of file solver_relax.h.

Referenced by LbmFsgrSolver::mainLoop().

#define EQEB   ( DFL3*(rho + (+ux-uz)*(4.5*(+ux-uz) + 3.0) - usqr))

Definition at line 87 of file solver_relax.h.

Referenced by LbmFsgrSolver::mainLoop().

#define EQET   ( DFL3*(rho + (+ux+uz)*(4.5*(+ux+uz) + 3.0) - usqr))

Definition at line 86 of file solver_relax.h.

Referenced by LbmFsgrSolver::mainLoop().

#define EQN   ( DFL2*(rho + uy*(4.5*uy + 3.0) - usqr))

Definition at line 71 of file solver_relax.h.

Referenced by LbmFsgrSolver::mainLoop().

#define EQNB   ( DFL3*(rho + (+uy-uz)*(4.5*(+uy-uz) + 3.0) - usqr))

Definition at line 83 of file solver_relax.h.

Referenced by LbmFsgrSolver::mainLoop().

#define EQNE   ( DFL3*(rho + (+ux+uy)*(4.5*(+ux+uy) + 3.0) - usqr))

Definition at line 78 of file solver_relax.h.

Referenced by LbmFsgrSolver::mainLoop().

#define EQNT   ( DFL3*(rho + (+uy+uz)*(4.5*(+uy+uz) + 3.0) - usqr))

Definition at line 82 of file solver_relax.h.

Referenced by LbmFsgrSolver::mainLoop().

#define EQNW   ( DFL3*(rho + (-ux+uy)*(4.5*(-ux+uy) + 3.0) - usqr))

Definition at line 79 of file solver_relax.h.

Referenced by LbmFsgrSolver::mainLoop().

#define EQS   ( DFL2*(rho + uy*(4.5*uy - 3.0) - usqr))

Definition at line 72 of file solver_relax.h.

Referenced by LbmFsgrSolver::mainLoop().

#define EQSB   ( DFL3*(rho + (-uy-uz)*(4.5*(-uy-uz) + 3.0) - usqr))

Definition at line 85 of file solver_relax.h.

Referenced by LbmFsgrSolver::mainLoop().

#define EQSE   ( DFL3*(rho + (+ux-uy)*(4.5*(+ux-uy) + 3.0) - usqr))

Definition at line 80 of file solver_relax.h.

Referenced by LbmFsgrSolver::mainLoop().

#define EQST   ( DFL3*(rho + (-uy+uz)*(4.5*(-uy+uz) + 3.0) - usqr))

Definition at line 84 of file solver_relax.h.

Referenced by LbmFsgrSolver::mainLoop().

#define EQSW   ( DFL3*(rho + (-ux-uy)*(4.5*(-ux-uy) + 3.0) - usqr))

Definition at line 81 of file solver_relax.h.

Referenced by LbmFsgrSolver::mainLoop().

#define EQT   ( DFL2*(rho + uz*(4.5*uz + 3.0) - usqr))

Definition at line 75 of file solver_relax.h.

Referenced by LbmFsgrSolver::mainLoop().

#define EQW   ( DFL2*(rho + ux*(4.5*ux - 3.0) - usqr))

Definition at line 74 of file solver_relax.h.

Referenced by LbmFsgrSolver::mainLoop().

#define EQWB   ( DFL3*(rho + (-ux-uz)*(4.5*(-ux-uz) + 3.0) - usqr))

Definition at line 89 of file solver_relax.h.

Referenced by LbmFsgrSolver::mainLoop().

#define EQWT   ( DFL3*(rho + (-ux+uz)*(4.5*(-ux+uz) + 3.0) - usqr))

Definition at line 88 of file solver_relax.h.

Referenced by LbmFsgrSolver::mainLoop().

#define IDF_WRITEBACK
Value:
LbmFloat omegaDst, omegaSrc;\
        /* smago new */ \
        LbmFloat feq[LBM_DFNUM]; \
        LbmFloat dfScale = mDfScaleDown; \
        FORDF0{ \
            feq[l] = getCollideEq(l, rho,ux,uy,uz); \
        } \
        if(mLevel[lev  ].lcsmago>0.0) {\
            LbmFloat Qo = this->getLesNoneqTensorCoeff(intDf,feq); \
            omegaDst  = this->getLesOmega(mLevel[lev+0].omega,mLevel[lev+0].lcsmago,Qo); \
            omegaSrc = this->getLesOmega(mLevel[lev-1].omega,mLevel[lev-1].lcsmago,Qo); \
        } else {\
            omegaDst = mLevel[lev+0].omega; \
            omegaSrc = mLevel[lev-1].omega;\
        } \
         \
        dfScale   = (mLevel[lev+0].timestep/mLevel[lev-1].timestep)* (1.0/omegaDst-1.0)/ (1.0/omegaSrc-1.0);  \
        FORDF0{ \
            /*errMsg("SMAGO"," org"<<mDfScaleDown<<" n"<<dfScale<<" qc"<< QCELL(lev,i,j,k, dstSet, l)<<" idf"<<intDf[l]<<" eq"<<feq[l] ); */ \
            QCELL(lev,i,j,k, dstSet, l) = (feq[l]+ (intDf[l]-feq[l])*dfScale);\
        } \
        /* check that all values are ok */ \
        INTDEBOUT

Definition at line 899 of file solver_relax.h.

Referenced by LbmFsgrSolver::interpolateCellFromCoarse().

#define IDF_WRITEBACK_
Value:
FORDF0{ \
            LbmFloat eq = getCollideEq(l, rho,ux,uy,uz);\
            QCELL(lev,i,j,k, dstSet, l) = (eq+ (intDf[l]-eq)*mDfScaleDown);\
        } \
        /* check that all values are ok */ \
        INTDEBOUT

Definition at line 892 of file solver_relax.h.

#define INTDEBOUT

Definition at line 868 of file solver_relax.h.

#define INTSTRICTCHECK (   x,
  y,
  z,
  caseId 
)

Definition at line 840 of file solver_relax.h.

#define INTUNUTCHECK (   ix,
  iy,
  iz 
)

Definition at line 841 of file solver_relax.h.

#define LBMDS_ADDMOV (   linv,
 
)
Value:
\
     \
    if((nbflag[linv]&CFBndMoving)&&(!(nbflag[l]&CFBnd))){ \
     \
    m[l]+=QCELL_NBINV(lev, i, j, k, SRCS(lev), l,l); \
    } \

Definition at line 195 of file solver_relax.h.

#define MSRC_B   m[dB ]

Definition at line 99 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define MSRC_C   m[dC ]

Definition at line 93 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define MSRC_E   m[dE ]

Definition at line 96 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define MSRC_EB   m[dEB]

Definition at line 109 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define MSRC_ET   m[dET]

Definition at line 108 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define MSRC_N   m[dN ]

Definition at line 94 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define MSRC_NB   m[dNB]

Definition at line 105 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define MSRC_NE   m[dNE]

Definition at line 100 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define MSRC_NT   m[dNT]

Definition at line 104 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define MSRC_NW   m[dNW]

Definition at line 101 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define MSRC_S   m[dS ]

Definition at line 95 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define MSRC_SB   m[dSB]

Definition at line 107 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define MSRC_SE   m[dSE]

Definition at line 102 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define MSRC_ST   m[dST]

Definition at line 106 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define MSRC_SW   m[dSW]

Definition at line 103 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define MSRC_T   m[dT ]

Definition at line 98 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define MSRC_W   m[dW ]

Definition at line 97 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define MSRC_WB   m[dWB]

Definition at line 111 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define MSRC_WT   m[dWT]

Definition at line 110 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseRestrictCell().

#define OMEGA (   l)    mLevel[(l)].omega

Definition at line 68 of file solver_relax.h.

#define OPTIMIZED_STREAMCOLLIDE
Value:
m[0] = RAC(ccel,0); \
    FORDF1 { \
     \
    if(RFLAG_NBINV(lev, i,j,k,SRCS(lev),l)&CFBnd) { errMsg("???", "bnd-err-nobndfl"); CAUSE_PANIC; \
    } else { m[l] = QCELL_NBINV(lev, i, j, k, SRCS(lev), l, l); } \
    STREAMCHECK(8, i+this->dfVecX[this->dfInv[l]], j+this->dfVecY[this->dfInv[l]],k+this->dfVecZ[this->dfInv[l]], l); \
    } \
    rho=m[0]; \
    DEFAULT_COLLIDEG(mLevel[lev].gravity) \

Definition at line 363 of file solver_relax.h.

Referenced by LbmFsgrSolver::coarseAdvance(), LbmFsgrSolver::mainLoop(), and LbmFsgrSolver::standingFluidPreinit().

#define OPTIMIZED_STREAMCOLLIDE___UNUSED
Value:
\
    this->collideArrays(lev, i,j,k, m, rho,ux,uy,uz, OMEGA(lev), mLevel[lev].gravity, mLevel[lev].lcsmago , &mDebugOmegaRet, &lcsmqo   ); \
    CSMOMEGA_STATS(lev,mDebugOmegaRet); \
    FORDF0 { RAC(tcel,l) = m[l]; } \
    usqr = 1.5 * (ux*ux + uy*uy + uz*uz); \
    COLLCHECK; \

Definition at line 376 of file solver_relax.h.

#define SRCS (   l)    mLevel[(l)].setCurr
#define STREAMCHECK (   id,
  ni,
  nj,
  nk,
  nl 
)

Definition at line 305 of file solver_relax.h.

#define TSET (   l)    mLevel[(l)].setOther
#define USQRMAXCHECK (   Cusqr,
  Cux,
  Cuy,
  Cuz,
  CmMaxVlen,
  CmMxvx,
  CmMxvy,
  CmMxvz 
)
Value:
if(Cusqr>CmMaxVlen) { \
                CmMxvx = Cux; CmMxvy = Cuy; CmMxvz = Cuz; CmMaxVlen = Cusqr; \
            }

Definition at line 767 of file solver_relax.h.

Referenced by LbmFsgrSolver::initMovingObstacles().

#define WO1D1   ( 1.0/ 2.0)

Definition at line 780 of file solver_relax.h.

Referenced by LbmFsgrSolver::interpolateCellFromCoarse().

#define WO1D2   ( 1.0/ 4.0)

Definition at line 781 of file solver_relax.h.

Referenced by LbmFsgrSolver::interpolateCellFromCoarse().

#define WO1D3   ( 1.0/ 8.0)

Definition at line 782 of file solver_relax.h.

Referenced by LbmFsgrSolver::interpolateCellFromCoarse().

#define WO2D1_1   (-1.0/16.0)

Definition at line 784 of file solver_relax.h.

#define WO2D1_9   ( 9.0/16.0)

Definition at line 785 of file solver_relax.h.

#define WO2D2_11   (WO2D1_1 * WO2D1_1)

Definition at line 787 of file solver_relax.h.

#define WO2D2_19   (WO2D1_9 * WO2D1_1)

Definition at line 788 of file solver_relax.h.

#define WO2D2_91   (WO2D1_9 * WO2D1_1)

Definition at line 789 of file solver_relax.h.

#define WO2D2_99   (WO2D1_9 * WO2D1_9)

Definition at line 790 of file solver_relax.h.

#define WO2D3_111   (WO2D1_1 * WO2D1_1 * WO2D1_1)

Definition at line 792 of file solver_relax.h.

#define WO2D3_119   (WO2D1_1 * WO2D1_1 * WO2D1_9)

Definition at line 796 of file solver_relax.h.

#define WO2D3_191   (WO2D1_9 * WO2D1_1 * WO2D1_1)

Definition at line 793 of file solver_relax.h.

#define WO2D3_199   (WO2D1_9 * WO2D1_1 * WO2D1_9)

Definition at line 797 of file solver_relax.h.

#define WO2D3_911   (WO2D1_9 * WO2D1_1 * WO2D1_1)

Definition at line 794 of file solver_relax.h.

#define WO2D3_919   (WO2D1_9 * WO2D1_1 * WO2D1_9)

Definition at line 798 of file solver_relax.h.

#define WO2D3_991   (WO2D1_9 * WO2D1_9 * WO2D1_1)

Definition at line 795 of file solver_relax.h.

#define WO2D3_999   (WO2D1_9 * WO2D1_9 * WO2D1_9)

Definition at line 799 of file solver_relax.h.

#define XSRC_B (   x)    RAC(ccel + ((x)+mLevel[lev].lOffsy) *dTotalNum, dB )

Definition at line 56 of file solver_relax.h.

#define XSRC_C (   x)    RAC(ccel + (x) *dTotalNum, dC )

Definition at line 46 of file solver_relax.h.

#define XSRC_E (   x)    RAC(ccel + ((x)-1) *dTotalNum, dE )

Definition at line 47 of file solver_relax.h.

#define XSRC_EB (   x)    RAC(ccel + ((x)+mLevel[lev].lOffsy-1) *dTotalNum, dEB)

Definition at line 58 of file solver_relax.h.

#define XSRC_ET (   x)    RAC(ccel + ((x)-mLevel[lev].lOffsy-1) *dTotalNum, dET)

Definition at line 57 of file solver_relax.h.

#define XSRC_N (   x)    RAC(ccel + ((x)-mLevel[lev].lOffsx) *dTotalNum, dN )

Definition at line 49 of file solver_relax.h.

#define XSRC_NB (   x)    RAC(ccel + ((x)+mLevel[lev].lOffsy-mLevel[lev].lOffsx) *dTotalNum, dNB)

Definition at line 62 of file solver_relax.h.

#define XSRC_NE (   x)    RAC(ccel + ((x)-mLevel[lev].lOffsx-1) *dTotalNum, dNE)

Definition at line 51 of file solver_relax.h.

#define XSRC_NT (   x)    RAC(ccel + ((x)-mLevel[lev].lOffsy-mLevel[lev].lOffsx) *dTotalNum, dNT)

Definition at line 61 of file solver_relax.h.

#define XSRC_NW (   x)    RAC(ccel + ((x)-mLevel[lev].lOffsx+1) *dTotalNum, dNW)

Definition at line 52 of file solver_relax.h.

#define XSRC_S (   x)    RAC(ccel + ((x)+mLevel[lev].lOffsx) *dTotalNum, dS )

Definition at line 50 of file solver_relax.h.

#define XSRC_SB (   x)    RAC(ccel + ((x)+mLevel[lev].lOffsy+mLevel[lev].lOffsx) *dTotalNum, dSB)

Definition at line 64 of file solver_relax.h.

#define XSRC_SE (   x)    RAC(ccel + ((x)+mLevel[lev].lOffsx-1) *dTotalNum, dSE)

Definition at line 53 of file solver_relax.h.

#define XSRC_ST (   x)    RAC(ccel + ((x)-mLevel[lev].lOffsy+mLevel[lev].lOffsx) *dTotalNum, dST)

Definition at line 63 of file solver_relax.h.

#define XSRC_SW (   x)    RAC(ccel + ((x)+mLevel[lev].lOffsx+1) *dTotalNum, dSW)

Definition at line 54 of file solver_relax.h.

#define XSRC_T (   x)    RAC(ccel + ((x)-mLevel[lev].lOffsy) *dTotalNum, dT )

Definition at line 55 of file solver_relax.h.

#define XSRC_W (   x)    RAC(ccel + ((x)+1) *dTotalNum, dW )

Definition at line 48 of file solver_relax.h.

#define XSRC_WB (   x)    RAC(ccel + ((x)+mLevel[lev].lOffsy+1) *dTotalNum, dWB)

Definition at line 60 of file solver_relax.h.

#define XSRC_WT (   x)    RAC(ccel + ((x)-mLevel[lev].lOffsy+1) *dTotalNum, dWT)

Definition at line 59 of file solver_relax.h.