Searched refs:VecDeltaMu (Results 1 – 2 of 2) sorted by relevance
| /petsc/src/ts/impls/explicit/rk/ |
| H A D | rk.h | 30 …Vec VecDeltaMu; /* Increment of the adjoint sensitivity w.r.t P at stage … member
|
| H A D | rk.c | 858 if (ts->vecs_sensip) PetscCall(VecDuplicate(ts->vecs_sensip[0], &rk->VecDeltaMu)); in TSAdjointSetUp_RK()
880 …Vec *Y = rk->Y, *VecsDeltaLam = rk->VecsDeltaLam, VecDeltaMu = rk->VecDeltaMu, *VecsSe… in TSAdjointStep_RK() local
936 PetscCall(MatMultTranspose(ts->Jacprhs, VecsSensiTemp[nadj], VecDeltaMu)); in TSAdjointStep_RK()
938 PetscCall(VecScale(VecDeltaMu, -h * b[i])); in TSAdjointStep_RK()
942 PetscCall(VecAXPY(VecDeltaMu, -h * b[i], VecDRDPTransCol)); in TSAdjointStep_RK()
947 PetscCall(VecScale(VecDeltaMu, -h)); in TSAdjointStep_RK()
949 … PetscCall(VecAXPY(ts->vecs_sensip[nadj], 1., VecDeltaMu)); /* update sensip for each stage */ in TSAdjointStep_RK()
1039 PetscCall(VecDestroy(&rk->VecDeltaMu)); in TSAdjointReset_RK()
|