Is there a work-around to use Solve with Quaternions?

I just noticed in the documentation that one can’t use “dom = Quaternions” as the last parameter of Solve. Is there a work-around? Developmental plans?

Solve works with polynomial systems primarily, and I doubt there is any plan to change it to understand quaternions explicitly. As for workaround, it is difficult to respond without having a concrete example at hand. That is a general statement that should be considered whenever posting a question to the forum.

Thanks. I’m curious though. One can certainly have polynomials of quaternions just like one can for complexes. Is it the non-commuting nature that is the problem?

Noncommutativity would be an issue, yes. Might be possible to embed in matrix algebra though.

Hmmm, - matrice representations …that may be the “work-around.” Thanks.

You can simulate quaternions with a linear combination of the 3 Pauli spin matrices and the 2x2 identity matrix.

Right - so why doesn’t MMA do that with “Solve?” and allow “dom=Quaternion?” Or maybe Solve doesn’t solve matrix equations either? Haven’t checked.