#include "InsParameters.h"
#include "TridiagonalSolver.h"
#include "insFactors.h"
#include "CompositeGrid.h"
#include "GridFunction.h"
#include "ParallelUtility.h"
#include "ParallelGridUtility.h"
#include "PlotIt.h"
#include "AdamsPCData.h"
#include "kkcdefs.h"

Include dependency graph for insFactors.C:

Macros
#define	ins_rfactor EXTERN_C_NAME(ins_rfactor)

#define	ins_rrfactor EXTERN_C_NAME(ins_rrfactor)

#define	ins_diagfactor EXTERN_C_NAME(ins_diagfactor)

#define	ins_mfactor EXTERN_C_NAME(ins_mfactor)

#define	ins_mfactor_opt EXTERN_C_NAME(ins_mfactor_opt)

#define	ins_fscoeff EXTERN_C_NAME(ins_fscoeff)

#define	ins_evalux EXTERN_C_NAME(ins_evalux)

#define	DEFINE_AF_FACTOR_SUBROUTINE(NAME)

#define	SETUP_PARAMETER_ARRAYS

#define	CALL_INS_FACTOR_ROUTINE(NAME)

#define	PRINT_MATRIX_ARRAYS(NAME)

#define	IS_DIRICHLET(SIDE, AXIS) ( bcLocal(SIDE,AXIS)==Parameters::dirichletBoundaryCondition \|\| bcLocal(SIDE,AXIS)==Parameters::noSlipWall \|\| bcLocal(SIDE,AXIS)==InsParameters::inflowWithVelocityGiven \|\| bcLocal(SIDE,AXIS)==InsParameters::outflow)

#define	GET_BOUNDS_INFO(ARRAY, LOCAL_ARRAY, EXTRA_FACTOR)

#define	ZERO_GHOST(LOCAL_ARRAY)

#define	REBUILD_AND_REFACTOR(CC) (((CC)==uc) \|\| (tc>wc && CC>=tc))

#define	USE_MFACTOR_OPT

#define	ug_starLocal ugLocal

#define	ug_starLocal ugLocal

#define	ug_starLocal ugLocal

#define	ug_starLocal fgLocal

#define	ug_starLocal ugLocal

#define	ug_starLocal uxy

#define	ugLocal uxy

Functions
	DEFINE_AF_FACTOR_SUBROUTINE (ins_rfactor)

	DEFINE_AF_FACTOR_SUBROUTINE (ins_rrfactor)

	DEFINE_AF_FACTOR_SUBROUTINE (ins_diagfactor)

	DEFINE_AF_FACTOR_SUBROUTINE (ins_mfactor)

	DEFINE_AF_FACTOR_SUBROUTINE (ins_mfactor_opt)

	DEFINE_AF_FACTOR_SUBROUTINE (ins_fscoeff)

	DEFINE_AF_FACTOR_SUBROUTINE (ins_evalux)

typedef	DEFINE_AF_FACTOR_SUBROUTINE ((*INSFactorSubroutine))

	getGhostIndex (gidLocal, side, axis, Ig1, Ig2, Ig3, gl)

	LOCAL_ARRAY (Ig1, Ig2, Ig3)=0.0

Variables
Index	Ig1

Index	Ig2

Index	Ig3

Macro Definition Documentation

#define CALL_INS_FACTOR_ROUTINE ( NAME )

Value:

NAME(numberOfDimensions, \
                               ug_starLocal.getBase(0), ug_starLocal.getBound(0),\
                               ug_starLocal.getBase(1), ug_starLocal.getBound(1),\
                               ug_starLocal.getBase(2), ug_starLocal.getBound(2),\
                               ug_starLocal.getBase(3), ug_starLocal.getBound(3),\
                               *maskLocal.getDataPointer(),\
                               *rsxyLocal.getDataPointer(), *ugLocal.getDataPointer(), *ug_starLocal.getDataPointer(), \
                               *gridVelocityLocal.getDataPointer(),\
                               bcLocal(0,0), mg.boundaryCondition(0,0), \
                               bcData.getBase(0),bcData.getBound(0),\
                               bcData.getBase(1),bcData.getBound(1),\
                               bcData.getBase(2),bcData.getBound(2),\
                               bcData.getBase(3),bcData.getBound(3),\
                               *bcData.getDataPointer(),\
                               ipar[0],rpar[0],pdb,\
                               mode, dir, component_to_solve_for,\
                               *dl2.getDataPointer(), *dl1.getDataPointer(), *d.getDataPointer(), \
                               *du1.getDataPointer(), *du2.getDataPointer(), \
                               *rhsLocal.getDataPointer(), ierr);

Referenced by CGINS_ApproximateFactorization::INS_Factor::addExplicitContribution(), CGINS_ApproximateFactorization::INS_Factor::solveLeftHandSide(), and CGINS_ApproximateFactorization::INS_Factor::solveRightHandSide().

#define DEFINE_AF_FACTOR_SUBROUTINE ( NAME )

Value:

void NAME(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 &mask, \
                                                    const real &rsxy, const real &u, const real &ul, const real &gv,  \
                                                    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, DataBase *pdb, \
                                                    const int &mode, const int&dir, const int &component,\
                                                    real &a, real &b, real &c, real &d, real &e, real &rhs, int &ierr )

#define GET_BOUNDS_INFO	(	ARRAY,
		LOCAL_ARRAY,
		EXTRA_FACTOR
	)

Referenced by CGINS_ApproximateFactorization::INS_Factor::addExplicitContribution(), CGINS_ApproximateFactorization::INS_Factor::solveLeftHandSide(), and CGINS_ApproximateFactorization::INS_Factor::solveRightHandSide().

#define ins_diagfactor EXTERN_C_NAME(ins_diagfactor)

#define ins_evalux EXTERN_C_NAME(ins_evalux)

Referenced by CGINS_ApproximateFactorization::INS_Factor::addExplicitContribution().

#define ins_fscoeff EXTERN_C_NAME(ins_fscoeff)

Referenced by CGINS_ApproximateFactorization::INS_Factor::addExplicitContribution(), CGINS_ApproximateFactorization::INS_Factor::solveLeftHandSide(), and CGINS_ApproximateFactorization::INS_Factor::solveRightHandSide().

#define ins_mfactor EXTERN_C_NAME(ins_mfactor)

#define ins_mfactor_opt EXTERN_C_NAME(ins_mfactor_opt)

#define ins_rfactor EXTERN_C_NAME(ins_rfactor)

#define ins_rrfactor EXTERN_C_NAME(ins_rrfactor)

#define IS_DIRICHLET	(	SIDE,
		AXIS
	)	( bcLocal(SIDE,AXIS)==Parameters::dirichletBoundaryCondition \|\| bcLocal(SIDE,AXIS)==Parameters::noSlipWall \|\| bcLocal(SIDE,AXIS)==InsParameters::inflowWithVelocityGiven \|\| bcLocal(SIDE,AXIS)==InsParameters::outflow)

Referenced by CGINS_ApproximateFactorization::INS_Factor::addExplicitContribution(), CGINS_ApproximateFactorization::INS_Factor::solveLeftHandSide(), and CGINS_ApproximateFactorization::INS_Factor::solveRightHandSide().

#define PRINT_MATRIX_ARRAYS ( NAME )

Value:

if ( parameters->dbase.get<int >("debug")>=4 ) \
                { \
                  cout<<"INS_FACTORS " #NAME " debug :: after ins_factor["<<type<<"], dir="<<getDirection()<<", cc="<<component_to_solve_for<<" : "<<endl; \
                  if ( isPenta) dl2.display("dl2"); \
                  dl1.display("dl1"); \
                  d.display("d"); \
                  du1.display("du1"); \
                  if ( isPenta ) du2.display("du2"); \
                  rhsLocal.display("rhs"); \
                } \