[Getfem-commits] r4850 - /trunk/getfem/src/getfem_models.cc

[logari81]

Author: logari81
Date: Tue Feb  3 09:48:29 2015
New Revision: 4850

URL: http://svn.gna.org/viewcvs/getfem?rev=4850&view=rev
Log:
fix bug in Dirichlet condition brick assembling

Modified:
    trunk/getfem/src/getfem_models.cc

Modified: trunk/getfem/src/getfem_models.cc
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/src/getfem_models.cc?rev=4850&r1=4849&r2=4850&view=diff
==============================================================================
--- trunk/getfem/src/getfem_models.cc   (original)
+++ trunk/getfem/src/getfem_models.cc   Tue Feb  3 09:48:29 2015
@@ -3941,9 +3941,9 @@
             asm_real_or_complex_1_param
               (*B, mim, mf_mult, *mf_H, *H, rg, (mf_u.get_qdim() == 1) ?
                "F=data(#2);"
-               "M(#1,#3)+=comp(Base(#1).Base(#3).Base(#2)(:,:,i).F(i))"
+               "M(#1,#3)+=comp(Base(#1).Base(#3).Base(#2))(:,:,i).F(i)"
                : "F=data(qdim(#1),qdim(#1),#2);"
-               
"M(#1,#3)+=comp(vBase(#1).vBase(#3).Base(#2)(:,i,:,j,k).F(i,j,k));", &mf_u);
+               
"M(#1,#3)+=comp(vBase(#1).vBase(#3).Base(#2))(:,i,:,j,k).F(i,j,k);", &mf_u);
           else {
             asm_real_or_complex_1_param
               (*B, mim, mf_mult, mf_u, *H, rg, (mf_u.get_qdim() == 1) ?