Restrict main MM routines to distance threshold

[urbananalyst]

Restrict main MM routines to distance threshold

[mpadge]

library (m4ra)
packageVersion ("m4ra")
#> [1] '0.1.0.193'

net_sc <- "/<path>/<to>/mannheim-sc.Rds"
gtfs <- "/<path>/<to>/mannheim-gtfs-with-transfers.Rds"
city <- city_name <- "mannheim"
from <- c ("1216123389", "4187247393", "1178997830", "4502052112", "659273998",
           "1386088239", "273373473", "269713455", "1742645215", "5045799953")
from_initial <- from
initial_mode <- final_mode <- "bicycle"
quiet <- FALSE
day <- start_time_limits <- NULL
fast <- FALSE

# Initial call to set up caching
x <- m4ra_times_multi_mode (net_sc = net_sc, gtfs = gtfs,
        city_name = city_name, initial_mode = initial_mode,
        final_mode = final_mode, from = from, quiet = TRUE)

system.time ({
    x1 <- m4ra_times_multi_mode (
        net_sc = net_sc,
        gtfs = gtfs,
        city_name = city_name,
        initial_mode = initial_mode,
        final_mode = final_mode,
        from = from
    )
})
#> ✔ Total time for data preparation: 00:00:00.37
#> ℹ Calculating single-mode travel times for bicycle
#> ✔ Calculated single-mode travel times for bicycle
#> ℹ Calculating multi-mode travel times
#> ✔ Calculated multi-mode travel times
#>    user  system elapsed 
#>  45.724   0.807   8.221
dim (x1)
#> [1]    10 62403
system.time ({
    x2 <- m4ra_times_multi_mode (
        net_sc = net_sc,
        gtfs = gtfs,
        city_name = city_name,
        initial_mode = initial_mode,
        final_mode = final_mode,
        from = from,
        duration_max = 900 # 15 minutes
    )
})
#> ✔ Total time for data preparation: 00:00:00.28
#> ℹ Calculating single-mode travel times for bicycle
#> ✔ Calculated single-mode travel times for bicycle
#> ℹ Calculating multi-mode travel times
#> ✔ Calculated multi-mode travel times
#>    user  system elapsed 
#>  31.249   0.473   5.646
dim (x2)
#> [1]    10 62354

^{Created on 2023-04-24 with reprex v2.0.2}

And passing a value for duration_max trims a few destination points off the final matrix which are reached with the full routine, but that example amounts to < 0.08% of all points, so can easily be ignored for general use.