Mercurial > repos > public > sbplib

view +util/get_convergence.m @ 678:06676c40e77f feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add scheme for complete elastic operator. Traction BC working with MMS variable coefficient.
author Martin Almquist <malmquist@stanford.edu>
date Mon, 05 Feb 2018 11:06:15 -0800
parents 48b6fb693025
children
line wrap: on
line source

function q = get_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