|
157 | 157 | return hopping rate of species at site |
158 | 158 | """ |
159 | 159 | function evalhoprate(hop_rates::HopRatesGraphDs, u, species, site, spatial_system) |
160 | | - @inbounds u[species, site] * hop_rates.hopping_constants[species] * |
161 | | - outdegree(spatial_system, site) |
| 160 | + @inbounds u[species, site] * hop_rates.hopping_constants[species] * outdegree(spatial_system, site) |
162 | 161 | end |
163 | 162 |
|
164 | 163 | ############## hopping rates of form D_{s,i} ################ |
|
197 | 196 | return hopping rate of species at site |
198 | 197 | """ |
199 | 198 | function evalhoprate(hop_rates::HopRatesGraphDsi, u, species, site, spatial_system) |
200 | | - @inbounds u[species, site] * hop_rates.hopping_constants[species, site] * |
201 | | - outdegree(spatial_system, site) |
| 199 | + @inbounds u[species, site] * hop_rates.hopping_constants[species, site] * outdegree(spatial_system, site) |
202 | 200 | end |
203 | 201 |
|
204 | 202 | ############## hopping rates of form D_{s,i,j} ################ |
@@ -344,8 +342,7 @@ function sample_target_site(hop_rates::HopRatesGraphDsLij, site, species, rng, |
344 | 342 | end |
345 | 343 |
|
346 | 344 | function evalhoprate(hop_rates::HopRatesGraphDsLij, u, species, site, spatial_system) |
347 | | - @inbounds u[species, site] * hop_rates.species_hop_constants[species] * |
348 | | - hop_rates.hop_const_cumulative_sums[site][end] |
| 345 | + @inbounds u[species, site] * hop_rates.species_hop_constants[species] * hop_rates.hop_const_cumulative_sums[site][end] |
349 | 346 | end |
350 | 347 |
|
351 | 348 | ############## hopping rates of form D_s * L_{i,j} optimized for cartesian grid ################ |
@@ -398,8 +395,7 @@ function sample_target_site(hop_rates::HopRatesGridDsLij, site, species, rng, gr |
398 | 395 | end |
399 | 396 |
|
400 | 397 | function evalhoprate(hop_rates::HopRatesGridDsLij, u, species, site, grid) |
401 | | - @inbounds u[species, site] * hop_rates.species_hop_constants[species] * |
402 | | - hop_rates.hop_const_cumulative_sums[end, site] |
| 398 | + @inbounds u[species, site] * hop_rates.species_hop_constants[species] * hop_rates.hop_const_cumulative_sums[end, site] |
403 | 399 | end |
404 | 400 |
|
405 | 401 | ############## hopping rates of form D_{s,i} * L_{i,j} ################ |
@@ -441,8 +437,7 @@ function sample_target_site(hop_rates::HopRatesGraphDsiLij, site, species, rng, |
441 | 437 | end |
442 | 438 |
|
443 | 439 | function evalhoprate(hop_rates::HopRatesGraphDsiLij, u, species, site, spatial_system) |
444 | | - @inbounds u[species, site] * hop_rates.species_hop_constants[species, site] * |
445 | | - hop_rates.hop_const_cumulative_sums[site][end] |
| 440 | + @inbounds u[species, site] * hop_rates.species_hop_constants[species, site] * hop_rates.hop_const_cumulative_sums[site][end] |
446 | 441 | end |
447 | 442 |
|
448 | 443 | ############## hopping rates of form D_{s,i} * L_{i,j} optimized for cartesian grid ################ |
@@ -497,6 +492,5 @@ function sample_target_site(hop_rates::HopRatesGridDsiLij, site, species, rng, g |
497 | 492 | end |
498 | 493 |
|
499 | 494 | function evalhoprate(hop_rates::HopRatesGridDsiLij, u, species, site, grid) |
500 | | - @inbounds u[species, site] * hop_rates.species_hop_constants[species, site] * |
501 | | - hop_rates.hop_const_cumulative_sums[end, site] |
| 495 | + @inbounds u[species, site] * hop_rates.species_hop_constants[species, site] * hop_rates.hop_const_cumulative_sums[end, site] |
502 | 496 | end |
0 commit comments