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