122   PetscScalar  Id, Iq;      /* Generator dq axis currents */  in SetInitialGuess()  local
154     Id = IGr * PetscCosScalar(theta) - IGi * PetscSinScalar(theta); /* d-axis stator current */  in SetInitialGuess()
160     Edp = Vd + Rs[i] * Id - Xqp[i] * Iq; /* d-axis transient EMF */  in SetInitialGuess()
161     Eqp = Vq + Rs[i] * Iq + Xdp[i] * Id; /* q-axis transient EMF */  in SetInitialGuess()
173     xgen[idx]     = Id;  in SetInitialGuess()
179     Efd = Eqp + (Xd[i] - Xdp[i]) * Id;  in SetInitialGuess()
211   PetscScalar  Id, Iq;             /* Generator dq axis currents */  in ResidualFunction()  local
246     Id    = xgen[idx + 4];  in ResidualFunction()
253     fgen[idx]     = (Eqp + (Xd[i] - Xdp[i]) * Id - Efd) / Td0p[i];  in ResidualFunction()
256 …fgen[idx + 3] = (-TM[i] + Edp * Id + Eqp * Iq + (Xqp[i] - Xdp[i]) * Id * Iq + D[i] * (w - w_s)) / …  in ResidualFunction()
271     fgen[idx + 4] = Zdq_inv[0] * (-Edp + Vd) + Zdq_inv[1] * (-Eqp + Vq) + Id;  in ResidualFunction()
275     PetscCall(dq2ri(Id, Iq, delta, &IGr, &IGi));  in ResidualFunction()
437   PetscScalar        Id, Iq;          /* Generator dq axis currents */  in ResidualJacobian()  local
468     Id    = xgen[idx + 4];  in ResidualJacobian()
507     val[1] = Id / M[i];  in ResidualJacobian()
510     val[4] = (Eqp + (Xqp[i] - Xdp[i]) * Id) / M[i];  in ResidualJacobian()
563     dIGr_ddelta = Id * PetscCosScalar(delta) - Iq * PetscSinScalar(delta);  in ResidualJacobian()
564     dIGi_ddelta = Id * PetscSinScalar(delta) + Iq * PetscCosScalar(delta);  in ResidualJacobian()