A closest point method library for PDEs on surfaces with parallel domain decomposition solvers and preconditioners