Benchmarking Cyclotomic Linear Algebra

Systems:

Operating Systems / Hardware:

Benchmarks:

  1. Matrix addition
  2. Matrix vector product
  3. Matrix multiplication
  4. Determinant
  5. Rank
  6. Characteristic Polynomial [:../charpoly: some benchmarks here]
  7. Reduced row echelon form
  8. Sparse reduced row echelon form
  9. Solving A*x = b and A*X = B
  10. Kernel
  11. Rational decomposition

A good source of matrices that we care about:

Below, replace 23 by any of 23, 67, 89, 199, and p=2 in the last line below by any positive integer:

ModularSymbols_clear_cache()
eps = DirichletGroup(23*3, CyclotomicField(11)).1^2
M = ModularSymbols(eps); M
t = M.hecke_matrix(2)