348   PetscScalar        Id, Iq;      /* Generator dq axis currents */  in SetInitialGuess()  local
381     Iq = IGr * PetscSinScalar(theta) + IGi * PetscCosScalar(theta); /* q-axis stator current */  in SetInitialGuess()
386     Edp = Vd + Rs[i] * Id - Xqp[i] * Iq; /* d-axis transient EMF */  in SetInitialGuess()
387     Eqp = Vq + Rs[i] * Iq + Xdp[i] * Id; /* q-axis transient EMF */  in SetInitialGuess()
400     xgen[idx + 1] = Iq;  in SetInitialGuess()
439   PetscScalar        Id, Iq;             /* Generator dq axis currents */  in ResidualFunction()  local
475     Iq    = xgen[idx + 5];  in ResidualFunction()
482     fgen[idx + 1] = (-Edp + (Xq[i] - Xqp[i]) * Iq) / Tq0p[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()
499     fgen[idx + 5] = Zdq_inv[2] * (-Edp + Vd) + Zdq_inv[3] * (-Eqp + Vq) + Iq;  in ResidualFunction()
502     PetscCall(dq2ri(Id, Iq, delta, &IGr, &IGi));  in ResidualFunction()
695   PetscScalar        Id, Iq; /* Generator dq axis currents */  in ResidualJacobian()  local
728     Iq    = xgen[idx + 5];  in ResidualJacobian()
765     val[0] = -Iq / M[i];  in ResidualJacobian()
768     val[3] = (-Edp - (Xqp[i] - Xdp[i]) * Iq) / M[i];  in ResidualJacobian()
822     dIGr_ddelta = Id * PetscCosScalar(delta) - Iq * PetscSinScalar(delta);  in ResidualJacobian()
823     dIGi_ddelta = Id * PetscSinScalar(delta) + Iq * PetscCosScalar(delta);  in ResidualJacobian()