Torkel Loman

Results 439 comments of Torkel Loman

Do we actually need to cache them? För the current ones it makes sense, since they depend on each others, and we do not want the user to have to...

The reason is that Catalyst internally represents its models as ModelingToolkit.jl NonlinearSystems (or other systems) (https://docs.sciml.ai/ModelingToolkit/stable/tutorials/nonlinear/). Next, ModeingToolkit gave us some quick function which essentially interprets a NonlinearSystem as a...

> Also @TorkelE, could Catalyst start providing jacobian sparsity? I am all for a 2x speed-up. This isn't something that MTK can/should determine automatically though? I.e. Catalyst jus tegenrates normal...

I mean, this entierly depends on what you are using MTK for right? For people who does DAEs with very few differential equations and many many algebraic equations, yes. But...

So there are cases where one option makes sense and another doesn't? In Sys bio the main reason for observables is pretty much the opposite (I have 3,000 differential equations,...

One could even imagine an heuristic where: - If fewer than 25 unknowns + observables, plot everything. - If fewer than 25 unknowns, but more than 1000 observables, plot only...

A potential solution here (instead of having MTK count variables/observables to figure out the best behaviour) could be to use the new system metadata interface, i.e. we could have a...

60 hours doesn't seem too outlandish for the BCR one.

It is the same with the `make_ReactionSystem_internal((Catalyst.Catalyst)`, which I think I just copied the former of. Maybe here: ```julia Catalyst.remake_ReactionSystem_internal( Catalyst.make_ReactionSystem_internal( $rxexprs, $tiv, setdiff(union($spssym, $varssym, $compssym), $obs_syms), $pssym; name =...

In a brief test, removing the `Catalys.` to: ```julia remake_ReactionSystem_internal( make_ReactionSystem_internal( $rxexprs, $tiv, setdiff(union($spssym, $varssym, $compssym), $obs_syms), $pssym; name = $name, spatial_ivs = $sivs, observed = $observed_eqs, continuous_events = $continuous_events_expr,...