215   PetscScalar  PD, QD, Vm0, *v0;  in ResidualFunction()  local
299     PD = QD = 0.0;  in ResidualFunction()
301 …for (k = 0; k < ld_nsegsq[i]; k++) QD += ld_alphaq[k] * QD0[i] * PetscPowScalar(Vm / Vm0, ld_betaq…  in ResidualFunction()
304     IDr = (PD * Vr + QD * Vi) / Vm2;  in ResidualFunction()
305     IDi = (-QD * Vr + PD * Vi) / Vm2;  in ResidualFunction()
451   PetscScalar        PD, QD, Vm0, *v0, Vm4;  in ResidualJacobian()  local
665     PD = QD = 0.0;  in ResidualJacobian()
673       QD += ld_alphaq[k] * QD0[i] * PetscPowScalar(Vm / Vm0, ld_betaq[k]);  in ResidualJacobian()
681     dIDr_dVr = (dPD_dVr * Vr + dQD_dVr * Vi + PD) / Vm2 - ((PD * Vr + QD * Vi) * 2 * Vr) / Vm4;  in ResidualJacobian()
682     dIDr_dVi = (dPD_dVi * Vr + dQD_dVi * Vi + QD) / Vm2 - ((PD * Vr + QD * Vi) * 2 * Vi) / Vm4;  in ResidualJacobian()
684     dIDi_dVr = (-dQD_dVr * Vr + dPD_dVr * Vi - QD) / Vm2 - ((-QD * Vr + PD * Vi) * 2 * Vr) / Vm4;  in ResidualJacobian()
685     dIDi_dVi = (-dQD_dVi * Vr + dPD_dVi * Vi + PD) / Vm2 - ((-QD * Vr + PD * Vi) * 2 * Vi) / Vm4;  in ResidualJacobian()