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 + vxBC;
vy = b*Bvy;

% Interior residuals

debug = 0;

if debug
   figure(4)
   subplot(2,2,1)
   plot(vx), title('vx')
   subplot(2,2,2)
   plot(Dx*vx), title('Dx*vx')
   subplot(2,2,3)
   plot(Dy2*vy), title('Dy2*vy')
   subplot(2,2,4)
   plot(Dx*vx + Dy2*vy), title('Dx*vx + Dy2*vy')
   end


Rp = wip( Dx*vx + Dy2*vy, Bp );
% Rp(1,1) = d(1,1)-0; % fix a pressure field value for a unique solution.

mean(mean(abs(Rp)));

% No boundary condition residuals

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