ensure type consistency in compute_c_coeffs

src/time_integration/paired_explicit_runge_kutta/methods_PERK3.jl

@@ -9,12 +9,12 @@ using DelimitedFiles: readdlm
 
 # Initialize Butcher array abscissae c for PairedExplicitRK3 based on SSPRK33 base method
 function compute_c_coeffs(num_stages, cS2)
-    c = zeros(num_stages)
+    c = zeros(eltype(cS2), num_stages)
 
     # Last timesteps as for SSPRK33, see motivation in Section 3.3 of
     # https://doi.org/10.1016/j.jcp.2024.113223
-    c[num_stages - 1] = 1.0f0
-    c[num_stages] = 0.5f0
+    c[num_stages - 1] = one(cS2)
+    c[num_stages] = convert(eltype(cS2), 0.5)
 
     # Linear increasing timestep for remainder
     for i in 2:(num_stages - 2)