Home
last modified time | relevance | path

Searched refs:VecsSensi2Temp (Results 1 – 3 of 3) sorted by relevance

/petsc/src/ts/impls/implicit/theta/
H A Dtheta.c38 …Vec *VecsSensi2Temp; /* Working vectors that holds the residual for the second-order ad… member
276 …taLam2 = th->VecsDeltaLam2, *VecsDeltaMu2 = th->VecsDeltaMu2, *VecsSensi2Temp = th->VecsSensi2Temp; in TSAdjointStepBEuler_Private() local
343 PetscCall(VecCopy(ts->vecs_sensi2[nadj], VecsSensi2Temp[nadj])); in TSAdjointStepBEuler_Private()
344 PetscCall(VecScale(VecsSensi2Temp[nadj], 1. / adjoint_time_step)); in TSAdjointStepBEuler_Private()
345 PetscCall(VecAXPY(VecsSensi2Temp[nadj], -1., ts->vecs_fuu[nadj])); in TSAdjointStepBEuler_Private()
346 if (ts->vecs_fup) PetscCall(VecAXPY(VecsSensi2Temp[nadj], -1., ts->vecs_fup[nadj])); in TSAdjointStepBEuler_Private()
351 PetscCall(KSPSolveTranspose(ksp, VecsSensi2Temp[nadj], VecsDeltaLam2[nadj])); in TSAdjointStepBEuler_Private()
372 …PetscCall(MatMultTransposeAdd(J, VecsDeltaLam2[nadj], VecsSensi2Temp[nadj], VecsSensi2Temp[nadj])); in TSAdjointStepBEuler_Private()
373 PetscCall(VecScale(VecsSensi2Temp[nadj], -adjoint_time_step)); in TSAdjointStepBEuler_Private()
374 PetscCall(VecCopy(VecsSensi2Temp[nadj], ts->vecs_sensi2[nadj])); in TSAdjointStepBEuler_Private()
[all …]
/petsc/src/ts/impls/explicit/rk/
H A Drk.h33 Vec *VecsSensi2Temp; member
H A Drk.c861 PetscCall(VecDuplicateVecs(ts->vecs_sensi2[0], ts->numcost, &rk->VecsSensi2Temp)); in TSAdjointSetUp_RK()
881 …DeltaLam2 = rk->VecsDeltaLam2, VecDeltaMu2 = rk->VecDeltaMu2, *VecsSensi2Temp = rk->VecsSensi2Temp; in TSAdjointStep_RK() local
990 PetscCall(VecCopy(ts->vecs_sensi2[nadj], VecsSensi2Temp[nadj])); in TSAdjointStep_RK()
991 PetscCall(VecMAXPY(VecsSensi2Temp[nadj], s - i - 1, w, &VecsDeltaLam2[nadj * s + i + 1])); in TSAdjointStep_RK()
992 PetscCall(MatMultTranspose(J, VecsSensi2Temp[nadj], VecsDeltaLam2[nadj * s + i])); in TSAdjointStep_RK()
999 PetscCall(VecMAXPY(VecsSensi2Temp[nadj], s - i - 1, w, &VecsDeltaLam2[nadj * s + i + 1])); in TSAdjointStep_RK()
1000 PetscCall(MatMultTranspose(J, VecsSensi2Temp[nadj], VecsDeltaLam2[nadj * s + i])); in TSAdjointStep_RK()
1006 PetscCall(MatMultTranspose(ts->Jacprhs, VecsSensi2Temp[nadj], VecDeltaMu2)); in TSAdjointStep_RK()
1042 PetscCall(VecDestroyVecs(ts->numcost, &rk->VecsSensi2Temp)); in TSAdjointReset_RK()