mini-EM fails when you run higher order
Created by: bathmatt
@egphill should the higher order work? It is failing here
#4 0x000000000063f55a in Intrepid2::Orientation::getEdgeOrientation (this=0x1af0740,
edgeOrt=0x1226e80, numEdge=54)
at /home/mbetten/Trilinos/Trilinos/packages/intrepid2/src/Orientation/Intrepid2_OrientationDef.hpp:480
#5 0x0000000000628efc in addDiscreteGradientToRequestHandler (linObjFactory=..., reqHandler=...)
at /home/mbetten/Trilinos/Trilinos/packages/panzer/mini-em/src/solvers/MiniEM_DiscreteGradient.hpp:111
#6 0x000000000062c81c in main (argc=2, argv=0x7fffffff9608)
at /home/mbetten/Trilinos/Trilinos/packages/panzer/mini-em/example/BlockPrec/main.cpp:362
It says that you have 54 edges, vs between 3 and 12. (tri and hex). It looks like there is an issue with the way you are using the orientations.
@kyungjoo-kim can you help look at this as well?
To reproduce, simply build mini-EM and run
./PanzerMiniEM_BlockPrec.exe --basis-order=2
@rppawlo @eric-c-cyr any ideas?
This is a modestly low priority as there are more issues to solve with mini-EM but don't want it off the radar.