I started with this
N.<c,d,a1,a2,a3,a4,b1,b2,b3,b4,s> = QQ[]
H.<i,j,k> = QuaternionAlgebra(c,d)
a = a1 + a2 * i + a3 * j + a4 * k
b = b1 + b2 * i + b3 * j + b4 * k
a+b
a*b
but the next step does not work
gravity = vector([0, 0, 1])
a * gravity * a.conjugate()
Not sure but maybe you need to convert to matrix beforehand? If I append zero to the gravity vector I get
sage: gravity = vector([0, 0, 1, 0])
sage: a.matrix() * gravity * a.conjugate().matrix()
(-c^2*d*a2*a4 - c^2*a2*a4 - d*a1*a3 + a1*a3, c*d*a2*a3 + c*a1*a4 + d*a1*a4 - a2*a3, -c^2*a2^2 - c^2*a4^2 + a1^2 - a3^2, -c*a1*a2 - c*a3*a4 + a1*a2 - a3*a4)
Asked: 2015-07-31 03:33:17 +0100
Seen: 314 times
Last updated: Jul 31 '15
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.