An importance sampling iterative algorithm for diagonalizing large matrices is upgraded and adopted for large scale nuclear shell model calculations using a spin uncoupled basis. Its numerical implementation shows that the iterative procedure converges rapidly to the exact eigensolutions achieving an effective drastic cut of the sizes of the Hamiltonian matrix.