zxss
November 5, 2025, 4:13pm
1
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.
prj
November 5, 2025, 7:17pm
2
Switch from PETSc to MUMPS for doing the factorization and it will work even with a single process.
zxss
November 5, 2025, 10:09pm
3
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?
prj
November 6, 2025, 5:06am
4
No, you just need to set the proper sparams, it’s in many examples from examples/hpddm.