Motivated by recent research on the iterative approach proposed for the smoothed particle hydrodynamics (ISPH) method, some ideas to improve the process are introduced. The standard procedure is enhanced iterating on the residuals preserving the matrix-free nature of the process. The method is appealing providing reasonable results with disordered data distribution too and no kernel variations are needed in the approximation. This work moves forward with a novel formulation requiring a lower number of iterations to reach a desired accuracy. The computational procedure is described and some results are introduced to appreciate the proposed formulation.