Lines Matching refs:ctx
30 AppCtx ctx; in main() local
55 PetscCall(MatCreate(PETSC_COMM_WORLD, &ctx.Jac)); in main()
56 PetscCall(MatSetSizes(ctx.Jac, n, n, PETSC_DETERMINE, PETSC_DETERMINE)); in main()
57 PetscCall(MatSetType(ctx.Jac, MATDENSE)); in main()
58 PetscCall(MatSetFromOptions(ctx.Jac)); in main()
59 PetscCall(MatSetUp(ctx.Jac)); in main()
60 PetscCall(MatCreateVecs(ctx.Jac, &U, NULL)); in main()
61 PetscCall(MatCreate(PETSC_COMM_WORLD, &ctx.Jacp)); in main()
62 PetscCall(MatSetSizes(ctx.Jacp, PETSC_DECIDE, PETSC_DECIDE, 2, 1)); in main()
63 PetscCall(MatSetFromOptions(ctx.Jacp)); in main()
64 PetscCall(MatSetUp(ctx.Jacp)); in main()
65 PetscCall(MatCreateDense(PETSC_COMM_WORLD, PETSC_DECIDE, PETSC_DECIDE, 1, 1, NULL, &ctx.DRDP)); in main()
66 PetscCall(MatSetUp(ctx.DRDP)); in main()
67 PetscCall(MatCreateDense(PETSC_COMM_WORLD, PETSC_DECIDE, PETSC_DECIDE, 2, 1, NULL, &ctx.DRDU)); in main()
68 PetscCall(MatSetUp(ctx.DRDU)); in main()
75 ctx.beta = 2; in main()
76 ctx.c = 10000.0; in main()
77 ctx.u_s = 1.0; in main()
78 ctx.omega_s = 1.0; in main()
79 ctx.omega_b = 120.0 * PETSC_PI; in main()
80 ctx.H = 5.0; in main()
81 PetscCall(PetscOptionsScalar("-Inertia", "", "", ctx.H, &ctx.H, NULL)); in main()
82 ctx.D = 5.0; in main()
83 PetscCall(PetscOptionsScalar("-D", "", "", ctx.D, &ctx.D, NULL)); in main()
84 ctx.E = 1.1378; in main()
85 ctx.V = 1.0; in main()
86 ctx.X = 0.545; in main()
87 ctx.Pmax = ctx.E * ctx.V / ctx.X; in main()
88 ctx.Pmax_ini = ctx.Pmax; in main()
89 PetscCall(PetscOptionsScalar("-Pmax", "", "", ctx.Pmax, &ctx.Pmax, NULL)); in main()
90 ctx.Pm = 1.1; in main()
91 PetscCall(PetscOptionsScalar("-Pm", "", "", ctx.Pm, &ctx.Pm, NULL)); in main()
92 ctx.tf = 0.1; in main()
93 ctx.tcl = 0.2; in main()
94 PetscCall(PetscOptionsReal("-tf", "Time to start fault", "", ctx.tf, &ctx.tf, NULL)); in main()
95 PetscCall(PetscOptionsReal("-tcl", "Time to end fault", "", ctx.tcl, &ctx.tcl, NULL)); in main()
98 ctx.tf = -1; in main()
99 ctx.tcl = -1; in main()
103 u[0] = PetscAsinScalar(ctx.Pm / ctx.Pmax); in main()
112 ctx.tf = -1; in main()
113 ctx.tcl = -1; in main()
126 PetscCall(TSSetRHSFunction(ts, NULL, (TSRHSFunctionFn *)RHSFunction, &ctx)); in main()
127 PetscCall(TSSetRHSJacobian(ts, ctx.Jac, ctx.Jac, (TSRHSJacobianFn *)RHSJacobian, &ctx)); in main()
135 PetscCall(TSSetRHSJacobianP(ts, ctx.Jacp, RHSJacobianP, &ctx)); in main()
138 PetscCall(TSSetRHSFunction(quadts, NULL, (TSRHSFunctionFn *)CostIntegrand, &ctx)); in main()
139 …PetscCall(TSSetRHSJacobian(quadts, ctx.DRDU, ctx.DRDU, (TSRHSJacobianFn *)DRDUJacobianTranspose, &… in main()
140 PetscCall(TSSetRHSJacobianP(quadts, ctx.DRDP, DRDPJacobianTranspose, &ctx)); in main()
146 PetscCall(MatCreateVecs(ctx.Jac, &lambda[0], NULL)); in main()
147 PetscCall(MatCreateVecs(ctx.Jacp, &mu[0], NULL)); in main()
159 val[0] = 1. / PetscSqrtScalar(1. - (ctx.Pm / ctx.Pmax) * (ctx.Pm / ctx.Pmax)) / ctx.Pmax; in main()
177 …ll(TSSetEventHandler(ts, 2, direction, terminate, EventFunction, PostEventFunction, (void *)&ctx)); in main()
185 u[0] = PetscAsinScalar(ctx.Pm / ctx.Pmax); in main()
186 u[1] = ctx.omega_s; in main()
221 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "\n cost function=%g\n", (double)(x_ptr[0] - ctx.Pm))); in main()
223 PetscCall(ComputeSensiP(lambda[0], mu[0], &ctx)); in main()
235 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "\n cost function=%g\n", (double)(s_ptr[0] - ctx.Pm))); in main()
246 PetscCall(MatDestroy(&ctx.Jac)); in main()
247 PetscCall(MatDestroy(&ctx.Jacp)); in main()
248 PetscCall(MatDestroy(&ctx.DRDU)); in main()
249 PetscCall(MatDestroy(&ctx.DRDP)); in main()