Lines Matching refs:tab
363 RKTableau t = &link->tab; in TSRKRegisterDestroy()
449 t = &link->tab; in TSRKRegister()
488 RKTableau tab = rk->tableau; in TSRKGetTableau_RK() local
491 if (s) *s = tab->s; in TSRKGetTableau_RK()
492 if (A) *A = tab->A; in TSRKGetTableau_RK()
493 if (b) *b = tab->b; in TSRKGetTableau_RK()
494 if (c) *c = tab->c; in TSRKGetTableau_RK()
495 if (bembed) *bembed = tab->bembed; in TSRKGetTableau_RK()
496 if (p) *p = tab->p; in TSRKGetTableau_RK()
497 if (binterp) *binterp = tab->binterp; in TSRKGetTableau_RK()
498 if (FSAL) *FSAL = tab->FSAL; in TSRKGetTableau_RK()
551 RKTableau tab = rk->tableau; in TSEvaluateStep_RK() local
554 PetscInt s = tab->s, j; in TSEvaluateStep_RK()
568 if (order == tab->order) { in TSEvaluateStep_RK()
571 for (j = 0; j < s; j++) w[j] = h * tab->b[j] / rk->dtratio; in TSEvaluateStep_RK()
575 } else if (order == tab->order - 1) { in TSEvaluateStep_RK()
576 if (!tab->bembed) goto unavailable; in TSEvaluateStep_RK()
579 for (j = 0; j < s; j++) w[j] = h * tab->bembed[j]; in TSEvaluateStep_RK()
583 for (j = 0; j < s; j++) w[j] = h * (tab->bembed[j] - tab->b[j]); in TSEvaluateStep_RK()
590 …onsider using -ts_adapt_type none or a different method that has an embedded estimate.", tab->name, in TSEvaluateStep_RK()
591 tab->order, order); in TSEvaluateStep_RK()
600 RKTableau tab = rk->tableau; in TSForwardCostIntegral_RK() local
601 const PetscInt s = tab->s; in TSForwardCostIntegral_RK()
602 const PetscReal *b = tab->b, *c = tab->c; in TSForwardCostIntegral_RK()
619 RKTableau tab = rk->tableau; in TSAdjointCostIntegral_RK() local
621 const PetscInt s = tab->s; in TSAdjointCostIntegral_RK()
622 const PetscReal *b = tab->b, *c = tab->c; in TSAdjointCostIntegral_RK()
639 RKTableau tab = rk->tableau; in TSRollBack_RK() local
640 const PetscInt s = tab->s; in TSRollBack_RK()
641 const PetscReal *b = tab->b, *c = tab->c; in TSRollBack_RK()
674 RKTableau tab = rk->tableau; in TSForwardStep_RK() local
676 const PetscInt s = tab->s; in TSForwardStep_RK()
677 const PetscReal *A = tab->A, *c = tab->c, *b = tab->b; in TSForwardStep_RK()
727 RKTableau tab = rk->tableau; in TSForwardGetStages_RK() local
730 if (ns) *ns = tab->s; in TSForwardGetStages_RK()
738 RKTableau tab = rk->tableau; in TSForwardSetUp_RK() local
745 PetscCall(PetscMalloc1(tab->s, &rk->MatsFwdStageSensip)); in TSForwardSetUp_RK()
746 PetscCall(PetscMalloc1(tab->s, &rk->MatsFwdSensipTemp)); in TSForwardSetUp_RK()
747 for (i = 0; i < tab->s; i++) { in TSForwardSetUp_RK()
758 RKTableau tab = rk->tableau; in TSForwardReset_RK() local
764 for (i = 0; i < tab->s; i++) PetscCall(MatDestroy(&rk->MatsFwdStageSensip[i])); in TSForwardReset_RK()
768 for (i = 0; i < tab->s; i++) PetscCall(MatDestroy(&rk->MatsFwdSensipTemp[i])); in TSForwardReset_RK()
778 RKTableau tab = rk->tableau; in TSStep_RK() local
779 const PetscInt s = tab->s; in TSStep_RK()
780 const PetscReal *A = tab->A, *c = tab->c; in TSStep_RK()
783 PetscBool FSAL = (PetscBool)(tab->FSAL && !rk->newtableau); in TSStep_RK()
814 PetscCall(TSEvaluateStep(ts, tab->order, ts->vec_sol, NULL)); in TSStep_RK()
818 …PetscCall(TSAdaptCandidateAdd(adapt, tab->name, tab->order, 1, tab->ccfl, (PetscReal)tab->s, PETSC… in TSStep_RK()
850 RKTableau tab = rk->tableau; in TSAdjointSetUp_RK() local
851 PetscInt s = tab->s; in TSAdjointSetUp_RK()
875 RKTableau tab = rk->tableau; in TSAdjointStep_RK() local
877 const PetscInt s = tab->s; in TSAdjointStep_RK()
878 const PetscReal *A = tab->A, *b = tab->b, *c = tab->c; in TSAdjointStep_RK()
893 if (tab->FSAL && i == s - 1) { in TSAdjointStep_RK()
1034 RKTableau tab = rk->tableau; in TSAdjointReset_RK() local
1037 PetscCall(VecDestroyVecs(tab->s * ts->numcost, &rk->VecsDeltaLam)); in TSAdjointReset_RK()
1040 PetscCall(VecDestroyVecs(tab->s * ts->numcost, &rk->VecsDeltaLam2)); in TSAdjointReset_RK()
1085 RKTableau tab = rk->tableau; in TSRKTableauReset() local
1088 if (!tab) PetscFunctionReturn(PETSC_SUCCESS); in TSRKTableauReset()
1090 PetscCall(VecDestroyVecs(tab->s, &rk->Y)); in TSRKTableauReset()
1091 PetscCall(VecDestroyVecs(tab->s, &rk->YdotRHS)); in TSRKTableauReset()
1134 RKTableau tab = rk->tableau; in TSRKTableauSetUp() local
1137 PetscCall(PetscMalloc1(tab->s, &rk->work)); in TSRKTableauSetUp()
1138 PetscCall(VecDuplicateVecs(ts->vec_sol, tab->s, &rk->Y)); in TSRKTableauSetUp()
1139 PetscCall(VecDuplicateVecs(ts->vec_sol, tab->s, &rk->YdotRHS)); in TSRKTableauSetUp()
1181 …ink = RKTableauList, count = 0; link; link = link->next, count++) namelist[count] = link->tab.name; in TSSetFromOptions_RK()
1203 RKTableau tab = rk->tableau; in TSView_RK() local
1213 PetscCall(PetscViewerASCIIPrintf(viewer, " Order: %" PetscInt_FMT "\n", tab->order)); in TSView_RK()
1333 PetscCall(PetscStrcmp(link->tab.name, rktype, &match)); in TSRKSetType_RK()
1336 rk->tableau = &link->tab; in TSRKSetType_RK()