On 10/16/24 09:53, Davit Gyulamiryan wrote:
Error on the initial grid is ~1.6e-3
Error on the second grid is ~9.2
Error on the third grid is ~11.5
My function:
1.578e-03
3.956e-04
1.021e-04
My implementation is rather long and is not really relevant to the
question itself. It seems to be working so can just use it instead of
the built-in functions, but I was interested if I have a bug somewhere
in the code.
Has anyone ever had this or similar problem, or know how to use the
built-in funcionality of deal.ii for my case?
It's of course hard to tell without seeing the complete code, but just
in case here's a guess:
* On the cells where you use the FE_Nothing, your finite element
solution is zero.
* On those same cells, your exact solution is probably *not* zero.
* VectorTools::integrate_difference() integrates over *all* cells.
* As a consequence, on the cells with FE_Nothing, you don't compute
|| u_h - u ||
but instead
|| -u ||
and that would probably not be what you want.
Best
W.
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/9cfe254f-3bcc-4cc8-8452-30b37c68013a%40colostate.edu.
