142   PetscScalar        Efd, RF, VR, Vr, Vi, Vm;  in EventFunction()  local
161     Vr = xnet[2 * gbus[i]];     /* Real part of generator terminal voltage */  in EventFunction()
163     Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in EventFunction()
194   PetscScalar  Vr, Vi, Vm;  in PostEventFunction()  local
245         Vr = xnet[2 * gbus[i]];     /* Real part of generator terminal voltage */  in PostEventFunction()
247         Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in PostEventFunction()
345   PetscScalar        Vr, Vi, IGr, IGi, Vm, Vm2;  in SetInitialGuess()  local
369     Vr  = xnet[2 * gbus[i]];     /* Real part of generator terminal voltage */  in SetInitialGuess()
371     Vm  = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in SetInitialGuess()
373     IGr = (Vr * PG[i] + Vi * QG[i]) / Vm2;  in SetInitialGuess()
374     IGi = (Vi * PG[i] - Vr * QG[i]) / Vm2;  in SetInitialGuess()
376     delta = PetscAtan2Real(Vi + Xq[i] * IGr, Vr - Xq[i] * IGi); /* Machine angle */  in SetInitialGuess()
383     Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta);  in SetInitialGuess()
384     Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta);  in SetInitialGuess()
436   PetscScalar        Vr, Vi, Vm, Vm2;  in ResidualFunction()  local
486     Vr = xnet[2 * gbus[i]];     /* Real part of generator terminal voltage */  in ResidualFunction()
489     PetscCall(ri2dq(Vr, Vi, delta, &Vd, &Vq));  in ResidualFunction()
523     Vr  = xnet[2 * lbus[i]];     /* Real part of load bus voltage */  in ResidualFunction()
525     Vm  = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in ResidualFunction()
533     IDr = (PD * Vr + QD * Vi) / Vm2;  in ResidualFunction()
534     IDi = (-QD * Vr + PD * Vi) / Vm2;  in ResidualFunction()
692   PetscScalar        Vr, Vi, Vm, Vm2;  in ResidualJacobian()  local
772     Vr = xnet[2 * gbus[i]];     /* Real part of generator terminal voltage */  in ResidualJacobian()
774     PetscCall(ri2dq(Vr, Vi, delta, &Vd, &Vq));  in ResidualJacobian()
787     dVd_ddelta = Vr * PetscCosScalar(delta) + Vi * PetscSinScalar(delta);  in ResidualJacobian()
788     dVq_ddelta = -Vr * PetscSinScalar(delta) + Vi * PetscCosScalar(delta);  in ResidualJacobian()
931     Vr  = xnet[2 * lbus[i]];     /* Real part of load bus voltage */  in ResidualJacobian()
933     Vm  = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in ResidualJacobian()
941 …dPD_dVr += ld_alphap[k] * ld_betap[k] * PD0[i] * PetscPowScalar(1 / Vm0, ld_betap[k]) * Vr * Petsc…  in ResidualJacobian()
946 …dQD_dVr += ld_alphaq[k] * ld_betaq[k] * QD0[i] * PetscPowScalar(1 / Vm0, ld_betaq[k]) * Vr * Petsc…  in ResidualJacobian()
953     dIDr_dVr = (dPD_dVr * Vr + dQD_dVr * Vi + PD) / Vm2 - ((PD * Vr + QD * Vi) * 2 * Vr) / Vm4;  in ResidualJacobian()
954     dIDr_dVi = (dPD_dVi * Vr + dQD_dVi * Vi + QD) / Vm2 - ((PD * Vr + QD * Vi) * 2 * Vi) / Vm4;  in ResidualJacobian()
956     dIDi_dVr = (-dQD_dVr * Vr + dPD_dVr * Vi - QD) / Vm2 - ((-QD * Vr + PD * Vi) * 2 * Vr) / Vm4;  in ResidualJacobian()
957     dIDi_dVi = (-dQD_dVi * Vr + dPD_dVi * Vi + PD) / Vm2 - ((-QD * Vr + PD * Vi) * 2 * Vi) / Vm4;  in ResidualJacobian()