From d108ce8311c197c15bf57ea6d31bf7a27d25e805 Mon Sep 17 00:00:00 2001 From: Maien Hamed Date: Mon, 5 Oct 2015 22:17:06 +0200 Subject: [PATCH] push forward one coordinate at a time in MappingFEField jacobian pushed-forward derivatives --- source/fe/mapping_fe_field.cc | 168 +++++++++++++++++++++++---------- source/fe/mapping_q_generic.cc | 40 ++++---- 2 files changed, 139 insertions(+), 69 deletions(-) diff --git a/source/fe/mapping_fe_field.cc b/source/fe/mapping_fe_field.cc index 4bf8092abe..bc6157a737 100644 --- a/source/fe/mapping_fe_field.cc +++ b/source/fe/mapping_fe_field.cc @@ -685,6 +685,7 @@ namespace internal if (cell_similarity != CellSimilarity::translation) { + double tmp[spacedim][spacedim][spacedim]; for (unsigned int point=0; point *second = @@ -700,25 +701,36 @@ namespace internal for (unsigned int l=0; l *third = @@ -819,28 +832,49 @@ namespace internal * data.local_dof_values[k]); } - // pushing forward the derivative coordinates + // push forward the j-coordinate + for (unsigned int i=0; i *fourth = @@ -948,31 +983,66 @@ namespace internal * data.local_dof_values[k]); } - // pushing forward the derivative coordinates + // push-forward the j-coordinate + for (unsigned int i=0; i