142   PetscScalar        Efd, RF, VR, Vr, Vi, Vm;  in EventFunction()  local
163     Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in EventFunction()
168 … = PetscRealPart((VR - KA[i] * RF + KA[i] * KF[i] * Efd / TF[i] - KA[i] * (Vref[i] - Vm)) / TA[i]);  in EventFunction()
173 … = PetscRealPart((VR - KA[i] * RF + KA[i] * KF[i] * Efd / TF[i] - KA[i] * (Vref[i] - Vm)) / TA[i]);  in EventFunction()
194   PetscScalar  Vr, Vi, Vm;  in PostEventFunction()  local
247         Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in PostEventFunction()
250           fvalue = (VR - KA[i] * RF + KA[i] * KF[i] * Efd / TF[i] - KA[i] * (Vref[i] - Vm)) / TA[i];  in PostEventFunction()
257           fvalue = (VR - KA[i] * RF + KA[i] * KF[i] * Efd / TF[i] - KA[i] * (Vref[i] - Vm)) / TA[i];  in PostEventFunction()
345   PetscScalar        Vr, Vi, IGr, IGi, Vm, Vm2;  in SetInitialGuess()  local
371     Vm  = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in SetInitialGuess()
372     Vm2 = Vm * Vm;  in SetInitialGuess()
414     Vref[i] = Vm + (VR / KA[i]);  in SetInitialGuess()
436   PetscScalar        Vr, Vi, Vm, Vm2;  in ResidualFunction()  local
507     Vm = PetscSqrtScalar(Vd * Vd + Vq * Vq);  in ResidualFunction()
516 …else fgen[idx + 8] = (-VR + KA[i] * RF - KA[i] * KF[i] * Efd / TF[i] + KA[i] * (Vref[i] - Vm)) / T…  in ResidualFunction()
525     Vm  = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in ResidualFunction()
526     Vm2 = Vm * Vm;  in ResidualFunction()
529 …for (k = 0; k < ld_nsegsp[i]; k++) PD += ld_alphap[k] * PD0[i] * PetscPowScalar(Vm / Vm0, ld_betap…  in ResidualFunction()
530 …for (k = 0; k < ld_nsegsq[i]; k++) QD += ld_alphaq[k] * QD0[i] * PetscPowScalar(Vm / Vm0, ld_betaq…  in ResidualFunction()
692   PetscScalar        Vr, Vi, Vm, Vm2;  in ResidualJacobian()  local
849     Vm = PetscSqrtScalar(Vd * Vd + Vq * Vq);  in ResidualJacobian()
886       dVm_dVd = Vd / Vm;  in ResidualJacobian()
887       dVm_dVq = Vq / Vm;  in ResidualJacobian()
933     Vm  = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in ResidualJacobian()
934     Vm2 = Vm * Vm;  in ResidualJacobian()
940       PD += ld_alphap[k] * PD0[i] * PetscPowScalar(Vm / Vm0, ld_betap[k]);  in ResidualJacobian()
941 …etap[k] * PD0[i] * PetscPowScalar(1 / Vm0, ld_betap[k]) * Vr * PetscPowScalar(Vm, ld_betap[k] - 2);  in ResidualJacobian()
942 …etap[k] * PD0[i] * PetscPowScalar(1 / Vm0, ld_betap[k]) * Vi * PetscPowScalar(Vm, ld_betap[k] - 2);  in ResidualJacobian()
945       QD += ld_alphaq[k] * QD0[i] * PetscPowScalar(Vm / Vm0, ld_betaq[k]);  in ResidualJacobian()
946 …etaq[k] * QD0[i] * PetscPowScalar(1 / Vm0, ld_betaq[k]) * Vr * PetscPowScalar(Vm, ld_betaq[k] - 2);  in ResidualJacobian()
947 …etaq[k] * QD0[i] * PetscPowScalar(1 / Vm0, ld_betaq[k]) * Vi * PetscPowScalar(Vm, ld_betaq[k] - 2);  in ResidualJacobian()