208   PetscScalar *xgen, *xnet, *fgen, *fnet;  in ResidualFunction()  local
239   PetscCall(VecGetArray(Fgen, &fgen));  in ResidualFunction()
255     fgen[idx]     = (Eqp + (Xd[i] - Xdp[i]) * Id - Efd) / Td0p[i];  in ResidualFunction()
256     fgen[idx + 1] = (Edp - (Xq[i] - Xqp[i]) * Iq) / Tq0p[i];  in ResidualFunction()
257     fgen[idx + 2] = -w + w_s;  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()
286     fgen[idx + 6] = (KE[i] * Efd + SE - VR) / TE[i];  in ResidualFunction()
287     fgen[idx + 7] = (RF - KF[i] * Efd / TF[i]) / TF[i];  in ResidualFunction()
288 …  fgen[idx + 8] = (VR - KA[i] * RF + KA[i] * KF[i] * Efd / TF[i] - KA[i] * (Vref[i] - Vm)) / TA[i];  in ResidualFunction()
315   PetscCall(VecRestoreArray(Fgen, &fgen));  in ResidualFunction()