Mercurial > repos > public > sbplib

view subsSymfun.m @ 1118:07d0caf915e4 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Introduce optFlag so that one can choose not to build the optimization-related cell of matrices called B. It is too computationally costly and should probably be done in a different way.
author Martin Almquist <malmquist@stanford.edu>
date Sun, 05 May 2019 19:05:31 -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