Euler-Maclaurin based summation for 2D problems