From: Ian Jackson Date: Fri, 17 Nov 2017 21:35:58 +0000 (+0000) Subject: symbolic.py: before new rotations X-Git-Url: https://www.chiark.greenend.org.uk/ucgi/~ian/git?a=commitdiff_plain;h=5c1ba2487d9dcf377a95bff694612a8ff8c0afdf;p=moebius3.git symbolic.py: before new rotations Signed-off-by: Ian Jackson --- diff --git a/symbolic.py b/symbolic.py index 29ca679..d4069ea 100755 --- a/symbolic.py +++ b/symbolic.py @@ -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')