We investigate the critical properties of the d = 3 random-field Ising model with a Gaussian field distribution at zero temperature. By implementing suitable graph-theoretical algorithms, we perform a large-scale numerical simulation of the model for a vast range of values of the disorder strength h and system sizesV = L× L× L, with L 156. Using the sample-to-sample fluctuations of various quan...