M6 m6Qy(double y, double ydot) { M6 m6; m6SetPP(m6, m3Ry(y)); m6SetPV(m6, m3O()); m6SetVP(m6, m3RyDot(y, ydot)); m6SetVV(m6, m3Ry(y)); return(m6); }
V6 equ2gal(V6 v6) { /* subtract e-terms */ v6 = ellab(B1950, v6, -1); v6 = m3v6(m3Rz(d2r(GAL_RA)), v6); v6 = m3v6(m3Ry(d2r(90-GAL_DEC)), v6); v6 = m3v6(m3Rz(d2r(90-GAL_LON)), v6); return(v6); }