- Timestamp:
- Oct 26, 2017 12:03:38 PM (7 years ago)
- Branches:
- master, core_shell_microgels, magnetic_model, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests
- Children:
- b9578fc
- Parents:
- 3bfd924
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
explore/angles.py
r20fe0cd rf39d4a3 161 161 162 162 # various vectors 163 xyz = Matrix([[x], [y], [z]]) 164 x_hat = Matrix([[x], [0], [0]]) 165 y_hat = Matrix([[0], [y], [0]]) 166 z_hat = Matrix([[0], [0], [z]]) 163 167 q_xy = Matrix([[qx], [qy], [0]]) 164 168 q_abc = Matrix([[qa], [qb], [qc]]) … … 173 177 if 0: # forward calculations 174 178 vprint(q_xyz, jitter*q_abc, "apply jitter") 179 #vprint(xyz, jitter*z_hat, "r") 180 #mprint(J, jitter, "forward jitter") 175 181 vprint(dq_xyz, view*q_xyz, "apply view after jitter") 182 #mprint(V, view, "forward view") 176 183 177 184 #vprint(dq_xyz, view*jitter*q_abc, "combine view and jitter") 178 mprint( M, view*jitter, "forward matrix")185 mprint(R, view*jitter, "forward matrix") 179 186 180 187 if 1: # reverse calculations … … 196 203 mprint(q_abc[index,:], (R*q_xy)[index,:], "applied reverse matrix") 197 204 #mprint(q_abc, J*V*q_xy, "applied reverse matrix") 198 #mprint( M, jitter_inv*view_inv, "reverse matrix direct")205 #mprint(R[index,:2], jitter_inv*view_inv, "reverse matrix direct") 199 206 200 207 #vprint(q_abc, M*q_xy, "matrix application") … … 210 217 ) 211 218 219 if 1: 212 220 comment("\n\n==== symmetric ====") 213 221 print_steps(
Note: See TracChangeset
for help on using the changeset viewer.