New version of mafast based on data.table (5s for 450000 obs aggr=3)