158   PetscScalar  Vr, Vi, IGr, IGi, Vm, Vm2;  in SetInitialGuess()  local
187     IGi = (Vi * PG[i] - Vr * 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
275     IGi = (Vi * PGv[i] - Vr * 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()
423     fnet[2 * gbus[i]] -= IGi;  in ResidualFunction()