Euler-Maclaurin based summation for 1D problems