swh:1:snp:a4c99a50dc49f82b591f268001b320f8c3ca0041
Tip revision: 1c016a1d9d1dc7690f7c256204f737661d261cb0 authored by jmc on 20 September 2007, 00:00:00 UTC
version 0.96
version 0.96
Tip revision: 1c016a1
trackSpeed.R
trackSpeed <- function(coords, time) {
dt <- diff(as.numeric(time))
d2 <- diff(as.matrix(coords))^2
speed <- d2[,1]
ndim <- ncol(d2)
if(ndim > 1)
speed <- d2 %*% rep(1,ndim)
else
speed <- d2
sqrt(as.vector(speed))/dt
}