Mercurial > repos > public > sbplib

view +noname/calcSol.m @ 1228:ef08adea56c4 feature/volcano

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Utilize grid module in Hypsyst3d. Pass cartesian grid and replace some of the previous functionality in Hypsyst3d with properties/methods for cartesian grids.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Mon, 18 Nov 2019 10:31:58 -0800
parents 127e562bd6d3
children
line wrap: on
line source

% Calculates the solution of a discr at a given time using aligned timesteps.
% Returns the solution as a grid function as defined in +grid
function gf = calculateSolution(discr, T, tsOpt)
    k_max = discr.getTimestep(tsOpt);
    [k,N] = alignedTimestep(k_max,T);
    tsOpt.k = k;
    ts = discr.getTimestepper(tsOpt);

    gf = ts.stepN(N-ts.n);
end