[Getfem-commits] (no subject)

[Markus Bürg]

branch: mb-Use_rtree_in_poly_composite
commit 22cd69d3739ee1238782cd0176487a5411e950c6
Author: mb <address@hidden>
Date:   Thu May 23 16:23:45 2019 +0200

    Create matrix only, if really needed.
---
 src/getfem/getfem_interpolation.h | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/src/getfem/getfem_interpolation.h 
b/src/getfem/getfem_interpolation.h
index 8dadccd..9a1aaf5 100644
--- a/src/getfem/getfem_interpolation.h
+++ b/src/getfem/getfem_interpolation.h
@@ -421,8 +421,8 @@ namespace getfem {
     dim_type qqdim = dim_type(gmm::vect_size(UU)/mf_source.nb_dof());
 
     std::vector<T> U(mf_source.nb_basic_dof()*qqdim);
-    gmm::row_matrix<gmm::rsvector<scalar_type> >
-      M(gmm::mat_nrows(MM), mf_source.nb_basic_dof());
+    gmm::row_matrix<gmm::rsvector<scalar_type> > M;
+    if (version != 0) M.resize(gmm::mat_nrows(MM), mf_source.nb_basic_dof());
 
     if (version == 0) mf_source.extend_vector(UU, U);
 
@@ -565,8 +565,8 @@ namespace getfem {
     size_type qqdimt = qqdim * mf_source.get_qdim()/mf_target.get_qdim();
     std::vector<T> V(mf_target.nb_basic_dof()*qqdimt);
     mf_target.extend_vector(VV,V);
-    gmm::row_matrix<gmm::rsvector<scalar_type> >
-      M(mf_target.nb_basic_dof(), mf_source.nb_dof());
+    gmm::row_matrix<gmm::rsvector<scalar_type> > M;
+    if (version != 0) M.resize(mf_target.nb_basic_dof(), mf_source.nb_dof());
 
     const mesh &msh(mf_source.linked_mesh());
     getfem::mesh_trans_inv mti(msh, EPS);