q_dirn_maincoords = p2q_rotate * q_dirn_owncoords;
q_maincoords = p2q_rotate * q_owncoords + p2q_translate
-dbg('q_maincoords','q_dirn_maincoords')
-dbg('q_maincoords.replace(t,0)','q_dirn_maincoords.replace(t,0)')
-
dbg('diff(p_dirn_rightvars,s)')
dbg('diff(q_dirn_maincoords,t)')
dbg('diff(q_dirn_maincoords,t).replace(t,0)')
assert(Eq(p2q_rotate * Matrix([0,1,mu]), p_dirn_rightvars))
-print('\n eye3 subs etc.\n')
-dbg('''Eq(eye(3) * Matrix([1,0,mu]),
- p_dirn_rightvars .cross(Matrix([0,0,1]) .subs(s,0)))''')
+#for v in 's','t','la','mu':
+# dbg('diff(q_maincoords,%s)' % v)
+
+#print('\n eye3 subs etc.\n')
+#dbg('''Eq(eye(3) * Matrix([1,0,mu]),
+# p_dirn_rightvars .cross(Matrix([0,0,1]) .subs(s,0)))''')
-dbg('''Eq(p2q_rotate * Matrix([1,0,mu]),
- p_dirn_rightvars .cross(Matrix([0,0,1])))''')
+#dbg('''Eq(p2q_rotate * Matrix([1,0,mu]),
+# p_dirn_rightvars .cross(Matrix([0,0,1])))''')
#eq = Eq(qmat * q_dirn_owncoords_0, p_dirn_rightvars)
#print
#pprint(eq)
#solve(eq, Q)
+
+dbg('q_maincoords.replace(t,0)','q_dirn_maincoords.replace(t,0)')
+
+dbg('q_maincoords','q_dirn_maincoords')