From 33aedfd4fe1071aeb3aa92a04286d4f86bc90268 Mon Sep 17 00:00:00 2001 From: Maien Hamed Date: Mon, 5 Oct 2015 18:24:39 +0200 Subject: [PATCH] push forward one coordinate at a time for jacobian derivatives --- source/fe/mapping_q_generic.cc | 163 +++++++++++++++++++++++---------- 1 file changed, 117 insertions(+), 46 deletions(-) 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