An incremental Gaussian-process sampler and adaptive ray marcher accelerate ray intersections with Gaussian process implicit surfaces. Step sizes are chosen to control the probability of missing a surface crossing, reducing rendering error by up to 46 times at equal time compared with the previous method.