Mercurial > repos > public > sbplib

view gridDerivatives.m @ 1136:eee71789f13b feature/laplace_curvilinear_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add LaplaceCurvilinear schemes where the minimum check will be implemented. The Virta scheme will be used for comparison only.
author Martin Almquist <malmquist@stanford.edu>
date Mon, 10 Jun 2019 10:43:12 +0200
parents 48b6fb693025
children
line wrap: on
line source

% Calculates derivatives in all directions of function.
% Uses multi-dim vector form.
%   gridDerivatives(F,Dx,Dy)
%   gridDerivatives(F,Dx,Dy,Dz)
function varargout = gridDerivatives(F, varargin)
    assert(length(varargin) == ndims(F));

    switch ndims(F)
        case 2
            varargout{1} = (varargin{1}*F')';
            varargout{2} = varargin{2}*F;
        otherwise
            error('Not implemented for ndims(F) = %d',ndims(F));
    end
end