From: Simon Sticko Date: Wed, 6 Apr 2022 16:09:36 +0000 (+0200) Subject: Template class ImmersedSurfaceQuadrature<0, 1> X-Git-Tag: v9.4.0-rc1~327^2~1 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=6327c8cb9e2bb9bb82e4e309a20f4c7ce2f69a3a;p=dealii.git Template class ImmersedSurfaceQuadrature<0, 1> --- diff --git a/source/non_matching/immersed_surface_quadrature.cc b/source/non_matching/immersed_surface_quadrature.cc index b5529864af..85f1c882f5 100644 --- 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 7c6e8a9311..c3f59e56f6 100644 --- 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 efb667c09a..f1d0c12ac9 100644 --- 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