Compute random projection statistics (parallelized over projections)

Usage

rp_stat(proj_x_ord, residuals, n_proj)