              p1 = drv.particles(i);
                
                % Update maximum velocity
                if (this.auto_freq)
                    vparticle = vecnorm(p1.veloc_trl);
                    if (vparticle > vmax)
                        vmax = vparticle;
                    end
                    if (p1.radius > rmax)
                        rmax = p1.radius;
                    end
                end
                
                % Reset verlet lists
                p1.verlet_p = Parti