residualold.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,'var')
unpack(A,'param')

vx    = a*Bvx;
vy    = b*Bvy + vyBC;
press = d*Bp;

% Interior residuals

Rvx = wip( DDx*vx + DDy*vx - Dx*press, Bvx );
Rvy = wip( DDx*vy + DDy*vy - Dy*press, Bvy );
Rp  = wip( Dx*vx + Dy*vy, Bp );
Rp(1,1) = d(1,1)-0; % set mean pressure = 0;

% No boundary condition residuals

A = set(A,'resid',{Rvx Rvy Rp});