I know that we can get the residual if use 'sweep' to solve the equation, but what is the formula for computing residual in fipy? Is it the difference between the two sides of the equation? And why the residual of my equation decrease very slowly? Can anyone help me? Thanks !
Asked
Active
Viewed 45 times
1 Answers
0
The residual returned by sweep() is the L2-norm of L*x - b (where L is the matrix, x is the solution vector, and b is the right-hand-side vector for the system of linear equations obtained from discretizing your equation(s)). I have filed issue #734 to add documentation of this.
As to why the residual is decreasing very slowly, that's impossible to say without knowing what equations you are solving and how you are solving them.
jeguyer
- 2,379
- 1
- 11
- 15
-
Thank you from the bottom of my heart. Wish you have a nice day. – Xiukai Wang Jun 19 '20 at 00:31