From a25274e81c5eeec8f45a51e88c4e3d5e1516c9b2 Mon Sep 17 00:00:00 2001 From: Luca Heltai Date: Thu, 3 Mar 2016 20:28:10 +0100 Subject: [PATCH] Fixed Jacobian for square mesh. --- source/fe/mapping_manifold.cc | 2 +- tests/fe/mapping_manifold_03.cc | 45 ++++++++++++++--------------- tests/fe/mapping_manifold_03.output | 13 +++++++++ 3 files changed, 35 insertions(+), 25 deletions(-) diff --git a/source/fe/mapping_manifold.cc b/source/fe/mapping_manifold.cc index c017d0b4a3..3a69a9ed25 100644 --- a/source/fe/mapping_manifold.cc +++ b/source/fe/mapping_manifold.cc @@ -441,7 +441,7 @@ namespace internal Tensor<1,spacedim> T = data.cell->get_manifold().get_tangent_vector(P, NP); for (unsigned int d=0; d triangulation; - Point center; - for (unsigned int i=0; i manifold(center); - - triangulation.set_all_manifold_ids(0); - triangulation.set_manifold (0, manifold); - - const QGauss quad(3); + const QGauss quad(5); FE_Q fe(1); @@ -82,15 +69,19 @@ void test() for (unsigned int q=0; q 1e-10) + { + deallog << "Jacobian from mapping manifold at point "<< q << std::endl; + for (unsigned int d=0; d(); + test<1,2>(); + test<1,3>(); + test<2,2>(); + test<2,3>(); + test<3,3>(); return 0; diff --git a/tests/fe/mapping_manifold_03.output b/tests/fe/mapping_manifold_03.output index e69de29bb2..040c04d08b 100644 --- a/tests/fe/mapping_manifold_03.output +++ b/tests/fe/mapping_manifold_03.output @@ -0,0 +1,13 @@ + +DEAL::dim=1, spacedim=1 +DEAL::OK +DEAL::dim=1, spacedim=2 +DEAL::OK +DEAL::dim=1, spacedim=3 +DEAL::OK +DEAL::dim=2, spacedim=2 +DEAL::OK +DEAL::dim=2, spacedim=3 +DEAL::OK +DEAL::dim=3, spacedim=3 +DEAL::OK -- 2.39.5