Calling createMesh externally in v24

[aprilnovak]

In Cardinal, we need to be able to call createMesh external to the NekRS code base. I think I've updated the function call correctly, to:

      auto [meshT, meshV] = createMesh(platform->comm.mpiComm, _order + 1, 0, nrs->cht, platform->kernelInfo);

but I get an error message at runtime due to a kernel not being found. Here's the whole error:

*** ERROR ***
Aborting in /home/anovak/nekRS/src/core/kernelRequestManager.cpp:load

Cannot find request <mesh-geometricFactorsHex3D_2>
Available:
	</home/anovak/cardinal/install/kernels/core/ogs/gather.okl>
	</home/anovak/cardinal/install/kernels/core/ogs/gatherMany.okl>
	</home/anovak/cardinal/install/kernels/core/ogs/gatherScatter.okl>
	</home/anovak/cardinal/install/kernels/core/ogs/gatherScatterMany.okl>
	</home/anovak/cardinal/install/kernels/core/ogs/gatherScatterNew.okl>
	</home/anovak/cardinal/install/kernels/core/ogs/oogs.okl>
	</home/anovak/cardinal/install/kernels/core/ogs/scatter.okl>
	</home/anovak/cardinal/install/kernels/core/ogs/scatterMany.okl>
	</home/anovak/cardinal/install/kernels/nrs/subCycleStrongCubatureVolumeHex3D.c::8c2199f8389c4bd0>
	<EToLdfloat>
	<EToLpfloat>
	<LToEdfloat>
	<LToEpfloat>
	<absoluteError>
	<add>
	<ady>
	<adyMany>
	<adyz>
	<amax>
	<amaxMany>
	<axdy>
	<axdyz>
	<axmy>
	<axmyMany>
	<axmyVector>
	<axmyz>
	<axmyzMany>
	<axpby>
	<axpbyMany>
	<axpbyz>
	<axpbyzMany>
	<aydx>
	<aydxMany>
	<cds-advectMeshVelocityHex3D>
	<cds-dirichletBC>
	<cds-maskCopy>
	<cds-maskCopy2>
	<cds-neumannBCHex3D>
	<cds-strongAdvectionCubatureVolumeHex3D>
	<cds-strongAdvectionVolumeHex3D>
	<cds-subCycleStrongVolumeHex3D>
	<computeFlux>
	<computeMaxVisc>
	<copyNekNekPoints>
	<core-copyDfloatToDouble>
	<core-copyDfloatToFloat>
	<core-copyDfloatToPfloat>
	<core-copyDoubleToDfloat>
	<core-copyFloatToDfloat>
	<core-copyPfloatToDfloat>
	<core-curlHex3D>
	<core-divergenceVolumeHex3D>
	<core-filterRTHex3D>
	<core-gradientVolumeHex3D>
	<core-nStagesSum3>
	<core-tensorProduct1DHex3D>
	<core-vectorFilterRTHex3D>
	<core-wDivergenceVolumeHex3D>
	<core-wGradientVolumeHex3D>
	<core-weakLaplacianHex3D>
	<crossProduct>
	<entrywiseMag>
	<extrapolateBoundary>
	<fill>
	<findptsLocalEvalMask_1>
	<findptsLocalEvalMask_10>
	<findptsLocalEvalMask_11>
	<findptsLocalEvalMask_12>
	<findptsLocalEvalMask_13>
	<findptsLocalEvalMask_14>
	<findptsLocalEvalMask_15>
	<findptsLocalEvalMask_2>
	<findptsLocalEvalMask_3>
	<findptsLocalEvalMask_4>
	<findptsLocalEvalMask_5>
	<findptsLocalEvalMask_6>
	<findptsLocalEvalMask_7>
	<findptsLocalEvalMask_8>
	<findptsLocalEvalMask_9>
	<findptsLocalEval_1>
	<findptsLocalEval_10>
	<findptsLocalEval_11>
	<findptsLocalEval_12>
	<findptsLocalEval_13>
	<findptsLocalEval_14>
	<findptsLocalEval_15>
	<findptsLocalEval_2>
	<findptsLocalEval_3>
	<findptsLocalEval_4>
	<findptsLocalEval_5>
	<findptsLocalEval_6>
	<findptsLocalEval_7>
	<findptsLocalEval_8>
	<findptsLocalEval_9>
	<findptsLocal_1>
	<findptsLocal_10>
	<findptsLocal_11>
	<findptsLocal_12>
	<findptsLocal_13>
	<findptsLocal_14>
	<findptsLocal_15>
	<findptsLocal_2>
	<findptsLocal_3>
	<findptsLocal_4>
	<findptsLocal_5>
	<findptsLocal_6>
	<findptsLocal_7>
	<findptsLocal_8>
	<findptsLocal_9>
	<fixSurfaceFlux>
	<gatherPlanarValuesXY>
	<gatherPlanarValuesXZ>
	<gatherPlanarValuesYZ>
	<hlong-mesh-sum>
	<innerProd>
	<interpolateP1>
	<linearCombination>
	<magSqrSymTensor>
	<magSqrSymTensorDiag>
	<magSqrTensor>
	<magSqrVector>
	<mapScalar>
	<mask>
	<max>
	<mesh-coarsenHex3D_Nf_10_Nc_5>
	<mesh-coarsenHex3D_Nf_11_Nc_5>
	<mesh-coarsenHex3D_Nf_12_Nc_5>
	<mesh-coarsenHex3D_Nf_13_Nc_5>
	<mesh-coarsenHex3D_Nf_14_Nc_5>
	<mesh-coarsenHex3D_Nf_15_Nc_5>
	<mesh-coarsenHex3D_Nf_5_Nc_1>
	<mesh-coarsenHex3D_Nf_5_Nc_2>
	<mesh-coarsenHex3D_Nf_5_Nc_3>
	<mesh-coarsenHex3D_Nf_5_Nc_4>
	<mesh-coarsenHex3D_Nf_5_Nc_5>
	<mesh-coarsenHex3D_Nf_6_Nc_5>
	<mesh-coarsenHex3D_Nf_7_Nc_5>
	<mesh-coarsenHex3D_Nf_8_Nc_5>
	<mesh-coarsenHex3D_Nf_9_Nc_5>
	<mesh-cubatureGeometricFactorsHex3D_5>
	<mesh-distanceHex3D>
	<mesh-geometricFactorsHex3D_5>
	<mesh-nStagesSumVector>
	<mesh-prolongateHex3D_Nf_10_Nc_5>
	<mesh-prolongateHex3D_Nf_11_Nc_5>
	<mesh-prolongateHex3D_Nf_12_Nc_5>
	<mesh-prolongateHex3D_Nf_13_Nc_5>
	<mesh-prolongateHex3D_Nf_14_Nc_5>
	<mesh-prolongateHex3D_Nf_15_Nc_5>
	<mesh-prolongateHex3D_Nf_5_Nc_1>
	<mesh-prolongateHex3D_Nf_5_Nc_2>
	<mesh-prolongateHex3D_Nf_5_Nc_3>
	<mesh-prolongateHex3D_Nf_5_Nc_4>
	<mesh-prolongateHex3D_Nf_5_Nc_5>
	<mesh-prolongateHex3D_Nf_6_Nc_5>
	<mesh-prolongateHex3D_Nf_7_N

*** ERROR ***

Do I need to change something in Cardinal's makefile to be able to link to/find the missing kernel? Or, if I'm misinterpreting something here, would you please help me to diagnose what is going wrong?

Thank you!

[stgeke]

What's the usecase for calling createMesh? It looks like the polynomial order is different?

[aprilnovak]

Yes, we need it in order to create a MOOSE version of NekRS's mesh to facilitate data transfers. It is used to create 27-node brick elements in libMesh corresponding to the (N+1)^3 brick elements in NekRS. We need it in order to figure out the face-center node positions on the 6 hex faces in the NekRS bricks.