Searched refs:vecs_sensip (Results 1 – 5 of 5) sorted by relevance
| /petsc/src/ts/interface/sensitivity/ |
| H A D | tssen.c | 883 ts->vecs_sensip = mu; in TSSetCostGradients()
912 if (mu) *mu = ts->vecs_sensip; in TSGetCostGradients()
1080 …PetscCheck(!ts->vecs_sensip || ts->Jacp || ts->Jacprhs, PetscObjectComm((PetscObject)ts), PETSC_ER… in TSAdjointSetUp()
1092 if (ts->vecs_sensip) PetscCall(VecDuplicate(ts->vecs_sensip[0], &ts->vec_drdp_col)); in TSAdjointSetUp()
1119 if (ts->vecs_sensip) PetscCall(VecDestroy(&ts->vec_drdp_col)); in TSAdjointReset()
1122 ts->vecs_sensip = NULL; in TSAdjointReset()
1591 …jointMonitor(ts, ts->steps, ts->ptime, ts->vec_sol, ts->numcost, ts->vecs_sensi, ts->vecs_sensip)); in TSAdjointSolve()
1598 …jointMonitor(ts, ts->steps, ts->ptime, ts->vec_sol, ts->numcost, ts->vecs_sensi, ts->vecs_sensip)); in TSAdjointSolve()
|
| /petsc/src/ts/impls/implicit/theta/ |
| H A D | theta.c | 378 if (ts->vecs_sensip) { in TSAdjointStepBEuler_Private()
391 PetscCall(VecAXPY(ts->vecs_sensip[nadj], -adjoint_time_step, VecsDeltaMu[nadj])); in TSAdjointStepBEuler_Private()
395 PetscCall(VecAXPY(ts->vecs_sensip[nadj], adjoint_time_step, ts->vec_drdp_col)); in TSAdjointStepBEuler_Private()
565 if (ts->vecs_sensip) { /* sensitivities wrt parameters */ in TSAdjointStep_Theta()
573 … PetscCall(VecAXPY(ts->vecs_sensip[nadj], -adjoint_time_step * th->Theta, VecsDeltaMu[nadj])); in TSAdjointStep_Theta()
577 … PetscCall(VecAXPY(ts->vecs_sensip[nadj], adjoint_time_step * th->Theta, ts->vec_drdp_col)); in TSAdjointStep_Theta()
608 …PetscCall(VecAXPY(ts->vecs_sensip[nadj], -adjoint_time_step * (1.0 - th->Theta), VecsDeltaMu[nadj]… in TSAdjointStep_Theta()
612 …PetscCall(VecAXPY(ts->vecs_sensip[nadj], adjoint_time_step * (1.0 - th->Theta), ts->vec_drdp_col)); in TSAdjointStep_Theta()
652 if (ts->vecs_sensip) { in TSAdjointStep_Theta()
659 PetscCall(VecAXPY(ts->vecs_sensip[nadj], -adjoint_time_step, VecsDeltaMu[nadj])); in TSAdjointStep_Theta()
[all …]
|
| /petsc/src/ts/impls/explicit/rk/ |
| H A D | rk.c | 858 if (ts->vecs_sensip) PetscCall(VecDuplicate(ts->vecs_sensip[0], &rk->VecDeltaMu)); in TSAdjointSetUp_RK()
900 if (ts->vecs_sensip) { in TSAdjointStep_RK()
935 if (ts->vecs_sensip) { in TSAdjointStep_RK()
949 … PetscCall(VecAXPY(ts->vecs_sensip[nadj], 1., VecDeltaMu)); /* update sensip for each stage */ in TSAdjointStep_RK()
963 if (ts->vecs_sensip) { in TSAdjointStep_RK()
995 …if (ts->vecs_sensip) PetscCall(VecAXPY(VecsDeltaLam2[nadj * s + i], -h * b[i], ts->vecs_gup[nadj])… in TSAdjointStep_RK()
1003 … if (ts->vecs_sensip) PetscCall(VecAXPY(VecsDeltaLam2[nadj * s + i], -h, ts->vecs_gup[nadj])); in TSAdjointStep_RK()
|
| /petsc/src/ts/impls/arkimex/ |
| H A D | arkimex.c | 1590 if (ts->vecs_sensip) { in TSAdjointStep_ARKIMEX()
1610 if (ts->vecs_sensip) { in TSAdjointStep_ARKIMEX()
1614 PetscCall(VecAXPY(ts->vecs_sensip[nadj], adjoint_time_step, VecsSensiPTemp[nadj])); in TSAdjointStep_ARKIMEX()
1631 if (ts->vecs_sensip) { in TSAdjointStep_ARKIMEX()
1635 PetscCall(VecAXPY(ts->vecs_sensip[nadj], adjoint_time_step, VecsSensiPTemp[nadj])); in TSAdjointStep_ARKIMEX()
1653 if (ts->vecs_sensip) { in TSAdjointStep_ARKIMEX()
1657 …PetscCall(VecAXPY(ts->vecs_sensip[nadj], -At[i * s + i] * adjoint_time_step, VecsSensiPTemp[nadj])… in TSAdjointStep_ARKIMEX()
2056 …if (ts->vecs_sensip) PetscCall(VecDuplicateVecs(ts->vecs_sensip[0], ts->numcost, &ark->VecsSensiPT… in TSAdjointSetUp_ARKIMEX()
|
| /petsc/include/petsc/private/ |
| H A D | tsimpl.h | 200 Vec *vecs_sensip; member
|