121   PetscScalar  Vr, Vi, IGr, IGi, Vm, Vm2;  in SetInitialGuess()  local
146     Vi  = xnet[2 * gbus[i] + 1]; /* Imaginary part of the generator terminal voltage */  in SetInitialGuess()
147     Vm  = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in SetInitialGuess()
149     IGr = (Vr * PG[i] + Vi * QG[i]) / Vm2;  in SetInitialGuess()
150     IGi = (Vi * PG[i] - Vr * QG[i]) / Vm2;  in SetInitialGuess()
152     delta = PetscAtan2Real(Vi + Xq[i] * IGr, Vr - Xq[i] * IGi); /* Machine angle */  in SetInitialGuess()
159     Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta);  in SetInitialGuess()
160     Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta);  in SetInitialGuess()
210   PetscScalar  Vr, Vi, Vm, Vm2;  in ResidualFunction()  local
261     Vi = xnet[2 * gbus[i] + 1]; /* Imaginary part of the generator terminal voltage */  in ResidualFunction()
263     PetscCall(ri2dq(Vr, Vi, delta, &Vd, &Vq));  in ResidualFunction()
296     Vi  = xnet[2 * lbus[i] + 1]; /* Imaginary part of the load bus voltage */  in ResidualFunction()
297     Vm  = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in ResidualFunction()
305     IDr = (PD * Vr + QD * Vi) / Vm2;  in ResidualFunction()
306     IDi = (-QD * Vr + PD * Vi) / Vm2;  in ResidualFunction()
435   PetscScalar        Vr, Vi, Vm, Vm2;  in ResidualJacobian()  local
515     Vi = xnet[2 * gbus[i] + 1]; /* Imaginary part of the generator terminal voltage */  in ResidualJacobian()
516     PetscCall(ri2dq(Vr, Vi, delta, &Vd, &Vq));  in ResidualJacobian()
529     dVd_ddelta = Vr * PetscCosScalar(delta) + Vi * PetscSinScalar(delta);  in ResidualJacobian()
530     dVq_ddelta = -Vr * PetscSinScalar(delta) + Vi * PetscCosScalar(delta);  in ResidualJacobian()
662     Vi  = xnet[2 * lbus[i] + 1]; /* Imaginary part of the load bus voltage */  in ResidualJacobian()
663     Vm  = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in ResidualJacobian()
672 …dPD_dVi += ld_alphap[k] * ld_betap[k] * PD0[i] * PetscPowScalar(1 / Vm0, ld_betap[k]) * Vi * Petsc…  in ResidualJacobian()
677 …dQD_dVi += ld_alphaq[k] * ld_betaq[k] * QD0[i] * PetscPowScalar(1 / Vm0, ld_betaq[k]) * Vi * Petsc…  in ResidualJacobian()
683     dIDr_dVr = (dPD_dVr * Vr + dQD_dVr * Vi + PD) / Vm2 - ((PD * Vr + QD * Vi) * 2 * Vr) / Vm4;  in ResidualJacobian()
684     dIDr_dVi = (dPD_dVi * Vr + dQD_dVi * Vi + QD) / Vm2 - ((PD * Vr + QD * Vi) * 2 * Vi) / Vm4;  in ResidualJacobian()
686     dIDi_dVr = (-dQD_dVr * Vr + dPD_dVr * Vi - QD) / Vm2 - ((-QD * Vr + PD * Vi) * 2 * Vr) / Vm4;  in ResidualJacobian()
687     dIDi_dVi = (-dQD_dVi * Vr + dPD_dVi * Vi + PD) / Vm2 - ((-QD * Vr + PD * Vi) * 2 * Vi) / Vm4;  in ResidualJacobian()