with parallel::export

[epclust.git] / epclust / R / utils.R

toInteger <- function(x, condition)
{
	if (!is.integer(x))
		tryCatch(
			{x = as.integer(x)[1]},
			error = function(e) paste("cannot convert argument",substitute(x),"to integer")
		)
	if (!condition(x))
		stop(paste("argument",substitute(x),"does not verify condition",body(condition)))
	x
}

writeCoeffs = function(coeffs)
{
	file = ".coeffs"
	#.........
	#C function (from data.frame, type of IDs ??! force integers ? [yes])
	#return raw vector
	#take raw vector, append it (binary mode) to a file
#TODO: appendCoeffs() en C --> serialize et append to file
}

readCoeffs = function(indices)
{
	#......
	file = ".coeffs"
	#C function (from file name)
}

getSeries(data, rank=NULL, id=NULL)
{
	#TODO:
}

curvesToCoeffs = function(series, wf)
{
	L = length(series[1,])
	D = ceiling( log2(L) )
	nb_sample_points = 2^D
	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) ) ) ) )
	})
}