Where communities thrive


  • Join over 1.5M+ people
  • Join over 100K+ communities
  • Free without limits
  • Create your own community
People
Repo info
Activity
  • Aug 05 15:55
    tknopp commented #503
  • Aug 05 15:50
    ChrisRackauckas commented #503
  • Aug 05 15:47
    tknopp commented #503
  • Aug 05 15:28
    ChrisRackauckas closed #884
  • Aug 05 15:28
    ChrisRackauckas commented #884
  • Aug 05 15:26
    ChrisRackauckas closed #800
  • Aug 05 15:25
    ChrisRackauckas closed #687
  • Aug 05 15:25
    ChrisRackauckas closed #636
  • Aug 05 15:25
    ChrisRackauckas commented #636
  • Aug 05 15:24
    ChrisRackauckas closed #607
  • Aug 05 15:24
    ChrisRackauckas commented #607
  • Aug 05 15:24
    ChrisRackauckas closed #503
  • Aug 05 15:24
    ChrisRackauckas commented #503
  • Aug 05 15:23
    ChrisRackauckas closed #472
  • Aug 05 15:23
    ChrisRackauckas commented #472
  • Aug 05 15:22
    ChrisRackauckas closed #453
  • Aug 05 15:22
    ChrisRackauckas commented #453
  • Aug 05 15:22
    ChrisRackauckas closed #312
  • Aug 05 15:22
    ChrisRackauckas commented #371
  • Aug 05 15:21
    ChrisRackauckas closed #371
Paweł Biernat
@pwl
that's why it takes so long, it's a theoretical discussion for the most part
and there are only two people working on it (me and @mauro3)
Christopher Rackauckas
@ChrisRackauckas
I see. You wrote DASSL? I suppose you plan to wrap that BDF solver too?
Paweł Biernat
@pwl
that's on our list
Christopher Rackauckas
@ChrisRackauckas
How do you guys plan on handling the dependencies?
Mauro
@mauro3
Not sure yet whether we'll go the Plots.jl-eval or the MathProgBase route.
Christopher Rackauckas
@ChrisRackauckas
I think I'm going to be doing the Plots.jl approach. I think it makes more sense given how I already set things up.
So are there any more things I should know about? Any other solvers (or "integrators") you're implementing? (Waveform relaxation???)
Oh and do all of your methods work for arbitrary number types?
Paweł Biernat
@pwl
they do
the only limitation now is that the state vector has to be indexible
Also, you should know that the new solvers in ODE.jl are not conforming to the iterator interface yet. They will be ported at some point (perhaps once we merge back into the original repo).
Christopher Rackauckas
@ChrisRackauckas
Okay.
Is ODE.jl going to still be part of JuliaLang?
Mauro
@mauro3
I guess that is something to be discussed.
Paweł Biernat
@pwl
That's not up to us to decide, we are only developing a rogue fork and we have no write access to the main package.
Christopher Rackauckas
@ChrisRackauckas
Oh I see.
I would've thought you guys had write access.
Mauro
@mauro3
One day ...
Christopher Rackauckas
@ChrisRackauckas
Who actually has write access to that REPO?
Paweł Biernat
@pwl
viral for sure, I don't know about others
Christopher Rackauckas
@ChrisRackauckas
Do they know about all of the changes you're doing and will pretty much accept? Or is there a chance that when you decide to pull it will hang for awhile?
Paweł Biernat
@pwl
I talked to viral and he is pretty eager to accept almost anything that would make ODE.jl more attractive.
and we plan to be backward compatible
Christopher Rackauckas
@ChrisRackauckas
Well I think I am going to go ahead with wrapping only ODE.jl, and clean up my code to conform better with Julia standard (I realized that my camelCase isn't what's used here, so I'll fix that all up).
I will use that to set benchmarks, and when that's all together I'll contact you guys again.
Paweł Biernat
@pwl
cool, let's keep in touch!
Mauro
@mauro3
Tnx!
Christopher Rackauckas
@ChrisRackauckas
I am thinking that we shouldn't both need to write wrappers, and I would be willing to change my API in many ways to match what you guys think is reasonable.
Paweł Biernat
@pwl
that's great, we will let you know once we fix our docs, you could write some wrappers yourself if you are up for it:-)
Christopher Rackauckas
@ChrisRackauckas
So I'll get to there so you guys can take a look at it. If you think it's something that could be used for doing the rest of the wrapping (and we agree on whatever API changes), then we should merge into some org and get the rest of it together.
Just let me know if that "new" interface sol = ODE.solve(ode,stepper;opts...) ends up changing at all.
I'll be watching that dev branch.
Paweł Biernat
@pwl
@ChrisRackauckas you mentioned wave relaxation, do you have any good books/overview papers on that?
Christopher Rackauckas
@ChrisRackauckas
I know Burrage's book covers it.
I was talking with Kevin Burrage about it not too long ago. Let me find the reference.
I guess that tome is kind of expensive and out of print.
Paweł Biernat
@pwl
that is some obscure source
Christopher Rackauckas
@ChrisRackauckas
I'll probably pick up a copy myself, but I'll send Kevin an email and see if I can get a PDF, or what the "latest" reference is (I know this is from early 2000)
(Late 1990s)
Paweł Biernat
@pwl
cool, let me know if you get it
Christopher Rackauckas
@ChrisRackauckas
Well it's not obscure, you'd find it if you were reading literature and looked at citations.
But it's just out of print.
Paweł Biernat
@pwl
what I mean by it is that this doesn't seem like a popular/well known method if the only book about it is out of print
Christopher Rackauckas
@ChrisRackauckas
None of them really are.
I mean, you can find a bunch of articles on it.
Paweł Biernat
@pwl
have you tested it?
Christopher Rackauckas
@ChrisRackauckas
There are a bunch of articles which I think tend to show that Waveform Relaxation "kind of" came out on top of the other parallel methods (parallel extrapolation, special Runge-Kutta methods)