162   PetscScalar  theta, Vd, Vq, SE;  in SetInitialGuess()  local
197     Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta);  in SetInitialGuess()
200     Eqp = Vq + Rs[i] * Iq + Xdp[i] * Id; /* q-axis transient EMF */  in SetInitialGuess()
250   PetscScalar  theta, Vd, Vq, SE;  in InitialGuess()  local
285     Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta);  in InitialGuess()
288     Eqp = Vq + Rs[i] * Iq + Xdp[i] * Id; /* q-axis transient EMF */  in InitialGuess()
359   PetscScalar  Vd, Vq, SE;  in ResidualFunction()  local
408     PetscCall(ri2dq(Vr, Vi, delta, &Vd, &Vq));  in ResidualFunction()
417     fgen[idx + 4] = Zdq_inv[0] * (-Edp + Vd) + Zdq_inv[1] * (-Eqp + Vq) + Id;  in ResidualFunction()
418     fgen[idx + 5] = Zdq_inv[2] * (-Edp + Vd) + Zdq_inv[3] * (-Eqp + Vq) + Iq;  in ResidualFunction()
426     Vm = PetscSqrtScalar(Vd * Vd + Vq * Vq);  in ResidualFunction()
582   PetscScalar        Vd, Vq;  in ResidualJacobian()  local
659     PetscCall(ri2dq(Vr, Vi, delta, &Vd, &Vq));  in ResidualJacobian()
734     Vm = PetscSqrtScalar(Vd * Vd + Vq * Vq);  in ResidualJacobian()
760     dVm_dVq = Vq / Vm;  in ResidualJacobian()