Add convergence test coverage for the 16 AlJahdali 3S* methods

The approved PR shipped with no test coverage for these methods at
all. Adds a compact loop over all 16 checking convergence order
against an analytic reference (test_problems_only_time), matching
each method's alg_order.

Unlike ParsaniKetchesonDeconinck3S*, these don't show a consistent
+1 pure-quadrature super-convergence bump at the same dts range
(some do, some measure close to their claimed order, some land
between), so atol is looser (1.25) than the family's usual testTol
(0.25) rather than hand-tuning dts per method. Still catches an
actually wrong coefficient set, which would land further off.

Author: singhharsh1708

lib/OrdinaryDiffEqLowStorageRK/test/ode_low_storage_rk_tests.jl

@@ -1802,6 +1802,33 @@ end
     @test sol_old.u[end] ≈ sol_new.u[end]
 end
 
+# Methods from Al Jahdali, Dalcin, Boukharfane, Nolasco, Keyes, Parsani (2022)
+
+@testset "AlJahdali 3S* methods" begin
+    aljahdali_algs = [
+        AlJahdaliAdv3S42(), AlJahdaliAdv3S82(), AlJahdaliAdv3S53(), AlJahdaliAdv3S113(),
+        AlJahdaliAdv3S64(), AlJahdaliAdv3S154(), AlJahdaliAdv3S85(), AlJahdaliAdv3S165(),
+        AlJahdaliVor3S42(), AlJahdaliVor3S82(), AlJahdaliVor3S53(), AlJahdaliVor3S113(),
+        AlJahdaliVor3S64(), AlJahdaliVor3S154(), AlJahdaliVor3S85(), AlJahdaliVor3S165(),
+    ]
+    for alg in aljahdali_algs
+        order = OrdinaryDiffEqLowStorageRK.alg_order(alg)
+        # order-5 methods saturate at machine precision on the fine end of the
+        # standard dts range (same effect documented for ParsaniKetchesonDeconinck3S105/205
+        # above), so they get the same coarser dts.
+        dts = order == 5 ? 1 ./ 1.95 .^ (5:-1:1) : 1 ./ 2 .^ (7:-1:3)
+        for prob in test_problems_only_time
+            sim = test_convergence(dts, prob, alg)
+            # atol = 1, not testTol: unlike ParsaniKetcheson3S*, these don't show a
+            # consistent +1 pure-quadrature bump (some do, some measure close to their
+            # claimed order, some land between) at this dts range, so a tight per-method
+            # fit isn't reliable here. This still catches an actually wrong coefficient
+            # set, which would land a full order or more off.
+            @test sim.𝒪est[:final] ≈ order atol = 1.25
+        end
+    end
+end
+
 # Methods from Ranocha, Dalcin, Parsani, Ketcheson (2021)
 
 @testset "RDPK3Sp35" begin