122   PetscScalar  Eqp, Edp, delta;  in SetInitialGuess()  local
163     Eqp = Vq + Rs[i] * Iq + Xdp[i] * Id; /* q-axis transient EMF */  in SetInitialGuess()
168     xgen[idx]     = Eqp;  in SetInitialGuess()
181     Efd = Eqp + (Xd[i] - Xdp[i]) * Id;  in SetInitialGuess()
211   PetscScalar  Eqp, Edp, delta, w; /* Generator variables */  in ResidualFunction()  local
244     Eqp   = xgen[idx];  in ResidualFunction()
255     fgen[idx]     = (Eqp + (Xd[i] - Xdp[i]) * Id - Efd) / Td0p[i];  in ResidualFunction()
258 …fgen[idx + 3] = (-TM[i] + Edp * Id + Eqp * Iq + (Xqp[i] - Xdp[i]) * Id * Iq + D[i] * (w - w_s)) / …  in ResidualFunction()
272     fgen[idx + 4] = Zdq_inv[0] * (-Edp + Vd) + Zdq_inv[1] * (-Eqp + Vq) + Id;  in ResidualFunction()
273     fgen[idx + 5] = Zdq_inv[2] * (-Edp + Vd) + Zdq_inv[3] * (-Eqp + Vq) + Iq;  in ResidualFunction()
436   PetscScalar        Eqp, Edp, delta; /* Generator variables */  in ResidualJacobian()  local
466     Eqp   = xgen[idx];  in ResidualJacobian()
511     val[4] = (Eqp + (Xqp[i] - Xdp[i]) * Id) / M[i];  in ResidualJacobian()