ceed / remhos Goto Github PK

Dear contributors,

I have been using the modified version of Laghos by implementing RZ (Cylindrical) geometry as well as multi-material formulations. Presently looking forward to implement a remap scheme using REMHOS miniapp. As far as I can understand, the RZ capability is not implemented with the present REMHOS version. I would like to know is there any future version planned to include this capability ? What are the changes that might required to make it compatible with RZ geometry (like changes in Convection Integrator, Mass Integrator, Trace Integrators, ....) ? Is there any documentation available ?

Thank you in advance
-- Sijoy

I noticed sometimes we run into floating point issues using inf(). Could we consider using max() instead?

Hi all,

I'm incorporating Remhos into my code to remap the L2 fields after remeshing.
I believe the remapping is functioning correctly with the conforming mesh, but I'm encountering some issues with a non-conforming mesh.

In the case of a conforming mesh, the mesh simply moves, and the field does not move because they are advected together. However, with a non-conforming mesh, the shape of the field is altered, leading to incorrect mapping (see attached image).

This is obvious when I use the example of Remhos.
I typed mpirun -np 8 remhos -m ./data/cube01_hex.mesh -p 10 -rs 1 -o 2 -dt 0.02 -tf 0.8 -ho 1 -lo 4 -fct 2 but I modified element attributes of 1st and 2nd element to two and Remhos like this way. You can see non-conforming mesh and final remapped field is highly distorted. I think advection is prohibited between coarse- and fine-mesh regions.

mesh->EnsureNCMesh(true);
Array<int> refs;

for (int i = 0; i < mesh->GetNE(); i++)
{
if(mesh->GetAttribute(i) >= 2)
{
refs.Append(i);
}
}

mesh->GeneralRefinement(refs, 1);
refs.DeleteAll();
mesh->Finalize(true);

Any comments would be greatly appreciated.

Sungho

I tried building with mfem develop, and ran into the following error
cannot find -lmfem-common. I think we may need to update our makefile.

Hi all,

I have a quick question
Can we make Remhos work for low-order elements like (order<2)?

Best,

Sungho