#4753: [with patch, needs work] make sparse * sparse = dense mod p like 50
frickin' times faster.
----------------------------+-----------------------------------------------
Reporter: was | Owner: was
Type: defect | Status: new
Priority: critical | Milestone: sage-3.2.4
Component: linear algebra | Resolution:
Keywords: |
----------------------------+-----------------------------------------------
Changes (by mabshoff):
* priority: major => critical
* summary: [with patch, needs review] make sparse * sparse = dense mod p
like 50 frickin' times faster. => [with patch,
needs work] make sparse * sparse = dense mod p
like 50 frickin' times faster.
* milestone: sage-3.4 => sage-3.2.4
Comment:
This patch set causes one doctest failure:
{{{
sage -t "devel/sage/sage/modular/modsym/ambient.py"
**********************************************************************
File
"/scratch/mabshoff/sage-3.2.4.alpha0/devel/sage/sage/modular/modsym/ambient.py",
line 20:
sage: M.T(2).matrix().fcp('x')
Exception raised:
Traceback (most recent call last):
File
"/home/mabshoff/build/sage-3.2.4-cycle/sage-3.2.4.alpha0/local/bin/ncadoctest.py",
line 1231, in run_one_test
self.run_one_example(test, example, filename, compileflags)
File
"/home/mabshoff/build/sage-3.2.4-cycle/sage-3.2.4.alpha0/local/bin/sagedoctest.py",
line 38, in run_one_example
OrigDocTestRunner.run_one_example(self, test, example, filename,
compileflags)
File
"/home/mabshoff/build/sage-3.2.4-cycle/sage-3.2.4.alpha0/local/bin/ncadoctest.py",
line 1172, in run_one_example
compileflags, 1) in test.globs
File "<doctest __main__.example_0[6]>", line 1, in <module>
M.T(Integer(2)).matrix().fcp('x')###line 20:
sage: M.T(2).matrix().fcp('x')
File "/scratch/mabshoff/sage-3.2.4.alpha0/local/lib/python2.5/site-
packages/sage/modular/hecke/hecke_operator.py", line 459, in matrix
self.__matrix = self.parent().hecke_matrix(self.__n)
File "/scratch/mabshoff/sage-3.2.4.alpha0/local/lib/python2.5/site-
packages/sage/modular/hecke/algebra.py", line 201, in hecke_matrix
return self.__M.hecke_matrix(n)
File "/scratch/mabshoff/sage-3.2.4.alpha0/local/lib/python2.5/site-
packages/sage/modular/hecke/module.py", line 917, in hecke_matrix
T = self._compute_hecke_matrix(n)
File "/scratch/mabshoff/sage-3.2.4.alpha0/local/lib/python2.5/site-
packages/sage/modular/hecke/module.py", line 117, in _compute_hecke_matrix
return self._compute_hecke_matrix_prime(n)
File "/scratch/mabshoff/sage-3.2.4.alpha0/local/lib/python2.5/site-
packages/sage/modular/modsym/ambient.py", line 653, in
_compute_hecke_matrix_prime
Tp = W._matrix_times_matrix_dense(R)
File "matrix_modn_sparse.pyx", line 370, in
sage.matrix.matrix_modn_sparse.Matrix_modn_sparse._matrix_times_matrix_dense
(sage/matrix/matrix_modn_sparse.c:6236)
TypeError: Cannot convert
sage.matrix.matrix_mod2_dense.Matrix_mod2_dense to
sage.matrix.matrix_modn_dense.Matrix_modn_dense
**********************************************************************
}}}
Cheers,
Michael
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/4753#comment:10>
Sage <http://sagemath.org/>
Sage - Open Source Mathematical Software: Building the Car Instead of
Reinventing the Wheel
--~--~---------~--~----~------------~-------~--~----~
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en
-~----------~----~----~----~------~----~------~--~---
