curvesToCoeffs = function(series)
{
	#... return wavelets coeffs : compute in parallel !
	#TODO: always keep ID in first column
}

coeffsToCurves = function(coeffs)
{
	#re-expand on wavelet basis
}

getClusters = function(data, K)
{
	pam_output = pam(data, K)
	return ( list( clusts=pam_output$clustering, medoids=pam_output$medoids,
		ranks=pam_output$id.med ) )
}