Mercurial > repos > public > sbplib
annotate convergence.m @ 971:e54c2f54dbfe feature/getBoundaryOperator
Merge with feature/poroelastic. Use only the changes made to multiblock.DiffOp and scheme.Elastic2dVariable. DiffOp.getBoundaryOperator/Quadrature now use scheme methods instead of propeties.
author | Martin Almquist <malmquist@stanford.edu> |
---|---|
date | Tue, 25 Dec 2018 07:50:07 +0100 |
parents | ad1a41d135a1 |
children |
rev | line source |
---|---|
23
ad1a41d135a1
Added function to calculate convergence given errors. Added function to print convergence table (not finished.)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
1 function q = convergence(e,h) |
ad1a41d135a1
Added function to calculate convergence given errors. Added function to print convergence table (not finished.)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
2 for i = 1:length(e)-1 |
ad1a41d135a1
Added function to calculate convergence given errors. Added function to print convergence table (not finished.)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
3 q(i) = log(e(i)/e(i+1))/log(h(i+1)/h(i)); |
ad1a41d135a1
Added function to calculate convergence given errors. Added function to print convergence table (not finished.)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
4 end |
ad1a41d135a1
Added function to calculate convergence given errors. Added function to print convergence table (not finished.)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
5 end |