Merge pull request #461 from SciML/reshape

fix reshaped ldiv

Author: ChrisRackauckas

src/operations.jl

@@ -66,6 +66,7 @@ Base.convert(::Type{Operation}, x::Number) = Operation(identity, Expression[Cons
 Base.convert(::Type{Operation}, x::Operation) = x
 Base.convert(::Type{Operation}, x::Expression) = Operation(identity, Expression[x])
 Operation(x) = convert(Operation, x)
+Operation(x::Operation) = x
 Base.Symbol(O::Operation) = Symbol(convert(Variable,O))
 Base.convert(::Type{Symbol},O::Operation) = Symbol(convert(Variable,O))
 

test/operation_overloads.jl

@@ -46,3 +46,14 @@ J = [Dx(eqs[1].rhs) Dy(eqs[1].rhs) Dz(eqs[1].rhs)
 J = expand_derivatives.(J)
 using LinearAlgebra
 luJ = lu(J,Val(false))
+
+using ModelingToolkit
+@variables M[1:2,1:2]
+inv(M)
+
+@variables b[1:2]
+M = [1 0; 0 2]
+M \ b
+M \ reshape(b,2,1)
+M = [1 1; 0 2]
+M \ reshape(b,2,1)