Merge pull request #696 from JuliaSymbolics/s/ref-broadcast

RefValue handling in symtype inference

Author: shashi

src/array-lib.jl

@@ -159,11 +159,12 @@ function Broadcast.materialize(bc::Broadcast.Broadcasted{SymBroadcast})
             subs = map(i-> extruded[i] && isonedim(x, i) ?
                        1 : subscripts[i], 1:ndims(x))
             x[subs...]
+        elseif x isa Base.RefValue
+            x[]
         else
             x
         end
     end
-
     expr = term(bc.f, expr_args′...) # Imagine x .=> y -- if you don't have a term
                                      # then you get pairs, and index matcher cannot
                                      # recurse into pairs

src/arrays.jl

@@ -558,21 +558,6 @@ function replace_by_scalarizing(ex, dict)
                  scalarize(~x, (map(j->substitute(j, dict), ~~i)...,)))
 
     simterm = (x, f, args; kws...) -> begin
-        if f === Base.literal_pow && length(args) == 3
-            #=
-            julia> @variables u[1:3]
-            1-element Vector{Symbolics.Arr{Num, 1}}:
-             u[1:3]
-
-            julia> u.^1
-            (broadcast(literal_pow, Base.RefValue{typeof(^)}(^), u, Base.RefValue{Val{1}}(Val{1}())))[1:3]
-            =#
-            base = args[2]
-            exp = val2num(only(args[3]))
-            f = only(args[1])
-            args = [base, exp]
-        end
-
         if metadata(x) !== nothing
             similarterm(x, f, args; metadata=metadata(x))
         else

test/arrays.jl

@@ -113,6 +113,12 @@ getdef(v) = getmetadata(v, Symbolics.VariableDefaultValue)
 
     ##653
     Symbolics.scalarize(inv(A)[1,1])
+
+    # ModelingToolkit.jl#1736
+    #
+    @variables t F(t)[1:1]
+
+    @test isequal(collect(F ./ t), [F[1] / t])
 end
 
 n = 2