--- /dev/null
+context("computeWerDists")
+
+test_that("computeWerDists output correct results",
+{
+ nbytes <- 8
+ endian <- "little"
+
+ # On two identical series
+ serie <- rnorm(212, sd=5)
+ series <- cbind(serie, serie)
+ getSeries <- function(indices) as.matrix(series[,indices])
+ dists <- computeWerDists(1:2, getSeries, 50, 3, 4, nbytes, endian,
+ verbose=TRUE)
+ expect_equal(dists, matrix(0.,nrow=2,ncol=2))
+
+ # On two constant series
+ # TODO: ref results. Ask Jairo to check function.
+})