Mercurial > repos > public > sbplib
view gridDerivatives.m @ 1065:c2bd7f15da48 feature/laplace_curvilinear_test
Add scheme where only corner points are multiplied by dim. Does not seem to decrease spectral radius much, but is stable.
author | Martin Almquist <malmquist@stanford.edu> |
---|---|
date | Thu, 17 Jan 2019 18:57:54 -0800 |
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