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