443   PetscScalar        PD, QD, Vm0, *v0;  in ResidualFunction()  local
528     PD = QD = 0.0;  in ResidualFunction()
529 …for (k = 0; k < ld_nsegsp[i]; k++) PD += ld_alphap[k] * PD0[i] * PetscPowScalar(Vm / Vm0, ld_betap…  in ResidualFunction()
533     IDr = (PD * Vr + QD * Vi) / Vm2;  in ResidualFunction()
534     IDi = (-QD * Vr + PD * Vi) / Vm2;  in ResidualFunction()
709   PetscScalar        PD, QD, Vm0, Vm4;  in ResidualJacobian()  local
937     PD = QD = 0.0;  in ResidualJacobian()
940       PD += ld_alphap[k] * PD0[i] * PetscPowScalar(Vm / Vm0, ld_betap[k]);  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()