346   PetscScalar        Eqp, Edp, delta;  in SetInitialGuess()  local
387     Eqp = Vq + Rs[i] * Iq + Xdp[i] * Id; /* q-axis transient EMF */  in SetInitialGuess()
392     xgen[idx]     = Eqp;  in SetInitialGuess()
405     Efd = Eqp + (Xd[i] - Xdp[i]) * Id;  in SetInitialGuess()
437   PetscScalar        Eqp, Edp, delta, w; /* Generator variables */  in ResidualFunction()  local
470     Eqp   = xgen[idx];  in ResidualFunction()
481     fgen[idx]     = (-Eqp - (Xd[i] - Xdp[i]) * Id + Efd) / Td0p[i];  in ResidualFunction()
484 …fgen[idx + 3] = (TM[i] - Edp * Id - Eqp * Iq - (Xqp[i] - Xdp[i]) * Id * Iq - D[i] * (w - w_s)) / M…  in ResidualFunction()
498     fgen[idx + 4] = Zdq_inv[0] * (-Edp + Vd) + Zdq_inv[1] * (-Eqp + Vq) + Id;  in ResidualFunction()
499     fgen[idx + 5] = Zdq_inv[2] * (-Edp + Vd) + Zdq_inv[3] * (-Eqp + Vq) + Iq;  in ResidualFunction()
693   PetscScalar        Eqp, Edp, delta; /* Generator variables */  in ResidualJacobian()  local
724     Eqp   = xgen[idx];  in ResidualJacobian()
769     val[4] = (-Eqp - (Xqp[i] - Xdp[i]) * Id) / M[i];  in ResidualJacobian()