Searched refs:adjoint_time_step (Results 1 – 2 of 2) sorted by relevance
| /petsc/src/ts/impls/implicit/theta/ |
| H A D | theta.c | 283 PetscReal adjoint_time_step; in TSAdjointStepBEuler_Private() local 301 adjoint_time_step = -ts->time_step; /* always positive since time_step is negative */ in TSAdjointStepBEuler_Private() 308 PetscCall(VecScale(VecsSensiTemp[nadj], 1. / adjoint_time_step)); /* lambda_{n+1}/h */ in TSAdjointStepBEuler_Private() 319 th->shift = 1. / adjoint_time_step; in TSAdjointStepBEuler_Private() 344 PetscCall(VecScale(VecsSensi2Temp[nadj], 1. / adjoint_time_step)); in TSAdjointStepBEuler_Private() 369 PetscCall(VecScale(VecsSensiTemp[nadj], -adjoint_time_step)); in TSAdjointStepBEuler_Private() 373 PetscCall(VecScale(VecsSensi2Temp[nadj], -adjoint_time_step)); in TSAdjointStepBEuler_Private() 379 …PetscCall(VecAXPBYPCZ(th->Xdot, -1. / adjoint_time_step, 1.0 / adjoint_time_step, 0, th->X0, ts->v… in TSAdjointStepBEuler_Private() 380 …PetscCall(TSComputeIJacobianP(ts, th->stage_time, ts->vec_sol, th->Xdot, 1. / adjoint_time_step, t… in TSAdjointStepBEuler_Private() 391 PetscCall(VecAXPY(ts->vecs_sensip[nadj], -adjoint_time_step, VecsDeltaMu[nadj])); in TSAdjointStepBEuler_Private() [all …]
|
| /petsc/src/ts/impls/arkimex/ |
| H A D | arkimex.c | 1571 …PetscReal adjoint_time_step = -ts->time_step; /* always positive since ts->time_step is neg… in TSAdjointStep_ARKIMEX() local 1581 ark->stage_time = t - adjoint_time_step * (1.0 - ct[i]); in TSAdjointStep_ARKIMEX() 1582 stage_time_ex = t - adjoint_time_step * (1.0 - c[i]); in TSAdjointStep_ARKIMEX() 1591 …(TSComputeIJacobianP(ts, ark->stage_time, Y[i], Ydot, ark->scoeff / adjoint_time_step, ts->Jacp, P… in TSAdjointStep_ARKIMEX() 1609 …PetscCall(VecAXPY(VecsDeltaLam[nadj * s + i], ark->scoeff / adjoint_time_step, VecsSensiTemp[nadj]… in TSAdjointStep_ARKIMEX() 1614 PetscCall(VecAXPY(ts->vecs_sensip[nadj], adjoint_time_step, VecsSensiPTemp[nadj])); in TSAdjointStep_ARKIMEX() 1635 PetscCall(VecAXPY(ts->vecs_sensip[nadj], adjoint_time_step, VecsSensiPTemp[nadj])); in TSAdjointStep_ARKIMEX() 1640 PetscCall(VecScale(VecsDeltaLam[nadj * s + i], adjoint_time_step)); in TSAdjointStep_ARKIMEX() 1657 …PetscCall(VecAXPY(ts->vecs_sensip[nadj], -At[i * s + i] * adjoint_time_step, VecsSensiPTemp[nadj])… in TSAdjointStep_ARKIMEX()
|