message: remove interior face integrators from bilinearforms where they do not exist;
diff: @@ -126,7 +126,7 @@ void MFBilinearFormExtension::FormLinearSystem(const Array<int> &ess_tdof_list,
  void MFBilinearFormExtension::Mult(const Vector &x, Vector &y) const
  {
  - MFEM_VERIFY(interior_face_integs.Size() == 0,
  + MFEM_VERIFY(a->GetFBFI()->Size() == 0,
  "the case of interior face integrators is not"
  " implemented");
  @@ -188,7 +188,7 @@ void MFBilinearFormExtension::Mult(const Vector &x, Vector &y) const
  void MFBilinearFormExtension::MultTranspose(const Vector &x, Vector &y) const
  {
  - MFEM_VERIFY(interior_face_integs.Size() == 0,
  + MFEM_VERIFY(a->GetFBFI()->Size() == 0,
  "the case of interior face integrators is not"
  " implemented");
  @@ -1213,9 +1213,6 @@ void PAMixedBilinearFormExtension::SetupMultInputs(
  void PAMixedBilinearFormExtension::Mult(const Vector &x, Vector &y) const
  {
  - MFEM_VERIFY(interior_face_integs.Size() == 0,
  - "the case of interior face integrators is not"
  - " implemented");
  y = 0.0;
  AddMult(x, y);
  }
  @@ -1223,9 +1220,6 @@ void PAMixedBilinearFormExtension::Mult(const Vector &x, Vector &y) const
  void PAMixedBilinearFormExtension::AddMult(const Vector &x, Vector &y,
  const double c) const
  {
  - MFEM_VERIFY(interior_face_integs.Size() == 0,
  - "the case of interior face integrators is not"
  - " implemented");
  Array<BilinearFormIntegrator*> &integrators = *a->GetDBFI();
  const int iSz = integrators.Size();
  @@ -1251,9 +1245,6 @@ void PAMixedBilinearFormExtension::AddMult(const Vector &x, Vector &y,
  void PAMixedBilinearFormExtension::MultTranspose(const Vector &x,
  Vector &y) const
  {
  - MFEM_VERIFY(interior_face_integs.Size() == 0,
  - "the case of interior face integrators is not"
  - " implemented");
  y = 0.0;
  AddMultTranspose(x, y);
  }
  @@ -1261,9 +1252,6 @@ void PAMixedBilinearFormExtension::MultTranspose(const Vector &x,
  void PAMixedBilinearFormExtension::AddMultTranspose(const Vector &x, Vector &y,
  const double c) const
  {
  - MFEM_VERIFY(interior_face_integs.Size() == 0,
  - "the case of interior face integrators is not"
  - " implemented");
  Array<BilinearFormIntegrator*> &integrators = *a->GetDBFI();
  const int iSz = integrators.Size();
  