Ask Your Question
0

Polynomials over Quaternions via Determinants

asked 2020-03-20 08:22:43 +0100

anonymous user

Anonymous

updated 2020-03-23 12:29:55 +0100

I currently trying to get polynomials over quaternions out of determinants of matrices. Example giving below:

Q.<i,j,k> = QuaternionAlgebra(SR, -1, -1)

A=matrix(4,4,[0,j,0,j,j,0,0,j,0,0,0,0,j,j,0,0])

I=matrix.identity(4)

P=(j*A-x*I)

P.det()

But it doesn't work. Any idea?

edit retag flag offensive close merge delete

Comments

1

What does "determinant" mean if the matrix is defined over a noncommutative ring? In your case you are lucky that P.change_ring(SR).det() makes sense, but that doesn't work in general.

rburing gravatar imagerburing ( 2020-03-23 12:51:23 +0100 )edit

1 Answer

Sort by » oldest newest most voted
0

answered 2020-03-26 09:44:58 +0100

Sébastien gravatar image

I do not know much about quaternion algebra, but I do know it is better to avoid involving the symbolic ring:

sage: Q.<i,j,k> = QuaternionAlgebra(QQ, -1, -1)
sage: K.<x> = Q[]
sage: K
Univariate Polynomial Ring in x over Quaternion Algebra (-1, -1) with base ring Rational Field

The variable x now lives in K:

sage: x.parent()
Univariate Polynomial Ring in x over Quaternion Algebra (-1, -1) with base ring Rational Field

and then:

sage: A = matrix(4,4,[0,j,0,j,j,0,0,j,0,0,0,0,j,j,0,0])
sage: I = matrix.identity(4)
sage: P = j*A - x*I
sage: P
[-x -1  0 -1]
[-1 -x  0 -1]
[ 0  0 -x  0]
[-1 -1  0 -x]
sage: P.parent()
Full MatrixSpace of 4 by 4 dense matrices over Univariate Polynomial Ring in x over Quaternion Algebra (-1, -1) with base ring Rational Field
sage: P.det()
x^4 - 3*x^2 + 2*x
edit flag offensive delete link more

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

1 follower

Stats

Asked: 2020-03-20 08:22:43 +0100

Seen: 492 times

Last updated: Mar 26 '20