Mercurial > repos > public > sbplib

view convergence.m @ 1083:a22b23098021 feature/laplace1dUpdate

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Divide interface penalty strength by 4, to make it sharp. Stability verified by eigenvalue computations.
author Martin Almquist <malmquist@stanford.edu>
date Tue, 05 Mar 2019 12:41:50 -0800
parents ad1a41d135a1
children
line wrap: on
line source

function q = convergence(e,h)
    for i = 1:length(e)-1
        q(i) = log(e(i)/e(i+1))/log(h(i+1)/h(i));
    end
end