@@ -781,7 +781,7 @@ using Aqua: Aqua
781
781
[rootedtree ([1 ]), rootedtree ([1 ]), rootedtree ([1 , 2 ])],
782
782
[rootedtree ([1 ]), rootedtree ([1 ]), rootedtree ([1 , 2 ])],
783
783
[rootedtree ([1 ]), rootedtree ([1 ]), rootedtree ([1 , 2 ])],
784
- [rootedtree ([1 ]), rootedtree ([1 ]), rootedtree ([1 ]), rootedtree ([1 ])],
784
+ [rootedtree ([1 ]), rootedtree ([1 ]), rootedtree ([1 ]), rootedtree ([1 ])]
785
785
]
786
786
reference_skeletons = [
787
787
rootedtree ([1 ]),
@@ -791,7 +791,7 @@ using Aqua: Aqua
791
791
rootedtree ([1 , 2 , 2 ]),
792
792
rootedtree ([1 , 2 , 3 ]),
793
793
rootedtree ([1 , 2 , 3 ]),
794
- rootedtree ([1 , 2 , 3 , 3 ]),
794
+ rootedtree ([1 , 2 , 3 , 3 ])
795
795
]
796
796
for forest in reference_forests
797
797
sort! (forest)
@@ -855,7 +855,7 @@ using Aqua: Aqua
855
855
([rootedtree ([1 ]), rootedtree ([1 ])], rootedtree ([1 , 2 , 2 ])),
856
856
([rootedtree ([1 ]), rootedtree ([1 ]), rootedtree ([1 ])], rootedtree ([1 , 2 ])),
857
857
([rootedtree ([1 , 2 ]), rootedtree ([1 ]), rootedtree ([1 ])], rootedtree ([1 ])),
858
- ([rootedtree ([1 , 2 , 3 , 2 , 2 ])], rootedtree (Int[])),
858
+ ([rootedtree ([1 , 2 , 3 , 2 , 2 ])], rootedtree (Int[]))
859
859
]
860
860
sort! (reference_forests_and_subtrees)
861
861
@@ -1502,24 +1502,24 @@ using Aqua: Aqua
1502
1502
[
1503
1503
rootedtree ([1 ], Bool[1 ]),
1504
1504
rootedtree ([1 ], Bool[1 ]),
1505
- rootedtree ([1 , 2 ], Bool[0 , 0 ]),
1505
+ rootedtree ([1 , 2 ], Bool[0 , 0 ])
1506
1506
],
1507
1507
[
1508
1508
rootedtree ([1 ], Bool[0 ]),
1509
1509
rootedtree ([1 ], Bool[1 ]),
1510
- rootedtree ([1 , 2 ], Bool[1 , 0 ]),
1510
+ rootedtree ([1 , 2 ], Bool[1 , 0 ])
1511
1511
],
1512
1512
[
1513
1513
rootedtree ([1 ], Bool[0 ]),
1514
1514
rootedtree ([1 ], Bool[1 ]),
1515
- rootedtree ([1 , 2 ], Bool[0 , 1 ]),
1515
+ rootedtree ([1 , 2 ], Bool[0 , 1 ])
1516
1516
],
1517
1517
[
1518
1518
rootedtree ([1 ], Bool[0 ]),
1519
1519
rootedtree ([1 ], Bool[0 ]),
1520
1520
rootedtree ([1 ], Bool[1 ]),
1521
- rootedtree ([1 ], Bool[1 ]),
1522
- ],
1521
+ rootedtree ([1 ], Bool[1 ])
1522
+ ]
1523
1523
]
1524
1524
reference_skeletons = [
1525
1525
rootedtree ([1 ], Bool[1 ]),
@@ -1529,7 +1529,7 @@ using Aqua: Aqua
1529
1529
rootedtree ([1 , 2 , 2 ], Bool[1 , 1 , 0 ]),
1530
1530
rootedtree ([1 , 2 , 3 ], Bool[1 , 0 , 0 ]),
1531
1531
rootedtree ([1 , 2 , 3 ], Bool[1 , 0 , 1 ]),
1532
- rootedtree ([1 , 2 , 3 , 3 ], Bool[1 , 0 , 1 , 0 ]),
1532
+ rootedtree ([1 , 2 , 3 , 3 ], Bool[1 , 0 , 1 , 0 ])
1533
1533
]
1534
1534
for forest in reference_forests
1535
1535
sort! (forest)
@@ -1594,7 +1594,8 @@ using Aqua: Aqua
1594
1594
let order = 4
1595
1595
for t in RootedTreeIterator (order)
1596
1596
@test elementary_weight (t, rk. A, rk. b, rk. c) ≈ elementary_weight (t, rk)
1597
- @test derivative_weight (t, RungeKuttaMethod (rk. A, rk. b, rk. c)) ≈ derivative_weight (t, rk)
1597
+ @test derivative_weight (t, RungeKuttaMethod (rk. A, rk. b, rk. c)) ≈
1598
+ derivative_weight (t, rk)
1598
1599
@test residual_order_condition (t, RungeKuttaMethod (rk. A, rk. b, rk. c)) ≈
1599
1600
residual_order_condition (t, rk)
1600
1601
end
@@ -1631,11 +1632,11 @@ using Aqua: Aqua
1631
1632
# Table II.2.1
1632
1633
As = [
1633
1634
[0 0 ; 1 // 2 1 // 2 ],
1634
- [1 // 2 0 ; 1 // 2 0 ],
1635
+ [1 // 2 0 ; 1 // 2 0 ]
1635
1636
]
1636
1637
bs = [
1637
1638
[1 // 2 , 1 // 2 ],
1638
- [1 // 2 , 1 // 2 ],
1639
+ [1 // 2 , 1 // 2 ]
1639
1640
]
1640
1641
ark = AdditiveRungeKuttaMethod (As, bs)
1641
1642
@@ -1660,11 +1661,11 @@ using Aqua: Aqua
1660
1661
# Table II.2.2
1661
1662
As = [
1662
1663
[0 0 0 ; 5 // 24 1 // 3 - 1 // 24 ; 1 // 6 2 // 3 1 // 6 ],
1663
- [1 // 6 - 1 // 6 0 ; 1 // 6 1 // 3 0 ; 1 // 6 5 // 6 0 ],
1664
+ [1 // 6 - 1 // 6 0 ; 1 // 6 1 // 3 0 ; 1 // 6 5 // 6 0 ]
1664
1665
]
1665
1666
bs = [
1666
1667
[1 // 6 , 2 // 3 , 1 // 6 ],
1667
- [1 // 6 , 2 // 3 , 1 // 6 ],
1668
+ [1 // 6 , 2 // 3 , 1 // 6 ]
1668
1669
]
1669
1670
ark = AdditiveRungeKuttaMethod (As, bs)
1670
1671
@@ -1725,7 +1726,7 @@ using Aqua: Aqua
1725
1726
1471266399579 / 7840856788654 ,
1726
1727
- 4482444167858 / 7529755066697 ,
1727
1728
11266239266428 / 11593286722821 ,
1728
- 1767732205903 / 4055673282236 ,
1729
+ 1767732205903 / 4055673282236
1729
1730
]
1730
1731
rk_explicit = RungeKuttaMethod (A_explicit, b_explicit)
1731
1732
A_implicit = @SArray [0 0 0 0
@@ -1736,7 +1737,7 @@ using Aqua: Aqua
1736
1737
1471266399579 / 7840856788654 ,
1737
1738
- 4482444167858 / 7529755066697 ,
1738
1739
11266239266428 / 11593286722821 ,
1739
- 1767732205903 / 4055673282236 ,
1740
+ 1767732205903 / 4055673282236
1740
1741
]
1741
1742
rk_implicit = RungeKuttaMethod (A_implicit, b_implicit)
1742
1743
ark = AdditiveRungeKuttaMethod ([rk_explicit, rk_implicit])
0 commit comments