[Getfem-commits] (no subject)

[Markus Bürg]

branch: mb-transInversion
commit 1872d5be4b6e2d232673852549e08f9b3f9f8a3e
Author: mb <address@hidden>
Date:   Tue Sep 5 16:10:54 2017 +0200

    Fix sizes for lower-dimensional elements.
---
 src/bgeot_geotrans_inv.cc       | 5 +++--
 src/getfem/bgeot_geotrans_inv.h | 4 ++--
 2 files changed, 5 insertions(+), 4 deletions(-)

diff --git a/src/bgeot_geotrans_inv.cc b/src/bgeot_geotrans_inv.cc
index 1ce0de3..feac928 100644
--- a/src/bgeot_geotrans_inv.cc
+++ b/src/bgeot_geotrans_inv.cc
@@ -125,10 +125,11 @@ 
nonlinear_storage::linearised_structure::linearised_structure(
     direct_points_ref.push_back(reference_nodes[i]);
   }
 
-  auto N = direct_points_ref.begin()->size();
+  auto N = direct_points.begin()->size();
+  auto N_ref = direct_points_ref.begin()->size();
   base_matrix K_linear(N, n_points - 1);
   B_linear.base_resize(N, n_points - 1);
-  K_ref_linear.base_resize(N, n_points - 1);
+  K_ref_linear.base_resize(N_ref, n_points - 1);
   P_linear = direct_points[0];
   P_ref_linear = direct_points_ref[0];
 
diff --git a/src/getfem/bgeot_geotrans_inv.h b/src/getfem/bgeot_geotrans_inv.h
index 10db1d3..881b978 100644
--- a/src/getfem/bgeot_geotrans_inv.h
+++ b/src/getfem/bgeot_geotrans_inv.h
@@ -190,8 +190,8 @@ namespace bgeot {
       // computation of the pseudo inverse
       update_B();
     } else {
-      nonlinear_storage.diff.resize(P);
-      nonlinear_storage.x_real.resize(P);
+      nonlinear_storage.diff.resize(N);
+      nonlinear_storage.x_real.resize(N);
       nonlinear_storage.x_ref.resize(P);
 
       if (pgt->complexity() > 1) {