Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Package not registered? #19

Open
pagnani opened this issue Jul 31, 2020 · 4 comments
Open

Package not registered? #19

pagnani opened this issue Jul 31, 2020 · 4 comments

Comments

@pagnani
Copy link

pagnani commented Jul 31, 2020

julia> using Pkg
julia> Pkg.add("Gillespie")
ERROR: The following package names could not be resolved:
 * Gillespie (not found in project, manifest or registry)

I installed directly with no problem via https. Just to let you know.

Well done and very efficient. Much faster than DifferentialEquations btw, at least in my 7 chemical species example.

I'm on

Julia Version 1.4.2
Commit 44fa15b150* (2020-05-23 18:35 UTC)
Platform Info:
  OS: Linux (x86_64-pc-linux-gnu)
  CPU: Intel(R) Xeon(R) Gold 6132 CPU @ 2.60GHz
  WORD_SIZE: 64
  LIBM: libopenlibm
  LLVM: libLLVM-8.0.1 (ORCJIT, skylake)
@ChrisRackauckas
Copy link
Contributor

Interesting, could you share the code you used to test that? Was that with SSAStepper and MassActionJump?

Indeed it isn't registered and we should probably get it re-registered.

@pagnani
Copy link
Author

pagnani commented Aug 1, 2020

The estimate was unfair indeed. I was using FunctionMap(). With SSAStepper the time is comparable. Didn't try with MassActionJump though

In case of interest, the model is here...

rab_model = @reaction_network RAB begin
    alpha5, 0 --> Rab5i
    gamma5*Rab5i, Rab5i --> 0
    kGEF5+kMAX5*(Rab5aE^n5/(kM5^n5+Rab5aE^n5)), Rab5i --> Rab5a
    kGAP5*Rab5a, Rab5a --> Rab5i
    kC5*Rab5a*E, Rab5a+E --> Rab5aE
    kD5*Rab5aE, Rab5aE --> Rab5a+E
    alpha11, 0 --> Rab11i
    gamma11*Rab11i, Rab11i --> 0
    kGEF11+kMAX11*(Rab11aE^n11/(kM11^n11+Rab11aE^n11)), Rab11i --> Rab11a
    kGAP11*Rab11a, Rab11a --> Rab11i
    kC11*Rab11a*E, Rab11a+E --> Rab11aE
    kD11*Rab11aE, Rab11aE --> Rab11a+E
    
end alpha5 gamma5 kGEF5 kGAP5 kMAX5 kM5 n5 kC5 kD5 alpha11 gamma11 kGEF11 kGAP11 kMAX11 kM11 n11 kC11 kD11

Thanks for the quick reply and the implicit suggestion.

Edit ... just for completeness
reasonable parameters (for us)
parms = (0.5,0.1,0.02,0.01,0.05,50.0,2.0,0.01,0.005,0.5,0.1,0.02,0.01,0.05,50.0,2.0,0.01,0.005)
tested initial condition
x0 = [20,0,0,50,20,0,0]

@ChrisRackauckas
Copy link
Contributor

@reaction_network will make use of MassActionJump, so that should be fine. Could you share the Gillespie code as well? I'd like to do some profiling on this and see if in general whether there's something the dense format specializes on here.

@JuliaRegistrator register()

@ChrisRackauckas
Copy link
Contributor

@sdwfrost you need to give JuliaRegistrator permissions to get it registered. Since it never bumped a v1.0 version it was dropped when everything changed registries for TOML files. But in my PRs I set it up with a Project.toml, so it should be good to register.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants