Description: Celestial-to-intermediate matrix given b-p-n.
The source code for the Fortran routine C2IBPN is available: