159   PetscScalar  Eqp, Edp, delta;  in SetInitialGuess()  local
199     Edp = Vd + Rs[i] * Id - Xqp[i] * Iq; /* d-axis transient EMF */  in SetInitialGuess()
206     xgen[idx + 1] = Edp;  in SetInitialGuess()
247   PetscScalar  Eqp, Edp, delta;  in InitialGuess()  local
287     Edp = Vd + Rs[i] * Id - Xqp[i] * Iq; /* d-axis transient EMF */  in InitialGuess()
292     xgen[idx + 1] = Edp;  in InitialGuess()
356   PetscScalar  Eqp, Edp, delta, w; /* Generator variables */  in ResidualFunction()  local
390     Edp   = xgen[idx + 1];  in ResidualFunction()
401     fgen[idx + 1] = (Edp - (Xq[i] - Xqp[i]) * Iq) / Tq0p[i];  in ResidualFunction()
403 …fgen[idx + 3] = (-TM[i] + Edp * Id + Eqp * Iq + (Xqp[i] - Xdp[i]) * Id * Iq + D[i] * (w - w_s)) / …  in ResidualFunction()
417     fgen[idx + 4] = Zdq_inv[0] * (-Edp + Vd) + Zdq_inv[1] * (-Eqp + Vq) + Id;  in ResidualFunction()
418     fgen[idx + 5] = Zdq_inv[2] * (-Edp + Vd) + Zdq_inv[3] * (-Eqp + Vq) + Iq;  in ResidualFunction()
579   PetscScalar        Eqp, Edp, delta; /* Generator variables */  in ResidualJacobian()  local
610     Edp   = xgen[idx + 1];  in ResidualJacobian()
653     val[3] = (Edp + (Xqp[i] - Xdp[i]) * Iq) / M[i];  in ResidualJacobian()