We provided a numerical solution to the inverse inference problem for close-packed structures by using a simulated annealing procedure which we call numerical stacking reconstruction (NSR). The method does not rely on any extra assumption besides the PCPs themselves. In this sense, it yields results as close as the numerical procedure is capable of. A key aspect of the developed procedure is the educated guess of the initial sequences used as input to the simulated annealing stage. The convergence and reliability of the result depended heavily on it. NSR incorporates the knowledge of previous models into its framework to advance initial guesses. This has the important implication that even in the worst scenario, NSR performs as well as previous approaches. This also implies that it is amenable to further extension with other models as previous guesses. Yet in all the numerical test performed, there was no need for additional models than the ones considered. In all cases, the convergence of the NSR the procedure was below seven minutes for a general purpose desktop computer with 8M of RAM and a Dual core Pentium processor.