A multigrid continuation algorithm for multiple peak solutions of the Gross-Pitaevskii equation in a periodic potential