view sortd.m @ 851:ab2e5a24ddde feature/burgers1d

- Fix bug when constructing closure for narrow stencils - Update the residual outside of the RK time steps. At least for the inner convergence rate, updating the residual inside does not seem to be required.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 21 Sep 2018 15:33:15 +0200
parents 48b6fb693025
children
line wrap: on
line source

% Sorts x with respect to the function d
% d is a function that accepts an element of x and returns a scalar.
function y = sortd(d,x)
    v = [];
    for i = 1:length(x)
        v(i) = d(x(i));
    end

    [~,I] = sort(v);
    y = x(I);
end