Merge pull request #874 from SciML/myb/tab_info

Add `stability_size`, `fsal`, and dense output fields

Author: ChrisRackauckas

src/tableaus.jl

@@ -3,20 +3,26 @@ $(TYPEDEF)
 
 Holds a tableau which defines an explicit Runge-Kutta method.
 """
-mutable struct ExplicitRKTableau{MType <: AbstractMatrix, VType <: AbstractVector, fsal} <:
+mutable struct ExplicitRKTableau{MType <: AbstractMatrix, VType <: AbstractVector, S} <:
                ODERKTableau
     A::MType
     c::VType
     α::VType
     αEEst::VType
+    d::VType # dense output coefficients
     stages::Int
     order::Int
     adaptiveorder::Int #The lower order of the pair. Only used for adaptivity.
+    fsal::Bool
+    stability_size::S
 end
 function ExplicitRKTableau(A::MType, c::VType, α::VType, order;
-                           adaptiveorder = 0, αEEst = VType(),
-                           fsal = false) where {MType, VType}
-    ExplicitRKTableau{MType, VType, fsal}(A, c, α, αEEst, length(α), order, adaptiveorder)
+                           adaptiveorder = 0, αEEst = similar(α, 0),
+                           fsal = false, stability_size = 0.0,
+                           d = similar(α, 0)) where {MType, VType}
+    S = typeof(stability_size)
+    ExplicitRKTableau{MType, VType, S}(A, c, α, αEEst, d, length(α), order, adaptiveorder,
+                                       fsal, stability_size)
 end
 
 """