Fix some bugs and renames I just introduced

src/caches.jl

@@ -92,8 +92,8 @@ function OrdinaryDiffEq.alg_cache(
     d = is_secondorder_ode ? length(u[1, :]) : length(u)
     D = d * (q + 1)
 
-    FAC = get_covariance_factorization(alg)
-    if FAC isa KroneckerCovariance && !(f.mass_matrix isa UniformScaling)
+    FAC = get_covariance_structure(alg)
+    if FAC isa IsometricKroneckerCovariance && !(f.mass_matrix isa UniformScaling)
         error(
             "The selected algorithm uses an efficient Kronecker-factorized implementation which is incompatible with the provided mass matrix. Try using the `EK1` instead.",
         )

src/covariance_structure.jl

@@ -11,7 +11,7 @@ function get_covariance_structure(alg)
         ) &&
         alg.prior isa IWP
     )
-        return IsometrciKroneckerCovariance()
+        return IsometricKroneckerCovariance()
     else
         return DenseCovariance()
     end

src/kronecker.jl

@@ -1,4 +1,4 @@
-"""
+@doc raw"""
     IsometricKroneckerProduct(left_factor_dim::Int64, right_factor::AbstractMatrix)
 
 Kronecker product of an identity and a generic matrix:

src/priors/common.jl

@@ -45,7 +45,7 @@ function initialize_transition_matrices(
     return A, Q, Ah, Qh, P, PI
 end
 function initialize_transition_matrices(
-    fac::CovarianceFactorization,
+    fac::CovarianceStructure,
     p::AbstractODEFilterPrior,
     dt,
 )

src/priors/iwp.jl

@@ -105,15 +105,20 @@ function discretize(p::IWP, dt::Real)
     return A, Q
 end
 
-function initialize_transition_matrices(::KroneckerCovariance, p::IWP{T}, dt) where {T}
+function initialize_transition_matrices(
+    ::IsometricKroneckerCovariance,
+    p::IWP{T},
+    dt,
+) where {T}
     A, Q = preconditioned_discretize(p)
     P, PI = initialize_preconditioner(p, dt)
     Ah = PI * A * P
     Qh = PSDMatrix(Q.R * PI)
     return A, Q, Ah, Qh, P, PI
 end
 function initialize_transition_matrices(::DenseCovariance, p::IWP{T}, dt) where {T}
-    A, Q, Ah, Qh, P, PI = initialize_transition_matrices(KroneckerCovariance(), p, dt)
+    A, Q, Ah, Qh, P, PI =
+        initialize_transition_matrices(IsometricKroneckerCovariance(), p, dt)
     P, PI = Diagonal(P), Diagonal(PI)
     A, Q, Ah, Qh = Matrix(A), PSDMatrix(Matrix(Q.R)), Matrix(Ah), PSDMatrix(Matrix(Qh.R))
     return A, Q, Ah, Qh, P, PI

src/projection.jl

@@ -23,7 +23,7 @@ function projection(
 end
 
 function projection(
-    ::KroneckerCovariance,
+    ::IsometricKroneckerCovariance,
     d::Integer,
     q::Integer,
     ::Type{elType}=typeof(1.0),