| /petsc/src/ts/tutorials/power_grid/ |
| H A D | ex3sa.c | 90 ctx.Pm = 1.1; in main()
91 PetscCall(PetscOptionsScalar("-Pm", "", "", ctx.Pm, &ctx.Pm, NULL)); in main()
103 u[0] = PetscAsinScalar(ctx.Pm / ctx.Pmax); in main()
159 val[0] = 1. / PetscSqrtScalar(1. - (ctx.Pm / ctx.Pmax) * (ctx.Pm / ctx.Pmax)) / ctx.Pmax; in main()
185 u[0] = PetscAsinScalar(ctx.Pm / ctx.Pmax); in main()
221 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "\n cost function=%g\n", (double)(x_ptr[0] - ctx.Pm))); in main()
235 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "\n cost function=%g\n", (double)(s_ptr[0] - ctx.Pm))); in main()
|
| H A D | ex9.c | 34 PetscScalar H, D, omega_b, omega_s, Pmax, Pm, E, V, X; member
54 …f[1] = (-Pmax * PetscSinScalar(u[0]) - ctx->D * (u[1] - ctx->omega_s) + ctx->Pm) * ctx->omega_s / … in RHSFunction()
139 ctx.Pm = 0.9; in main()
140 PetscCall(PetscOptionsScalar("-Pm", "", "", ctx.Pm, &ctx.Pm, NULL)); in main()
152 u[0] = PetscAsinScalar(ctx.Pm / ctx.Pmax); in main()
195 u[0] = PetscAsinScalar(ctx.Pm / ctx.Pmax); in main()
|
| H A D | ex3.c | 84 ctx.Pm = 0.9; in main()
85 PetscCall(PetscOptionsScalar("-Pm", "", "", ctx.Pm, &ctx.Pm, NULL)); in main()
97 u[0] = PetscAsinScalar(ctx.Pm / ctx.Pmax); in main()
147 u[0] = PetscAsinScalar(ctx.Pm / ctx.Pmax); in main()
|
| H A D | ex2.c | 34 PetscScalar H, D, omega_s, Pmax, Pm, E, V, X; member
55 f[1] = 2.0 * ctx->H * udot[1] + Pmax * PetscSinScalar(u[0]) + ctx->D * (u[1] - 1.0) - ctx->Pm; in IFunction()
159 ctx.Pm = 0.9; in main()
160 PetscCall(PetscOptionsScalar("-Pm", "", "", ctx.Pm, &ctx.Pm, NULL)); in main()
172 u[0] = PetscAsinScalar(ctx.Pm / ctx.Pmax); in main()
216 u[0] = PetscAsinScalar(ctx.Pm / ctx.Pmax); in main()
|
| H A D | ex9opt.c | 36 PetscScalar H, D, omega_b, omega_s, Pmax, Pm, E, V, X, u_s, c; member
60 …f[1] = (-Pmax * PetscSinScalar(u[0]) - ctx->D * (u[1] - ctx->omega_s) + ctx->Pm) * ctx->omega_s / … in RHSFunction()
160 …sensip = 1. / PetscSqrtScalar(1. - (ctx->Pm / ctx->Pmax) * (ctx->Pm / ctx->Pmax)) / ctx->Pmax * x[… in ComputeSensiP()
211 ctx.Pm = PetscRealConstant(1.0194); in main()
212 PetscCall(PetscOptionsScalar("-Pm", "", "", ctx.Pm, &ctx.Pm, NULL)); in main()
281 x_ptr[0] = ctx.Pm; in main()
353 ctx->Pm = x_ptr[0]; in FormFunction()
371 u[0] = PetscAsinScalar(ctx->Pm / ctx->Pmax); in FormFunction()
381 *f = -ctx->Pm + x_ptr[0]; in FormFunction()
399 ctx->Pm = x_ptr[0]; in FormGradient()
[all …]
|
| H A D | ex3opt_fd.c | 84 ctx.Pm = 1.06; in main()
85 PetscCall(PetscOptionsScalar("-Pm", "", "", ctx.Pm, &ctx.Pm, NULL)); in main()
106 x_ptr[0] = ctx.Pm; in main()
176 ctx->Pm = x_ptr[0]; in FormFunction()
202 u[0] = PetscAsinScalar(ctx->Pm / ctx->Pmax); in FormFunction()
232 *f = -ctx->Pm + qx_ptr[0]; in FormFunction()
|
| H A D | ex9adj.c | 34 PetscScalar H, D, omega_b, omega_s, Pmax, Pm, E, V, X, u_s, c; member
70 …f[1] = (-Pmax * PetscSinScalar(u[0]) - ctx->D * (u[1] - ctx->omega_s) + ctx->Pm) * ctx->omega_s / … in RHSFunction()
170 …sensip = 1. / PetscSqrtScalar(1. - (ctx->Pm / ctx->Pmax) * (ctx->Pm / ctx->Pmax)) / ctx->Pmax * x[… in ComputeSensiP()
243 ctx.Pm = 1.1; in main()
244 PetscCall(PetscOptionsScalar("-Pm", "", "", ctx.Pm, &ctx.Pm, NULL)); in main()
256 u[0] = PetscAsinScalar(ctx.Pm / ctx.Pmax); in main()
323 u[0] = PetscAsinScalar(ctx.Pm / ctx.Pmax); in main()
359 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "\n cost function=%g\n", (double)(x_ptr[0] - ctx.Pm))); in main()
|
| H A D | ex3opt.c | 84 ctx.Pm = 1.06; in main()
85 PetscCall(PetscOptionsScalar("-Pm", "", "", ctx.Pm, &ctx.Pm, NULL)); in main()
168 x_ptr[0] = ctx.Pm; in main()
250 ctx->Pm = x_ptr[0]; in FormFunctionGradient()
255 u[0] = PetscAsinScalar(ctx->Pm / ctx->Pmax); in FormFunctionGradient()
284 val[0] = 1. / PetscSqrtScalar(1. - (ctx->Pm / ctx->Pmax) * (ctx->Pm / ctx->Pmax)) / ctx->Pmax; in FormFunctionGradient()
325 *f = -ctx->Pm + x_ptr[0]; in FormFunctionGradient()
|
| H A D | ex3.h | 10 PetscScalar H, D, omega_b, omega_s, Pmax, Pmax_ini, Pm, E, V, X, u_s, c; member
65 …f[1] = ctx->omega_s / (2.0 * ctx->H) * (ctx->Pm - Pmax * PetscSinScalar(u[0]) - ctx->D * (u[1] - c… in RHSFunction()
115 … / ctx->omega_s * udot[1] + Pmax * PetscSinScalar(u[0]) + ctx->D * (u[1] - ctx->omega_s) - ctx->Pm; in IFunction()
220 …sensip = 1. / PetscSqrtScalar(1. - (ctx->Pm / ctx->Pmax) * (ctx->Pm / ctx->Pmax)) / ctx->Pmax * x[… in ComputeSensiP()
|