Description: Initialize r-matrix to identity.
The source code for the Fortran routine IR is available: