Mercurial > repos > public > sbplib

view +parametrization/old/triang_interp_pts.m @ 845:1e057b0f2fed feature/burgers1d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add RK6 with residual viscosity update and reduce computational effort of spatial scheme - Add RK6 with residual updates - Change the D2 operator for upwind schemes to one less computationally expensive.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Wed, 19 Sep 2018 16:32:05 +0200
parents 81e0ead29431
children
line wrap: on
line source

% Creates a transfinite interpolation from connecting the four points wiht straight lines.
function [S, g1, g2, g3] = triang_interp_pts(p1,p2,p3)
    if size(p1) ~= [2 1]
        error('p1 is strange!');
    end

    g1 = @(t)(p1 + t*(p2-p1));
    g2 = @(t)(p2 + t*(p3-p2));
    g3 = @(t)(p3 + t*(p1-p3));

    S = parametrization.triang_interp(g1,g2,g3);
end