162   PetscScalar  theta, Vd, Vq, SE;  in SetInitialGuess()  local
196     Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta);  in SetInitialGuess()
199     Edp = Vd + Rs[i] * Id - Xqp[i] * Iq; /* d-axis transient EMF */  in SetInitialGuess()
250   PetscScalar  theta, Vd, Vq, SE;  in InitialGuess()  local
284     Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta);  in InitialGuess()
287     Edp = Vd + Rs[i] * Id - Xqp[i] * Iq; /* d-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()
759     dVm_dVd = Vd / Vm;  in ResidualJacobian()