Simplify to_binary_operation

Author: asterycs

src/simplify.jl

@@ -96,6 +96,15 @@ function get_last_letter(indices::IndexList)
     return current_last
 end
 
+
+function to_binary_operation(op::Op, terms::AbstractArray) where {Op}
+    if length(terms) == 1
+        return first(terms)
+    end
+
+    return BinaryOperation{Op}(to_binary_operation(op, terms[1:(end-1)]), terms[end])
+end
+
 function simplify(::Mult, arg1::BinaryOperation{Mult}, arg2::KrD)
     if is_diag(arg1)
         d = get_diag_delta(arg1)

src/std.jl

@@ -601,25 +601,6 @@ function to_standard(arg::Real)
     return arg
 end
 
-function to_binary_operation(op::Op, terms::AbstractArray) where {Op}
-    binop = nothing
-
-    for t ∈ terms
-        if isnothing(binop)
-            binop = t
-            continue
-        end
-
-        binop = BinaryOperation{Op}(binop, t)
-    end
-
-    return binop
-end
-
-function to_binary_operation(op::Op, term) where {Op}
-    return term
-end
-
 # Recursive adjoint
 function radjoint(arg::T) where {T<:UnaryOperation}
     return T(arg.arg')