getCoeffs = function(series)
{
	#... return wavelets coeffs : compute in parallel !
}

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