On 4/7/20 4:53 PM, 孙翔 wrote:
I am still unclear about what is happening. That zero is just a typo. Even if
I remove the zero, these two methods can not get the same solution. What I did
is just using two methods to calculate the right-hand side generated by time
discretization. One way is to form a mass matrix then time the old solution,
the other way is to form the local right-hand side in element level then
assemble it into system right-hand side array.
Right, these two approaches should result in the same solution. So the
question is: how do they differ? Do they result in different right hand side
vectors, or in different solutions? If they differ in the solutions, once you
visualize them, is it obvious that the problem is at the boundary? At hanging
nodes? These are useful tools to figure out where the bug may be.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.colostate.edu/~bangerth/
--
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/4f4e2b7f-11bc-7d89-69d4-938d672b20f6%40colostate.edu.
