[Getfem-commits] r5452 - /trunk/getfem/src/getfem_generic_assembly.cc

[Yves . Renard]

Author: renard
Date: Fri Oct 28 13:51:20 2016
New Revision: 5452

URL: http://svn.gna.org/viewcvs/getfem?rev=5452&view=rev
Log:
small fix

Modified:
    trunk/getfem/src/getfem_generic_assembly.cc

Modified: trunk/getfem/src/getfem_generic_assembly.cc
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/src/getfem_generic_assembly.cc?rev=5452&r1=5451&r2=5452&view=diff
==============================================================================
--- trunk/getfem/src/getfem_generic_assembly.cc (original)
+++ trunk/getfem/src/getfem_generic_assembly.cc Fri Oct 28 13:51:20 2016
@@ -12091,8 +12091,8 @@
           // iterations on interpolation points
           gis.nbpt = pspt->size();
           for (size_type ii = 0; ii < ind.size(); ++ii) {
-            gis.ipt = ind[ii];
-            if (gis.ctx.have_pgp()) gis.ctx.set_ii(gis.ipt);
+            gis.ipt = ii;
+            if (gis.ctx.have_pgp()) gis.ctx.set_ii(ind[ii]);
             else gis.ctx.set_xref((*pspt)[gis.ipt]);
 
             if (ii == 0 || !(pgt->is_linear())) {
@@ -12111,7 +12111,7 @@
             if (ii == 0)
               for (size_type j = 0; j < gilb.size(); ++j) j += gilb[j]->exec();
             for (size_type j = 0; j < gil.size(); ++j) j += gil[j]->exec();
-            gic.store_result(v.cv(), gis.ipt, workspace.assembled_tensor());
+            gic.store_result(v.cv(), ind[ii], workspace.assembled_tensor());
           }
         }
       }