Laghos fails with message -- unknown kernel -- higher order test with PA and higher order runs with -ok >= 8 fails with FA

[cjoycd]

Dear all,

I was trying to test three-material, triple point problem using the following command.

mpirun -np 1 ./laghos -p 3 -m data/rectangle01_quad.mesh -rs 1 -tf 5.0 -pa -ok 5 -ot 4

and it shows the following error message and stops.

Unknown kernel 0x2a

MFEM abort: Unknown kernel
... in function: void mfem::hydrodynamics::QUpdate::UpdateQuadratureData(cont mfem::Vector&, mfem::hydrodynamics::QuadratureData&)
... in file: laghos_solver.cpp:1329

This happens when I choose PA and -ok greater than 4 (order of the kinematic field).

Similarly, with FA, for the same problem with -ok 8 -ot 7, it shows the following error.

Verification failed: (Q1D <= MQ) is false.
--> Quadrature rules more than 14 1D points are not supported!
... in function: void mfem::internal::quadrature_interpolator::TensorValues(int, int, const mfem::DofToQuad&, const mfem::Vector&, mfem::Vector&) [with mfem::QVectorLayout::byVDIM]
... in file: fem/qinterp/eval_by_vdim.cpp:74

Please suggest..

Thank you in advance..

Sijoy C.D.

[vladotomov]

Hi Sijoy,

The first error happens because we don't have that kernel instance. You can get it running with this #180. You can add higher orders the same way.

The second verification is an MFEM restriction, see here. You can create an MFEM issue and request more 1D quad points. I guess nobody has tried to run so high yet.

[vladotomov]

@cjoycd you can also try with mfem/mfem#3722

[cjoycd]

Dear Vladimir, Okay, I will follow the suggestions, thanks a lot…

…

On Wed, 6 Sep 2023 at 12:32 AM, Vladimir Tomov ***@***.***> wrote: You can also try with mfem/mfem#3722 <mfem/mfem#3722> — Reply to this email directly, view it on GitHub <#179 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AIUNJIEW4GBQGVLAOLSOLRTXY5ZLFANCNFSM6AAAAAA4J4CKLM> . You are receiving this because you authored the thread.Message ID: ***@***.***>