Mercurial > repos > public > sbplib
annotate +scheme/error1d.m @ 774:66eb4a2bbb72 feature/grids
Remove default scaling of the system.
The scaling doens't seem to help actual solutions. One example that fails in the flexural code.
With large timesteps the solutions seems to blow up. One particular example is profilePresentation
on the tdb_presentation_figures branch with k = 0.0005
author | Jonatan Werpers <jonatan@werpers.com> |
---|---|
date | Wed, 18 Jul 2018 15:42:52 -0700 |
parents | ed6a704b028d |
children |
rev | line source |
---|---|
26
ed6a704b028d
Made some changes to error functions and comparison functions before finalizing convergence script.
Jonatan Werpers <jonatan@werpers.com>
parents:
15
diff
changeset
|
1 function e = error1d(discr, v1, v2) |
ed6a704b028d
Made some changes to error functions and comparison functions before finalizing convergence script.
Jonatan Werpers <jonatan@werpers.com>
parents:
15
diff
changeset
|
2 h = discr.h; |
15
16bad7c459da
Added abstract methods to allow for comparing of solutions. Added a bunch of error functions.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
3 e = sqrt(h*sum((v1-v2).^2)); |
16bad7c459da
Added abstract methods to allow for comparing of solutions. Added a bunch of error functions.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
4 end |