Mercurial > repos > public > sbplib

annotate convergence.m @ 255:df3cc9c5dffc operator_remake

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Added ordinary 12th order accurate, with D1*D1 as 2nd derivative.
author Martin Almquist <martin.almquist@it.uu.se>
date Wed, 07 Sep 2016 15:54:41 +0200
parents ad1a41d135a1
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
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