[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Getfem-commits] (no subject)
From: |
Markus Bürg |
Subject: |
[Getfem-commits] (no subject) |
Date: |
Tue, 5 Sep 2017 10:10:58 -0400 (EDT) |
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) {