Skip to content

Commit 2d72836

Browse files
committed
fix tableaus
1 parent 6cb8122 commit 2d72836

File tree

1 file changed

+14
-14
lines changed

1 file changed

+14
-14
lines changed

lib/OrdinaryDiffEqRosenbrock/src/rosenbrock_tableaus.jl

Lines changed: 14 additions & 14 deletions
Original file line numberDiff line numberDiff line change
@@ -253,8 +253,8 @@ function Rodas4Tableau(T, T2)
253253
7.496443313967647 -10.24680431464352 -33.99990352819905 11.70890893206160 0
254254
8.083246795921522 -7.981132988064893 -31.52159432874371 16.31930543123136 -6.058818238834054
255255
]
256-
c = T2[0, 0.386, 0.21, 0.63, 1, 1, 1]
257-
d = T[0.25, -0.1043, 0.1035, -0.0362, 0, 0, 0]
256+
c = T2[0, 0.386, 0.21, 0.63, 1, 1]
257+
d = T[0.25, -0.1043, 0.1035, -0.0362, 0, 0]
258258
H = T[10.12623508344586 -7.487995877610167 -34.80091861555747 -7.992771707568823 1.025137723295662 0
259259
-0.6762803392801253 6.087714651680015 16.43084320892478 24.76722511418386 -6.594389125716872 0]
260260
RodasTableau(A, C, gamma, c, d, H)
@@ -263,12 +263,12 @@ end
263263
function Rodas42Tableau(T, T2)
264264
gamma = convert(T, 1 // 4)
265265
A = T[
266-
0.0 0 0 0 0 0
267-
1.4028884 0 0 0 0 0
268-
0.6581213 -1.320936088384301 0 0 0 0
269-
7.131197445744498 16.02964143958207 -5.561572550509766 0 0 0
270-
22.73885722420363 67.38147284535289 -31.21877493038560 0.7285641833203814 0 0
271-
22.73885722420363 67.38147284535289 -31.21877493038560 0.7285641833203814 1 0
266+
0.0 0 0 0 0 0
267+
1.4028884 0 0 0 0 0
268+
0.6581212688557198 -1.320936088384301 0 0 0 0
269+
7.131197445744498 16.02964143958207 -5.561572550509766 0 0 0
270+
22.73885722420363 67.38147284535289 -31.21877493038560 0.7285641833203814 0 0
271+
22.73885722420363 67.38147284535289 -31.21877493038560 0.7285641833203814 1 0
272272
]
273273
C = T[
274274
0 0 0 0 0
@@ -278,8 +278,8 @@ function Rodas42Tableau(T, T2)
278278
-76.46023087151691 -278.5942120829058 153.9294840910643 10.97101866258358 0
279279
-76.29701586804983 -294.2795630511232 162.0029695867566 23.65166903095270 -7.652977706771382
280280
]
281-
c = T2[0, 0.3507221, 0.2557041, 0.681779, 1, 1, 1]
282-
d = T[0.25, -0.0690221, -0.0009672, -0.087979, 0, 0, 0]
281+
c = T2[0, 0.3507221, 0.2557041, 0.681779, 1, 1]
282+
d = T[0.25, -0.0690221, -0.0009672, -0.087979, 0, 0]
283283
H = T[-38.71940424117216 -135.8025833007622 64.51068857505875 -4.192663174613162 -2.531932050335060 0
284284
-14.99268484949843 -76.30242396627033 58.65928432851416 16.61359034616402 -0.6758691794084156 0]
285285
RodasTableau(A, C, gamma, c, d, H)
@@ -306,8 +306,8 @@ function Rodas4PTableau(T, T2)
306306
34.19095006749676 15.49671153725963 54.74760875964130 14.16005392148534 0
307307
34.62605830930532 15.30084976114473 56.99955578662667 18.40807009793095 -5.714285714285717
308308
]
309-
c = T2[0, 0.75, 0.21, 0.63, 1, 1, 1]
310-
d = T[0.25, -0.5, -0.023504, -0.0362, 0, 0, 0]
309+
c = T2[0, 0.75, 0.21, 0.63, 1, 1]
310+
d = T[0.25, -0.5, -0.023504, -0.0362, 0, 0]
311311
H = T[25.09876703708589 11.62013104361867 28.49148307714626 -5.664021568594133 0 0
312312
1.638054557396973 -0.7373619806678748 8.477918219238990 15.99253148779520 -1.882352941176471 0]
313313
RodasTableau(A, C, gamma, c, d, H)
@@ -331,8 +331,8 @@ function Rodas4P2Tableau(T, T2)
331331
-5.888975457523102 -8.157396617841821 24.805546872612922 12.790401512796979 0
332332
-4.408651676063871 -6.692003137674639 24.625568527593117 16.627521966636085 -5.714285714285718
333333
]
334-
c = T2[0, 0.75, 0.321448134013046, 0.519745732277726, 1, 1, 1]
335-
d = T[0.25, -0.5, -0.189532918363016, 0.085612108792769, 0, 0, 0]
334+
c = T2[0, 0.75, 0.321448134013046, 0.519745732277726, 1, 1]
335+
d = T[0.25, -0.5, -0.189532918363016, 0.085612108792769, 0, 0]
336336
H = [-5.323528268423303 -10.042123754867493 17.175254928256965 -5.079931171878093 -0.016185991706112 0
337337
6.984505741529879 6.914061169603662 -0.849178943070653 18.104410789349338 -3.516963011559032 0]
338338
RodasTableau(A, C, gamma, c, d, H)

0 commit comments

Comments
 (0)