chiark / gitweb /
symbolic.py: before new rotations
authorIan Jackson <ijackson@chiark.greenend.org.uk>
Fri, 17 Nov 2017 21:35:58 +0000 (21:35 +0000)
committerIan Jackson <ijackson@chiark.greenend.org.uk>
Fri, 17 Nov 2017 21:35:58 +0000 (21:35 +0000)
Signed-off-by: Ian Jackson <ijackson@chiark.greenend.org.uk>
symbolic.py

index 29ca6794e50461ef065be99b4e9571d3ddc4edb6..d4069ea0f28db7bfb20fd62573915aa64ab9d371 100755 (executable)
@@ -72,23 +72,27 @@ dbg('p2q_rotate')
 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')