Ask Your Question

List (or set) of matrices

asked 2013-03-23 18:08:40 -0500

sjcavazos2 gravatar image

updated 2014-11-10 09:41:11 -0500

tmonteil gravatar image

Hi, I am trying to write a function that will put all the 2x2 matrices over a finite field Z_p into a list or set. I am just switching from Mathematica to Sage, so this might be way off, but could you please help me go in the right direction?

def ListM(p): 
    for i1 in range(0,p): 
       for i2 in range(0,p): 
          for i3 in range(0,p):
             for i4 in range(0,p):

I am trying to create a dummy list UT, and recursively add to it. I want it to serve as a variable though. The output gives me some error as well. I've been on it for a few hours and cant seem to figure it out.

edit retag flag offensive close merge delete

2 answers

Sort by » oldest newest most voted

answered 2013-03-24 00:43:20 -0500

ppurka gravatar image

You should mention what error you got. For creating all the 2x2 matrices over some finite field use the MatrixSpace class.

sage: M = MatrixSpace(GF(3), 2, 2)
sage: M
Full MatrixSpace of 2 by 2 dense matrices over Finite Field of size 3
sage: Mlist = M.list() # List containing all the matrices
sage: len(Mlist)

As for your function, it is also not correct, unless you have copy-pasted incorrectly.

  1. First, the indentation in the function is wrong.
  2. Secondly, you need not define UT to be a symbolic variable because you are redefining UT to be a list.
  3. Thirdly, you are not returning anything from the function, so I don't know whether your function is incomplete, or this is what you intended.
edit flag offensive delete link more

answered 2013-03-24 02:41:21 -0500

updated 2013-03-24 03:40:36 -0500

To complete ppurka's answer, here is your code with corrections:

def ListM(p): 
    for i1 in range(p): 
       for i2 in range(p): 
          for i3 in range(p):
             for i4 in range(p):
    return UT

Note that output of ListM(3) looks ugly, but this is a formatting "bug" when you display such a list, output is correct.

Python has very nice syntax for "list comprehensions", you may write:

def rows(p):
    return [[x,y] for x in range(p) for y in range(p)]
def ListM(p): 
    return [matrix([r1,r2]) for r1 in rows(p) for r2 in rows(p)]

Of course in this case using MatrixSpace(GF(p),2,2) is definitely the right solution. Good luck with Sage !

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


Asked: 2013-03-23 18:08:40 -0500

Seen: 287 times

Last updated: Mar 24 '13