| /petsc/src/ts/tutorials/power_grid/ |
| H A D | ex3.h | 10 PetscScalar H, D, omega_b, omega_s, Pmax, Pmax_ini, Pm, E, V, X, u_s, c; member
64 f[0] = ctx->omega_b * (u[1] - ctx->omega_s); 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()
86 J[1][0] = -ctx->omega_s / (2.0 * ctx->H) * Pmax * PetscCosScalar(u[0]); in RHSJacobian()
87 J[1][1] = -ctx->omega_s / (2.0 * ctx->H) * ctx->D; in RHSJacobian()
114 f[0] = udot[0] - ctx->omega_b * (u[1] - ctx->omega_s); 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()
138 J[1][1] = 2.0 * ctx->H / ctx->omega_s * a + ctx->D; in IJacobian()
163 J[1][0] = ctx->omega_s / (2.0 * ctx->H); in RHSJacobianP()
|
| H A D | ex1.c | 25 PetscScalar H, omega_s, E, V, X; member
47 …f[0] = 2.0 * ctx->H * udot[0] / ctx->omega_s + ctx->E * ctx->V * PetscSinScalar(u[1]) / ctx->X - R; in IFunction()
48 f[1] = udot[1] - u[0] + ctx->omega_s; in IFunction()
68 J[0][0] = 2.0 * ctx->H * a / ctx->omega_s; in IJacobian()
118 ctx.omega_s = 1.0; in main()
119 PetscCall(PetscOptionsScalar("-omega_s", "", "", ctx.omega_s, &ctx.omega_s, NULL)); in main()
|
| H A D | ex9.c | 34 PetscScalar H, D, omega_b, omega_s, Pmax, Pm, E, V, X; member
53 f[0] = ctx->omega_b * (u[1] - ctx->omega_s); in RHSFunction()
54 …f[1] = (-Pmax * PetscSinScalar(u[0]) - ctx->D * (u[1] - ctx->omega_s) + ctx->Pm) * ctx->omega_s / … in RHSFunction()
77 J[1][1] = -ctx->D * ctx->omega_s / (2.0 * ctx->H); in RHSJacobian()
78 J[1][0] = -Pmax * PetscCosScalar(u[0]) * ctx->omega_s / (2.0 * ctx->H); in RHSJacobian()
129 ctx.omega_s = 2.0 * PETSC_PI * 60.0; in main()
196 u[1] = ctx.omega_s; in main()
|
| H A D | ex9adj.c | 34 PetscScalar H, D, omega_b, omega_s, Pmax, Pm, E, V, X, u_s, c; member
69 f[0] = ctx->omega_b * (u[1] - ctx->omega_s); in RHSFunction()
70 …f[1] = (-Pmax * PetscSinScalar(u[0]) - ctx->D * (u[1] - ctx->omega_s) + ctx->Pm) * ctx->omega_s / … in RHSFunction()
93 J[1][1] = -ctx->D * ctx->omega_s / (2.0 * ctx->H); in RHSJacobian()
94 J[1][0] = -Pmax * PetscCosScalar(u[0]) * ctx->omega_s / (2.0 * ctx->H); in RHSJacobian()
116 J[1][0] = ctx->omega_s / (2.0 * ctx->H); in RHSJacobianP()
232 ctx.omega_s = 1.0; in main()
324 u[1] = ctx.omega_s; in main()
|
| H A D | ex2.c | 34 PetscScalar H, D, omega_s, Pmax, Pm, E, V, X; member
54 f[0] = udot[0] - ctx->omega_s * (u[1] - 1.0); in IFunction()
80 J[0][1] = -ctx->omega_s; in IJacobian()
149 ctx.omega_s = 2.0 * PETSC_PI * 60.0; in main()
217 u[1] = ctx.omega_s; in main()
|
| H A D | ex9opt.c | 36 PetscScalar H, D, omega_b, omega_s, Pmax, Pm, E, V, X, u_s, c; member
59 f[0] = ctx->omega_b * (u[1] - ctx->omega_s); in RHSFunction()
60 …f[1] = (-Pmax * PetscSinScalar(u[0]) - ctx->D * (u[1] - ctx->omega_s) + ctx->Pm) * ctx->omega_s / … in RHSFunction()
83 J[1][1] = -ctx->D * ctx->omega_s / (2.0 * ctx->H); in RHSJacobian()
84 J[1][0] = -Pmax * PetscCosScalar(u[0]) * ctx->omega_s / (2.0 * ctx->H); in RHSJacobian()
106 J[1][0] = ctx->omega_s / (2.0 * ctx->H); in RHSJacobianP()
200 ctx.omega_s = PetscRealConstant(1.0); in main()
|
| H A D | ex3.c | 73 ctx.omega_s = 2.0 * PETSC_PI * 60.0; in main()
148 u[1] = ctx.omega_s; in main()
|
| H A D | ex3sa.c | 78 ctx.omega_s = 1.0; in main()
186 u[1] = ctx.omega_s; in main()
|
| H A D | ex3opt_fd.c | 72 ctx.omega_s = 1.0; in main()
|
| H A D | ex3opt.c | 72 ctx.omega_s = 1.0; in main()
|