Lines Matching refs:tab
1030 ARKTableau t = &link->tab; in TSARKIMEXRegisterDestroy()
1116 PetscCall(PetscStrcmp(link->tab.name, name, &match)); in TSARKIMEXRegister()
1120 t = &link->tab; in TSARKIMEXRegister()
1228 ARKTableau tab = ark->tableau; in TSEvaluateStep_ARKIMEX() local
1231 PetscInt s = tab->s, j; in TSEvaluateStep_ARKIMEX()
1246 if (order == tab->order) { in TSEvaluateStep_ARKIMEX()
1248 …if (!ark->imex && tab->stiffly_accurate) { /* Only the stiffly accurate implicit formula is used */ in TSEvaluateStep_ARKIMEX()
1252 for (j = 0; j < s; j++) w[j] = h * tab->bt[j]; in TSEvaluateStep_ARKIMEX()
1254 if (tab->additive && ark->imex) { /* Method is IMEX, complete the explicit formula */ in TSEvaluateStep_ARKIMEX()
1257 for (j = 0; j < s; j++) w[j] = h * tab->b[j]; in TSEvaluateStep_ARKIMEX()
1265 } else if (order == tab->order - 1) { in TSEvaluateStep_ARKIMEX()
1266 if (!tab->bembedt) goto unavailable; in TSEvaluateStep_ARKIMEX()
1269 for (j = 0; j < s; j++) w[j] = h * tab->bembedt[j]; in TSEvaluateStep_ARKIMEX()
1271 if (tab->additive) { in TSEvaluateStep_ARKIMEX()
1274 for (j = 0; j < s; j++) w[j] = h * tab->bembed[j]; in TSEvaluateStep_ARKIMEX()
1280 for (j = 0; j < s; j++) w[j] = h * (tab->bembedt[j] - tab->bt[j]); in TSEvaluateStep_ARKIMEX()
1281 PetscCall(VecMAXPY(X, tab->s, w, ark->YdotI)); in TSEvaluateStep_ARKIMEX()
1282 if (tab->additive) { in TSEvaluateStep_ARKIMEX()
1285 for (j = 0; j < s; j++) w[j] = h * (tab->bembed[j] - tab->b[j]); in TSEvaluateStep_ARKIMEX()
1295 tab->name, tab->order, order); in TSEvaluateStep_ARKIMEX()
1333 ARKTableau tab = ark->tableau; in TSStep_ARKIMEX() local
1334 const PetscInt s = tab->s; in TSStep_ARKIMEX()
1335 const PetscReal *At = tab->At, *A = tab->A, *ct = tab->ct, *c = tab->c; in TSStep_ARKIMEX()
1343 …PetscBool hasE = PETSC_FALSE, dirk = (PetscBool)(!tab->additive), stageok, accept = PETSC_T… in TSStep_ARKIMEX()
1348 PetscCall(VecDuplicateVecs(ts->vec_sol, tab->s, &ark->Y_prev)); in TSStep_ARKIMEX()
1349 PetscCall(VecDuplicateVecs(ts->vec_sol, tab->s, &ark->YdotI_prev)); in TSStep_ARKIMEX()
1350 if (tab->additive) PetscCall(VecDuplicateVecs(ts->vec_sol, tab->s, &ark->YdotRHS_prev)); in TSStep_ARKIMEX()
1364 if (tab->additive && hasE) PetscCall(VecCopy(YdotRHS[i], ark->YdotRHS_prev[i])); in TSStep_ARKIMEX()
1379 if (dirk && tab->explicit_first_stage && (ts->steprestart || ts->stepresize)) { in TSStep_ARKIMEX()
1394 …if (!dirk && ts->equation_type >= TS_EQ_IMPLICIT && tab->explicit_first_stage && ts->steprestart) { in TSStep_ARKIMEX()
1444 if (tab->additive && hasE) { in TSStep_ARKIMEX()
1456 if (tab->additive && hasE) { in TSStep_ARKIMEX()
1483 if (i == 0 && tab->explicit_first_stage) { in TSStep_ARKIMEX()
1484 …PetscCheck(tab->stiffly_accurate, PetscObjectComm((PetscObject)ts), PETSC_ERR_SUP, "%s %s is not s… in TSStep_ARKIMEX()
1491 if (i == 0 && tab->explicit_first_stage) { in TSStep_ARKIMEX()
1510 PetscCall(TSEvaluateStep_ARKIMEX(ts, tab->order, ts->vec_sol, NULL)); in TSStep_ARKIMEX()
1514 …PetscCall(TSAdaptCandidateAdd(adapt, tab->name, tab->order, 1, tab->ccfl, (PetscReal)tab->s, PETSC… in TSStep_ARKIMEX()
1563 ARKTableau tab = ark->tableau; in TSAdjointStep_ARKIMEX() local
1564 const PetscInt s = tab->s; in TSAdjointStep_ARKIMEX()
1565 …const PetscReal *At = tab->At, *A = tab->A, *ct = tab->ct, *c = tab->c, *b = tab->b, *bt = tab->bt; in TSAdjointStep_ARKIMEX()
1672 ARKTableau tab = ark->tableau; in TSInterpolate_ARKIMEX() local
1673 PetscInt s = tab->s, pinterp = tab->pinterp, i, j; in TSInterpolate_ARKIMEX()
1677 const PetscReal *Bt = tab->binterpt, *B = tab->binterp; in TSInterpolate_ARKIMEX()
1703 if (tab->additive) { in TSInterpolate_ARKIMEX()
1714 ARKTableau tab = ark->tableau; in TSExtrapolate_ARKIMEX() local
1715 PetscInt s = tab->s, pinterp = tab->pinterp, i, j; in TSExtrapolate_ARKIMEX()
1718 const PetscReal *Bt = tab->binterpt, *B = tab->binterp; in TSExtrapolate_ARKIMEX()
1735 if (tab->additive) { in TSExtrapolate_ARKIMEX()
1746 ARKTableau tab = ark->tableau; in TSARKIMEXTableauReset() local
1749 if (!tab) PetscFunctionReturn(PETSC_SUCCESS); in TSARKIMEXTableauReset()
1751 PetscCall(VecDestroyVecs(tab->s, &ark->Y)); in TSARKIMEXTableauReset()
1752 PetscCall(VecDestroyVecs(tab->s, &ark->YdotI)); in TSARKIMEXTableauReset()
1753 PetscCall(VecDestroyVecs(tab->s, &ark->YdotRHS)); in TSARKIMEXTableauReset()
1754 PetscCall(VecDestroyVecs(tab->s, &ark->Y_prev)); in TSARKIMEXTableauReset()
1755 PetscCall(VecDestroyVecs(tab->s, &ark->YdotI_prev)); in TSARKIMEXTableauReset()
1756 PetscCall(VecDestroyVecs(tab->s, &ark->YdotRHS_prev)); in TSARKIMEXTableauReset()
1780 ARKTableau tab = ark->tableau; in TSAdjointReset_ARKIMEX() local
1783 PetscCall(VecDestroyVecs(tab->s * ts->numcost, &ark->VecsDeltaLam)); in TSAdjointReset_ARKIMEX()
2010 ARKTableau tab = ark->tableau; in TSARKIMEXTableauSetUp() local
2013 PetscCall(PetscMalloc1(2 * tab->s, &ark->work)); in TSARKIMEXTableauSetUp()
2014 PetscCall(VecDuplicateVecs(ts->vec_sol, tab->s, &ark->Y)); in TSARKIMEXTableauSetUp()
2015 PetscCall(VecDuplicateVecs(ts->vec_sol, tab->s, &ark->YdotI)); in TSARKIMEXTableauSetUp()
2016 if (tab->additive) PetscCall(VecDuplicateVecs(ts->vec_sol, tab->s, &ark->YdotRHS)); in TSARKIMEXTableauSetUp()
2018 PetscCall(VecDuplicateVecs(ts->vec_sol, tab->s, &ark->Y_prev)); in TSARKIMEXTableauSetUp()
2019 PetscCall(VecDuplicateVecs(ts->vec_sol, tab->s, &ark->YdotI_prev)); in TSARKIMEXTableauSetUp()
2020 if (tab->additive) PetscCall(VecDuplicateVecs(ts->vec_sol, tab->s, &ark->YdotRHS_prev)); in TSARKIMEXTableauSetUp()
2051 ARKTableau tab = ark->tableau; in TSAdjointSetUp_ARKIMEX() local
2054 PetscCall(VecDuplicateVecs(ts->vecs_sensi[0], tab->s * ts->numcost, &ark->VecsDeltaLam)); in TSAdjointSetUp_ARKIMEX()
2079 if (!dirk && link->tab.additive) count++; in TSSetFromOptions_ARKIMEX()
2080 if (dirk && !link->tab.additive) count++; in TSSetFromOptions_ARKIMEX()
2084 if (!dirk && link->tab.additive) namelist[count++] = link->tab.name; in TSSetFromOptions_ARKIMEX()
2085 if (dirk && !link->tab.additive) namelist[count++] = link->tab.name; in TSSetFromOptions_ARKIMEX()
2124 ARKTableau tab = ark->tableau; in TSView_ARKIMEX() local
2132 PetscCall(PetscFormatRealArray(buf, sizeof(buf), "% 8.6f", tab->s, tab->ct)); in TSView_ARKIMEX()
2137 for (PetscInt i = 0; i < tab->s; i++) { in TSView_ARKIMEX()
2138 PetscCall(PetscFormatRealArray(buf, sizeof(buf), "% 8.6f", tab->s, tab->At + i * tab->s)); in TSView_ARKIMEX()
2141 PetscCall(PetscFormatRealArray(buf, sizeof(buf), "% 8.6f", tab->s, tab->bt)); in TSView_ARKIMEX()
2143 PetscCall(PetscFormatRealArray(buf, sizeof(buf), "% 8.6f", tab->s, tab->bembedt)); in TSView_ARKIMEX()
2147 …PetscCall(PetscViewerASCIIPrintf(viewer, "Stiffly accurate: %s\n", tab->stiffly_accurate ? "yes" :… in TSView_ARKIMEX()
2148 …PetscCall(PetscViewerASCIIPrintf(viewer, "Explicit first stage: %s\n", tab->explicit_first_stage ?… in TSView_ARKIMEX()
2149 …PetscCall(PetscViewerASCIIPrintf(viewer, "FSAL property: %s\n", tab->FSAL_implicit ? "yes" : "no")… in TSView_ARKIMEX()
2151 PetscCall(PetscFormatRealArray(buf, sizeof(buf), "% 8.6f", tab->s, tab->c)); in TSView_ARKIMEX()
2293 PetscCall(PetscStrcmp(link->tab.name, arktype, &match)); in TSARKIMEXSetType_ARKIMEX()
2296 ark->tableau = &link->tab; in TSARKIMEXSetType_ARKIMEX()