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