Home
last modified time | relevance | path

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

/petsc/src/ts/impls/implicit/theta/
H A Dtheta.c29 Vec *VecsSensiTemp; /* Vector to be multiplied with Jacobian transpose */ member
275 …ecsDeltaLam = th->VecsDeltaLam, *VecsDeltaMu = th->VecsDeltaMu, *VecsSensiTemp = th->VecsSensiTemp; in TSAdjointStepBEuler_Private() local
307 PetscCall(VecCopy(ts->vecs_sensi[nadj], VecsSensiTemp[nadj])); in TSAdjointStepBEuler_Private()
308 PetscCall(VecScale(VecsSensiTemp[nadj], 1. / adjoint_time_step)); /* lambda_{n+1}/h */ in TSAdjointStepBEuler_Private()
312 PetscCall(VecAXPY(VecsSensiTemp[nadj], 1., ts->vec_drdu_col)); in TSAdjointStepBEuler_Private()
326 PetscCall(KSPSolveTranspose(ksp, VecsSensiTemp[nadj], VecsDeltaLam[nadj])); in TSAdjointStepBEuler_Private()
367 PetscCall(VecScale(VecsSensiTemp[nadj], -1.)); in TSAdjointStepBEuler_Private()
368 … PetscCall(MatMultTransposeAdd(J, VecsDeltaLam[nadj], VecsSensiTemp[nadj], VecsSensiTemp[nadj])); in TSAdjointStepBEuler_Private()
369 PetscCall(VecScale(VecsSensiTemp[nadj], -adjoint_time_step)); in TSAdjointStepBEuler_Private()
370 PetscCall(VecCopy(VecsSensiTemp[nadj], ts->vecs_sensi[nadj])); in TSAdjointStepBEuler_Private()
[all …]
/petsc/src/ts/impls/explicit/rk/
H A Drk.h29 Vec *VecsSensiTemp; member
H A Drk.c857 PetscCall(VecDuplicateVecs(ts->vecs_sensi[0], ts->numcost, &rk->VecsSensiTemp)); in TSAdjointSetUp_RK()
880 … *VecsDeltaLam = rk->VecsDeltaLam, VecDeltaMu = rk->VecDeltaMu, *VecsSensiTemp = rk->VecsSensiTemp; in TSAdjointStep_RK() local
915 …PetscCall(VecCopy(ts->vecs_sensi[nadj], VecsSensiTemp[nadj])); /* VecDeltaLam is an vec array of s… in TSAdjointStep_RK()
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()
928 PetscCall(VecSet(VecsSensiTemp[nadj], 0)); 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()
936 PetscCall(MatMultTranspose(ts->Jacprhs, VecsSensiTemp[nadj], VecDeltaMu)); in TSAdjointStep_RK()
958 …PetscCall(TSComputeRHSHessianProductFunctionUU(ts, rk->stage_time, Y[i], VecsSensiTemp, ts->vec_se… in TSAdjointStep_RK()
[all …]
/petsc/src/ts/impls/arkimex/
H A Darkimex.h53 Vec *VecsSensiTemp; /* Vectors to be multiplied with Jacobian transpose */ member
H A Darkimex.c1567 …->Y, Ydot = ark->Ydot, *VecsDeltaLam = ark->VecsDeltaLam, *VecsSensiTemp = ark->VecsSensiTemp, *Ve… in TSAdjointStep_ARKIMEX() local
1597 PetscCall(VecSet(VecsSensiTemp[nadj], 0)); in TSAdjointStep_ARKIMEX()
1601 PetscCall(VecMAXPY(VecsSensiTemp[nadj], s - i - 1, w, &VecsDeltaLam[nadj * s + i + 1])); in TSAdjointStep_ARKIMEX()
1604 PetscCall(VecAXPY(VecsSensiTemp[nadj], -bt[i], ts->vecs_sensi[nadj])); in TSAdjointStep_ARKIMEX()
1607 PetscCall(MatMultTranspose(Jim, VecsSensiTemp[nadj], VecsDeltaLam[nadj * s + i])); in TSAdjointStep_ARKIMEX()
1609 …tscCall(VecAXPY(VecsDeltaLam[nadj * s + i], ark->scoeff / adjoint_time_step, VecsSensiTemp[nadj])); in TSAdjointStep_ARKIMEX()
1612 PetscCall(MatMultTranspose(ts->Jacp, VecsSensiTemp[nadj], VecsSensiPTemp[nadj])); in TSAdjointStep_ARKIMEX()
1620 PetscCall(VecSet(VecsSensiTemp[nadj], 0)); in TSAdjointStep_ARKIMEX()
1624 PetscCall(VecMAXPY(VecsSensiTemp[nadj], s - i - 1, w, &VecsDeltaLam[nadj * s + i + 1])); in TSAdjointStep_ARKIMEX()
1627 PetscCall(VecAXPY(VecsSensiTemp[nadj], b[i], ts->vecs_sensi[nadj])); in TSAdjointStep_ARKIMEX()
[all …]