Mercurial > repos > public > sbplib

view +util/semidiscrete2function.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 48b6fb693025
children
line wrap: on
line source

function F = semidiscrete2function(D,S)
    % TODO: Find threshhold when sparse makes sense.
    % TODO: Why does having the spase block slow down convergence.m (or does it?)
    if size(D,1) > 500
        D = sparse(D);
    end

    if ~isa(S,'function_handle')
        F = @(w,t)(D*w + S);
    else
        F = @(w,t)(D*w + S(t));
    end
end