Home
last modified time | relevance | path

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

/petsc/src/ts/impls/explicit/rk/
H A Drk.c860 PetscCall(VecDuplicateVecs(ts->vecs_sensi[0], s * ts->numcost, &rk->VecsDeltaLam2)); in TSAdjointSetUp_RK()
881 …Vec *VecsDeltaLam2 = rk->VecsDeltaLam2, VecDeltaMu2 = rk->VecDeltaMu2, *VecsSensi2Temp… in TSAdjointStep_RK() local
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()
993 PetscCall(VecScale(VecsDeltaLam2[nadj * s + i], -h * b[i])); in TSAdjointStep_RK()
994 PetscCall(VecAXPY(VecsDeltaLam2[nadj * s + i], -h * b[i], ts->vecs_guu[nadj])); in TSAdjointStep_RK()
995 …if (ts->vecs_sensip) PetscCall(VecAXPY(VecsDeltaLam2[nadj * s + i], -h * b[i], ts->vecs_gup[nadj])… in TSAdjointStep_RK()
998 PetscCall(VecSet(VecsDeltaLam2[nadj * s + i], 0)); 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()
[all …]
H A Drk.h31 Vec *VecsDeltaLam2; /* Increment of the 2nd-order adjoint sensitivity w.r.t IC at stage */ member
/petsc/src/ts/impls/implicit/theta/
H A Dtheta.c36 …Vec *VecsDeltaLam2; /* Increment of the 2nd-order adjoint sensitivity w.r.t IC at stag… member
276 …Vec *VecsDeltaLam2 = th->VecsDeltaLam2, *VecsDeltaMu2 = th->VecsDeltaMu2, *VecsSensi2Tem… in TSAdjointStepBEuler_Private() local
290 VecsDeltaLam2 = ts->vecs_sensi2; 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()
400 PetscCall(MatMultTranspose(ts->Jacp, VecsDeltaLam2[nadj], VecsDeltaMu2[nadj])); in TSAdjointStepBEuler_Private()
421 …Vec *VecsDeltaLam2 = th->VecsDeltaLam2, *VecsDeltaMu2 = th->VecsDeltaMu2, *VecsSensi2Temp … in TSAdjointStep_Theta() local
514 PetscCall(KSPSolveTranspose(ksp, VecsSensi2Temp[nadj], VecsDeltaLam2[nadj])); in TSAdjointStep_Theta()
556 PetscCall(MatMultTranspose(J, VecsDeltaLam2[nadj], ts->vecs_sensi2[nadj])); in TSAdjointStep_Theta()
595 PetscCall(MatMultTranspose(ts->Jacp, VecsDeltaLam2[nadj], VecsDeltaMu2[nadj])); in TSAdjointStep_Theta()
[all …]