Home
last modified time | relevance | path

Searched refs:VecDeltaMu2 (Results 1 – 2 of 2) sorted by relevance

/petsc/src/ts/impls/explicit/rk/
H A Drk.h32 Vec VecDeltaMu2; /* Increment of the 2nd-order adjoint sensitivity w.r.t P at stage */ member
H A Drk.c863 if (ts->vecs_sensi2p) PetscCall(VecDuplicate(ts->vecs_sensi2p[0], &rk->VecDeltaMu2)); in TSAdjointSetUp_RK()
881 …Vec *VecsDeltaLam2 = rk->VecsDeltaLam2, VecDeltaMu2 = rk->VecDeltaMu2, *VecsSensi2Temp… in TSAdjointStep_RK() local
1006 PetscCall(MatMultTranspose(ts->Jacprhs, VecsSensi2Temp[nadj], VecDeltaMu2)); in TSAdjointStep_RK()
1008 PetscCall(VecScale(VecDeltaMu2, -h * b[i])); in TSAdjointStep_RK()
1009 PetscCall(VecAXPY(VecDeltaMu2, -h * b[i], ts->vecs_gpu[nadj])); in TSAdjointStep_RK()
1010 PetscCall(VecAXPY(VecDeltaMu2, -h * b[i], ts->vecs_gpp[nadj])); in TSAdjointStep_RK()
1012 PetscCall(VecScale(VecDeltaMu2, -h)); in TSAdjointStep_RK()
1013 PetscCall(VecAXPY(VecDeltaMu2, -h, ts->vecs_gpu[nadj])); in TSAdjointStep_RK()
1014 PetscCall(VecAXPY(VecDeltaMu2, -h, ts->vecs_gpp[nadj])); in TSAdjointStep_RK()
1016 … PetscCall(VecAXPY(ts->vecs_sensi2p[nadj], 1, VecDeltaMu2)); /* update sensi2p for each stage */ in TSAdjointStep_RK()
[all …]