Template class ImmersedSurfaceQuadrature<0, 1>

diff --git a/source/non_matching/immersed_surface_quadrature.cc b/source/non_matching/immersed_surface_quadrature.cc
index b5529864afe00153e650f567941065a4c7afa74f..85f1c882f5e053fbea24d5a1e157405e3332bcea 100644 (file)
--- a/source/non_matching/immersed_surface_quadrature.cc
+++ b/source/non_matching/immersed_surface_quadrature.cc
@@ -77,6 +77,7 @@ namespace NonMatching
   template class ImmersedSurfaceQuadrature<1, 1>;
   template class ImmersedSurfaceQuadrature<2, 2>;
   template class ImmersedSurfaceQuadrature<3, 3>;
+  template class ImmersedSurfaceQuadrature<0, 1>;
   template class ImmersedSurfaceQuadrature<1, 2>;
   template class ImmersedSurfaceQuadrature<2, 3>;
 

diff --git a/tests/non_matching/immersed_surface_quadrature.cc b/tests/non_matching/immersed_surface_quadrature.cc
index 7c6e8a9311a1550d22d9d41c9ec5982287b5dc9c..c3f59e56f6aa9e3c9e624820b1e0ec6a50be5097 100644 (file)
--- a/tests/non_matching/immersed_surface_quadrature.cc
+++ b/tests/non_matching/immersed_surface_quadrature.cc
@@ -33,7 +33,9 @@ print_quadrature(
 {
   for (unsigned int i = 0; i < quadrature.size(); ++i)
     {
-      deallog << quadrature.point(i) << ", " << quadrature.weight(i) << ", "
+      if (dim > 0) // Can't print 0-dim points.
+        deallog << quadrature.point(i);
+      deallog << ", " << quadrature.weight(i) << ", "
               << quadrature.normal_vector(i) << std::endl;
     }
 }
@@ -108,6 +110,7 @@ main()
   construct_quadrature_and_print_points<3, 3>();
 
   // Face quadrature
+  construct_quadrature_and_print_points<0, 1>();
   construct_quadrature_and_print_points<1, 2>();
   construct_quadrature_and_print_points<2, 3>();
 }

diff --git a/tests/non_matching/immersed_surface_quadrature.output b/tests/non_matching/immersed_surface_quadrature.output
index efb667c09afb9fec94987b5ce09b8a17742d6657..f1d0c12ac90f0416b31902a189013e136274413b 100644 (file)
--- a/tests/non_matching/immersed_surface_quadrature.output
+++ b/tests/non_matching/immersed_surface_quadrature.output
@@ -12,6 +12,10 @@ DEAL::0.00000 0.00000 0.00000, 1.00000, 0.00000 0.00000 1.00000
 DEAL::Using constructor
 DEAL::0.00000 0.00000 0.00000, 1.00000, 0.00000 0.00000 1.00000
 DEAL::Using push_back
+DEAL::, 1.00000, 1.00000
+DEAL::Using constructor
+DEAL::, 1.00000, 1.00000
+DEAL::Using push_back
 DEAL::0.00000, 1.00000, 0.00000 1.00000
 DEAL::Using constructor
 DEAL::0.00000, 1.00000, 0.00000 1.00000