Should v+=s go after draw(u,round(v))?

This algorithm is a lot better than having to check if each square is part of the line or not and also obeys the 'diamond' drawing principle if you are sampling at the midpoint of each square.