[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: benchmark octave with atlas and vecLib
From: |
David Bateman |
Subject: |
Re: benchmark octave with atlas and vecLib |
Date: |
Thu, 6 Nov 2003 10:52:12 +0100 |
User-agent: |
Mutt/1.3.28i |
It's not a problem with ATLAS... See liboctave/Changelog
2003-02-18 David Bateman <address@hidden>
* dMatrix.cc (Matrix::inverse, Matrix::determinant, Matrix::solve):
Use Lapack instead of Linpack.
* Cmatrix.cc (ComplexMatrix::inverse, ComplexMatrix::determinant,
ComplexMatrix::solve): Likewise.
* dMatrix.cc (Matrix::determinant, Matrix::inverse): New arg,
calc_cond. If 0, skip condition number calculation.
* CMatrix.cc (ComplexMatrix::determinant, ComplexMatrix::inverse):
Likewise.
* CmplxLU.cc (ComplexLU::ComplexLU): Allow non-square matrices.
* dbleLU.cc (LU::LU): Likewise.
* base-lu.cc (base_lu::L), base_lu::U, base_lu::P): Likewise.
Basically 2.1.36 uses LINPACK for many of its matrix operations and so
doesn't have any advantage of using ATLAS. I removed all LINPACK functions
earlier this year, and thus the speed up you are seeing in the specific
functions you observed. "chol" was using ATLAS in 2.1.36, but was one of
the only functions that was.
As for random, I don't think you have octave-forge install, where the faster
"rand" and "randn"
> Sorting of 2,000,000 random values__________________ (sec): 8.45625
versus
> Sorting of 2,000,000 random values__________________ (sec): 9.244
Try installing octave-forge and rerunning your tests.
Regards
D.
--
David Bateman address@hidden
Motorola CRM +33 1 69 35 48 04 (Ph)
Parc Les Algorithmes, Commune de St Aubin +33 1 69 35 77 01 (Fax)
91193 Gif-Sur-Yvette FRANCE
The information contained in this communication has been classified as:
[x] General Business Information
[ ] Motorola Internal Use Only
[ ] Motorola Confidential Proprietary
-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.
Octave's home on the web: http://www.octave.org
How to fund new projects: http://www.octave.org/funding.html
Subscription information: http://www.octave.org/archive.html
-------------------------------------------------------------