Lines Matching refs:Ydot

43     else *Ydotstage = gl->Ydot[gl->stage];  in TSGLLEGetVecs()
69 Vec Ydot, Ydot_c; in DMRestrictHook_TSGLLE() local
72 PetscCall(TSGLLEGetVecs(ts, fine, NULL, &Ydot)); in DMRestrictHook_TSGLLE()
74 PetscCall(MatRestrict(restrct, Ydot, Ydot_c)); in DMRestrictHook_TSGLLE()
76 PetscCall(TSGLLERestoreVecs(ts, fine, NULL, &Ydot)); in DMRestrictHook_TSGLLE()
90 Vec Ydot, Ydot_s; in DMSubDomainRestrictHook_TSGLLE() local
93 PetscCall(TSGLLEGetVecs(ts, dm, NULL, &Ydot)); in DMSubDomainRestrictHook_TSGLLE()
96 PetscCall(VecScatterBegin(gscat, Ydot, Ydot_s, INSERT_VALUES, SCATTER_FORWARD)); in DMSubDomainRestrictHook_TSGLLE()
97 PetscCall(VecScatterEnd(gscat, Ydot, Ydot_s, INSERT_VALUES, SCATTER_FORWARD)); in DMSubDomainRestrictHook_TSGLLE()
99 PetscCall(TSGLLERestoreVecs(ts, dm, NULL, &Ydot)); in DMSubDomainRestrictHook_TSGLLE()
358 static PetscErrorCode TSGLLEEstimateHigherMoments_Default(TSGLLEScheme sc, PetscReal h, Vec Ydot[],… in TSGLLEEstimateHigherMoments_Default() argument
370 PetscCall(VecMAXPY(hm[i], sc->s, phih, Ydot)); in TSGLLEEstimateHigherMoments_Default()
376 …TSGLLEScheme sc, PetscReal h, TSGLLEScheme next_sc, PetscReal next_h, Vec Ydot[], Vec Xold[], Vec … in TSGLLECompleteStep_Rescale() argument
388 PetscCall(VecMAXPY(X[i], s, brow, Ydot)); in TSGLLECompleteStep_Rescale()
395 …TSGLLEScheme sc, PetscReal h, TSGLLEScheme next_sc, PetscReal next_h, Vec Ydot[], Vec Xold[], Vec … in TSGLLECompleteStep_RescaleAndModify() argument
412 PetscCall(VecMAXPY(X[i], s, brow, Ydot)); in TSGLLECompleteStep_RescaleAndModify()
422 PetscCall(VecMAXPY(X[r], s, brow, Ydot)); in TSGLLECompleteStep_RescaleAndModify()
875 Vec *X, *Ydot, Y; in TSSolve_GLLE() local
885 Ydot = gl->Ydot; in TSSolve_GLLE()
913 for (j = 0; j < i; j++) PetscCall(VecAXPY(gl->Z, -shift * h * a[i * s + j], Ydot[j])); in TSSolve_GLLE()
923 PetscCall(VecAXPY(Y, (c[i] - c[i - 1]) * h, Ydot[i - 1])); in TSSolve_GLLE()
945 PetscCall((*gl->EstimateHigherMoments)(scheme, h, Ydot, gl->X, gl->himom)); in TSSolve_GLLE()
987 …PetscCall((*gl->CompleteStep)(scheme, h, gl->schemes[next_scheme], next_h, Ydot, gl->Xold, gl->X)); in TSSolve_GLLE()
1016 PetscCall(VecDestroyVecs(max_s, &gl->Ydot)); in TSReset_GLLE()
1052 Vec Z, Ydot; in SNESTSFormFunction_GLLE() local
1057 PetscCall(TSGLLEGetVecs(ts, dm, &Z, &Ydot)); in SNESTSFormFunction_GLLE()
1058 PetscCall(VecWAXPY(Ydot, gl->scoeff / ts->time_step, x, Z)); in SNESTSFormFunction_GLLE()
1061 PetscCall(TSComputeIFunction(ts, gl->stage_time, x, Ydot, f, PETSC_FALSE)); in SNESTSFormFunction_GLLE()
1063 PetscCall(TSGLLERestoreVecs(ts, dm, &Z, &Ydot)); in SNESTSFormFunction_GLLE()
1070 Vec Z, Ydot; in SNESTSFormJacobian_GLLE() local
1075 PetscCall(TSGLLEGetVecs(ts, dm, &Z, &Ydot)); in SNESTSFormJacobian_GLLE()
1079 …PetscCall(TSComputeIJacobian(ts, gl->stage_time, x, gl->Ydot[gl->stage], gl->scoeff / ts->time_ste… in SNESTSFormJacobian_GLLE()
1081 PetscCall(TSGLLERestoreVecs(ts, dm, &Z, &Ydot)); in SNESTSFormJacobian_GLLE()
1096 PetscCall(VecDuplicateVecs(ts->vec_sol, max_s, &gl->Ydot)); in TSSetUp_GLLE()