DArray: Implement in-place matrix-matrix multiply

Co-authored-by: Julian P Samaroo <[email protected]>

Author: Rabab53

src/Dagger.jl

@@ -72,6 +72,7 @@ include("array/matrix.jl")
 include("array/sparse_partition.jl")
 include("array/sort.jl")
 include("array/linalg.jl")
+include("array/mul.jl")
 include("array/cholesky.jl")
 
 # Visualization

src/array/matrix.jl

@@ -35,10 +35,9 @@ end
 size(x::MatMul) = mul_size(x.a, x.b)
 MatMul(a,b) =
     MatMul{promote_type(eltype(a), eltype(b)), length(mul_size(a,b))}(a,b)
-(*)(a::ArrayOp, b::ArrayOp) = _to_darray(MatMul(a,b))
+
 # Bonus method for matrix-vector multiplication
 (*)(a::ArrayOp, b::Vector) = _to_darray(MatMul(a,PromotePartition(b)))
-(*)(a::AbstractArray, b::ArrayOp) = _to_darray(MatMul(PromotePartition(a), b))
 
 function (*)(a::ArrayDomain{2}, b::ArrayDomain{2})
     if size(a, 2) != size(b, 1)