© 2004
Canonical Ltd.
•
Terms of use
•
Data privacy
•
Contact Launchpad Support
•
Blog
•
Careers
•
System status
•
0e1f616
(Get the code!)
The printouts from your expression, show that the expression is evaluated at the midpoints of each edge for each cell, e.g. these lines for one particular cell:
x: 1.000000, y: 0.250000, n.x: 0.000000, n.y: -1.000000
x: 0.750000, y: 0.250000, n.x: 0.000000, n.y: -1.000000
x: 0.750000, y: 0.000000, n.x: 0.000000, n.y: -1.000000
What happens is that:
1) Your expression is interpolated into the local fem basis on a cell, which involves the above evaluations on all edges of the cell.
2) Integration over the facet is performed, using the dofs on the cell computed in step 1.
So the prints look fine, and I'm closing this bug as invalid. If you get wrong results in the integration, you need to justify that otherwise and reopen the bug.