Generalized relativistic pseudopotentials (GRPP) of atomic cores implying the use of different potentials for atomic electronic shells with different principal quantum numbers give rise to accurate and reliable relativistic electronic structure models of atoms, molecules, clusters, and solids. These models readily incorporate the effects of Breit electron-electron interactions and one-loop quantum electrodynamics effects. Here we report the computational procedure for evaluating one-electron integrals of GRPP over contracted Gaussian functions. This procedure was implemented in a library of routines named LIBGRPP, which can be integrated into existing quantum chemistry software, thus enabling the application of various methods to solve the many-electron problem with GRPPs. Pilot applications to electronic transitions in the ThO and UO2 molecules using the new library and intermediate-Hamiltonian Fock space relativistic coupled cluster method are presented. The results clearly demonstrate that rather economical tiny-core GRPP models can exceed in accuracy relativistic all-electron models defined by Dirac-Coulomb and Dirac-Coulomb-Gaunt Hamiltonians.