10   PetscScalar H, D, omega_b, omega_s, Pmax, Pmax_ini, Pm, E, V, X, u_s, c;  member
41     if (forwardsolve) user->Pmax = 0.0; /* Apply disturbance - this is done by setting Pmax = 0 */  in PostEventFunction()
42     else user->Pmax = user->Pmax_ini;   /* Going backward, reversal of event */  in PostEventFunction()
44 …if (forwardsolve) user->Pmax = user->Pmax_ini; /* Remove the fault  - this is done by setting Pmax…  in PostEventFunction()
45     else user->Pmax = 0.0;                         /* Going backward, reversal of event */  in PostEventFunction()
56   PetscScalar       *f, Pmax;  in RHSFunction()  local
63   Pmax = ctx->Pmax;  in RHSFunction()
65 …f[1] = ctx->omega_s / (2.0 * ctx->H) * (ctx->Pm - Pmax * PetscSinScalar(u[0]) - ctx->D * (u[1] - c…  in RHSFunction()
78   PetscScalar        J[2][2], Pmax;  in RHSJacobian()  local
83   Pmax    = ctx->Pmax;  in RHSJacobian()
86   J[1][0] = -ctx->omega_s / (2.0 * ctx->H) * Pmax * PetscCosScalar(u[0]);  in RHSJacobian()
105   PetscScalar       *f, Pmax;  in IFunction()  local
113   Pmax = ctx->Pmax;  in IFunction()
115 …f[1] = 2.0 * ctx->H / ctx->omega_s * udot[1] + Pmax * PetscSinScalar(u[0]) + ctx->D * (u[1] - ctx-…  in IFunction()
129   PetscScalar        J[2][2], Pmax;  in IJacobian()  local
135   Pmax    = ctx->Pmax;  in IJacobian()
139   J[1][0] = Pmax * PetscCosScalar(u[0]);  in IJacobian()
220 …sensip = 1. / PetscSqrtScalar(1. - (ctx->Pm / ctx->Pmax) * (ctx->Pm / ctx->Pmax)) / ctx->Pmax * x[…  in ComputeSensiP()