decimate_tsd() now handles multidimensional data

fixes #15