Robbie Rosati

I think this reproduces it ```julia using AutoOptimize using DifferentialEquations function odefunc1(du,u,p,t) du[:] .= u end prob = ODEProblem(odefunc1,[0.1,0.2],(0.0,10.0)) prob,_ = auto_optimize(prob) sol = solve(prob) numfunc(t) = sol(t)[end] function odefunc2(du,u,p,t)...

That's a good idea, and if possible could make this a lot easier. I'm not quite sure how you could enforce indices matching, etc in Symbolics.jl, but maybe I just...

Wow, awesome work. Yeah I think a PR should work well, as far as I know nobody is using this repo in production so we can just commit to master....

This is a syntax issue, this seems to work just fine: ```julia replace_with_arrays(A(i)*A(j)*field.metric(-i,-j), Dict(A(-i) => [0, 1], field => [1 0; 0 1]) ) ``` Still, should we modify `replace_with_arrays`...

I haven't tested very much in Pluto. This is probably an issue with my override of `Base.show`

Thanks for opening all of these. The package is not very user friendly yet, sorry about that. I never implemented `epsilon` from sympy.tensor.tensor, mostly because I don't like how it...

To quote arrays the way you describe, I've been doing this: ```julia julia> x = [:a , :b] 2-element Array{Symbol,1}: :a :b julia> arr = :([ $(x...) ]) :([a, b])...