434   PetscScalar       *fgen, *fnet;  in ResidualFunction()  local
465   PetscCall(VecGetArrayWrite(Fgen, &fgen));  in ResidualFunction()
481     fgen[idx]     = (-Eqp - (Xd[i] - Xdp[i]) * Id + Efd) / Td0p[i];  in ResidualFunction()
482     fgen[idx + 1] = (-Edp + (Xq[i] - Xqp[i]) * Iq) / Tq0p[i];  in ResidualFunction()
483     fgen[idx + 2] = w - w_s;  in ResidualFunction()
484 …fgen[idx + 3] = (TM[i] - Edp * Id - Eqp * Iq - (Xqp[i] - Xdp[i]) * Id * Iq - D[i] * (w - w_s)) / M…  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()
512     fgen[idx + 6] = (-KE[i] * Efd - SE + VR) / TE[i];  in ResidualFunction()
513     fgen[idx + 7] = (-RF + KF[i] * Efd / TF[i]) / TF[i];  in ResidualFunction()
514     if (VRatmax[i]) fgen[idx + 8] = VR - VRMAX[i];  in ResidualFunction()
515     else if (VRatmin[i]) fgen[idx + 8] = VRMIN[i] - VR;  in ResidualFunction()
516 …else fgen[idx + 8] = (-VR + KA[i] * RF - KA[i] * KF[i] * Efd / TF[i] + KA[i] * (Vref[i] - Vm)) / T…  in ResidualFunction()
543   PetscCall(VecRestoreArrayWrite(Fgen, &fgen));  in ResidualFunction()