On The Efficient Computation Of Multidimensional Singular Sums