Hi Jean-Paul,
Thanks for your response. I do think am doing the scaling of the area
correctly, which I based on your answer to my previous question about the
normal vectors. I just did so in a way that is not so clear in the code. To
make it more clear, I have written the code to be
Tensor<1,
Hi Alex,
Yes, at first glance it does look to me like the integration is a problem. You
are using the same FEFaceValues object to perform the integration in both the
material and spatial configurations, and this is not quite correct with the way
that you have things set up. The following
