Differences between revisions 2 and 3
Revision 2 as of 2007-02-15 19:24:53
Size: 407
Editor: wstein
Comment:
Revision 3 as of 2007-02-16 07:41:44
Size: 1004
Editor: wstein
Comment:
Deletions are marked like this. Additions are marked like this.
Line 2: Line 2:

== Issues to address ==

 * Robert Bradshaw implemented multimodular matrix multiply over ZZ.
    1. This seems to work fine on 32-bit machines, but is totally broken on 64-bit machines, so is currently disabled (in sage-2.1.2).
    2. It is interesting to fine tune the algorithm, and decide when to switch over to a multimodular method.
    3. Why is this so much '''slower''' than linbox?
    4. Why is it slower than MAGMA? (How much slower?)

 * Linbox: Charpoly and minpoly over ZZ.
    1. They hang on 32-bit Debian Linux and on sage.math, so are disabled in the sage-2.1.2.

William Stein: Optimized exact sparse and dense linear algebra (especially for computing modular forms)

Issues to address

  • Robert Bradshaw implemented multimodular matrix multiply over ZZ.
    1. This seems to work fine on 32-bit machines, but is totally broken on 64-bit machines, so is currently disabled (in sage-2.1.2).
    2. It is interesting to fine tune the algorithm, and decide when to switch over to a multimodular method.
    3. Why is this so much slower than linbox?

    4. Why is it slower than MAGMA? (How much slower?)
  • Linbox: Charpoly and minpoly over ZZ.
    1. They hang on 32-bit Debian Linux and on sage.math, so are disabled in the sage-2.1.2.

Benchmarks to help direct our work

New code written from scratch in C/SageX/Python

IML -- Integer Matrix Library

Linbox

The SmithForm (or invariant factors) algorithm, which gives the invariant factors, is literally a hundred times slower than MAGMA.

days3/sprints/exactlinalg (last edited 2008-11-14 13:41:58 by anonymous)