explicitly zero caches

Author: ChrisRackauckas

src/composite_operators.jl

@@ -53,6 +53,7 @@ struct DiffEqOperatorCombination{T,O<:Tuple{Vararg{AbstractDiffEqLinearOperator{
     T = eltype(ops[1])
     if cache == nothing
       cache = Vector{T}(undef, size(ops[1], 1))
+      fill!(cache,0)
     end
     # TODO: safecheck dimensions
     new{T,typeof(ops),typeof(cache)}(ops, cache)
@@ -96,7 +97,9 @@ struct DiffEqOperatorComposition{T,O<:Tuple{Vararg{AbstractDiffEqLinearOperator{
       # Construct a list of caches to be used by mul! and ldiv!
       caches = []
       for op in ops[1:end-1]
-        push!(caches, Vector{T}(undef, size(op, 1)))
+        tmp = Vector{T}(undef, size(op, 1))
+        fill!(tmp,0)
+        push!(caches,tmp)
       end
       caches = tuple(caches...)
     end

src/derivative_operators/abstract_operator_functions.jl

@@ -166,7 +166,7 @@ end
 ################################################################################
 
 function *(coeff_func::Function, A::DerivativeOperator{T,N,Wind}) where {T,N,Wind}
-    coefficients = A.coefficients === nothing ? Vector{T}(undef, A.len) : A.coefficients
+    coefficients = A.coefficients === nothing ? zeros(T,A.len) : A.coefficients
 
     DerivativeOperator{T,N,Wind,typeof(A.dx),typeof(A.stencil_coefs),
                        typeof(A.low_boundary_coefs),typeof(coefficients),

src/derivative_operators/derivative_operator.jl

@@ -103,7 +103,7 @@ function CenteredDifference{N}(derivative_order::Int,
                                                                   calculate_weights(derivative_order, high_boundary_x[end-i], high_boundary_x)) for i in boundary_point_count:-1:1]
     high_boundary_coefs      = convert(SVector{boundary_point_count},_high_boundary_coefs)
 
-    coefficients            = coeff_func isa Nothing ? nothing : Vector{T}(undef,len)
+    coefficients            = coeff_func isa Nothing ? nothing : zeros(T,len)
 
     DerivativeOperator{T,N,false,typeof(dx),typeof(stencil_coefs),
         typeof(low_boundary_coefs),typeof(coefficients),
@@ -169,7 +169,7 @@ function UpwindDifference{N}(derivative_order::Int,
     _high_boundary_coefs     = SVector{boundary_stencil_length, T}[convert(SVector{boundary_stencil_length, T}, ((-1/dx)^derivative_order) * calculate_weights(derivative_order, oneunit(T)*x0, high_boundary_x)) for x0 in R_boundary_deriv_spots]
     high_boundary_coefs = convert(SVector{boundary_point_count},_high_boundary_coefs)
 
-    coefficients = Vector{T}(undef,len)
+    coefficients = zeros(T,len)
     if coeff_func != nothing
         compute_coeffs!(coeff_func, coefficients)
     end
@@ -215,7 +215,7 @@ function UpwindDifference{N}(derivative_order::Int,
     high_boundary_coefs = [_upwind_coefs ; _downwind_coefs]
 
     # Compute coefficients
-    coefficients = Vector{T}(undef,len)
+    coefficients = zeros(T,len)
     if coeff_func != nothing
         compute_coeffs!(coeff_func, coefficients)
     end

test/upwind_operators_interface.jl

@@ -586,6 +586,7 @@ end
       L = UpwindDifference(1,2, 0.1, N, 1.0)
       analyticL = 10.0*analyticOneTwoPos()
       x = rand(7)
+      @show x
 
       # Test that multiplication agrees with analytic multiplication
       @test L*x ≈ analyticL*x