+ sgMakeRotMat4( sgLOCAL,
+ f->get_Psi() * RAD_TO_DEG,
+ f->get_Phi() * RAD_TO_DEG,
+ f->get_Theta() * RAD_TO_DEG );
+ /*
+ cout << "FG derived LOCAL matrix using sg routines" << endl;
+ MAT3mat print;
+ int i;
+ int j;
+ for ( i = 0; i < 4; i++ ) {
+ for ( j = 0; j < 4; j++ ) {
+ print[i][j] = sgLOCAL[i][j];
+ }
+ }
+ MAT3print( print, stdout);
+
+ sgMat4 sgIDENT;
+ sgMakeIdentMat4( sgIDENT );
+ for ( i = 0; i < 4; i++ ) {
+ for ( j = 0; j < 4; j++ ) {
+ print[i][j] = sgIDENT[i][j];
+ }
+ }
+ MAT3print( print, stdout);
+ */