view +parametrization/old/triang_map.m @ 1314:58df4a35fe43 feature/poroelastic

Add traction operators that include the viscous terms and use them to clean up the imposition of traction BC.
author Martin Almquist <malmquist@stanford.edu>
date Sat, 25 Jul 2020 15:56:31 -0700
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