Home
last modified time | relevance | path

Searched refs:VecsDeltaLam (Results 1 – 5 of 5) sorted by relevance

/petsc/src/ts/impls/implicit/theta/
H A Dtheta.c27 Vec *VecsDeltaLam; /* Increment of the adjoint sensitivity w.r.t IC at stage */ member
275 …Vec *VecsDeltaLam = th->VecsDeltaLam, *VecsDeltaMu = th->VecsDeltaMu, *VecsSensiTemp = t… in TSAdjointStepBEuler_Private() local
289 VecsDeltaLam = ts->vecs_sensi; in TSAdjointStepBEuler_Private()
326 PetscCall(KSPSolveTranspose(ksp, VecsSensiTemp[nadj], VecsDeltaLam[nadj])); in TSAdjointStepBEuler_Private()
339 …PetscCall(TSComputeIHessianProductFunctionUU(ts, th->stage_time, ts->vec_sol, VecsDeltaLam, ts->ve… in TSAdjointStepBEuler_Private()
341 …PetscCall(TSComputeIHessianProductFunctionUP(ts, th->stage_time, ts->vec_sol, VecsDeltaLam, ts->ve… in TSAdjointStepBEuler_Private()
368 … PetscCall(MatMultTransposeAdd(J, VecsDeltaLam[nadj], VecsSensiTemp[nadj], VecsSensiTemp[nadj])); in TSAdjointStepBEuler_Private()
384 …PetscCall(TSComputeIHessianProductFunctionPU(ts, th->stage_time, ts->vec_sol, VecsDeltaLam, ts->ve… in TSAdjointStepBEuler_Private()
386 …PetscCall(TSComputeIHessianProductFunctionPP(ts, th->stage_time, ts->vec_sol, VecsDeltaLam, ts->ve… in TSAdjointStepBEuler_Private()
390 PetscCall(MatMultTranspose(ts->Jacp, VecsDeltaLam[nadj], VecsDeltaMu[nadj])); in TSAdjointStepBEuler_Private()
[all …]
/petsc/src/ts/impls/explicit/rk/
H A Drk.h28 …Vec *VecsDeltaLam; /* Increment of the adjoint sensitivity w.r.t IC at stage … member
H A Drk.c856 PetscCall(VecDuplicateVecs(ts->vecs_sensi[0], s * ts->numcost, &rk->VecsDeltaLam)); in TSAdjointSetUp_RK()
880 …Vec *Y = rk->Y, *VecsDeltaLam = rk->VecsDeltaLam, VecDeltaMu = rk->VecDeltaMu, *VecsSe… in TSAdjointStep_RK() local
916 PetscCall(VecMAXPY(VecsSensiTemp[nadj], s - i - 1, w, &VecsDeltaLam[nadj * s + i + 1])); in TSAdjointStep_RK()
917 …PetscCall(MatMultTranspose(J, VecsSensiTemp[nadj], VecsDeltaLam[nadj * s + i])); /* VecsSensiTemp … in TSAdjointStep_RK()
918 PetscCall(VecScale(VecsDeltaLam[nadj * s + i], -h * b[i])); in TSAdjointStep_RK()
922 PetscCall(VecAXPY(VecsDeltaLam[nadj * s + i], -h * b[i], VecDRDUTransCol)); in TSAdjointStep_RK()
929 PetscCall(VecMAXPY(VecsSensiTemp[nadj], s - i - 1, w, &VecsDeltaLam[nadj * s + i + 1])); in TSAdjointStep_RK()
930 PetscCall(MatMultTranspose(J, VecsSensiTemp[nadj], VecsDeltaLam[nadj * s + i])); in TSAdjointStep_RK()
931 PetscCall(VecScale(VecsDeltaLam[nadj * s + i], -h)); in TSAdjointStep_RK()
1024 PetscCall(VecMAXPY(ts->vecs_sensi[nadj], s, w, &VecsDeltaLam[nadj * s])); in TSAdjointStep_RK()
[all …]
/petsc/src/ts/impls/arkimex/
H A Darkimex.h52 Vec *VecsDeltaLam; /* Increment of the adjoint sensitivity w.r.t IC at stage */ member
H A Darkimex.c1567 …Vec *Y = ark->Y, Ydot = ark->Ydot, *VecsDeltaLam = ark->VecsDeltaLam, *VecsSensiTemp =… in TSAdjointStep_ARKIMEX() local
1601 PetscCall(VecMAXPY(VecsSensiTemp[nadj], s - i - 1, w, &VecsDeltaLam[nadj * s + i + 1])); in TSAdjointStep_ARKIMEX()
1607 PetscCall(MatMultTranspose(Jim, VecsSensiTemp[nadj], VecsDeltaLam[nadj * s + i])); in TSAdjointStep_ARKIMEX()
1609 …PetscCall(VecAXPY(VecsDeltaLam[nadj * s + i], ark->scoeff / adjoint_time_step, VecsSensiTemp[nadj]… in TSAdjointStep_ARKIMEX()
1617 PetscCall(VecSet(VecsDeltaLam[nadj * s + i], 0)); // make sure it is initialized in TSAdjointStep_ARKIMEX()
1624 PetscCall(VecMAXPY(VecsSensiTemp[nadj], s - i - 1, w, &VecsDeltaLam[nadj * s + i + 1])); in TSAdjointStep_ARKIMEX()
1630 …PetscCall(MatMultTransposeAdd(Jex, VecsSensiTemp[nadj], VecsDeltaLam[nadj * s + i], VecsDeltaLam[n… in TSAdjointStep_ARKIMEX()
1640 PetscCall(VecScale(VecsDeltaLam[nadj * s + i], adjoint_time_step)); in TSAdjointStep_ARKIMEX()
1646 PetscCall(VecScale(VecsDeltaLam[nadj * s + i], 1. / At[i * s + i])); in TSAdjointStep_ARKIMEX()
1647 PetscCall(KSPSolveTranspose(ksp, VecsDeltaLam[nadj * s + i], VecsDeltaLam[nadj * s + i])); in TSAdjointStep_ARKIMEX()
[all …]