residual.m

MDPSAS工具箱是马里兰大学开发的

function A = residual(A)

% residual.m Nonlinear equation model object residual update.
%
% INPUT/OUTPUT PARAMETERS
%        A : an object derived from the neqmodel class
%
% Copyright (c) by Raymond A. Adomaitis, 1998-2003

unpack(A)

% update parameter d in momobj with variable d (of continuity
% object A), then evaluate

momobj = set(momobj,'param',d,'d');
momobj = equation(momobj);
unpack(momobj)

vx = a*Bvx;
vy = b*Bvy + vyBC;

% Interior residuals

Rp = wip( Dx*vx + Dy*vy, Bp );
Rp(1,1) = d(1,1)-0; % fix mean pressure field value at zero

% No boundary condition residuals

A = set(A,'resid',{Rp});
A = set(A,'param',momobj,'momobj');