X-Git-Url: https://git.auder.net/css/rpsls.css?a=blobdiff_plain;f=epclust%2FR%2FcomputeCoeffs.R;fp=epclust%2FR%2FcomputeCoeffs.R;h=0000000000000000000000000000000000000000;hb=0e2dce80a3fddaca50c96c6c27a8b32468095d6c;hp=fca3b9188d3ceda45281070211e98adfc1adfb00;hpb=48108c3999d28d973443fa5e78f73a0a9f2bfc07;p=epclust.git

diff --git a/epclust/R/computeCoeffs.R b/epclust/R/computeCoeffs.R
deleted file mode 100644
index fca3b91..0000000
--- a/epclust/R/computeCoeffs.R
+++ /dev/null
@@ -1,43 +0,0 @@
-computeCoeffs = function(data, index, nb_series_per_chunk, wf)
-{
-	coeffs_chunk = NULL
-	if (is.data.frame(data) && index < nrow(data))
-	{
-		#full data matrix
-		coeffs_chunk = curvesToCoeffs(
-			data[index:(min(index+nb_series_per_chunk-1,nrow(data))),], wf)
-	}
-	else if (is.function(data))
-	{
-		#custom user function to retrieve next n curves, probably to read from DB
-		coeffs_chunk = curvesToCoeffs( data(rank=(index-1)+seq_len(nb_series_per_chunk)), wf )
-	}
-	else if (exists(data_con))
-	{
-		#incremental connection ; TODO: more efficient way to parse than using a temp file
-		ascii_lines = readLines(data_con, nb_series_per_chunk)
-		if (length(ascii_lines > 0))
-		{
-			series_chunk_file = ".series_chunk"
-			writeLines(ascii_lines, series_chunk_file)
-			coeffs_chunk = curvesToCoeffs( read.csv(series_chunk_file), wf )
-			unlink(series_chunk_file)
-		}
-	}
-	coeffs_chunk
-}
-
-curvesToCoeffs = function(series, wf)
-{
-	if (!require(wavelets, quietly=TRUE))
-		stop("Couldn't load wavelets library")
-	L = length(series[1,])
-	D = ceiling( log2(L) )
-	nb_sample_points = 2^D
-	#TODO: parallel::parApply() ?!
-	as.data.frame( apply(series, 1, function(x) {
-		interpolated_curve = spline(1:L, x, n=nb_sample_points)$y
-		W = wavelets::dwt(interpolated_curve, filter=wf, D)@W
-		rev( sapply( W, function(v) ( sqrt( sum(v^2) ) ) ) )
-	}) )
-}