incorporated RateSystem file

src/CriticalTransitions.jl

@@ -21,6 +21,7 @@ using StochasticDiffEq:
     u_modified!
 using DynamicalSystemsBase:
     DynamicalSystemsBase,
+    ContinuousTimeDynamicalSystem,
     CoupledSDEs,
     CoupledODEs,
     dynamic_rule,
@@ -55,6 +56,7 @@ include("extention_functions.jl")
 include("utils.jl")
 include("system_utils.jl")
 include("io.jl")
+include("RateSystem.jl")
 include("trajectories/simulation.jl")
 include("trajectories/transition.jl")
 include("trajectories/equib.jl")
@@ -69,6 +71,7 @@ using .CTLibrary
 # Core types
 export CoupledSDEs, CoupledODEs, noise_process, covariance_matrix, diffusion_matrix
 export dynamic_rule, current_state, set_state!, trajectory
+export RateSystem
 
 # Methods
 export drift, div_drift

src/RateSys.jl → src/RateSystem.jl

@@ -2,12 +2,9 @@
 # t_i    autonomous    t_ni   non-autonomous   t_nf   autonomous   t_f
 #  |                    |                        |                  |
 
-using DynamicalSystems
-using Plots
-
 # we write
 
-function RateSyst(tni,tnf,f,λ,p_λ,initvals)
+function RateSystem(tni,tnf,f,λ,p_λ,initvals)
     func(u,p,t) = combined_system(u,t,tni,tnf,f,λ,p_λ);
     return ContinuousDynamicalSystem(func, initvals, Float64[], t0=tni)
 end
@@ -21,26 +18,7 @@ end;
 ###############
 # user writes
 
-function f(u::SVector{1, Float64},p::SVector{1, Float64},t::Float64)
-    x = u[1]
-    λ = p[1]
-    dx = (x+λ)^2 - 1
-    return SVector{1}(dx)
-end;
-
-function λ(p::Vector{Float64},t::Float64)
-    r,λ_max = p
-    lambda = (λ_max/2)*(tanh(λ_max*r*t/2) +1)
-    return SVector{1}(lambda)
-end;
-
-tni=-10.; tnf=10.; p_λ = [1.0,3.0]; initvals = [-4.];
-sys = RateSyst(tni,tnf,f,λ,p_λ,initvals)
-
-traj=trajectory(sys,40.,t0=-20.)
-
-trajic=[traj[1][i][1] for i in 1:length(traj[1])]
-plot(collect(traj[2]),trajic)
+# ...moved to test/RateSystem.jl (Reyk)
 
 
 # further ideas

test/RateSystem.jl

@@ -0,0 +1,22 @@
+using Plots
+
+function f(u::SVector{1, Float64},p::SVector{1, Float64},t::Float64)
+    x = u[1]
+    λ = p[1]
+    dx = (x+λ)^2 - 1
+    return SVector{1}(dx)
+end;
+
+function λ(p::Vector{Float64},t::Float64)
+    r,λ_max = p
+    lambda = (λ_max/2)*(tanh(λ_max*r*t/2) +1)
+    return SVector{1}(lambda)
+end;
+
+tni=-10.; tnf=10.; p_λ = [1.0,3.0]; initvals = [-4.];
+sys = RateSystem(tni,tnf,f,λ,p_λ,initvals)
+
+traj=trajectory(sys,40.,t0=-20.)
+
+trajic=[traj[1][i][1] for i in 1:length(traj[1])]
+plot(collect(traj[2]),trajic)