fix compats and undefined memory

Author: ChrisRackauckas

Project.toml

@@ -22,8 +22,8 @@ BandedMatrices = "0.15.11"
 BlockBandedMatrices = "0.8.6, 0.9"
 DiffEqBase = "6.4.1"
 ForwardDiff = "0.10"
-LazyArrays = "0.16.12"
-LazyBandedMatrices = "0.2.11"
+LazyArrays = "0.17"
+LazyBandedMatrices = "0.3"
 ModelingToolkit = "0.10.0, 1.0, 2.0, 3.0"
 NNlib = "0.6, 0.7"
 StaticArrays = "0.10, 0.11, 0.12"

src/derivative_operators/multi_dim_bc_operators.jl

@@ -16,7 +16,7 @@ slice_rmul lets you multiply each vector like strip of an array `u` with a linea
 """
 function slice_rmul(A::AbstractDiffEqLinearOperator, u::AbstractArray{T,N}, dim::Int) where {T,N}
     @assert N != 1
-    u_temp = similar(u)
+    u_temp = zero(u)
 
     _slice_rmul!(u_temp, A, u, dim, CartesianIndices(axes(u)[1:dim-1]), CartesianIndices(axes(u)[(dim+1):end]))
 

src/matrixfree_operators.jl

@@ -24,7 +24,7 @@ end
   return n <= length(M.size) ? M.size[n] : 1
 end
 @inline ==(M1::MatrixFreeOperator, M2::MatrixFreeOperator) = M1.f == M2.f && M1.args == M2.args && M1.size == M2.size && M1.opnorm == M2.opnorm && M1.ishermitian == M2.ishermitian
-@inline Base.:*(A::MatrixFreeOperator, X::AbstractVecOrMat) = mul!(similar(X), A, X)
+@inline Base.:*(A::MatrixFreeOperator, X::AbstractVecOrMat) = mul!(zero(X), A, X)
 
 # Overloading LinearAlgebra functions
 LinearAlgebra.ishermitian(M::MatrixFreeOperator) = M.ishermitian
@@ -60,7 +60,7 @@ end
 
 function (M::MatrixFreeOperator{F,N})(u, p, t) where {F,N}
   update_coefficients!(M,u,p,t)
-  du = similar(u)
+  du = zero(u)
   if isconstant(M)
     M.f(du, u, p)
   else