Numerical study of black-hole formation in Painlevé-Gullstrand coordinates