- Join over
**1.5M+ people** - Join over
**100K+ communities** - Free
**without limits** - Create
**your own community**

- 10:02fredrikekre commented #529
- 10:01ChrisRackauckas commented #529
- 10:00fredrikekre commented #529
- 10:00fredrikekre commented #529
- 09:59ChrisRackauckas commented #529
- 09:57fredrikekre commented #529
- 09:56ChrisRackauckas commented #529
- 09:30mortenpi commented #529
- 09:29fredrikekre commented #529
- 09:28mortenpi commented #529
- 09:26
ChrisRackauckas on master

Update README.md (compare)

- 09:26ChrisRackauckas commented #529
- 09:26ChrisRackauckas closed #529
- 09:26ChrisRackauckas commented #529
- 09:25fredrikekre commented #529
- 09:24mortenpi commented #529
- 09:23ChrisRackauckas edited #529
- 09:23ChrisRackauckas commented #529
- 09:22
ChrisRackauckas on master

Update README.md (compare)

- 09:14ChrisRackauckas opened #529

But in the middle, 64-~500 bits, ArbReals is the way to go.

I kow Sanders and Benet have been on the emails with Jeffrey Sarnoff on ArbFloats/ArbReals, so I wonder if that's why.

Where's do I have the error?

:worried:

Thanks for letting me know!

I tried the simple code

'using DifferentialEquations

function vanDerPolExample(u₀=[0,sqrt(3)])

f1(u,t) = (1-u[2].^2)*u[1] - u[2]

f2(u,t) = u[1]

f(u,t) = [f1(u,t);f2(u,t)]

return(ODEProblem(f,u₀))

end

prob = vanDerPolExample()

Δt = 1//2^(4) #The initial timestepping size. It will automatically assigned if not given.

tspan = [0,20] # The timespan. This is the default if not given.

sol = solve(prob::ODEProblem,tspan,Δt=Δt,alg=:RK4)'

but it gives the error (line 12 is the solve):

ERROR: LoadError: BoundsError

in getindex at ./number.jl:21 [inlined]

in (::#f1#10)(::Rational{Int64}, ::Array{Float64,1}) at /Users/rveltz/Downloads/ode.jl:4

in (::#f#12{#f1#10,#f2#11})(::Rational{Int64}, ::Array{Float64,1}) at /Users/rveltz/Downloads/ode.jl:6

in (::DifferentialEquations.##500#507{DifferentialEquations.ODEProblem{Array{Float64,1},Float64}})(::Rational{Int64}, ::Array{Float64,1}, ::Array{Float64,1}) at /Users/rveltz/.julia/v0.5/DifferentialEquations/src/ode/ode_solve.jl:79

in ode_solve(::DifferentialEquations.ODEIntegrator{:RK4,Array{Float64,1},Float64,2,Rational{Int64}}) at /Users/rveltz/.julia/v0.5/DifferentialEquations/src/ode/ode_integrators.jl:376

in #solve#499(::Array{Any,1}, ::Function, ::DifferentialEquations.ODEProblem{Array{Float64,1},Float64}, ::Array{Int64,1}) at /Users/rveltz/.julia/v0.5/DifferentialEquations/src/ode/ode_solve.jl:134

in (::DifferentialEquations.#kw##solve)(::Array{Any,1}, ::DifferentialEquations.#solve, ::DifferentialEquations.ODEProblem{Array{Float64,1},Float64}, ::Array{Int64,1}) at ./<missing>:0

in include_from_node1(::String) at ./loading.jl:426

in include_from_node1(::String) at /Applications/Julia-0.5.app/Contents/Resources/julia/lib/julia/sys.dylib:?

while loading /Users/rveltz/Downloads/ode.jl, in expression starting on line 12

Do you have any idea?

Thank you,

'using DifferentialEquations

function vanDerPolExample(u₀=[0,sqrt(3)])

f1(u,t) = (1-u[2].^2)*u[1] - u[2]

f2(u,t) = u[1]

f(u,t) = [f1(u,t);f2(u,t)]

return(ODEProblem(f,u₀))

end

prob = vanDerPolExample()

Δt = 1//2^(4) #The initial timestepping size. It will automatically assigned if not given.

tspan = [0,20] # The timespan. This is the default if not given.

sol = solve(prob::ODEProblem,tspan,Δt=Δt,alg=:RK4)'

but it gives the error (line 12 is the solve):

ERROR: LoadError: BoundsError

in getindex at ./number.jl:21 [inlined]

in (::#f1#10)(::Rational{Int64}, ::Array{Float64,1}) at /Users/rveltz/Downloads/ode.jl:4

in (::#f#12{#f1#10,#f2#11})(::Rational{Int64}, ::Array{Float64,1}) at /Users/rveltz/Downloads/ode.jl:6

in (::DifferentialEquations.##500#507{DifferentialEquations.ODEProblem{Array{Float64,1},Float64}})(::Rational{Int64}, ::Array{Float64,1}, ::Array{Float64,1}) at /Users/rveltz/.julia/v0.5/DifferentialEquations/src/ode/ode_solve.jl:79

in ode_solve(::DifferentialEquations.ODEIntegrator{:RK4,Array{Float64,1},Float64,2,Rational{Int64}}) at /Users/rveltz/.julia/v0.5/DifferentialEquations/src/ode/ode_integrators.jl:376

in #solve#499(::Array{Any,1}, ::Function, ::DifferentialEquations.ODEProblem{Array{Float64,1},Float64}, ::Array{Int64,1}) at /Users/rveltz/.julia/v0.5/DifferentialEquations/src/ode/ode_solve.jl:134

in (::DifferentialEquations.#kw##solve)(::Array{Any,1}, ::DifferentialEquations.#solve, ::DifferentialEquations.ODEProblem{Array{Float64,1},Float64}, ::Array{Int64,1}) at ./<missing>:0

in include_from_node1(::String) at ./loading.jl:426

in include_from_node1(::String) at /Applications/Julia-0.5.app/Contents/Resources/julia/lib/julia/sys.dylib:?

while loading /Users/rveltz/Downloads/ode.jl, in expression starting on line 12

Do you have any idea?

Thank you,

Yes. It's (t,u,du). Julia usually has the in-place moded variable first, but all the Fortran codes that we have wrapped have them last. So it's slightly more efficient for the Fortran codes to just have everything be (t,u,du). ODE.jl for some reason already had that ordering, so that means all DE packages (that I know of) are standardized on it now.

So it's different from the Julia standard, but at least all DE packages are the same? I'm not sure about it still, but that's mostly because of the Fortran wrappers.

So it's now

`prob = prob_ode_linear`

The docs should all be updated for that.

An example for using Sundials is here: https://github.com/ChrisRackauckas/DifferentialEquations.jl/blob/master/test/ode/Sundials_tests.jl

That needs to be updated though.

One sec on that.

When using sundials, I always get:

ERROR: UndefVarError: cvode_fulloutput not defined

in #solve#499(::Array{Any,1}, ::Function, ::DifferentialEquations.ODEProblem{Float64,Float64}, ::Array{Int64,1}) at /Users/rveltz/.julia/v0.5/DifferentialEquations/src/ode/ode_solve.jl:255

in (::DifferentialEquations.#kw##solve)(::Array{Any,1}, ::DifferentialEquations.#solve, ::DifferentialEquations.ODEProblem{Float64,Float64}, ::Array{Int64,1}) at ./<missing>:0

Oh, did you check the conditional deps docs? There's a PR that needs to go through, so right now it's compatible with the branch for that PR (since it's right about to merge and change everything). Either use the branch `handles`

from `alyst/Sundials.jl`

, i.e.

```
Pkg.clone("https://github.com/alyst/Sundials.jl")
Pkg.checkout("Sundials","handles")
```

or just wait for the PR. It's here: JuliaDiffEq/Sundials.jl#67. It should've gone through, but Travis hasn't been working and we can't merge until it does...

It'll go through whenever Travis works again. Likely within a day or so? Lots of things are changing there, but it should settle down once that PR goes through.

Hi there. I'm trying to get this package to replace ODE.jl in my code, since I really like the design of it. My first attempt results in the following error:

```
ERROR: LoadError: LoadError: ArgumentError: argument is not a generic function
in methods at reflection.jl:180
in numparameters at /Users/cbinz/.julia/v0.4/DifferentialEquations/src/general/problems.jl:460
in call at /Users/cbinz/.julia/v0.4/DifferentialEquations/src/general/problems.jl:395
...
```

Here's the line that throws the error:

` p = ODEProblem( (t,y)->vopeom(t, y, sat, propOpt), x)`

On v0.4, you'd want to use "generic functions" anyways because they are a lot faster. Anonymous functions (in v0.4) have a lot of overhead (they are kept as code, and don't compile). So a quick fix would be:

```
f(t,y) = vopeom(t, y, sat, propOpt)
p = ODEProblem(f, x)
```

(In v0.5 you don't have to worry about this). But if you're on Julia v0.4 with DifferentialEquations v0.3.0 this issue shouldn't be there (according to Travis and AppVoyer?).

@/all Hello all. This chatroom will be deprecated. Please use the shared JuliaDiffEq chat for further discussions. Thanks! https://gitter.im/JuliaDiffEq/Lobby

Can DifferentialEquations.jl support Cudanative???

Hey, this is the wrong channel. This is our older channel. Our current channel is https://gitter.im/JuliaDiffEq/Lobby

Can DifferentialEquations.jl support Cudanative???

Yes it did in v0.6, but we do need to fix it.

Basically, the broadcast changes are going to be a bit complicated to unravel with some of the setups