158   PetscScalar  Vr, Vi, IGr, IGi, Vm, Vm2;  in SetInitialGuess()  local
184     Vm  = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in SetInitialGuess()
185     Vm2 = Vm * Vm;  in SetInitialGuess()
227     Vref[i] = Vm + (VR / KA[i]);  in SetInitialGuess()
246   PetscScalar  Vr, Vi, IGr, IGi, Vm, Vm2;  in InitialGuess()  local
272     Vm  = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in InitialGuess()
273     Vm2 = Vm * Vm;  in InitialGuess()
355   PetscScalar  Vr, Vi, Vm, Vm2;  in ResidualFunction()  local
426     Vm = PetscSqrtScalar(Vd * Vd + Vq * Vq);  in ResidualFunction()
433 …  fgen[idx + 8] = (VR - KA[i] * RF + KA[i] * KF[i] * Efd / TF[i] - KA[i] * (Vref[i] - Vm)) / TA[i];  in ResidualFunction()
442     Vm  = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in ResidualFunction()
443     Vm2 = Vm * Vm;  in ResidualFunction()
446 …for (k = 0; k < ld_nsegsp[i]; k++) PD += ld_alphap[k] * PD0[i] * PetscPowScalar(Vm / Vm0, ld_betap…  in ResidualFunction()
447 …for (k = 0; k < ld_nsegsq[i]; k++) QD += ld_alphaq[k] * QD0[i] * PetscPowScalar(Vm / Vm0, ld_betaq…  in ResidualFunction()
578   PetscScalar        Vr, Vi, Vm, Vm2;  in ResidualJacobian()  local
734     Vm = PetscSqrtScalar(Vd * Vd + Vq * Vq);  in ResidualJacobian()
759     dVm_dVd = Vd / Vm;  in ResidualJacobian()
760     dVm_dVq = Vq / Vm;  in ResidualJacobian()
805     Vm  = PetscSqrtScalar(Vr * Vr + Vi * Vi);  in ResidualJacobian()
806     Vm2 = Vm * Vm;  in ResidualJacobian()
812       PD += ld_alphap[k] * PD0[i] * PetscPowScalar(Vm / Vm0, ld_betap[k]);  in ResidualJacobian()
813 …etap[k] * PD0[i] * PetscPowScalar(1 / Vm0, ld_betap[k]) * Vr * PetscPowScalar(Vm, ld_betap[k] - 2);  in ResidualJacobian()
814 …etap[k] * PD0[i] * PetscPowScalar(1 / Vm0, ld_betap[k]) * Vi * PetscPowScalar(Vm, ld_betap[k] - 2);  in ResidualJacobian()
817       QD += ld_alphaq[k] * QD0[i] * PetscPowScalar(Vm / Vm0, ld_betaq[k]);  in ResidualJacobian()
818 …etaq[k] * QD0[i] * PetscPowScalar(1 / Vm0, ld_betaq[k]) * Vr * PetscPowScalar(Vm, ld_betaq[k] - 2);  in ResidualJacobian()
819 …etaq[k] * QD0[i] * PetscPowScalar(1 / Vm0, ld_betaq[k]) * Vi * PetscPowScalar(Vm, ld_betaq[k] - 2);  in ResidualJacobian()