def cross_p(u, v):
    """ returns the cross product of two vectors u x v - orthogonal to u and v """
    return (u[1]*v[2] - u[2]*v[1],
            u[2]*v[0] - u[0]*v[2],
            u[0]*v[1] - u[1]*v[0])