From: Luca Heltai Date: Thu, 3 Mar 2016 19:28:10 +0000 (+0100) Subject: Fixed Jacobian for square mesh. X-Git-Tag: v8.5.0-rc1~1131^2~18 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=a25274e81c5eeec8f45a51e88c4e3d5e1516c9b2;p=dealii.git Fixed Jacobian for square mesh. --- 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