residual.m

来自「MDPSAS工具箱是马里兰大学开发的」· M 代码 · 共 35 行

function [A, Jacobian, Rescol] = residual(A)

% residual.m Nonlinear equation model object residual update.

unpack(A,'var')
unpack(A,'param')

% residuals

Rescell{1} = DDx*vx + DDy*vx - Dx*p;
Rescell{1} = setbval(Rescell{1},vx,'x','min'); % BC residuals
Rescell{1} = setbval(Rescell{1},vx,'x','max');
Rescell{1} = setbval(Rescell{1},vx,'y','min');
Rescell{1} = setbval(Rescell{1},vx,'y','max');

Rescell{2} = DDx*vy + DDy*vy - Dy*p;
vyres = vy - vybc;
Rescell{2} = setbval(Rescell{2},vy,'x','min'); % BC residuals
Rescell{2} = setbval(Rescell{2},vy,'y','min');
Rescell{2} = setbval(Rescell{2},vy,'y','max');
Rescell{2} = setbval(Rescell{2},vyres,'x','max');

Rescell{3} = Dx*vx + Dy*vy;
Dxp = Dx*p;
Dyp = Dy*p;
Rescell{3} = setbval(Rescell{3},Dxp,'x','min'); % BC residuals
Rescell{3} = setbval(Rescell{3},Dxp,'x','max');
Rescell{3} = setbval(Rescell{3},Dyp,'y','min');
Rescell{3} = setbval(Rescell{3},Dyp,'y','max');

Rescell{3} = setval(Rescell{3},getval(p,2,2),2,2);

A = set(A,'resid',Rescell);