How do I apply the L norms?
I know this is a real newbie question.
I want to apply the L2-norm for a laplace equation. How should I proceed? Should I simply calculate the difference between the values of the current and the last iteration, then the sum of the squares of these values and finally the square root? Is this norm used to determine when to stop iterating? Thanks a lot! |
|
Thanks, that's very interesting!
If I still want to compare the solutions from the latest two iterations, do I have to treat the domensions sepparately or not? Obviusly, by applying a L2-norm for a 3D domain will naturally provide higher values than an L2-norm for a 1D domain (because will be more points in the grid). Further can you recommend me a good documentation on applying L-norms in CFD? |
If you divide by norm of rhs as I wrote, number of grid points does not matter.
|
I'm implementing the basic laplace equation, i.e. the RHS is zero... How should I proceed?
Thanks! |
If the right hand side is zero, the problem is pretty boring, isn't it? No boundary conditions, no source term?
|
Ohh...now I understand, my mistake. The RHS also includes the boundary conditions (which are specified in my problem)...
|
Hi
Hi
You can use the following norm for stopping your iterations. It just calculates the relative error in the iterations total_error=sqrt((\sun_{j=1}^N(\frac{u^(n+1) - u^(n)}{u^n})^2)_j) which corresponds to all the points in the grid. You can compare this to an epsilon value and if the condition is met the iteration stops. n in the above formula represents the iteration count. |
All times are GMT -4. The time now is 01:49. |