119   PetscScalar  Vr, Vi, IGr, IGi, Vm, Vm2;  in SetInitialGuess()  local
144     Vi  = xnet[2 * gbus[i] + 1]; /* Imaginary part of the generator terminal voltage */  in SetInitialGuess()
145     Vm  = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in SetInitialGuess()
147     IGr = (Vr * PG[i] + Vi * QG[i]) / Vm2;  in SetInitialGuess()
148     IGi = (Vi * PG[i] - Vr * QG[i]) / Vm2;  in SetInitialGuess()
150     delta = PetscAtan2Real(Vi + Xq[i] * IGr, Vr - Xq[i] * IGi); /* Machine angle */  in SetInitialGuess()
157     Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta);  in SetInitialGuess()
158     Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta);  in SetInitialGuess()
208   PetscScalar  Vr, Vi, Vm, Vm2;  in ResidualFunction()  local
259     Vi = xnet[2 * gbus[i] + 1]; /* Imaginary part of the generator terminal voltage */  in ResidualFunction()
261     PetscCall(ri2dq(Vr, Vi, delta, &Vd, &Vq));  in ResidualFunction()
295     Vi  = xnet[2 * lbus[i] + 1]; /* Imaginary part of the load bus voltage */  in ResidualFunction()
296     Vm  = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in ResidualFunction()
304     IDr = (PD * Vr + QD * Vi) / Vm2;  in ResidualFunction()
305     IDi = (-QD * Vr + PD * Vi) / Vm2;  in ResidualFunction()
434   PetscScalar        Vr, Vi, Vm, Vm2;  in ResidualJacobian()  local
514     Vi = xnet[2 * gbus[i] + 1]; /* Imaginary part of the generator terminal voltage */  in ResidualJacobian()
515     PetscCall(ri2dq(Vr, Vi, delta, &Vd, &Vq));  in ResidualJacobian()
528     dVd_ddelta = Vr * PetscCosScalar(delta) + Vi * PetscSinScalar(delta);  in ResidualJacobian()
529     dVq_ddelta = -Vr * PetscSinScalar(delta) + Vi * PetscCosScalar(delta);  in ResidualJacobian()
660     Vi  = xnet[2 * lbus[i] + 1]; /* Imaginary part of the load bus voltage */  in ResidualJacobian()
661     Vm  = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in ResidualJacobian()
670 …dPD_dVi += ld_alphap[k] * ld_betap[k] * PD0[i] * PetscPowScalar(1 / Vm0, ld_betap[k]) * Vi * Petsc…  in ResidualJacobian()
675 …dQD_dVi += ld_alphaq[k] * ld_betaq[k] * QD0[i] * PetscPowScalar(1 / Vm0, ld_betaq[k]) * Vi * Petsc…  in ResidualJacobian()
681     dIDr_dVr = (dPD_dVr * Vr + dQD_dVr * Vi + PD) / Vm2 - ((PD * Vr + QD * Vi) * 2 * Vr) / Vm4;  in ResidualJacobian()
682     dIDr_dVi = (dPD_dVi * Vr + dQD_dVi * Vi + QD) / Vm2 - ((PD * Vr + QD * Vi) * 2 * Vi) / Vm4;  in ResidualJacobian()
684     dIDi_dVr = (-dQD_dVr * Vr + dPD_dVr * Vi - QD) / Vm2 - ((-QD * Vr + PD * Vi) * 2 * Vr) / Vm4;  in ResidualJacobian()
685     dIDi_dVi = (-dQD_dVi * Vr + dPD_dVi * Vi + PD) / Vm2 - ((-QD * Vr + PD * Vi) * 2 * Vi) / Vm4;  in ResidualJacobian()