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) ?
_______________________________________________
Getfem-commits mailing list
[email protected]
https://mail.gna.org/listinfo/getfem-commits
