The lower limits needed an offset of +1 to avoid overwriting part of the previous results. Additionally, I have optimized the subroutine by fusing the three loops into one and computing the divisions only once.
Maybe @oscarmarino could also check if this also fixes #96 .
The lower limits needed an offset of +1 to avoid overwriting part of the previous results. Additionally, I have optimized the subroutine by fusing the three loops into one and computing the divisions only once.
Maybe @oscarmarino could also check if this also fixes #96 .