158   PetscScalar  Vr, Vi, IGr, IGi, Vm, Vm2;  in SetInitialGuess()  local
186     IGr = (Vr * PG[i] + Vi * QG[i]) / Vm2;  in SetInitialGuess()
189     delta = PetscAtan2Real(Vi + Xq[i] * IGr, Vr - Xq[i] * IGi); /* Machine angle */  in SetInitialGuess()
193     Id = IGr * PetscCosScalar(theta) - IGi * PetscSinScalar(theta); /* d-axis stator current */  in SetInitialGuess()
194     Iq = IGr * PetscSinScalar(theta) + IGi * PetscCosScalar(theta); /* q-axis stator current */  in SetInitialGuess()
246   PetscScalar  Vr, Vi, IGr, IGi, Vm, Vm2;  in InitialGuess()  local
274     IGr = (Vr * PGv[i] + Vi * QG[i]) / Vm2;  in InitialGuess()
277     delta = PetscAtan2Real(Vi + Xq[i] * IGr, Vr - Xq[i] * IGi); /* Machine angle */  in InitialGuess()
281     Id = IGr * PetscCosScalar(theta) - IGi * PetscSinScalar(theta); /* d-axis stator current */  in InitialGuess()
282     Iq = IGr * PetscSinScalar(theta) + IGi * PetscCosScalar(theta); /* q-axis stator current */  in InitialGuess()
360   PetscScalar  IGr, IGi, IDr, IDi;  in ResidualFunction()  local
421     PetscCall(dq2ri(Id, Iq, delta, &IGr, &IGi));  in ResidualFunction()
424     fnet[2 * gbus[i] + 1] -= IGr;  in ResidualFunction()