Mercurial > repos > public > sbplib

view convergence.m @ 679:247b58a4dbe8 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add support for Dirichlet and Traction BC on different components at the same boundary. Remove some unused variables and improve comments.
author Martin Almquist <malmquist@stanford.edu>
date Mon, 05 Feb 2018 14:45:26 -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