[Getfem-commits] [getfem-commits] branch devel-logari81 updated: Fix wrong point size in mesh method

[Konstantinos Poulios via Getfem-commits]

This is an automated email from the git hooks/post-receive script.

logari81 pushed a commit to branch devel-logari81
in repository getfem.

The following commit(s) were added to refs/heads/devel-logari81 by this push:
     new 60d7d962 Fix wrong point size in mesh method
60d7d962 is described below

commit 60d7d962e77f6d582ac4d01c3897212f687ea7bd
Author: Konstantinos Poulios <logar...@gmail.com>
AuthorDate: Fri Jun 24 14:45:06 2022 +0200

    Fix wrong point size in mesh method
---
 src/getfem_mesh.cc | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/src/getfem_mesh.cc b/src/getfem_mesh.cc
index 4c631a38..d6e25378 100644
--- a/src/getfem_mesh.cc
+++ b/src/getfem_mesh.cc
@@ -1150,7 +1150,7 @@ namespace getfem {
       pgp = bgeot::geotrans_precomp(pgt, pspt, 0);
     }
 
-    base_node pt(n);
+    base_node pt(dim());
     ipt.resize(pspt->size());
     for (size_type ip = 0; ip < pspt->size(); ++ip) {
       pgp->transform(points_of_convex(i), ip, pt);
@@ -1229,6 +1229,7 @@ namespace getfem {
       else return icc;
     }
     else if (ref) Bank_refine_normal_convex(i);
+
     return size_type(-1);
   }
 
@@ -1283,7 +1284,7 @@ namespace getfem {
     bgeot::pgeotrans_precomp pgp = bgeot::geotrans_precomp(pgt1, pspt1, 0);
 
     std::vector<size_type> ipt1(pspt1->size());
-    base_node pt(n);
+    base_node pt(dim());
     for (size_type i = 0; i < pspt1->size(); ++i) {
       pgp->transform(mesh2.points_of_convex(ic1), i, pt);
       ipt1[i] = mesh2.add_point(pt);