| /petsc/src/ts/impls/explicit/rk/ ! |
| H A D | rk.c | 197 const PetscReal bembed[2] = {RC(1.0), 0}; in TSRKRegisterAll() local
198 PetscCall(TSRKRegister(TSRK2A, 2, 2, &A[0][0], b, NULL, bembed, 0, NULL)); in TSRKRegisterAll()
225 …const PetscReal bembed[4] = {RC(7.0) / RC(24.0), RC(1.0) / RC(4.0), RC(1.0) / RC(3.0), RC(1.0) / R… in TSRKRegisterAll() local
226 PetscCall(TSRKRegister(TSRK3BS, 3, 4, &A[0][0], b, NULL, bembed, 0, NULL)); in TSRKRegisterAll()
248 …const PetscReal bembed[6] = {RC(25.0) / RC(216.0), 0, RC(1408.0) / RC(2565.0), RC(2197.0) / RC(410… in TSRKRegisterAll() local
249 PetscCall(TSRKRegister(TSRK5F, 5, 6, &A[0][0], b, NULL, bembed, 0, NULL)); in TSRKRegisterAll()
262 …const PetscReal bembed[7] = {RC(5179.0) / RC(57600.0), 0, RC(7571.0) / RC(16695.0), RC(393.0) … in TSRKRegisterAll() local
272 PetscCall(TSRKRegister(TSRK5DP, 5, 7, &A[0][0], b, NULL, bembed, 5, binterp[0])); in TSRKRegisterAll()
286 …const PetscReal bembed[8] = {RC(2479.0) / RC(34992.0), 0, RC(123.0) / RC(416.0), RC(612941.0) / RC… in TSRKRegisterAll() local
287 PetscCall(TSRKRegister(TSRK5BS, 5, 8, &A[0][0], b, NULL, bembed, 0, NULL)); in TSRKRegisterAll()
[all …]
|
| H A D | rk.h | 10 PetscReal *bembed; /* Embedded formula of order one less (order-1) */ member
|
| /petsc/src/ts/impls/arkimex/ ! |
| H A D | arkimex.c | 451 const PetscReal bembed[2] = {RC(1.0), RC(0.0)}; in TSARKIMEXRegisterAll() local
452 PetscCall(TSDIRKRegister(TSDIRKS212, 2, 2, &A[0][0], b, NULL, bembed, 1, b)); in TSARKIMEXRegisterAll()
462 const PetscReal bembed[2] = {RC(0.5), RC(0.5)}; in TSARKIMEXRegisterAll() local
463 PetscCall(TSDIRKRegister(TSDIRKES122SAL, 1, 2, &A[0][0], b, NULL, bembed, 1, b)); in TSARKIMEXRegisterAll()
476 …const PetscReal bembed[3] = {(RC(1.0) - s2 / RC(4.0)) / RC(3.0), (RC(3.0) * s2 / RC(4.0) + RC(1.0)… in TSARKIMEXRegisterAll() local
477 PetscCall(TSDIRKRegister(TSDIRKES213SAL, 2, 3, &A[0][0], b, NULL, bembed, 1, b)); in TSARKIMEXRegisterAll()
511 const PetscReal bembed[4] = {RC(1.0) - bh2 - bh3 - bh4, bh2, bh3, bh4}; in TSARKIMEXRegisterAll() local
512 PetscCall(TSDIRKRegister(TSDIRKES324SAL, 3, 4, &A[0][0], b, NULL, bembed, 1, b)); in TSARKIMEXRegisterAll()
535 …const PetscReal bembed[5] = {RC(23705.0) / RC(104328.0), RC(29720.0) / RC(91287.0), RC(4225.0) / R… in TSARKIMEXRegisterAll() local
536 PetscCall(TSDIRKRegister(TSDIRKES325SAL, 3, 5, &A[0][0], b, NULL, bembed, 1, b)); in TSARKIMEXRegisterAll()
[all …]
|
| H A D | arkimex.h | 13 PetscReal *bembedt, *bembed; /* Embedded formula of order one less (order-1) */ member
|
| H A D | fsarkimex.c | 166 for (j = 0; j < s; j++) w[j] = h * tab->bembed[j]; in TSEvaluateStep_ARKIMEX_FastSlowSplit()
181 for (j = 0; j < s; j++) w[j] = h * (tab->bembed[j] - tab->b[j]); in TSEvaluateStep_ARKIMEX_FastSlowSplit()
|
| /petsc/src/ts/tutorials/ ! |
| H A D | ex19.c | 122 …678118654752440, 0.50000000000000000000, 0.29289321881345247560}}, *bembedt = NULL, *bembed = NULL; in RegisterMyARK2() local
123 …ster("myark2", 2, 3, &At[0][0], NULL, NULL, &A[0][0], NULL, NULL, bembedt, bembed, 0, NULL, NULL)); in RegisterMyARK2()
|
| H A D | ex16.c | 149 …678118654752440, 0.50000000000000000000, 0.29289321881345247560}}, *bembedt = NULL, *bembed = NULL; in RegisterMyARK2() local
150 …ster("myark2", 2, 3, &At[0][0], NULL, NULL, &A[0][0], NULL, NULL, bembedt, bembed, 0, NULL, NULL)); in RegisterMyARK2()
|
| /petsc/src/ts/impls/rosw/ ! |
| H A D | rosw.c | 33 PetscReal *bembed; /* Step completion table for embedded method of order one less */ member
764 PetscCall(PetscFree2(t->bembed, t->bembedt)); in TSRosWRegisterDestroy()
830 …cReal A[], const PetscReal Gamma[], const PetscReal b[], const PetscReal bembed[], PetscInt pinter… in TSRosWRegister() argument
842 if (bembed) PetscAssertPointer(bembed, 7); in TSRosWRegister()
856 if (bembed) { in TSRosWRegister()
857 PetscCall(PetscMalloc2(s, &t->bembed, s, &t->bembedt)); in TSRosWRegister()
858 PetscCall(PetscArraycpy(t->bembed, bembed, s)); in TSRosWRegister()
924 if (bembed) { in TSRosWRegister()
926 for (j = 0; j < s; j++) t->bembedt[i] += t->bembed[j] * t->GammaInv[j * s + i]; in TSRosWRegister()
1675 ts->default_adapt_type = ros->tableau->bembed ? TSADAPTBASIC : TSADAPTNONE; in TSRosWSetType_RosW()
|