349   PetscScalar        theta, Vd, Vq, SE;  in SetInitialGuess()  local
384     Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta);  in SetInitialGuess()
387     Eqp = Vq + Rs[i] * Iq + Xdp[i] * Id; /* q-axis transient EMF */  in SetInitialGuess()
440   PetscScalar        Vd, Vq, SE;  in ResidualFunction()  local
489     PetscCall(ri2dq(Vr, Vi, delta, &Vd, &Vq));  in ResidualFunction()
498     fgen[idx + 4] = Zdq_inv[0] * (-Edp + Vd) + Zdq_inv[1] * (-Eqp + Vq) + Id;  in ResidualFunction()
499     fgen[idx + 5] = Zdq_inv[2] * (-Edp + Vd) + Zdq_inv[3] * (-Eqp + Vq) + Iq;  in ResidualFunction()
507     Vm = PetscSqrtScalar(Vd * Vd + Vq * Vq);  in ResidualFunction()
696   PetscScalar        Vd, Vq;  in ResidualJacobian()  local
774     PetscCall(ri2dq(Vr, Vi, delta, &Vd, &Vq));  in ResidualJacobian()
849     Vm = PetscSqrtScalar(Vd * Vd + Vq * Vq);  in ResidualJacobian()
887       dVm_dVq = Vq / Vm;  in ResidualJacobian()