Home
last modified time | relevance | path

Searched refs:Jacp (Results 1 – 21 of 21) sorted by relevance

/petsc/src/ts/tutorials/autodiff/
H A Dex16adj.cxx171 Mat Jacp; /* JacobianP matrix */ in main() local
213 PetscCall(MatCreate(PETSC_COMM_WORLD, &Jacp)); in main()
214 PetscCall(MatSetSizes(Jacp, PETSC_DECIDE, PETSC_DECIDE, 2, 1)); in main()
215 PetscCall(MatSetFromOptions(Jacp)); in main()
216 PetscCall(MatSetUp(Jacp)); in main()
286 PetscCall(MatCreateVecs(Jacp, &mu[0], NULL)); in main()
287 PetscCall(MatCreateVecs(Jacp, &mu[1], NULL)); in main()
297 PetscCall(TSSetRHSJacobianP(ts, Jacp, RHSJacobianP, &user)); in main()
311 PetscCall(MatDestroy(&Jacp)); in main()
H A Dex16adj_tl.cxx191 Mat Jacp; /* JacobianP matrix */ in main() local
233 PetscCall(MatCreate(PETSC_COMM_WORLD, &Jacp)); in main()
234 PetscCall(MatSetSizes(Jacp, PETSC_DECIDE, PETSC_DECIDE, 2, 1)); in main()
235 PetscCall(MatSetFromOptions(Jacp)); in main()
236 PetscCall(MatSetUp(Jacp)); in main()
296 PetscCall(MatCreateVecs(Jacp, &mu[0], NULL)); in main()
297 PetscCall(MatCreateVecs(Jacp, &mu[1], NULL)); in main()
307 PetscCall(TSSetRHSJacobianP(ts, Jacp, RHSJacobianP, &user)); in main()
321 PetscCall(MatDestroy(&Jacp)); in main()
/petsc/src/ts/tutorials/
H A Dex16fwd.c140 Mat Jacp; /* JacobianP matrix */ in main() local
175 PetscCall(MatCreate(PETSC_COMM_WORLD, &Jacp)); in main()
176 PetscCall(MatSetSizes(Jacp, PETSC_DECIDE, PETSC_DECIDE, 2, 3)); in main()
177 PetscCall(MatSetFromOptions(Jacp)); in main()
178 PetscCall(MatSetUp(Jacp)); in main()
196 PetscCall(TSSetRHSJacobianP(ts, Jacp, RHSJacobianP, &user)); in main()
230 PetscCall(MatDestroy(&Jacp)); in main()
H A Dex23fwdadj.c21 Mat Jacp; /* JacobianP matrix */ member
133 PetscCall(MatCreate(PETSC_COMM_WORLD, &user.Jacp)); in main()
134 PetscCall(MatSetSizes(user.Jacp, PETSC_DECIDE, PETSC_DECIDE, rows, cols)); in main()
135 PetscCall(MatSetFromOptions(user.Jacp)); in main()
136 PetscCall(MatSetUp(user.Jacp)); in main()
140 PetscCall(TSSetIJacobianP(ts, user.Jacp, IJacobianP, &user)); in main()
163 PetscCall(MatCreateVecs(user.Jacp, &user.mup[0], NULL)); in main()
176 PetscCall(MatDestroy(&user.Jacp)); in main()
H A Dex20adj.c63 Mat Jacp; /* IJacobianP matrix */ member
275 PetscCall(MatCreate(PETSC_COMM_WORLD, &user.Jacp)); in main()
276 PetscCall(MatSetSizes(user.Jacp, PETSC_DECIDE, PETSC_DECIDE, 2, 1)); in main()
277 PetscCall(MatSetFromOptions(user.Jacp)); in main()
278 PetscCall(MatSetUp(user.Jacp)); in main()
279 PetscCall(MatDuplicate(user.Jacp, MAT_DO_NOT_COPY_VALUES, &user.Jacprhs)); in main()
290 PetscCall(TSSetIJacobianP(ts, user.Jacp, IJacobianP, &user)); in main()
299 PetscCall(TSSetIJacobianP(ts, user.Jacp, IJacobianP, &user)); in main()
304 PetscCall(TSSetRHSJacobianP(ts, user.Jacp, RHSJacobianP, &user)); in main()
351 PetscCall(MatCreateVecs(user.Jacp, &user.mup[0], NULL)); in main()
[all …]
H A Dex20fwd.c27 Mat Jacp; /* JacobianP matrix */ member
184 PetscCall(MatCreate(PETSC_COMM_WORLD, &user.Jacp)); in main()
185 PetscCall(MatSetSizes(user.Jacp, PETSC_DECIDE, PETSC_DECIDE, rows, cols)); in main()
186 PetscCall(MatSetFromOptions(user.Jacp)); in main()
187 PetscCall(MatSetUp(user.Jacp)); in main()
196 PetscCall(TSSetRHSJacobianP(ts, user.Jacp, RHSJacobianP, &user)); in main()
223 PetscCall(MatDestroy(&user.Jacp)); in main()
H A Dex20td.c80 Mat Jacp; /* JacobianP matrix */ member
299 PetscCall(MatCreate(PETSC_COMM_WORLD, &user.Jacp)); in main()
300 if (sa == SA_TRACK) PetscCall(MatSetSizes(user.Jacp, PETSC_DECIDE, PETSC_DECIDE, 2, 2)); in main()
301 …if (sa == SA_GLOBAL) PetscCall(MatSetSizes(user.Jacp, PETSC_DECIDE, PETSC_DECIDE, 2, user.max_step… in main()
302 PetscCall(MatSetFromOptions(user.Jacp)); in main()
303 PetscCall(MatSetUp(user.Jacp)); in main()
314 if (sa == SA_TRACK) PetscCall(TSSetRHSJacobianP(ts, user.Jacp, RHSJacobianP_track, &user)); in main()
315 if (sa == SA_GLOBAL) PetscCall(TSSetRHSJacobianP(ts, user.Jacp, RHSJacobianP_global, &user)); in main()
354 PetscCall(MatCreateVecs(user.Jacp, &user.mup, NULL)); in main()
412 PetscCall(MatDestroy(&user.Jacp)); in main()
H A Dex20opt_p.c24 Mat Jacp; /* JacobianP matrix */ member
410 PetscCall(MatCreate(PETSC_COMM_WORLD, &user.Jacp)); in main()
411 PetscCall(MatSetSizes(user.Jacp, PETSC_DECIDE, PETSC_DECIDE, 2, 1)); in main()
412 PetscCall(MatSetFromOptions(user.Jacp)); in main()
413 PetscCall(MatSetUp(user.Jacp)); in main()
414 PetscCall(MatCreateVecs(user.Jacp, &user.Dir, NULL)); in main()
415 PetscCall(MatCreateVecs(user.Jacp, &user.Mup[0], NULL)); in main()
416 PetscCall(MatCreateVecs(user.Jacp, &user.Mup2[0], NULL)); in main()
417 PetscCall(MatCreateVecs(user.Jacp, &user.Ihp3[0], NULL)); in main()
418 PetscCall(MatCreateVecs(user.Jacp, &user.Ihp4[0], NULL)); in main()
[all …]
H A Dex20opt_ic.c24 Mat Jacp; /* JacobianP matrix for ODE*/ member
/petsc/src/ts/tutorials/power_grid/
H A Dex3sa.c61 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()
135 PetscCall(TSSetRHSJacobianP(ts, ctx.Jacp, RHSJacobianP, &ctx)); in main()
147 PetscCall(MatCreateVecs(ctx.Jacp, &mu[0], NULL)); in main()
247 PetscCall(MatDestroy(&ctx.Jacp)); in main()
H A Dex9adj.c182 Mat Jacp; /* Jacobian matrix */ in main() local
214 PetscCall(MatCreate(PETSC_COMM_WORLD, &Jacp)); in main()
215 PetscCall(MatSetSizes(Jacp, PETSC_DECIDE, PETSC_DECIDE, 2, 1)); in main()
216 PetscCall(MatSetFromOptions(Jacp)); in main()
217 PetscCall(MatSetUp(Jacp)); in main()
285 PetscCall(TSSetRHSJacobianP(ts, Jacp, RHSJacobianP, &ctx)); in main()
304 PetscCall(MatCreateVecs(Jacp, &mu[0], NULL)); in main()
368 PetscCall(MatDestroy(&Jacp)); in main()
H A Dex3opt.c105 PetscCall(MatCreate(PETSC_COMM_WORLD, &ctx.Jacp)); in main()
106 PetscCall(MatSetSizes(ctx.Jacp, PETSC_DECIDE, PETSC_DECIDE, 2, 1)); in main()
107 PetscCall(MatSetFromOptions(ctx.Jacp)); in main()
108 PetscCall(MatSetUp(ctx.Jacp)); in main()
123 PetscCall(TSSetRHSJacobianP(ctx.ts, ctx.Jacp, RHSJacobianP, &ctx)); in main()
127 PetscCall(MatCreateVecs(ctx.Jacp, &mu[0], NULL)); in main()
203 PetscCall(MatDestroy(&ctx.Jacp)); in main()
H A Dex9opt.c179 Mat Jacp; /* Jacobian matrix */ in main() local
231 PetscCall(MatCreate(PETSC_COMM_WORLD, &Jacp)); in main()
232 PetscCall(MatSetSizes(Jacp, PETSC_DECIDE, PETSC_DECIDE, 2, 1)); in main()
233 PetscCall(MatSetFromOptions(Jacp)); in main()
234 PetscCall(MatSetUp(Jacp)); in main()
252 PetscCall(MatCreateVecs(Jacp, &mu[0], NULL)); in main()
254 PetscCall(TSSetRHSJacobianP(ctx.ts, Jacp, RHSJacobianP, &ctx)); in main()
321 PetscCall(MatDestroy(&Jacp)); in main()
H A Dex3.h17 Mat Jacp; /* Jacobianp matrix */ member
/petsc/src/ts/tutorials/optimal_control/
H A Dex1.c11 Mat Jacp, DRDU, DRDP; member
270 PetscCall(MatCreate(PETSC_COMM_WORLD, &aircraft.Jacp)); in main()
271 PetscCall(MatSetSizes(aircraft.Jacp, PETSC_DECIDE, PETSC_DECIDE, 2, 2 * aircraft.nsteps)); in main()
272 PetscCall(MatSetFromOptions(aircraft.Jacp)); in main()
273 PetscCall(MatSetUp(aircraft.Jacp)); in main()
274 PetscCall(MatSetOption(aircraft.Jacp, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE)); in main()
285 PetscCall(TSSetRHSJacobianP(aircraft.ts, aircraft.Jacp, RHSJacobianP, &aircraft)); in main()
312 PetscCall(MatCreateVecs(aircraft.Jacp, &aircraft.Mup[0], NULL)); in main()
316 PetscCall(MatCreateVecs(aircraft.Jacp, &aircraft.rhshp3[0], NULL)); in main()
317 PetscCall(MatCreateVecs(aircraft.Jacp, &aircraft.rhshp4[0], NULL)); in main()
[all …]
/petsc/src/ts/impls/implicit/theta/
H A Dtheta.c380 …nP(ts, th->stage_time, ts->vec_sol, th->Xdot, 1. / adjoint_time_step, ts->Jacp, PETSC_FALSE)); /* … in TSAdjointStepBEuler_Private()
390 PetscCall(MatMultTranspose(ts->Jacp, VecsDeltaLam[nadj], VecsDeltaMu[nadj])); in TSAdjointStepBEuler_Private()
400 PetscCall(MatMultTranspose(ts->Jacp, VecsDeltaLam2[nadj], VecsDeltaMu2[nadj])); in TSAdjointStepBEuler_Private()
569 …>stage_time, ts->vec_sol, th->Xdot, -1. / (th->Theta * adjoint_time_step), ts->Jacp, PETSC_FALSE)); in TSAdjointStep_Theta()
572 PetscCall(MatMultTranspose(ts->Jacp, VecsDeltaLam[nadj], VecsDeltaMu[nadj])); in TSAdjointStep_Theta()
595 PetscCall(MatMultTranspose(ts->Jacp, VecsDeltaLam2[nadj], VecsDeltaMu2[nadj])); in TSAdjointStep_Theta()
604 …int_ptime, th->X0, th->Xdot, 1. / ((th->Theta - 1.0) * adjoint_time_step), ts->Jacp, PETSC_FALSE)); in TSAdjointStep_Theta()
607 PetscCall(MatMultTranspose(ts->Jacp, VecsDeltaLam[nadj], VecsDeltaMu[nadj])); in TSAdjointStep_Theta()
628 PetscCall(MatMultTranspose(ts->Jacp, VecsDeltaLam2[nadj], VecsDeltaMu2[nadj])); in TSAdjointStep_Theta()
655 …PetscCall(TSComputeIJacobianP(ts, th->stage_time, th->X, th->Xdot, th->shift, ts->Jacp, PETSC_FALS… in TSAdjointStep_Theta()
[all …]
/petsc/src/ts/interface/sensitivity/
H A Dtssen.c147 PetscCall(MatDestroy(&ts->Jacp)); in TSSetIJacobianP()
148 ts->Jacp = Amat; in TSSetIJacobianP()
189 if (Amat) *Amat = ts->Jacp; in TSGetIJacobianP()
1080 …PetscCheck(!ts->vecs_sensip || ts->Jacp || ts->Jacprhs, PetscObjectComm((PetscObject)ts), PETSC_ER… in TSAdjointSetUp()
1174 PetscCall(MatDestroy(&ts->Jacp)); in TSAdjointSetRHSJacobian()
1175 ts->Jacp = Amat; in TSAdjointSetRHSJacobian()
/petsc/src/ts/tutorials/power_grid/stability_9bus/
H A Dex9busopt.c106 Mat J, Jacp; member
1128 PetscCall(MatCreate(PETSC_COMM_WORLD, &user.Jacp)); in main()
1129 PetscCall(MatSetSizes(user.Jacp, PETSC_DECIDE, PETSC_DECIDE, user.neqs_pgrid, 3)); in main()
1130 PetscCall(MatSetFromOptions(user.Jacp)); in main()
1131 PetscCall(MatSetUp(user.Jacp)); in main()
1132 PetscCall(MatZeroEntries(user.Jacp)); /* initialize to zeros */ in main()
1194 PetscCall(MatDestroy(&user.Jacp)); in main()
1261 PetscCall(TSSetRHSJacobianP(ts, ctx->Jacp, RHSJacobianP, ctx)); in FormFunctionGradient()
1388 PetscCall(MatCreateVecs(ctx->Jacp, &mu[0], NULL)); in FormFunctionGradient()
/petsc/include/petsc/private/
H A Dtsimpl.h209 Mat Jacp, Jacprhs; member
/petsc/src/ts/impls/arkimex/
H A Darkimex.c1591 …anP(ts, ark->stage_time, Y[i], Ydot, ark->scoeff / adjoint_time_step, ts->Jacp, PETSC_TRUE)); // g… in TSAdjointStep_ARKIMEX()
1612 PetscCall(MatMultTranspose(ts->Jacp, VecsSensiTemp[nadj], VecsSensiPTemp[nadj])); in TSAdjointStep_ARKIMEX()
1655 PetscCall(MatMultTranspose(ts->Jacp, VecsDeltaLam[nadj * s + i], VecsSensiPTemp[nadj])); in TSAdjointStep_ARKIMEX()
/petsc/src/ts/interface/
H A Dts.c2508 …if (!ts->Jacp && ts->Jacprhs) { /* IJacobianP shares the same matrix with RHSJacobianP if only RHS… in TSSetUp()
2510 ts->Jacp = ts->Jacprhs; in TSSetUp()
2609 PetscCall(MatDestroy(&ts->Jacp)); in TSReset()