From 8f7969e7b64aa034c4ed96b00ede7060e79cc775 Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Thu, 1 Jun 2017 21:15:57 +0200 Subject: [PATCH] Fix bug introduced by transfinite manifold commit. --- source/fe/mapping_q_generic.cc | 5 +++++ tests/manifold/transfinite_manifold_03.output | 18 +++++++++--------- 2 files changed, 14 insertions(+), 9 deletions(-) diff --git a/source/fe/mapping_q_generic.cc b/source/fe/mapping_q_generic.cc index 0320705483..557d0a045d 100644 --- a/source/fe/mapping_q_generic.cc +++ b/source/fe/mapping_q_generic.cc @@ -3809,6 +3809,11 @@ compute_mapping_support_points(const typename Triangulation::cell_ all_manifold_ids_are_equal = false; } + if (all_manifold_ids_are_equal) + { + std::vector > vertices(a); + cell->get_manifold().add_new_points(vertices, support_point_weights_cell, a); + } else switch (dim) { diff --git a/tests/manifold/transfinite_manifold_03.output b/tests/manifold/transfinite_manifold_03.output index afd60ee953..276ab46d87 100644 --- a/tests/manifold/transfinite_manifold_03.output +++ b/tests/manifold/transfinite_manifold_03.output @@ -7,19 +7,19 @@ DEAL::80 cells, degree=4: 3.14159265350 DEAL:: DEAL::Testing with SphericalManifold dim=2, spacedim=2 DEAL::80 cells, degree=1: 3.06146745892 -DEAL::80 cells, degree=2: 3.14143771670 -DEAL::80 cells, degree=3: 3.14158706712 -DEAL::80 cells, degree=4: 3.14159263998 +DEAL::80 cells, degree=2: 3.14064643171 +DEAL::80 cells, degree=3: 3.14158597871 +DEAL::80 cells, degree=4: 3.14159257389 DEAL:: DEAL::Testing with SphericalManifold dim=3, spacedim=3 DEAL::56 cells, degree=1: 3.21034343105 -DEAL::56 cells, degree=2: 4.18080434784 -DEAL::56 cells, degree=3: 4.18759484733 -DEAL::56 cells, degree=4: 4.18877458327 +DEAL::56 cells, degree=2: 4.14326491350 +DEAL::56 cells, degree=3: 4.18822553211 +DEAL::56 cells, degree=4: 4.18868242879 DEAL:: DEAL::Testing with CylindricalManifold in 3d DEAL::80 cells, degree=1: 5.65685424949 -DEAL::80 cells, degree=2: 6.27829514062 -DEAL::80 cells, degree=3: 6.28283288022 -DEAL::80 cells, degree=4: 6.28318111226 +DEAL::80 cells, degree=2: 6.25775643968 +DEAL::80 cells, degree=3: 6.28283082764 +DEAL::80 cells, degree=4: 6.28312774156 DEAL:: -- 2.39.5