Macroscopic B0 magnetic field inhomogeneity adversely affects Gradient-Recalled-Echo (GRE) MRI images. A previously introduced voxel spread function (VSF) method has been proposed to correct for these adverse effects if data are acquired using multi-GRE (mGRE) sequences. The method accounts for both, through-slice and in-slice field inhomogeneities and accounts for between-voxel signal propagation effects. The direct voxel-wise VSF calculation requires long computation time. We propose a library-driven approach for the VSF calculation that significantly reduces computation time, thus providing a platform for the applications of mGRE-based quantitative techniques in clinical practices.