Re: [sage-support] Re: Matrix_generic_dense.det() is quite slow compared to Numpy

2023-01-09 Thread Thierry Dumont
When doing numérical computations, you must use RDF floats, except for very special purpose. RDF is your processor's float. The default Real Field, based on MPFR (which by default as the same precision as RDF) is slow as it is based on MPFR, a very good library, but MPFR is a software based

[sage-support] Re: Matrix_generic_dense.det() is quite slow compared to Numpy

2023-01-08 Thread John H Palmieri
I am not at all an expert in numerical linear algebra, but (a) certainly parts of Sage can be improved and (b) contributions are welcome. Others can speak more knowledgeably about this particular issue. On Sunday, January 8, 2023 at 11:45:00 AM UTC-8 ivana...@gmail.com wrote: > Is there any

[sage-support] Re: Matrix_generic_dense.det() is quite slow compared to Numpy

2023-01-08 Thread 陳孜穎
Is there any space to improve the computation of MPFR which is used by RealField? John H Palmieri 在 2023年1月6日 星期五凌晨3:01:37 [UTC+8] 的信中寫道: > One way to speed it up might be to work with RDF or QQ or RLF. The > documentation for the generic determinant method says, "Note that for > matrices

[sage-support] Re: Matrix_generic_dense.det() is quite slow compared to Numpy

2023-01-05 Thread John H Palmieri
One way to speed it up might be to work with RDF or QQ or RLF. The documentation for the generic determinant method says, "Note that for matrices over most rings, more sophisticated algorithms can be used." sage: %time ones_matrix(RDF, 600, 600).determinant() CPU times: user 78.6 ms, sys: 7.6