I'm trying to define quaternion-valued differential forms in Sage. For example, take e0, e1, e2, e3 to be the quaternion generators and dx0, dx1, dx2, dx3 be the generators of the exterior algebra. So for example I'd like to be able to compute dX∧dX†, where dX=e0dx0+e1dx1+e2dx2+e3dx3, dX†=e0dx0−e1dx1−e2dx2−e3dx3.