CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)
-   -   New Solver for a modified pUCoupledFoam (https://www.cfd-online.com/Forums/openfoam/175598-new-solver-modified-pucoupledfoam.html)

noiseMaker August 2, 2016 06:14

New Solver for a modified pUCoupledFoam
 
Hello :),

As part of my master thesis I am suppose to modify pUCoupledFoam in a way so that the matrix can be solved with an LU-SGS solver. However, when I try to run a case (pitzDaily) with what I have done (modified pUCoupledFoam where the matrix is solver by myDummy, the LU-SGS solver), I always get this error:

--> FOAM FATAL IO ERROR:
Unknown asymmetric matrix solver myDummy


"Valid asymmetric matrix solvers are :

4
(
AMG
BiCGStab
GMRES
GaussSeidel
)


file: /home/nelly/Desktop/pitzDaily/system/fvSolution::solvers::Up from line 21 to line 28.

From function BlockLduSolver<Type>::New
(
const word& solverName,
const word& fieldName,
BlockLduMatrix<Type>& matrix,
const dictionary& dict
)
in file /opt/foam/foam-extend-3.2/src/foam/lnInclude/BlockLduSolver.C at line 245.
"



even though the source file reads:
namespace Foam
{
#define makeSolver(type, Type, args...) \
typedef myDummySolver<type> block##Type##myDummySolver; \
makeBlockSolverTypeName(block##Type##myDummy); \
addSolverToBlockMatrix(Type, block##Type##myDummySolver, symMatrix); \
addSolverToBlockMatrix(Type, block##Type##myDummySolver, asymMatrix); \

// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //

makeBlockSolvers(myDummy);

// BiCGStab as a symmetric solver: use where the diagonal coefficient block
// is assymetric. HJ, 24/Jan/2008

addSymSolverToBlockMatrix(myDummy);

addAsymSolverToBlockMatrix(myDummy);

//-makeBlockSolvers(blockAmgSolver);
//+// Excluded tensor: needs reduced version. HJ, 28/May/2014


makeBlockSolverTypeName(myDummyScalar);
makeBlockSolverTypeName(myDummyVector);

addSolverToBlockMatrix(Scalar, myDummyScalar, symMatrix);
addSolverToBlockMatrix(Vector, myDummyVector, symMatrix);

addSolverToBlockMatrix(Scalar, myDummyScalar, asymMatrix);
addSolverToBlockMatrix(Vector, myDummyVector, asymMatrix);




Does anybody have clue?

Thanks

Noisemaker


All times are GMT -4. The time now is 19:27.