view subsSymfun.m @ 1270:9a7cac4202d0 feature/poroelastic

In CruvilinearAnisotropicUpwind: Make D1Upwind default opSet. Use D2Standard for metric derivatives. Bugfix global quadrature.
author Martin Almquist <malmquist@stanford.edu>
date Sun, 31 May 2020 20:18:45 -0700
parents 489bea7fc33f
children
line wrap: on
line source

% Subs for a symfun
% f remains a symbolic function. If any of it's arguments is eliminated
% it is removed from the argument list while preserving the order of the
% other arguments
function f = subsSymfun(f, old, new)
    args = argnames(f);

    newExpr = subs(f, old, new);
    vars = symvar(subs(args, old, new));

    newArgs = args(ismember(args,vars));

    f = symfun(newExpr, newArgs);
end