In this paper, three different uniqueness data are investigated to reconstruct the potential function in the Sturm-Liouville boundary value problem in the normal form. Taking account of R\"{o}hrl's objective function, the steepest descent method is used in the computation of potential functions. To decrease the volume of computation, we propose a theorem to precalculate the minimization parameter that is required in the optimization. Further, we propose a novel time-saving algorithm in which the obligation of using the asymptotics of eigenvalues and eigenfunctions and the appropriateness of selected boundary conditions are also eliminated. As partial data, we take two spectra, the set of the $j$th elements of the infinite numbers of spectra obtained by changing boundary conditions in the problem, and one spectrum with the set of terminal velocities. In order to show the efficiency of the proposed method, numerical results are given for three test potentials which are smooth, nonsmooth continuous, and noncontinuous, respectively.