+#' \dontrun{
+#' # WER distances computations are too long for CRAN (for now)
+#'
+#' # Random series around cos(x,2x,3x)/sin(x,2x,3x)
+#' x = seq(0,500,0.05)
+#' L = length(x) #10001
+#' ref_series = matrix( c(cos(x),cos(2*x),cos(3*x),sin(x),sin(2*x),sin(3*x)), ncol=6 )
+#' library(wmtsa)
+#' series = do.call( cbind, lapply( 1:6, function(i)
+#' do.call(cbind, wmtsa::wavBootstrap(ref_series[,i], n.realization=400)) ) )
+#' #dim(series) #c(2400,10001)
+#' res_ascii = claws(series, K1=60, K2=6, 200, verbose=TRUE)
+#'
+#' # Same example, from CSV file
+#' csv_file = "/tmp/epclust_series.csv"
+#' write.table(series, csv_file, sep=",", row.names=FALSE, col.names=FALSE)
+#' res_csv = claws(csv_file, K1=60, K2=6, 200)
+#'
+#' # Same example, from binary file
+#' bin_file <- "/tmp/epclust_series.bin"
+#' nbytes <- 8
+#' endian <- "little"
+#' binarize(csv_file, bin_file, 500, nbytes, endian)
+#' getSeries <- function(indices) getDataInFile(indices, bin_file, nbytes, endian)
+#' res_bin <- claws(getSeries, K1=60, K2=6, 200)
+#' unlink(csv_file)
+#' unlink(bin_file)
+#'
+#' # Same example, from SQLite database
+#' library(DBI)
+#' series_db <- dbConnect(RSQLite::SQLite(), "file::memory:")
+#' # Prepare data.frame in DB-format
+#' n <- nrow(series)
+#' time_values <- data.frame(
+#' id = rep(1:n,each=L),
+#' time = rep( as.POSIXct(1800*(0:n),"GMT",origin="2001-01-01"), L ),
+#' value = as.double(t(series)) )
+#' dbWriteTable(series_db, "times_values", times_values)
+#' # Fill associative array, map index to identifier
+#' indexToID_inDB <- as.character(
+#' dbGetQuery(series_db, 'SELECT DISTINCT id FROM time_values')[,"id"] )
+#' serie_length <- as.integer( dbGetQuery(series_db,
+#' paste("SELECT COUNT * FROM time_values WHERE id == ",indexToID_inDB[1],sep="")) )
+#' getSeries <- function(indices) {
+#' request <- "SELECT id,value FROM times_values WHERE id in ("
+#' for (i in indices)
+#' request <- paste(request, indexToID_inDB[i], ",", sep="")
+#' request <- paste(request, ")", sep="")
+#' df_series <- dbGetQuery(series_db, request)
+#' if (length(df_series) >= 1)
+#' as.matrix(df_series[,"value"], nrow=serie_length)
+#' else
+#' NULL
+#' }
+#' res_db = claws(getSeries, K1=60, K2=6, 200))
+#' dbDisconnect(series_db)
+#'
+#' # All results should be the same:
+#' library(digest)
+#' digest::sha1(res_ascii)
+#' digest::sha1(res_csv)
+#' digest::sha1(res_bin)
+#' digest::sha1(res_db)