Searched refs:Jacprhs (Results 1 – 6 of 6) sorted by relevance
| /petsc/src/ts/tutorials/ |
| H A D | ex20adj.c | 64 Mat Jacprhs; /* RHSJacobianP matrix */ member
279 PetscCall(MatDuplicate(user.Jacp, MAT_DO_NOT_COPY_VALUES, &user.Jacprhs)); in main()
280 PetscCall(MatZeroEntries(user.Jacprhs)); in main()
293 PetscCall(TSSetRHSJacobianP(ts, user.Jacprhs, NULL, &user)); in main()
390 PetscCall(MatDestroy(&user.Jacprhs)); in main()
|
| /petsc/src/ts/interface/sensitivity/ |
| H A D | tssen.c | 38 PetscCall(MatDestroy(&ts->Jacprhs)); in TSSetRHSJacobianP()
39 ts->Jacprhs = Amat; in TSSetRHSJacobianP()
69 if (Amat) *Amat = ts->Jacprhs; in TSGetRHSJacobianP()
234 if (ts->rhsjacobianp) PetscCall(TSComputeRHSJacobianP(ts, t, U, ts->Jacprhs)); in TSComputeIJacobianP()
235 if (ts->Jacprhs == Amat) { /* No IJacobian, so we only have the RHS matrix */ in TSComputeIJacobianP()
237 } else if (ts->Jacprhs) { /* Both IJacobian and RHSJacobian */ in TSComputeIJacobianP()
242 PetscCall(MatAXPY(Amat, -1, ts->Jacprhs, axpy)); in TSComputeIJacobianP()
1080 …PetscCheck(!ts->vecs_sensip || ts->Jacp || ts->Jacprhs, PetscObjectComm((PetscObject)ts), PETSC_ER… in TSAdjointSetUp()
|
| /petsc/src/ts/impls/explicit/rk/ |
| H A D | rk.c | 704 if (ts->Jacprhs) { in TSForwardStep_RK()
705 PetscCall(TSComputeRHSJacobianP(ts, stage_time, Y[i], ts->Jacprhs)); /* get f_p */ in TSForwardStep_RK()
710 …PetscCall(MatMultAdd(ts->Jacprhs, ts->vec_dir, rk->VecDeltaFwdSensipCol, rk->VecDeltaFwdSensipCol)… in TSForwardStep_RK()
714 PetscCall(MatAXPY(MatsFwdSensipTemp[i], 1., ts->Jacprhs, SUBSET_NONZERO_PATTERN)); in TSForwardStep_RK()
901 …PetscCall(TSComputeRHSJacobianP(ts, rk->stage_time, Y[i], ts->Jacprhs)); /* ge… in TSAdjointStep_RK()
902 …if (quadts) PetscCall(TSComputeRHSJacobianP(quadts, rk->stage_time, Y[i], quadts->Jacprhs)); /* ge… in TSAdjointStep_RK()
936 PetscCall(MatMultTranspose(ts->Jacprhs, VecsSensiTemp[nadj], VecDeltaMu)); in TSAdjointStep_RK()
940 PetscCall(MatDenseGetColumn(quadts->Jacprhs, nadj, &xarr)); in TSAdjointStep_RK()
944 PetscCall(MatDenseRestoreColumn(quadts->Jacprhs, &xarr)); in TSAdjointStep_RK()
1006 PetscCall(MatMultTranspose(ts->Jacprhs, VecsSensi2Temp[nadj], VecDeltaMu2)); in TSAdjointStep_RK()
|
| /petsc/include/petsc/private/ |
| H A D | tsimpl.h | 209 Mat Jacp, Jacprhs; member
|
| /petsc/src/ts/interface/ |
| H A D | ts.c | 2508 …if (!ts->Jacp && ts->Jacprhs) { /* IJacobianP shares the same matrix with RHSJacobianP if only RHS… in TSSetUp()
2509 PetscCall(PetscObjectReference((PetscObject)ts->Jacprhs)); in TSSetUp()
2510 ts->Jacp = ts->Jacprhs; in TSSetUp()
2608 PetscCall(MatDestroy(&ts->Jacprhs)); in TSReset()
|
| /petsc/src/ts/impls/arkimex/ |
| H A D | arkimex.c | 1592 …PetscCall(TSComputeRHSJacobianP(ts, stage_time_ex, Y[i], ts->Jacprhs)); … in TSAdjointStep_ARKIMEX()
1633 PetscCall(MatMultTranspose(ts->Jacprhs, VecsSensiTemp[nadj], VecsSensiPTemp[nadj])); in TSAdjointStep_ARKIMEX()
|