Preconditioning Large Scale Iterative Solution of Ax = b Using a Statistical Method with Application to Matrix-Free Spectral Solution of Helmholtz Equation