message: adding checks for input based on variable order space
diff: @@ -78,6 +78,8 @@ void FindPointsGSLIB::Setup(Mesh &m, const double bb_t, const double newt_tol,
  MFEM_VERIFY(m.GetNodes() != NULL, "Mesh nodes are required.");
  MFEM_VERIFY(m.GetNumGeometries(m.Dimension()) == 1,
  "Mixed meshes are not currently supported in FindPointsGSLIB.");
  + MFEM_VERIFY(!(m.GetNodes()->FESpace()->IsVariableOrder()),
  + "Variable order mesh is not currently supported.");
  // call FreeData if FindPointsGSLIB::Setup has been called already
  if (setupflag) { FreeData(); }
  @@ -590,7 +592,8 @@ void FindPointsGSLIB::Interpolate(const GridFunction &field_in,
  const L2_FECollection *fec_l2 = dynamic_cast<const L2_FECollection *>(fec_in);
  if (fec_h1 && gf_order == mesh_order &&
  - fec_h1->GetBasisType() == BasisType::GaussLobatto)
  + fec_h1->GetBasisType() == BasisType::GaussLobatto &&
  + !field_in.FESpace()->IsVariableOrder())
  {
  InterpolateH1(field_in, field_out);
  return;
  @@ -857,6 +860,8 @@ void OversetFindPointsGSLIB::Setup(Mesh &m, const int meshid,
  MFEM_VERIFY(m.GetNodes() != NULL, "Mesh nodes are required.");
  MFEM_VERIFY(m.GetNumGeometries(m.Dimension()) == 1,
  "Mixed meshes are not currently supported in FindPointsGSLIB.");
  + MFEM_VERIFY(!(m.GetNodes()->FESpace()->IsVariableOrder()),
  + "Variable order mesh is not currently supported.");
  // FreeData if OversetFindPointsGSLIB::Setup has been called already
  if (setupflag) { FreeData(); }
  