[noiseMaker]

Hello,

I am currently doing my master thesis in CFD using OpenFoam and one of the tasks I have is to implement and test a new linear solver (called myDummy) for a modified verison of pUCoupledFoam (i.e. solve the Up matrix with myDummy instead of BiCStab, for instance)

However, everytime I try to run a case I get the following fatal 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.

FOAM exiting
>>

This error appears because the linear solver myDummy is not being added to the hash table of asymmetric solvers.



The Solver header file myDummys.H reads:
<<
typedef BlockBiCGStabSolver<scalar> blockBiCGStabSolverScalar;
typedef BlockBiCGStabSolver<vector> blockBiCGStabSolverVector;
typedef BlockBiCGStabSolver<tensor> blockBiCGStabSolverTensor;
>>


And the souce file myDummys.C reads:

<<
namespace Foam
{

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

makeBlockSolvers(blockMyDummySolver);

addSymSolverToBlockMatrix(blockMyDummySolver);

addAsymSolverToBlockMatrix(blockMyDummySolver);


} // End namespace Foam


Both the solver and the linera solver compile. I just can't run the case.

1.Does anybody have an idea on what else I could do in order to add the linear solver myDummy to the hash table of solvers of asymmetric matrices?

2. Am I using the right macros here?

Thank you!!

Noisemaker

[Tobi]

Hi,

I made a test with a normal solver (that is very easy):

--> FOAM FATAL IO ERROR: 
Unknown asymmetric matrix solver PBiC
Valid asymmetric matrix solvers are :

13
(
BICCG
BiCG
BiCGStab
FPEAMG
GAMG
GMRES
MPEAMG
PBiCG
RREAMG
Shorty
amgSolver
deflation
smoothSolver
)

Then I tried with a blockCoupled one and had to figure out the place to manipulate / insert the new one. For the coupled stuff you need to make new matrix solver here:

cp -r foam/matrices/blockLduMatrix/BlockLduSolvers/BlockBiCGStab/ cp foam/matrices/blockLduMatrix/BlockLduSolvers/myNewGuy
cd foam/matrices/blockLduMatrix/BlockLduSolvers/myNewGuy

manipulate everything and compile

The result for the pUCoupledFoam is then:

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

Valid asymmetric matrix solvers are :

5
(
AMG
BiCGStab
MyNewGuy
GMRES
GaussSeidel
)

Hope this is the answer to your question. If you manipulate and insert the solver in coupledMatrix, it will not find the stuff. You have to figure out yourself how they are connected but that was not your question

[noiseMaker]

So, you're saying that I am suppose to copy myDummy folder in the solvers folders and compile it there?

[Tobi]

No, i am saying that the solvers that are available in the hash table for the coupled solvers are based on the matrix solvers. So you should copy one of them and build your own one based on that stuff.

I just made a test. So you have to check the hierarchy of the libraries yourself.

[noiseMaker]

Dear Tobi,

In order to understand what is going wrong I ran the test you showed me in the hope of obtaining similar results.

So I did the following:

1- Copy BlockBiCGStabSolver to Desktop
2- Change the name from BlockbiCGStabSolver to BlockNewGuySolver without changing the core of the code
3- Compile and test

The problem didn't go away.

So I am thinking that maybe copying the Folder of the solver into foam-extend and compiling the whole thing might be the solution (the reason I didn't do that is because I don't have root access at the institute, where I am doing my Master thesis). Could it be?

Sorry to bother you with something as trivial as this. And thank you so much for your help.

Best Regards,

Nelly

[Tobi]

Hi,

you (do not) have to compile in the foam-extend folder because you have to modify the Make/files file and add your new library (solver). Otherwise foam will not find it.

If you make a standalone lib, then you have to include that in the controlDict like

libs
(
    BlockNewGuySolver.so
);

I do not know out of the box how you should do it but it will also work (I think).

• If I understand you correct, you compile your new solver on the desctop because of root access.
• For that you made a new Make folder and files + options (or)?
• The name of the library should be included in your controlDict
• and of course, wmake libso (:

Hope it is clear and it will help you.


PS: I think in the user guide there is a section about "How to include own libraries".

[noiseMaker]

Yes Yes Yes. I did exactly that. Still not working

[Tobi]

Will try it

[noiseMaker]

Important update

there is a source file called blockVectorNSolvers.C. It contains macros for the definition and inclusion of linear solvers into the hash table. I think this is the one I need to manipulate in order to get myDummy added to the list at runtime.

But I can't do it today. I am compiling and OF extended version that I can manipulate and modify whichever way I want without root access. Seemed like the easiest option.

Let me know how it goes Tobi,

Liebe Grüße ,

Noisemaker

[Tobi]

Okay, if you use a local compilation, then do it first. If you have still problems, I will check it out (always busy and the best things are, if people solve their problems themselves )

Good luck.

[noiseMaker]

it's like that:

1- copy/paste/changeName of solver
2- manipulate blockVectorNSolvers.C (add one #include and one typedef)
3- In the files file: - blockVectorNSolvers.C
- BlockMyNewGuySolvers.C
4- wmake
5- wmake libso


Result

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

Valid asymmetric matrix solvers are :

5
(
AMG
BiCGStab
GMRES
GaussSeidel
MyNewGuy
)