Linear solve did not converge due to DIVERGED_NANORINF iterations

Hi,

I am now implementing 2D fluid structure interaction in parallel. I try to simulate a benchmark: oscillating flag. Here is the code

ParallelOscillatingFlag.edp (5.8 KB)

It works well when I use more than 2 processors (ff-mpirn -n 2 or more). However, there exist a problem when I choose only precossor, says

Does anyone have an idea? Here I directly use ‘lu‘ as my conditioner.

Switch from PETSc to MUMPS for doing the factorization and it will work even with a single process.

Thanks, professor.

I am not familiar with MUMPS. I have a few questions.

Do you mean replace “load ‘PETSc’ ‘‘ with “load ‘MUMPS’ ‘‘? Can the distributed mesh and Mat (in PETSc) be still remained in MUMPS. Is there a simple example to start with MUMPS?

No, you just need to set the proper sparams, it’s in many examples from examples/hpddm.