Mercurial > repos > public > sbplib
view +parametrization/old/triang_map.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 | 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