How to trunc in a parallel code

Dear all,

Is there a command to trunc in parallel code.
For example, when I solve fluid structure interaction problem, I need to trunc the solid part. It is easy to implement in sequential code like Ths = trunc(Th,region==solid) , however, in parallel code, the mesh is seperated and then I find it will report error if I directly copy from sequential code without change. So my question is how to trunc partial mesh in a parallel code?

You can trunc your global mesh before domain(mesh) decomposition.

Yes, I know that. But my mesh will move every time step, so I still need trunc every time to track the solid part. For example,



these are two different time configuration, I need to trunc the solid part every time step to calculate the solid part equation. If not, my code is wrong, so it’s inevitable to trunc the solid part after domain decomposition.

So is there method to achieve this?

It is not inevitable, why are you saying this?

Sorry, I’m being too absolute. Just in my code, it’s inevitable to trunc the solid part, because in my equation, the solid part is keeping changing, I need to track it every time step.

new_ale_method.edp (3.5 KB)
This is my sequential fsi code, now I want to parallelize it. And the problem is how I change this line in parallel version

Keep a backup of the global mesh, centralize your distributed solution, move your mesh, perform the domain decomposition. It’s done in multiple examples in the examples/hpddm folder.

Thank you very much for your precious advice, professor. Let me have a look.