Fixed un-initialization related bug in mul! for block matrices

gridap · Aug 21, 2020 · ed8d3bf · ed8d3bf
1 parent ee5bf7e
commit ed8d3bf
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/src/Arrays/BlockArraysCoo.jl b/src/Arrays/BlockArraysCoo.jl
@@ -318,10 +318,12 @@ function Base.:*(a::BlockMatrixCoo,b::BlockMatrixCoo)
 end
 
 function LinearAlgebra.mul!(c::BlockVectorCoo,a::BlockMatrixCoo,b::BlockVectorCoo)
+  fill!(c,zero(eltype(c)))
   mul!(c,a,b,1,0)
 end
 
 function LinearAlgebra.mul!(c::BlockMatrixCoo,a::BlockMatrixCoo,b::BlockMatrixCoo)
+  fill!(c,zero(eltype(c)))
   mul!(c,a,b,1,0)
 end