Transposing coord (i.e. let the column index correspond to x,y,z and the row index to the different particles) helps a little bit (but not so much). Julia uses column major order, Python is different I believe. How big is the difference with Numba?
Op vrijdag 29 januari 2016 21:26:07 UTC+1 schreef STAR0SS: > > A significant part of the time is spent in computing the norms (r_sq = > rx*rx + ry*ry + rz*rz) and the absolute value (abs(rz) > boxl2), I don't > know if there's more efficient ways to compute those. >
