by

Where communities thrive


  • Join over 1.5M+ people
  • Join over 100K+ communities
  • Free without limits
  • Create your own community
People
Activity
  • 21:47
    ChrisRackauckas commented #1193
  • 21:47
    ChrisRackauckas edited #1193
  • 21:46
    ChrisRackauckas reopened #1193
  • 21:44

    kanav99 on gh-pages

    build based on 4078de6 (compare)

  • 21:40
    YingboMa edited #508
  • 21:40
    YingboMa edited #508
  • 21:39
    YingboMa synchronize #508
  • 21:39

    YingboMa on skipzerofil

    Interpolate fill_array_with_zer… fill_zero -> fillzeros (compare)

  • 21:39

    kanav99 on gh-pages

    build based on ef3f5d4 (compare)

  • 21:34
    shashi opened #509
  • 21:33

    shashi on fix-consting

    fix when op is a Variable or Ex… (compare)

  • 21:33

    kanav99 on gh-pages

    build based on 8a0d245 (compare)

  • 21:32

    shashi on fix-505

    fix when op is a Variable or Ex… (compare)

  • 21:32

    YingboMa on skipzeros

    (compare)

  • 21:32
    YingboMa closed #504
  • 21:32
    YingboMa commented #504
  • 21:31
    YingboMa review_requested #508
  • 21:28
    YingboMa opened #508
  • 21:27

    shashi on skipzerofil

    Add fill_zero kwarg to build_fu… (compare)

  • 21:24

    kanav99 on gh-pages

    build based on e65ed6b (compare)

BridgingBot
@GitterIRCbot
[slack] <chrisrackauckas> interesting
[slack] <chrisrackauckas> that will need an issue
BridgingBot
@GitterIRCbot
[slack] <Wei> Hi, for the error estimates err_scaled = err/(abstol+max(uprev,u)reltol), what is the value *err here? Algorithm depended or a default value for all different algorithm?
[slack] <chrisrackauckas> algorithm dependent
[slack] <chrisrackauckas> it's something computed as part of the ODE solver, so it's not a parameter
[slack] <chrisrackauckas> it's a local error estimate
BridgingBot
@GitterIRCbot
[slack] <Wei> Thanks. For a SDEproblem, is it possible that if I set a larger g , the solver need more iterations and smaller dt to meet the abstol and reltol requirement? It means for the same abstol and reltol , the larger of g , the increasing of iterations and smaller adaptive dt ?
[slack] <chrisrackauckas> yes
[slack] <chrisrackauckas> that's quite common
[slack] <Wei> Ok, I got it. I don't know before. But I realized when I increasing the g value and keeping the other same, the solver needs more time to return a result.
BridgingBot
@GitterIRCbot
[slack] <chrisrackauckas> it's a stability thing
[slack] <Wei> Ok, I got it. Thanks.
BridgingBot
@GitterIRCbot
[slack] <simeondavidschaub99> Is ModelingToolkit.expand_derivatives supposed to recurse into a sum of two non-expanded derivatives?
[slack] <chrisrackauckas> it should

[slack] <simeondavidschaub99> The following example doesn't:
```using ModelingToolkit
@variables t x(t)
@derivatives ∂ₜ'~t
L = .5 ∂ₜ(x)^2 - .5 x^2
euler_lagrange(L, x) = expand_derivatives(Differential(x)(L) - ∂ₜ(Differential(∂ₜ(x))(L)))

julia> euler_lagrange(L, x)
derivative(0.5 derivative(x(t), t) ^ 2 + -0.5 x(t) ^ 2, x(t)) + -1 derivative(derivative(0.5 derivative(x(t), t) ^ 2 + -0.5 * x(t) ^ 2, derivative(x(t), t)), t)```

[slack] <chrisrackauckas> is that latest master?
BridgingBot
@GitterIRCbot
[slack] <chrisrackauckas> @shashi very recently put a fix in there.
[slack] <simeondavidschaub99> Yes:
``(@v1.5) pkg> st ModelingToolkit Status~/.julia/environments/v1.5/Project.toml[961ee093] ModelingToolkit v3.13.0https://github.com/SciML/ModelingToolkit.jl.git#master````
[slack] <simeondavidschaub99> Just checked out the latest master
BridgingBot
@GitterIRCbot
[slack] <shashi> @yingbo_ma take a look
BridgingBot
@GitterIRCbot
[slack] <Adam Gerlach> @AlCap23 In case you haven't seen this "Forecasting Sequential Data Using Consistent Koopman Autoencoders" may interest you. Includes Python code https://github.com/erichson/koopmanAE
BridgingBot
@GitterIRCbot
[slack] <isaacsas> Should hashing of Operations now be working? (i.e. are they safe to use as keys in a Dict?)
[slack] <chrisrackauckas> I think @shashi fixed that
[slack] <chrisrackauckas> hashing operations is good
[slack] <chrisrackauckas> hashing variables is weird IIRC
[slack] <isaacsas> That's the opposite of my experience.
[slack] <shashi> yeah they should be safe in dicts
[slack] <chrisrackauckas> we have a hash defined for operations and it's the one that's tested?
[slack] <isaacsas> OK, interesting. I think all the MT code I've used with Dicts used Variables (because previously I had issues with Operations). (I'm still trying to figure out this issetequal not working issue, for which it was suggested that maybe the problem is not having hash defined, but I guess it is something else.)
[slack] <isaacsas> Good to know that should be safe to do going forward at least.
BridgingBot
@GitterIRCbot
[slack] <isaacsas> @chrisrackauckas You had mentioned separating out the DEBio docs; is that something it makes sense to do now?
[slack] <chrisrackauckas> I think so
[slack] <chrisrackauckas> the DEBio docs probably are worth a few pages, at least for tutorials and such, right?
BridgingBot
@GitterIRCbot
[slack] <isaacsas> Probably. I was more thinking, if it is possible to do this, we should pull the relevant ReactionSystem docs from MT into them too. I wasn't sure if that is something you want to pull into the DiffEqDocs in any way.
BridgingBot
@GitterIRCbot
[slack] <chrisrackauckas> maybe, depends on what's easiest to maintain and what's easiest for the users
BridgingBot
@GitterIRCbot
[slack] <isaacsas> OK, let's leave in the DiffEqDocs for now and I'll just pull in the few MT defs we need there.
BridgingBot
@GitterIRCbot
[slack] <chrisrackauckas> I think it may even make sense to move ReactionSystem to DEBio / Catalyst (let's do the breaking release with the name change BTW)
[slack] <chrisrackauckas> JumpSystem and all of that is MTK, but I think ReactionSystem and the DSL for ReactionSystem might make sense as what Cataylst is?
[slack] <chrisrackauckas> I think that would make more sense then double docs for ReactionSystem
BridgingBot
@GitterIRCbot
[slack] <isaacsas> Yeah, that that would be fine with me, and seems a clean separation.
[slack] <isaacsas> It is a bit weird to have some stuff in MT and some in DEBio
BridgingBot
@GitterIRCbot
[slack] <chrisrackauckas> yeah, I wasn't sure how it would all pan out
[slack] <chrisrackauckas> I thought the split was, abstract systems in MTK
[slack] <chrisrackauckas> I think instead it's, abstract systems for symbolic numerics in MTK, and abstract systems for domains in a domain modeling package
[slack] <chrisrackauckas> so ODESystem MTK, JumpSystem MTK, ReactionSystem is it's own whole world
BridgingBot
@GitterIRCbot
[slack] <simeondavidschaub99> Should I open an issue on GitHub about this?
BridgingBot
@GitterIRCbot
[slack] <chrisrackauckas> is the new master good? I thought @shashi just handled this
BridgingBot
@GitterIRCbot
[slack] <simeondavidschaub99> I just updated to the latest master and I still get this.
BridgingBot
@GitterIRCbot
[slack] <shashi> Yeah do open an issue
BridgingBot
@GitterIRCbot
[slack] <simeondavidschaub99> SciML/ModelingToolkit.jl#505