Mercurial > repos > public > sbplib

view +util/get_grid.m @ 741:5a9acf282b34 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add scheme Heat2Dcurvilinear. Neumann and Dirichlet seem to work. Only tested for stretched Cartesian grids though.
author Martin Almquist <malmquist@stanford.edu>
date Wed, 09 May 2018 19:29:12 -0700
parents 7067bf8adbfa
children
line wrap: on
line source

% [x,h] = get_grid(a,b,m)
function [x,h] = get_grid(a,b,m)
    % TODO: allow the interval to be a vector
    x = linspace(a,b,m)';
    h = (b-a)/(m-1);
end