158   PetscScalar  Vr, Vi, IGr, IGi, Vm, Vm2;  in SetInitialGuess()  local
183     Vi  = xnet[2 * gbus[i] + 1]; /* Imaginary part of the generator terminal voltage */  in SetInitialGuess()
184     Vm  = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in SetInitialGuess()
186     IGr = (Vr * PG[i] + Vi * QG[i]) / Vm2;  in SetInitialGuess()
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()
196     Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta);  in SetInitialGuess()
197     Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta);  in SetInitialGuess()
246   PetscScalar  Vr, Vi, IGr, IGi, Vm, Vm2;  in InitialGuess()  local
271     Vi  = xnet[2 * gbus[i] + 1]; /* Imaginary part of the generator terminal voltage */  in InitialGuess()
272     Vm  = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in InitialGuess()
274     IGr = (Vr * PGv[i] + Vi * QG[i]) / Vm2;  in InitialGuess()
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()
284     Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta);  in InitialGuess()
285     Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta);  in InitialGuess()
355   PetscScalar  Vr, Vi, Vm, Vm2;  in ResidualFunction()  local
406     Vi = xnet[2 * gbus[i] + 1]; /* Imaginary part of the generator terminal voltage */  in ResidualFunction()
408     PetscCall(ri2dq(Vr, Vi, delta, &Vd, &Vq));  in ResidualFunction()
441     Vi  = xnet[2 * lbus[i] + 1]; /* Imaginary part of the load bus voltage */  in ResidualFunction()
442     Vm  = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in ResidualFunction()
450     IDr = (PD * Vr + QD * Vi) / Vm2;  in ResidualFunction()
451     IDi = (-QD * Vr + PD * Vi) / Vm2;  in ResidualFunction()
578   PetscScalar        Vr, Vi, Vm, Vm2;  in ResidualJacobian()  local
658     Vi = xnet[2 * gbus[i] + 1]; /* Imaginary part of the generator terminal voltage */  in ResidualJacobian()
659     PetscCall(ri2dq(Vr, Vi, delta, &Vd, &Vq));  in ResidualJacobian()
672     dVd_ddelta = Vr * PetscCosScalar(delta) + Vi * PetscSinScalar(delta);  in ResidualJacobian()
673     dVq_ddelta = -Vr * PetscSinScalar(delta) + Vi * PetscCosScalar(delta);  in ResidualJacobian()
804     Vi  = xnet[2 * lbus[i] + 1]; /* Imaginary part of the load bus voltage */  in ResidualJacobian()
805     Vm  = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in ResidualJacobian()
814 …dPD_dVi += ld_alphap[k] * ld_betap[k] * PD0[i] * PetscPowScalar(1 / Vm0, ld_betap[k]) * Vi * Petsc…  in ResidualJacobian()
819 …dQD_dVi += ld_alphaq[k] * ld_betaq[k] * QD0[i] * PetscPowScalar(1 / Vm0, ld_betaq[k]) * Vi * Petsc…  in ResidualJacobian()
825     dIDr_dVr = (dPD_dVr * Vr + dQD_dVr * Vi + PD) / Vm2 - ((PD * Vr + QD * Vi) * 2 * Vr) / Vm4;  in ResidualJacobian()
826     dIDr_dVi = (dPD_dVi * Vr + dQD_dVi * Vi + QD) / Vm2 - ((PD * Vr + QD * Vi) * 2 * Vi) / Vm4;  in ResidualJacobian()
828     dIDi_dVr = (-dQD_dVr * Vr + dPD_dVr * Vi - QD) / Vm2 - ((-QD * Vr + PD * Vi) * 2 * Vr) / Vm4;  in ResidualJacobian()
829     dIDi_dVi = (-dQD_dVi * Vr + dPD_dVi * Vi + PD) / Vm2 - ((-QD * Vr + PD * Vi) * 2 * Vi) / Vm4;  in ResidualJacobian()