drop fastmaths

Author: isaacsas

src/aggregators/direct.jl

@@ -58,16 +58,16 @@ end
 # calculate the next jump / jump time
 function generate_jumps!(p::DirectJumpAggregation, integrator, u, params, t)
     p.sum_rate, ttnj = time_to_next_jump(p, u, params, t)
-    @fastmath p.next_jump_time = t + ttnj
+    p.next_jump_time = t + ttnj
     @inbounds p.next_jump = searchsortedfirst(p.cur_rates, rand(p.rng) * p.sum_rate)
     nothing
 end
 
 ######################## SSA specific helper routines ########################
 
 # tuple-based constant jumps
-@fastmath function time_to_next_jump(p::DirectJumpAggregation{T, S, F1, F2, RNG}, u, params,
-                                     t) where {T, S, F1 <: Tuple, F2 <: Tuple, RNG}
+function time_to_next_jump(p::DirectJumpAggregation{T, S, F1, F2, RNG}, u, params,
+                           t) where {T, S, F1 <: Tuple, F2 <: Tuple, RNG}
     prev_rate = zero(t)
     new_rate = zero(t)
     cur_rates = p.cur_rates
@@ -108,9 +108,8 @@ end
 end
 
 # function wrapper-based constant jumps
-@fastmath function time_to_next_jump(p::DirectJumpAggregation{T, S, F1, F2, RNG}, u, params,
-                                     t) where {T, S, F1 <: AbstractArray,
-                                               F2 <: AbstractArray, RNG}
+function time_to_next_jump(p::DirectJumpAggregation{T, S, F1, F2, RNG}, u, params,
+                           t) where {T, S, F1 <: AbstractArray, F2 <: AbstractArray, RNG}
     prev_rate = zero(t)
     new_rate = zero(t)
     cur_rates = p.cur_rates

src/aggregators/frm.jl

@@ -75,7 +75,7 @@ end
 ######################## SSA specific helper routines ########################
 
 # mass action jumps
-@fastmath function next_ma_jump(p::FRMJumpAggregation, u, params, t)
+function next_ma_jump(p::FRMJumpAggregation, u, params, t)
     ttnj = typemax(typeof(t))
     nextrx = zero(Int)
     majumps = p.ma_jumps
@@ -91,9 +91,8 @@ end
 end
 
 # tuple-based constant jumps
-@fastmath function next_constant_rate_jump(p::FRMJumpAggregation{T, S, F1, F2, RNG}, u,
-                                           params,
-                                           t) where {T, S, F1 <: Tuple, F2 <: Tuple, RNG}
+function next_constant_rate_jump(p::FRMJumpAggregation{T, S, F1, F2, RNG}, u, params,
+                                 t) where {T, S, F1 <: Tuple, F2 <: Tuple, RNG}
     ttnj = typemax(typeof(t))
     nextrx = zero(Int)
     if !isempty(p.rates)
@@ -111,10 +110,9 @@ end
 end
 
 # function wrapper-based constant jumps
-@fastmath function next_constant_rate_jump(p::FRMJumpAggregation{T, S, F1, F2, RNG}, u,
-                                           params,
-                                           t) where {T, S, F1 <: AbstractArray,
-                                                     F2 <: AbstractArray, RNG}
+function next_constant_rate_jump(p::FRMJumpAggregation{T, S, F1, F2, RNG}, u, params,
+                                 t) where {T, S, F1 <: AbstractArray, F2 <: AbstractArray,
+                                           RNG}
     ttnj = typemax(typeof(t))
     nextrx = zero(Int)
     if !isempty(p.rates)

src/massaction_rates.jl

@@ -3,8 +3,8 @@
 # stochiometric coefficient.
 ###############################################################################
 
-@inline @fastmath function evalrxrate(speciesvec::AbstractVector{T}, rxidx::S,
-                                      majump::MassActionJump{U, V, W, X})::R where
+@inline function evalrxrate(speciesvec::AbstractVector{T}, rxidx::S,
+                            majump::MassActionJump{U, V, W, X})::R where
     {T, S, R, U <: AbstractVector{R}, V, W, X}
     val = one(T)
     @inbounds for specstoch in majump.reactant_stoch[rxidx]
@@ -19,17 +19,17 @@
     @inbounds return val * majump.scaled_rates[rxidx]
 end
 
-@inline @fastmath function executerx!(speciesvec::AbstractVector{T}, rxidx::S,
-                                      majump::M) where {T, S, M <: AbstractMassActionJump}
+@inline function executerx!(speciesvec::AbstractVector{T}, rxidx::S,
+                            majump::M) where {T, S, M <: AbstractMassActionJump}
     @inbounds net_stoch = majump.net_stoch[rxidx]
     @inbounds for specstoch in net_stoch
         speciesvec[specstoch[1]] += specstoch[2]
     end
     nothing
 end
 
-@inline @fastmath function executerx(speciesvec::SVector{T}, rxidx::S,
-                                     majump::M) where {T, S, M <: AbstractMassActionJump}
+@inline function executerx(speciesvec::SVector{T}, rxidx::S,
+                           majump::M) where {T, S, M <: AbstractMassActionJump}
     @inbounds net_stoch = majump.net_stoch[rxidx]
     @inbounds for specstoch in net_stoch
         speciesvec = setindex(speciesvec, speciesvec[specstoch[1]] + specstoch[2],