implements #39 with 2 levels

src/models.jl

@@ -270,3 +270,13 @@ function medianlarvasize(larvae::Vector{LarvaModel}; n::Int=100)
     end
     return median(sizes)
 end
+
+function simultaneouslarvae(larvae)
+    laststep = maximum([larva.alignedsteps[end] for larva in values(larvae)])
+    n = 0
+    for step in 1:20:laststep
+        n = max(n, count(larva -> larva.alignedsteps[1] <= step <= larva.alignedsteps[end],
+                         values(larvae)))
+    end
+    return n
+end

src/viewer.jl

@@ -67,7 +67,14 @@ function larvaviewer(controller;
 
     controller[:player] = player = timecontroller(times)
     controller[:larva] = larva = LarvaController(controller, model, tag_lut, player)
-    delayed_controller = slave(larva, PeriodicObservation(1.0))
+
+    n_simultaneous_larvae = simultaneouslarvae(model)
+    if 40 < n_simultaneous_larvae
+        @info "Refresh rate throttled to 1Hz" n_simultaneous_larvae
+        delayed_controller = slave(larva, PeriodicObservation(1.0))
+    else
+        delayed_controller = larva
+    end
 
     on(getactivelarva(controller)) do id
         @info isnothing(id) ? "No active larva" : "Activating larva #$(id)"

src/wgl.jl

@@ -129,10 +129,15 @@ function AssayPlot(ctrl, larvae::DecoratedLarvae; size=(1000, 750))
     AssayPlot(ctrl, fig, ax, larvae, button)
 end
 
-function AssayPlot(ctrl, larvae::Vector{LarvaModel}; kwargs...)
-    sample = downsampler()
+function AssayPlot(ctrl, larvae::Vector{LarvaModel}; downsampling_threshold=1e5, kwargs...)
+    npoints = sum([length(larva.path) for larva in larvae])
+    if downsampling_threshold <= npoints
+        @info "Trajectory down-sampling active" npoints
+        sample = downsampler()
+        larvae = sample.(larvae)
+    end
     plot = AssayPlot(ctrl,
-                     DecoratedLarvae(sample.(larvae),
+                     DecoratedLarvae(larvae,
                                      gettags(ctrl),
                                      gettimestep(ctrl));
                      kwargs...)