Searched refs:vecs_sensi (Results 1 – 5 of 5) sorted by relevance
| /petsc/src/ts/interface/sensitivity/ |
| H A D | tssen.c | 882 ts->vecs_sensi = lambda; in TSSetCostGradients()
911 if (lambda) *lambda = ts->vecs_sensi; in TSGetCostGradients()
1079 …PetscCheck(ts->vecs_sensi, PetscObjectComm((PetscObject)ts), PETSC_ERR_ARG_WRONGSTATE, "Must call … in TSAdjointSetUp()
1091 PetscCall(VecDuplicate(ts->vecs_sensi[0], &ts->vec_drdu_col)); in TSAdjointSetUp()
1121 ts->vecs_sensi = NULL; in TSAdjointReset()
1591 …PetscCall(TSAdjointMonitor(ts, ts->steps, ts->ptime, ts->vec_sol, ts->numcost, ts->vecs_sensi, ts-… in TSAdjointSolve()
1598 …PetscCall(TSAdjointMonitor(ts, ts->steps, ts->ptime, ts->vec_sol, ts->numcost, ts->vecs_sensi, ts-… in TSAdjointSolve()
1602 PetscCall(VecViewFromOptions(ts->vecs_sensi[0], (PetscObject)ts, "-ts_adjoint_view_solution")); in TSAdjointSolve()
|
| /petsc/src/ts/impls/implicit/theta/ |
| H A D | theta.c | 289 VecsDeltaLam = ts->vecs_sensi; in TSAdjointStepBEuler_Private()
307 PetscCall(VecCopy(ts->vecs_sensi[nadj], VecsSensiTemp[nadj])); in TSAdjointStepBEuler_Private()
370 PetscCall(VecCopy(VecsSensiTemp[nadj], ts->vecs_sensi[nadj])); in TSAdjointStepBEuler_Private()
462 PetscCall(VecCopy(ts->vecs_sensi[nadj], VecsSensiTemp[nadj])); in TSAdjointStep_Theta()
532 PetscCall(MatMultTranspose(J, VecsDeltaLam[nadj], ts->vecs_sensi[nadj])); in TSAdjointStep_Theta()
536 PetscCall(VecAXPY(ts->vecs_sensi[nadj], -1., ts->vec_drdu_col)); in TSAdjointStep_Theta()
540 PetscCall(VecScale(ts->vecs_sensi[nadj], 1. / th->shift)); in TSAdjointStep_Theta()
643 PetscCall(VecAXPY(ts->vecs_sensi[nadj], -adjoint_time_step, VecsSensiTemp[nadj])); in TSAdjointStep_Theta()
647 PetscCall(VecAXPY(ts->vecs_sensi[nadj], adjoint_time_step, ts->vec_drdu_col)); in TSAdjointStep_Theta()
1064 PetscCall(VecDuplicateVecs(ts->vecs_sensi[0], ts->numcost, &th->VecsDeltaLam)); in TSAdjointSetUp_Theta()
[all …]
|
| /petsc/src/ts/impls/explicit/rk/ |
| H A D | rk.c | 856 PetscCall(VecDuplicateVecs(ts->vecs_sensi[0], s * ts->numcost, &rk->VecsDeltaLam)); in TSAdjointSetUp_RK()
857 PetscCall(VecDuplicateVecs(ts->vecs_sensi[0], ts->numcost, &rk->VecsSensiTemp)); in TSAdjointSetUp_RK()
860 PetscCall(VecDuplicateVecs(ts->vecs_sensi[0], s * ts->numcost, &rk->VecsDeltaLam2)); in TSAdjointSetUp_RK()
915 …PetscCall(VecCopy(ts->vecs_sensi[nadj], VecsSensiTemp[nadj])); /* VecDeltaLam is an vec array of s… in TSAdjointStep_RK()
1024 PetscCall(VecMAXPY(ts->vecs_sensi[nadj], s, w, &VecsDeltaLam[nadj * s])); in TSAdjointStep_RK()
|
| /petsc/include/petsc/private/ |
| H A D | tsimpl.h | 199 Vec *vecs_sensi; /* one vector for each cost function */ member
|
| /petsc/src/ts/impls/arkimex/ |
| H A D | arkimex.c | 1604 PetscCall(VecAXPY(VecsSensiTemp[nadj], -bt[i], ts->vecs_sensi[nadj])); in TSAdjointStep_ARKIMEX()
1627 PetscCall(VecAXPY(VecsSensiTemp[nadj], b[i], ts->vecs_sensi[nadj])); in TSAdjointStep_ARKIMEX()
1664 PetscCall(VecMAXPY(ts->vecs_sensi[nadj], s, w, &VecsDeltaLam[nadj * s])); in TSAdjointStep_ARKIMEX()
2054 PetscCall(VecDuplicateVecs(ts->vecs_sensi[0], tab->s * ts->numcost, &ark->VecsDeltaLam)); in TSAdjointSetUp_ARKIMEX()
2055 PetscCall(VecDuplicateVecs(ts->vecs_sensi[0], ts->numcost, &ark->VecsSensiTemp)); in TSAdjointSetUp_ARKIMEX()
|