I do this in a loop now:
while:
A = random_matrix(...)
if(not A.is_singular()):
....But this is innefficient. How do I generate random matrices with more concrete properties in sage?
You could do it one row at a time by choosing random vectors:
def random_nonsingular_matrix(base_ring=QQ, size=3):
V = base_ring**size
vectors = []
for i in range(size):
v = V.random_element()
while v in V.span(vectors):
v = V.random_element()
vectors.append(v)
return(matrix(vectors))Examples:
sage: random_nonsingular_matrix(size=4).determinant()
1337/102
sage: random_nonsingular_matrix(size=4).determinant()
127/16
sage: random_nonsingular_matrix(size=4).determinant()
815/16
sage: random_nonsingular_matrix(size=4).determinant()
96533/326It's unfortunately much slower than random_matrix:
sage: %timeit random_nonsingular_matrix(size=20)
10 loops, best of 3: 84.9 ms per loop
sage: %timeit random_matrix(QQ, 20, 20)
10000 loops, best of 3: 58.2 µs per loopThe documentation helps with that.
It can be accessed using
sage: random_matrix?and in particular hints at
sage: random_matrix(algorithm='unimodular')which will return a matrix of determinant one.
What about constructing a matrix of determinant K?
One could start from a determinant one matrix and multiply one row or one column by K.
One could then multiply by another random determinant one matrix (on the left or on the right or both) to mix things up a bit more.
Asked: 6 years ago
Seen: 2,108 times
Last updated: Mar 01 '19
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.
