rthsort <- function(x) {
   require(Rcpp)
   dyn.load("rthsort.so")
   .Call("rthsort",x)
}