message: Setting up GeomToPerfGeomJac for pyramids
diff: #include "fem.hpp"
  #include "../mesh/wedge.hpp"
  +#include "../mesh/pyramid.hpp"
  namespace mfem
  {
  @@ -241,6 +242,14 @@ Geometry::Geometry()
  *GeomToPerfGeomJac[PRISM] = pri_T.Jacobian();
  CalcInverse(pri_T.Jacobian(), *PerfGeomToGeomJac[PRISM]);
  }
  + {
  + IsoparametricTransformation pyr_T;
  + pyr_T.SetFE(&PyramidFE);
  + GetPerfPointMat (PYRAMID, pyr_T.GetPointMat());
  + pyr_T.SetIntPoint(&GeomCenter[PYRAMID]);
  + *GeomToPerfGeomJac[PYRAMID] = pyr_T.Jacobian();
  + CalcInverse(pyr_T.Jacobian(), *PerfGeomToGeomJac[PYRAMID]);
  + }
  }
  Geometry::~Geometry()
  