Example #1
0
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);

}
Example #2
0
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);
}