Description: Celestial-to-intermediate matrix given CIP.
The source code for the ANSI C routine C2ixy is available: