It may seem that the construction of the library is a bottleneck of the method, as for the spectra with huge overlaps, peaks shifting or changing shape in time, an indication of proper areas corresponding to reagents can be infeasible. However, the proposed method is flexible and robust to such difficulties. Oftentimes, providing approximate regions where peaks from reagents are expected is sufficient for the algorithm to work properly. The approximate tolerance for peak shifting is a user-defined parameter (denoising penalty).