From: Maien Hamed Date: Mon, 5 Oct 2015 16:24:39 +0000 (+0200) Subject: push forward one coordinate at a time for jacobian derivatives X-Git-Tag: v8.4.0-rc2~336^2~1 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=33aedfd4fe1071aeb3aa92a04286d4f86bc90268;p=dealii.git push forward one coordinate at a time for jacobian derivatives --- diff --git a/source/fe/mapping_q_generic.cc b/source/fe/mapping_q_generic.cc index 3a1bbda100..f85c3f77ee 100644 --- a/source/fe/mapping_q_generic.cc +++ b/source/fe/mapping_q_generic.cc @@ -2263,6 +2263,7 @@ namespace internal if (cell_similarity != CellSimilarity::translation) { + double tmp[spacedim][spacedim][spacedim]; for (unsigned int point=0; point *second = @@ -2282,23 +2283,34 @@ namespace internal * data.mapping_support_points[k][i]); - // pushing forward the derivative coordinates + // first push forward the j-components for (unsigned int i=0; i *third = @@ -2397,27 +2410,49 @@ namespace internal * data.mapping_support_points[k][i]); + // push forward the j-coordinate + for (unsigned int i=0; i *fourth = @@ -2521,6 +2557,53 @@ namespace internal * data.mapping_support_points[k][i]); + // push-forward the j-coordinate + for (unsigned int i=0; i