Matrix Group over Symbolic Ring

asked 7 years ago

dimahphone
41 ●1 ●2 ●8

updated 1 year ago

FrédéricC

5426 ●4 ●44 ●121

I have problem on generate matrix group over symbolic ring.

First, I define

eta=I;
eta2=(1+I)*sqrt(2)/2;

Then, define a generator matrix

T=matrix(SR,4,[eta**(i*j)*eta2/2 for i in range(4) for j in range(4)]);

I try to make a matrix group by

G=MatrixGroup(T);

What I get is only very long computation that does not give any result. Can somebody help me? Thank you very much.

I have checked the order of T, I got 8 since T**8=I where I is identity matrix.

Comments

I've tried instead:

sage: K.<a> = QuadraticField( -1 )
sage: R.<X> = K[]
sage: L.<b> =K.extension( X^2 - 2 )
sage: eta = a
sage: eta2 = (1+a)*b/2
sage: T = matrix(L, 4, [eta**(k*kk)*eta2/2 for k in range(4) for kk in range(4)])
sage: T^8
[1 0 0 0]
[0 1 0 0]
[0 0 1 0]
[0 0 0 1]
sage: G = MatrixGroup(T)

running into NotImplementedError: Currently, only simple algebraic extensions are implemented in gap. Well... Perhaps gap cannot digest the input in both cases.

Which is the "real need"? Depending on it, the "other way to get the result" may be simpler or more complicated. In the above case one can use block matrices representing $\sqrt{-1}$ and $\sqrt 2$ as commuting $2\times 2$ matrices. If such cases are enough, i will write the code.

dan_fulea ( 7 years ago )

Do you mean "real need" is what I need on that code? Sorry if I am wrong understanding your comment.

dimahphone ( 7 years ago )

add a comment

2 Answers

Sort by » oldest newest most voted

answered 7 years ago

dimahphone
41 ●1 ●2 ●8

Case closed, I finally use CyclotomicField(8) and it works.

link

Comments

@dimahphone -- You can accept your answer by clicking on the accept (tick mark) button, so that this will appear as the accepted answer, and the question will appear as solved in the list of questions.

slelievre ( 7 years ago )

add a comment

answered 7 years ago

slelievre
17799 ●22 ●163 ●352 http://carva.org/samue...

Whenever possible, avoid working with Sage's symbolic ring.

Here, you can use the symbolic ring as an intermediate step but move to QQbar.

This will make your code work. So you could start as you did:

sage: eta = I
sage: eta2 = (1+I)*sqrt(2)/2

then instead of defining

sage: T = matrix(SR, 4, [eta**(i*j)*eta2/2 for i in range(4) for j in range(4)])

you could define

sage: T = matrix(QQbar, 4, [eta**(i*j)*eta2/2 for i in range(4) for j in range(4)])

And then things work.

sage: G = MatrixGroup(T)
sage: G
Matrix group over Algebraic Field with 1 generators (
[ 0.3535533905932738? + 0.3535533905932738?*I  0.3535533905932738? + 0.3535533905932738?*I  0.3535533905932738? + 0.3535533905932738?*I  0.3535533905932738? + 0.3535533905932738?*I]
[ 0.3535533905932738? + 0.3535533905932738?*I -0.3535533905932738? + 0.3535533905932738?*I -0.3535533905932738? - 0.3535533905932738?*I  0.3535533905932738? - 0.3535533905932738?*I]
[ 0.3535533905932738? + 0.3535533905932738?*I -0.3535533905932738? - 0.3535533905932738?*I  0.3535533905932738? + 0.3535533905932738?*I -0.3535533905932738? - 0.3535533905932738?*I]
[ 0.3535533905932738? + 0.3535533905932738?*I  0.3535533905932738? - 0.3535533905932738?*I -0.3535533905932738? - 0.3535533905932738?*I -0.3535533905932738? + 0.3535533905932738?*I]
)

Calculations in QQbar are exact and efficient.

If you need to see the radical expression of an element in QQbar (if it has one), use .radical_expression().

For example:

sage: T[0,0].radical_expression()
1/2*(-1)^(1/4)

link

Comments

Thank you for the suggestion. It works. But, I still have problem when I use QQbar. I can not do anything for groups I get. For example, I can not know the element order of the group. By hand, I get that the order is 8. Do you have any suggestion? Thank you very much. Sorry if it disturbs your time.

dimahphone ( 7 years ago )

add a comment

Your Answer

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

Add Answer

Matrix Group over Symbolic Ring

Comments

2 Answers

Comments

Comments

Your Answer

Question Tools

Stats

Related questions

Matrix Group over Symbolic Ring savecancel

Comments

2 Answers

Comments

Comments

Your Answer

Question Tools

Stats

Related questions

Matrix Group over Symbolic Ring