Diffeqpy, using ensembling

[ChrisRackauckas]

@Abolfazl_Ziaeemehr

julia_f = Main.eval("""
function f!(du, u, p, t)
    du .= 1.01 * u
end
""")

prob_func = Main.eval("""
function prob_func(prob, i, repeat)
    remake(prob, u0 = rand() * prob.u0)
    end
""")

# def prob_func(prob, i, repeat):
#     de.remake(prob, u0=rand() * prob.u0)


if __name__ == "__main__":

    tspan = (0.0, 10.0)
    prob = de.ODEProblem(julia_f, [0.5], tspan)
    ensemble_prob = de.EnsembleProblem(
        prob,
        prob_func=prob_func)
    
    sol = de.solve(ensemble_prob,
                   de.Tsit5(),
                   de.EnsembleThreads(),
                   trajectories=2)

Should be simpler and work. Using arrays instead of scalars for ODEs is probably a lot more common too.