gcall.m

来自「用于电力系统的一个很好的分析软件」· M 代码 · 共 10 行

function gcall(p)

global DAE Exc

if ~p.n, return, end

DAE.g = DAE.g - sparse(Exc.vref(p.exc),1,DAE.x(p.v),DAE.m,1);
DAE.g = DAE.g + sparse(p.If,1,ifield(p,1)-DAE.y(p.If),DAE.m,1);