v_cross = vect1.cross(vect2)
v_dot = vect1.dot(vect2)
if v_dot == 0:
v_angle = pi/2.
else:
v_angle =acos( v_dot / ( vect2.magnitude()*vect_base.magnitude() ) )R = N.orient_new_axis('R', v_angle,v_cross)
--
You received this message because you are subscribed to the Google Groups "sympy" group.
To unsubscribe from this group and stop receiving emails from it, send an email to sympy+un...@googlegroups.com.
To post to this group, send email to sy...@googlegroups.com.
Visit this group at https://groups.google.com/group/sympy.
To view this discussion on the web visit https://groups.google.com/d/msgid/sympy/082f4048-691d-4b29-83cb-97c3e03300b1%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.