[Getfem-commits] r4644 - /trunk/getfem/src/getfem/getfem_deformable_mesh.h

[andriy . andreykiv]

Author: andrico
Date: Mon May 12 16:45:33 2014
New Revision: 4644

URL: http://svn.gna.org/viewcvs/getfem?rev=4644&view=rev
Log:
extending the size of the vector to basic_dof in deformable_mesh::deform_mesh

Modified:
    trunk/getfem/src/getfem/getfem_deformable_mesh.h

Modified: trunk/getfem/src/getfem/getfem_deformable_mesh.h
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/src/getfem/getfem_deformable_mesh.h?rev=4644&r1=4643&r2=4644&view=diff
==============================================================================
--- trunk/getfem/src/getfem/getfem_deformable_mesh.h    (original)
+++ trunk/getfem/src/getfem/getfem_deformable_mesh.h    Mon May 12 16:45:33 2014
@@ -79,6 +79,8 @@
       //this vector will track if a point can be deformed
       std::vector<bool> deform_pt_flag(ppts.size(), true);
       size_type cv;
+      base_vector dU_basic(mf.nb_basic_dof());
+      mf.extend_vector(dU,dU_basic);
       for(cv << conv_indices; 
         cv!=bgeot::size_type(-1); cv << conv_indices) 
       {
@@ -92,13 +94,14 @@
           "mesh_fem should be isoparametric to the mesh, "
           "with nb_points() of convex * dim == size of 
ind_basic_dof_of_element");
 
+
         for(size_type pt = 0; pt < num_points; ++pt) 
         { 
           /** iterate through each components of point [pt]and deform the 
component*/
           if(deform_pt_flag[pt_index[pt]])
             for (size_type comp = 0; comp < ddim; ++comp)
               //move pts by dU;
-                ppts[pt_index[pt]][comp] += dU[dof[pt*ddim + comp]];
+                ppts[pt_index[pt]][comp] += dU_basic[dof[pt*ddim + comp]];
 
           //flag current [pt] to deformed
           deform_pt_flag[pt_index[pt]] = false;