CG
Version 25
|
#include "Cgins.h"
#include "TridiagonalSolver.h"
#include "MappedGridOperators.h"
#include "LineSolve.h"
#include "ParallelUtility.h"
#include "ParallelGridUtility.h"
Macros | |
#define | tridiagonalSolver(c, axis, grid) lineSolve.pTridiagonalSolvers[c+maxNumberOfSystems*(axis+numberOfDimensions*(grid))] |
#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 | insLineSetup EXTERN_C_NAME(inslinesetup) |
#define | insLineSetupNew EXTERN_C_NAME(inslinesetupnew) |
#define | insLineSolveBC EXTERN_C_NAME(inslinesolvebc) |
#define | computeResidual EXTERN_C_NAME(computeresidual) |
#define | computeResidualNew EXTERN_C_NAME(computeresidualnew) |
#define | FOR_3D(i1, i2, i3, I1, I2, I3) int I1Base =I1.getBase(), I2Base =I2.getBase(), I3Base =I3.getBase(); int I1Bound=I1.getBound(), I2Bound=I2.getBound(), I3Bound=I3.getBound(); for(i3=I3Base; i3<=I3Bound; i3++) for(i2=I2Base; i2<=I2Bound; i2++) for(i1=I1Base; i1<=I1Bound; i1++) |
#define | FOR_3(i1, i2, i3, I1, I2, I3) I1Base =I1.getBase(), I2Base =I2.getBase(), I3Base =I3.getBase(); I1Bound=I1.getBound(), I2Bound=I2.getBound(), I3Bound=I3.getBound(); for(i3=I3Base; i3<=I3Bound; i3++) for(i2=I2Base; i2<=I2Bound; i2++) for(i1=I1Base; i1<=I1Bound; i1++) |
#define | MASK(i0, i1, i2) maskp[i0+maskDim0*(i1+maskDim1*(i2))] |
#define | F(i0, i1, i2, i3) fp[i0+fDim0*(i1+fDim1*(i2+fDim2*(i3)))] |
#define | U(i0, i1, i2, i3) up[i0+uDim0*(i1+uDim1*(i2+uDim2*(i3)))] |
#define | RSXY(I1, I2, I3, m1, m2) rsxy(I1,I2,I3,(m1)+numberOfDimensions*(m2)) |
#define | RSXY(I1, I2, I3, m1, m2) rsxy(I1,I2,I3,(m1)+numberOfDimensions*(m2)) |
Functions | |
void | insLineSetup (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 &md1a, const int &md1b, const int &md2a, const int &md2b, const int &md3a, const int &md3b, const int &mask, const real &rx, const real &u, const real &gv, const real &dt, real &f, const real &dw, const int &dir, real &am, real &bm, real &cm, real &dm, real &em, const int &bc, const int &boundaryCondition, const int &ndbcd1a, const int &ndbcd1b, const int &ndbcd2a, const int &ndbcd2b, const int &ndbcd3a, const int &ndbcd3b, const int &ndbcd4a, const int &ndbcd4b, const real &bcData, const int &ipar, const real &rpar, const int &ierr) |
void | insLineSetupNew (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 &md1a, const int &md1b, const int &md2a, const int &md2b, const int &md3a, const int &md3b, const int &mask, const real &rx, const real &u, const real &gv, const real &dt, real &f, const real &dw, const int &dir, real &am, real &bm, real &cm, real &dm, real &em, const int &bc, const int &boundaryCondition, const int &ndbcd1a, const int &ndbcd1b, const int &ndbcd2a, const int &ndbcd2b, const int &ndbcd3a, const int &ndbcd3b, const int &ndbcd4a, const int &ndbcd4b, const real &bcData, const int &ipar, const real &rpar, const int &ierr) |
void | insLineSolveBC (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 &md1a, const int &md1b, const int &md2a, const int &md2b, const int &md3a, const int &md3b, const int &mask, const real &rx, const real &u, const real &gv, const real &dt, real &f, const real &dw, const int &dir, real &am, real &bm, real &cm, real &dm, real &em, const int &bc, const int &boundaryCondition, const int &ndbcd1a, const int &ndbcd1b, const int &ndbcd2a, const int &ndbcd2b, const int &ndbcd3a, const int &ndbcd3b, const int &ndbcd4a, const int &ndbcd4b, const real &bcData, const int &ipar, const real &rpar, const int &ierr) |
void | computeResidual (const int &nd, const int &n1a, const int &n1b, const int &n2a, const int &n2b, const int &n3a, const int &n3b, const int &nd1a, const int &nd1b, const int &nd2a, const int &nd2b, const int &nd3a, const int &nd3b, const int &nd4a, const int &nd4b, const int &mask, const real &rx, const real &u, const real &gv, const real &dt, const real &f, const real &dw, real &residual, const int &bc, const int &ipar, const real &rpar, const int &ierr) |
void | computeResidualNew (const int &nd, const int &n1a, const int &n1b, const int &n2a, const int &n2b, const int &n3a, const int &n3b, const int &nd1a, const int &nd1b, const int &nd2a, const int &nd2b, const int &nd3a, const int &nd3b, const int &nd4a, const int &nd4b, const int &mask, const real &rx, const real &u, const real &gv, const real &dt, const real &f, const real &dw, real &residual, const int &bc, const int &ipar, const real &rpar, const int &ierr) |
#define computeResidual EXTERN_C_NAME(computeresidual) |
#define computeResidualNew EXTERN_C_NAME(computeresidualnew) |
Referenced by Cgins::advanceLineSolveNew(), and Cgins::advanceLineSolveOld().
#define FOR_3 | ( | i1, | |
i2, | |||
i3, | |||
I1, | |||
I2, | |||
I3 | |||
) | I1Base =I1.getBase(), I2Base =I2.getBase(), I3Base =I3.getBase(); I1Bound=I1.getBound(), I2Bound=I2.getBound(), I3Bound=I3.getBound(); for(i3=I3Base; i3<=I3Bound; i3++) for(i2=I2Base; i2<=I2Bound; i2++) for(i1=I1Base; i1<=I1Bound; i1++) |
#define FOR_3D | ( | i1, | |
i2, | |||
i3, | |||
I1, | |||
I2, | |||
I3 | |||
) | int I1Base =I1.getBase(), I2Base =I2.getBase(), I3Base =I3.getBase(); int I1Bound=I1.getBound(), I2Bound=I2.getBound(), I3Bound=I3.getBound(); for(i3=I3Base; i3<=I3Bound; i3++) for(i2=I2Base; i2<=I2Bound; i2++) for(i1=I1Base; i1<=I1Bound; i1++) |
#define insLineSetup EXTERN_C_NAME(inslinesetup) |
Referenced by Cgins::advanceLineSolveOld().
#define insLineSetupNew EXTERN_C_NAME(inslinesetupnew) |
Referenced by Cgins::advanceLineSolveNew(), and Cgins::advanceLineSolveOld().
#define insLineSolveBC EXTERN_C_NAME(inslinesolvebc) |
Referenced by Cgins::advanceLineSolveNew(), and Cgins::advanceLineSolveOld().
#define mixedCoeff | ( | component, | |
side, | |||
axis, | |||
grid | |||
) | bcData(component+numberOfComponents*(1),side,axis,grid) |
Referenced by Cgins::assignLineSolverBoundaryConditions(), and Cgins::getLineSolverBoundaryConditions().
#define mixedNormalCoeff | ( | component, | |
side, | |||
axis, | |||
grid | |||
) | bcData(component+numberOfComponents*(2),side,axis,grid) |
Referenced by Cgins::assignLineSolverBoundaryConditions(), and Cgins::getLineSolverBoundaryConditions().
#define mixedRHS | ( | component, | |
side, | |||
axis, | |||
grid | |||
) | bcData(component+numberOfComponents*(0),side,axis,grid) |
Referenced by Cgins::assignLineSolverBoundaryConditions(), and Cgins::getLineSolverBoundaryConditions().
Referenced by Cgins::assignLineSolverBoundaryConditions(), and Cgins::lineSolverBoundaryConditions().
#define tridiagonalSolver | ( | c, | |
axis, | |||
grid | |||
) | lineSolve.pTridiagonalSolvers[c+maxNumberOfSystems*(axis+numberOfDimensions*(grid))] |
Referenced by Cgad::advanceADI(), and Cgins::advanceLineSolveNew().
void computeResidual | ( | const int & | nd, |
const int & | n1a, | ||
const int & | n1b, | ||
const int & | n2a, | ||
const int & | n2b, | ||
const int & | n3a, | ||
const int & | n3b, | ||
const int & | nd1a, | ||
const int & | nd1b, | ||
const int & | nd2a, | ||
const int & | nd2b, | ||
const int & | nd3a, | ||
const int & | nd3b, | ||
const int & | nd4a, | ||
const int & | nd4b, | ||
const int & | mask, | ||
const real & | rx, | ||
const real & | u, | ||
const real & | gv, | ||
const real & | dt, | ||
const real & | f, | ||
const real & | dw, | ||
real & | residual, | ||
const int & | bc, | ||
const int & | ipar, | ||
const real & | rpar, | ||
const int & | ierr | ||
) |
void computeResidualNew | ( | const int & | nd, |
const int & | n1a, | ||
const int & | n1b, | ||
const int & | n2a, | ||
const int & | n2b, | ||
const int & | n3a, | ||
const int & | n3b, | ||
const int & | nd1a, | ||
const int & | nd1b, | ||
const int & | nd2a, | ||
const int & | nd2b, | ||
const int & | nd3a, | ||
const int & | nd3b, | ||
const int & | nd4a, | ||
const int & | nd4b, | ||
const int & | mask, | ||
const real & | rx, | ||
const real & | u, | ||
const real & | gv, | ||
const real & | dt, | ||
const real & | f, | ||
const real & | dw, | ||
real & | residual, | ||
const int & | bc, | ||
const int & | ipar, | ||
const real & | rpar, | ||
const int & | ierr | ||
) |
void insLineSetup | ( | 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 & | md1a, | ||
const int & | md1b, | ||
const int & | md2a, | ||
const int & | md2b, | ||
const int & | md3a, | ||
const int & | md3b, | ||
const int & | mask, | ||
const real & | rx, | ||
const real & | u, | ||
const real & | gv, | ||
const real & | dt, | ||
real & | f, | ||
const real & | dw, | ||
const int & | dir, | ||
real & | am, | ||
real & | bm, | ||
real & | cm, | ||
real & | dm, | ||
real & | em, | ||
const int & | bc, | ||
const int & | boundaryCondition, | ||
const int & | ndbcd1a, | ||
const int & | ndbcd1b, | ||
const int & | ndbcd2a, | ||
const int & | ndbcd2b, | ||
const int & | ndbcd3a, | ||
const int & | ndbcd3b, | ||
const int & | ndbcd4a, | ||
const int & | ndbcd4b, | ||
const real & | bcData, | ||
const int & | ipar, | ||
const real & | rpar, | ||
const int & | ierr | ||
) |
void insLineSetupNew | ( | 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 & | md1a, | ||
const int & | md1b, | ||
const int & | md2a, | ||
const int & | md2b, | ||
const int & | md3a, | ||
const int & | md3b, | ||
const int & | mask, | ||
const real & | rx, | ||
const real & | u, | ||
const real & | gv, | ||
const real & | dt, | ||
real & | f, | ||
const real & | dw, | ||
const int & | dir, | ||
real & | am, | ||
real & | bm, | ||
real & | cm, | ||
real & | dm, | ||
real & | em, | ||
const int & | bc, | ||
const int & | boundaryCondition, | ||
const int & | ndbcd1a, | ||
const int & | ndbcd1b, | ||
const int & | ndbcd2a, | ||
const int & | ndbcd2b, | ||
const int & | ndbcd3a, | ||
const int & | ndbcd3b, | ||
const int & | ndbcd4a, | ||
const int & | ndbcd4b, | ||
const real & | bcData, | ||
const int & | ipar, | ||
const real & | rpar, | ||
const int & | ierr | ||
) |
void insLineSolveBC | ( | 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 & | md1a, | ||
const int & | md1b, | ||
const int & | md2a, | ||
const int & | md2b, | ||
const int & | md3a, | ||
const int & | md3b, | ||
const int & | mask, | ||
const real & | rx, | ||
const real & | u, | ||
const real & | gv, | ||
const real & | dt, | ||
real & | f, | ||
const real & | dw, | ||
const int & | dir, | ||
real & | am, | ||
real & | bm, | ||
real & | cm, | ||
real & | dm, | ||
real & | em, | ||
const int & | bc, | ||
const int & | boundaryCondition, | ||
const int & | ndbcd1a, | ||
const int & | ndbcd1b, | ||
const int & | ndbcd2a, | ||
const int & | ndbcd2b, | ||
const int & | ndbcd3a, | ||
const int & | ndbcd3b, | ||
const int & | ndbcd4a, | ||
const int & | ndbcd4b, | ||
const real & | bcData, | ||
const int & | ipar, | ||
const real & | rpar, | ||
const int & | ierr | ||
) |