| /petsc/src/ts/impls/implicit/theta/ |
| H A D | theta.c | 125 TS quadts = ts->quadraturets; in TSThetaEvaluateCostIntegral() local
131 PetscCall(TSComputeRHSFunction(quadts, th->ptime0, th->X0, ts->vec_costintegrand)); in TSThetaEvaluateCostIntegral()
132 … PetscCall(VecAXPY(quadts->vec_sol, th->time_step0 * (1.0 - th->Theta), ts->vec_costintegrand)); in TSThetaEvaluateCostIntegral()
134 PetscCall(TSComputeRHSFunction(quadts, ts->ptime, ts->vec_sol, ts->vec_costintegrand)); in TSThetaEvaluateCostIntegral()
135 PetscCall(VecAXPY(quadts->vec_sol, th->time_step0 * th->Theta, ts->vec_costintegrand)); in TSThetaEvaluateCostIntegral()
137 PetscCall(TSComputeRHSFunction(quadts, th->stage_time, th->X, ts->vec_costintegrand)); in TSThetaEvaluateCostIntegral()
138 PetscCall(VecAXPY(quadts->vec_sol, th->time_step0, ts->vec_costintegrand)); in TSThetaEvaluateCostIntegral()
146 TS quadts = ts->quadraturets; in TSForwardCostIntegral_Theta() local
150 PetscCall(VecCopy(quadts->vec_sol, th->VecCostIntegral0)); in TSForwardCostIntegral_Theta()
274 TS quadts = ts->quadraturets; in TSAdjointStepBEuler_Private() local
[all …]
|
| /petsc/src/ts/tutorials/power_grid/ |
| H A D | ex3opt.c | 130 PetscCall(TSCreateQuadratureTS(ctx.ts, PETSC_FALSE, &ctx.quadts)); in main()
131 PetscCall(TSSetRHSFunction(ctx.quadts, NULL, (TSRHSFunctionFn *)CostIntegrand, &ctx)); in main()
132 …PetscCall(TSSetRHSJacobian(ctx.quadts, ctx.DRDU, ctx.DRDU, (TSRHSJacobianFn *)DRDUJacobianTranspos… in main()
133 PetscCall(TSSetRHSJacobianP(ctx.quadts, ctx.DRDP, DRDPJacobianTranspose, &ctx)); in main()
139 PetscCall(TSCreateQuadratureTS(ctx.ts, PETSC_TRUE, &ctx.quadts)); in main()
140 PetscCall(TSForwardSetSensitivities(ctx.quadts, 1, qgrad)); in main()
141 PetscCall(TSSetRHSFunction(ctx.quadts, NULL, (TSRHSFunctionFn *)CostIntegrand, &ctx)); in main()
142 …PetscCall(TSSetRHSJacobian(ctx.quadts, ctx.DRDU, ctx.DRDU, (TSRHSJacobianFn *)DRDUJacobianTranspos… in main()
143 PetscCall(TSSetRHSJacobianP(ctx.quadts, ctx.DRDP, DRDPJacobianTranspose, &ctx)); in main()
270 PetscCall(TSGetQuadratureTS(ctx->ts, NULL, &ctx->quadts)); in FormFunctionGradient()
[all …]
|
| H A D | ex3sa.c | 26 TS ts, quadts; /* ODE integrator */ in main() local
137 PetscCall(TSCreateQuadratureTS(ts, PETSC_FALSE, &quadts)); 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()
158 PetscCall(TSForwardSetSensitivities(quadts, 1, qgrad)); in main()
|
| H A D | ex3opt_fd.c | 162 TS ts, quadts; in FormFunction() local
213 PetscCall(TSCreateQuadratureTS(ts, PETSC_TRUE, &quadts)); in FormFunction()
214 PetscCall(TSGetSolution(quadts, &q)); in FormFunction()
216 PetscCall(TSSetRHSFunction(quadts, NULL, (TSRHSFunctionFn *)CostIntegrand, ctx)); in FormFunction()
|
| H A D | ex9adj.c | 179 TS ts, quadts; /* ODE integrator */ in main() local
280 PetscCall(TSCreateQuadratureTS(ts, PETSC_TRUE, &quadts)); in main()
281 PetscCall(TSSetRHSFunction(quadts, NULL, (TSRHSFunctionFn *)CostIntegrand, &ctx)); in main()
282 PetscCall(TSSetRHSJacobian(quadts, DRDU, DRDU, (TSRHSJacobianFn *)DRDUJacobianTranspose, &ctx)); in main()
283 PetscCall(TSSetRHSJacobianP(quadts, DRDP, (TSRHSJacobianPFn *)DRDPJacobianTranspose, &ctx)); in main()
|
| H A D | ex9opt.c | 182 TS quadts; in main() local
265 PetscCall(TSCreateQuadratureTS(ctx.ts, PETSC_TRUE, &quadts)); in main()
266 PetscCall(TSSetRHSFunction(quadts, NULL, (TSRHSFunctionFn *)CostIntegrand, &ctx)); in main()
267 PetscCall(TSSetRHSJacobian(quadts, DRDU, DRDU, (TSRHSJacobianFn *)DRDUJacobianTranspose, &ctx)); in main()
268 PetscCall(TSSetRHSJacobianP(quadts, DRDP, (TSRHSJacobianPFn *)DRDPJacobianTranspose, &ctx)); in main()
|
| H A D | ex3.h | 14 TS ts, quadts; member
|
| /petsc/src/ts/impls/explicit/rk/ |
| H A D | rk.c | 599 TS quadts = ts->quadraturets; in TSForwardCostIntegral_RK() local
610 …PetscCall(TSComputeRHSFunction(quadts, rk->ptime + rk->time_step * c[i], Y[i], ts->vec_costintegra… in TSForwardCostIntegral_RK()
611 PetscCall(VecAXPY(quadts->vec_sol, rk->time_step * b[i], ts->vec_costintegrand)); in TSForwardCostIntegral_RK()
620 TS quadts = ts->quadraturets; in TSAdjointCostIntegral_RK() local
629 …PetscCall(TSComputeRHSFunction(quadts, ts->ptime + ts->time_step * (1.0 - c[i]), Y[i], ts->vec_cos… in TSAdjointCostIntegral_RK()
630 PetscCall(VecAXPY(quadts->vec_sol, -ts->time_step * b[i], ts->vec_costintegrand)); in TSAdjointCostIntegral_RK()
638 TS quadts = ts->quadraturets; in TSRollBack_RK() local
661 if (quadts && ts->costintegralfwd) { in TSRollBack_RK()
664 PetscCall(TSComputeRHSFunction(quadts, rk->ptime + h * c[j], Y[j], ts->vec_costintegrand)); in TSRollBack_RK()
665 PetscCall(VecAXPY(quadts->vec_sol, -h * b[j], ts->vec_costintegrand)); in TSRollBack_RK()
[all …]
|
| /petsc/src/ts/interface/sensitivity/ |
| H A D | tssen.c | 345 TS quadts; in TSGetCostIntegral() local
350 PetscCall(TSGetQuadratureTS(ts, NULL, &quadts)); in TSGetCostIntegral()
351 *v = quadts->vec_sol; in TSGetCostIntegral()
1708 TS quadts = ts->quadraturets; in TSForwardReset() local
1714 if (quadts) PetscCall(MatDestroy(&quadts->mat_sensip)); in TSForwardReset()
1940 PetscErrorCode TSCreateQuadratureTS(TS ts, PetscBool fwd, TS *quadts) in TSCreateQuadratureTS() argument
1946 PetscAssertPointer(quadts, 3); in TSCreateQuadratureTS()
1952 *quadts = ts->quadraturets; in TSCreateQuadratureTS()
1955 PetscCall(VecCreateSeq(PETSC_COMM_SELF, ts->numcost, &(*quadts)->vec_sol)); in TSCreateQuadratureTS()
1957 PetscCall(VecCreateSeq(PETSC_COMM_SELF, 1, &(*quadts)->vec_sol)); in TSCreateQuadratureTS()
[all …]
|
| /petsc/src/ts/tutorials/optimal_control/ |
| H A D | ex1.c | 6 TS ts, quadts; member
307 PetscCall(TSCreateQuadratureTS(aircraft.ts, PETSC_FALSE, &aircraft.quadts)); in main()
308 PetscCall(TSSetRHSFunction(aircraft.quadts, NULL, (TSRHSFunctionFn *)CostIntegrand, &aircraft)); in main()
309 …PetscCall(TSSetRHSJacobian(aircraft.quadts, aircraft.DRDU, aircraft.DRDU, (TSRHSJacobian)DRDUJacob… in main()
310 …PetscCall(TSSetRHSJacobianP(aircraft.quadts, aircraft.DRDP, (TSRHSJacobianPFn *)DRDPJacobianTransp… in main()
324 …PetscCall(TSSetRHSHessianProduct(aircraft.quadts, aircraft.inthp1, IntegrandHessianProductUU, airc… in main()
|
| /petsc/src/ts/tutorials/power_grid/stability_9bus/ |
| H A D | ex9busopt.c | 1222 TS ts, quadts; in FormFunctionGradient() local
1263 PetscCall(TSCreateQuadratureTS(ts, PETSC_FALSE, &quadts)); in FormFunctionGradient()
1264 PetscCall(TSSetRHSFunction(quadts, NULL, (TSRHSFunctionFn *)CostIntegrand, ctx)); in FormFunctionGradient()
1265 …PetscCall(TSSetRHSJacobian(quadts, ctx->DRDU, ctx->DRDU, (TSRHSJacobianFn *)DRDUJacobianTranspose,… in FormFunctionGradient()
1266 PetscCall(TSSetRHSJacobianP(quadts, ctx->DRDP, (TSRHSJacobianPFn *)DRDPJacobianTranspose, ctx)); in FormFunctionGradient()
1436 PetscCall(TSGetQuadratureTS(ts, NULL, &quadts)); in FormFunctionGradient()
1437 PetscCall(TSGetSolution(quadts, &q)); in FormFunctionGradient()
|
| /petsc/doc/manual/ |
| H A D | ts.md | 1193 TSCreateQuadratureTS(TS ts,PetscBool fwd,TS *quadts);
1200 TSSetRHSFunction(TS quadts,Vec R,PetscErrorCode (*rf)(TS,PetscReal,Vec,Vec,PetscCtx),PetscCtxctx)
1207 TSSetRHSJacobian(TS quadts,Vec DRDU,Vec DRDU,PetscErrorCode (*drdyf)(TS,PetscReal,Vec,Vec*,PetscCtx…
1208 TSSetRHSJacobianP(TS quadts,Vec DRDU,Vec DRDU,PetscErrorCode (*drdyp)(TS,PetscReal,Vec,Vec*,PetscCt…
1227 or accessing directly the solution vector used by `quadts`.
|