Lines Matching refs:shift
16 PetscReal shift; /* Shift parameter for SNES Jacobian, used by forward, TLM and adjoint */ member
212 th->shift = 1 / (th->Theta * ts->time_step); in TSStep_Theta()
215 if (th->extrapolate && !ts->steprestart) PetscCall(VecAXPY(th->X, 1 / th->shift, th->Xdot)); in TSStep_Theta()
231 …PetscCall(VecAXPBYPCZ(th->Xdot, -th->shift, th->shift, 0, th->X0, th->X)); /* th->Xdot is needed b… in TSStep_Theta()
319 th->shift = 1. / adjoint_time_step; in TSAdjointStepBEuler_Private()
362 th->shift = 0.0; in TSAdjointStepBEuler_Private()
474 th->shift = 1. / (th->Theta * adjoint_time_step); in TSAdjointStep_Theta()
507 PetscCall(VecScale(VecsSensi2Temp[nadj], th->shift)); in TSAdjointStep_Theta()
525 th->shift = 1. / ((th->Theta - 1.) * adjoint_time_step); in TSAdjointStep_Theta()
540 PetscCall(VecScale(ts->vecs_sensi[nadj], 1. / th->shift)); in TSAdjointStep_Theta()
559 PetscCall(VecScale(ts->vecs_sensi2[nadj], 1. / th->shift)); in TSAdjointStep_Theta()
567 th->shift = 1.0 / (adjoint_time_step * th->Theta); in TSAdjointStep_Theta()
568 PetscCall(VecAXPBYPCZ(th->Xdot, -th->shift, th->shift, 0, th->X0, ts->vec_sol)); in TSAdjointStep_Theta()
637 th->shift = 0.0; in TSAdjointStep_Theta()
653 th->shift = 1.0 / (adjoint_time_step * th->Theta); in TSAdjointStep_Theta()
654 PetscCall(VecAXPBYPCZ(th->Xdot, -th->shift, th->shift, 0, th->X0, th->X)); in TSAdjointStep_Theta()
655 …PetscCall(TSComputeIJacobianP(ts, th->stage_time, th->X, th->Xdot, th->shift, ts->Jacp, PETSC_FALS… in TSAdjointStep_Theta()
751 previous_shift = th->shift; in TSForwardStep_Theta()
764 th->shift = 1. / ((th->Theta - 1.) * th->time_step0); in TSForwardStep_Theta()
765 … PetscCall(TSComputeIJacobian(ts, th->ptime0, th->X0, th->Xdot, th->shift, J, Jpre, PETSC_FALSE)); in TSForwardStep_Theta()
772 …PetscCall(TSComputeIJacobianP(ts, th->ptime0, th->X0, th->Xdot, th->shift, ts->Jacp, PETSC_FALSE)); in TSForwardStep_Theta()
774 th->shift = previous_shift; in TSForwardStep_Theta()
775 PetscCall(VecAXPBYPCZ(th->Xdot, -th->shift, th->shift, 0, th->X0, ts->vec_sol)); in TSForwardStep_Theta()
776 …PetscCall(TSComputeIJacobianP(ts, th->stage_time, ts->vec_sol, th->Xdot, th->shift, ts->Jacp, PETS… in TSForwardStep_Theta()
780 th->shift = 0.0; in TSForwardStep_Theta()
781 …PetscCall(TSComputeIJacobian(ts, th->stage_time, th->X, th->Xdot, th->shift, J, Jpre, PETSC_FALSE)… in TSForwardStep_Theta()
787 th->shift = previous_shift; in TSForwardStep_Theta()
788 PetscCall(VecAXPBYPCZ(th->Xdot, -th->shift, th->shift, 0, th->X0, th->X)); in TSForwardStep_Theta()
789 …PetscCall(TSComputeIJacobianP(ts, th->stage_time, th->X, th->Xdot, th->shift, ts->Jacp, PETSC_FALS… in TSForwardStep_Theta()
795 th->shift = previous_shift; /* recover the previous shift used in TSStep_Theta() */ in TSForwardStep_Theta()
797 …PetscCall(TSComputeIJacobian(ts, th->stage_time, ts->vec_sol, th->Xdot, th->shift, J, Jpre, PETSC_… in TSForwardStep_Theta()
799 …PetscCall(TSComputeIJacobian(ts, th->stage_time, th->X, th->Xdot, th->shift, J, Jpre, PETSC_FALSE)… in TSForwardStep_Theta()
948 PetscReal shift = th->shift; in SNESTSFormFunction_Theta() local
955 PetscCall(VecAXPBYPCZ(Xdot, -shift, shift, 0, X0, x)); in SNESTSFormFunction_Theta()
973 PetscReal shift = th->shift; in SNESTSFormJacobian_Theta() local
982 PetscCall(TSComputeIJacobian(ts, th->stage_time, x, Xdot, shift, A, B, PETSC_FALSE)); in SNESTSFormJacobian_Theta()
1040 th->shift = 1 / (th->Theta * ts->time_step); in TSSetUp_Theta()