message: Add check for quadrature spaces that are invalid on mixed meshes
diff: @@ -82,6 +82,8 @@ QuadratureSpace::QuadratureSpace(Mesh *mesh_, std::istream &in)
  QuadratureSpace::QuadratureSpace(Mesh &mesh_, const IntegrationRule &ir)
  : QuadratureSpaceBase(mesh_, mesh_.GetElementGeometry(0), ir)
  {
  + MFEM_VERIFY(mesh.GetNumGeometries(mesh.Dimension()) == 1,
  + "Constructor not valid for mixed meshes");
  ConstructOffsets();
  }
  @@ -107,6 +109,8 @@ FaceQuadratureSpace::FaceQuadratureSpace(Mesh &mesh_, const IntegrationRule &ir,
  face_type(face_type_),
  num_faces(mesh.GetNFbyType(face_type))
  {
  + MFEM_VERIFY(mesh.GetNumGeometries(mesh.Dimension() - 1) == 1,
  + "Constructor not valid for mixed meshes");
  ConstructOffsets();
  }
  