Dear Torquil,
Your program give me exactly the same resluts in the three cases. Can you load
the last svn version of getfem. I think this problem has been already fixed.
Best regards,
Yves.
----- Original Message -----
From: "Torquil Macdonald Sørensen" <address@hidden>
To: address@hidden
Sent: Tuesday, May 5, 2015 5:35:50 PM
Subject: [Getfem-users] High-level assembly
Hi!
I'm experimenting with integration on faces using the high-level
assembly. I was not expecting that the following two expressions would
give different results when integrating on a single face in a
FEM_QK(2,1) mesh:
Case 1: "Test_u*Grad_u.Normal"
Case 2: "Test_u*Normal.Grad_u"
For all the details, please see the attached test program which shows
that these expressions give different results in the case of integrating
along face 0 in convex 0, in a 2x2 FEM_QK(2,1) mesh. I expect Case 1 to
be the correct behaviour, because its results coincide with the result
of selecting the x-component of Grad_u explicitly (since the outward
unit normal for the face in question is [1,0]):
Case 3: "Test_u*Grad_u(1)"
So unless I have misunderstood, something wrong happens in Case 2. I'm
using a very recent SVN-version of Getfem++ and the latest GCC 5 snapshot.
Best regards,
Torquil Sørensen
_______________________________________________
Getfem-users mailing list
address@hidden
https://mail.gna.org/listinfo/getfem-users