view +parametrization/old/triang_map.m @ 1065:c2bd7f15da48 feature/laplace_curvilinear_test

Add scheme where only corner points are multiplied by dim. Does not seem to decrease spectral radius much, but is stable.
author Martin Almquist <malmquist@stanford.edu>
date Thu, 17 Jan 2019 18:57:54 -0800
parents 3a3cf386bb7e
children
line wrap: on
line source

% Creates a grid [X,Y] from the mapping function S at points in vectors u,v
function [X, Y] = traing_map(S,u,v)
    error('not done')
    if nargin == 2
        v = u;
    end

    if isscalar(u)
        u = linspace(0,1,u);
    end

    if isscalar(v)
        v = linspace(0,1,v);
    end

    nu = length(u);
    nv = length(v);

    X = zeros(nu,nv);
    Y = zeros(nu,nv);

    for i = 1:nu
        for j = 1:nv
            [x,y] = S(u(i),v(j));
            X(i,j) = x;
            Y(i,j) = y;
        end
    end
end