Description: Celestial-to-intermediate matrix given B-P-N.
The source code for the ANSI C routine C2ibpn is available: