161   PetscScalar  Id, Iq;      /* Generator dq axis currents */  in SetInitialGuess()  local
194     Iq = IGr * PetscSinScalar(theta) + IGi * PetscCosScalar(theta); /* q-axis stator current */  in SetInitialGuess()
199     Edp = Vd + Rs[i] * Id - Xqp[i] * Iq; /* d-axis transient EMF */  in SetInitialGuess()
200     Eqp = Vq + Rs[i] * Iq + Xdp[i] * Id; /* q-axis transient EMF */  in SetInitialGuess()
213     xgen[idx + 1] = Iq;  in SetInitialGuess()
249   PetscScalar  Id, Iq;      /* Generator dq axis currents */  in InitialGuess()  local
282     Iq = IGr * PetscSinScalar(theta) + IGi * PetscCosScalar(theta); /* q-axis stator current */  in InitialGuess()
287     Edp = Vd + Rs[i] * Id - Xqp[i] * Iq; /* d-axis transient EMF */  in InitialGuess()
288     Eqp = Vq + Rs[i] * Iq + Xdp[i] * Id; /* q-axis transient EMF */  in InitialGuess()
299     xgen[idx + 1] = Iq;  in InitialGuess()
358   PetscScalar  Id, Iq;             /* Generator dq axis currents */  in ResidualFunction()  local
394     Iq    = xgen[idx + 5];  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()
418     fgen[idx + 5] = Zdq_inv[2] * (-Edp + Vd) + Zdq_inv[3] * (-Eqp + Vq) + Iq;  in ResidualFunction()
421     PetscCall(dq2ri(Id, Iq, delta, &IGr, &IGi));  in ResidualFunction()
581   PetscScalar        Id, Iq;          /* Generator dq axis currents */  in ResidualJacobian()  local
613     Iq    = xgen[idx + 5];  in ResidualJacobian()
650     val[0] = Iq / M[i];  in ResidualJacobian()
653     val[3] = (Edp + (Xqp[i] - Xdp[i]) * Iq) / M[i];  in ResidualJacobian()
707     dIGr_ddelta = Id * PetscCosScalar(delta) - Iq * PetscSinScalar(delta);  in ResidualJacobian()
708     dIGi_ddelta = Id * PetscSinScalar(delta) + Iq * PetscCosScalar(delta);  in ResidualJacobian()