equation.m

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

function A = equation(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)

press = d*Bp;

% Interior residuals

vxrhs   = wip( Dx*press, Bvx );
Rvx     = mprod(vxcoeffinv,vxrhs);

vyrhs   = wip( - (DDx*vyBC) - DDy*vyBC + Dy*press, Bvy );
Rvy     = mprod(vycoeffinv,vyrhs);

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