Opened 16 years ago
Closed 13 years ago
#623 closed task (remind)
Force CHOLMOD to use an optimiced BLAS library
Reported by: | Víctor de Buen Remiro | Owned by: | Jorge |
---|---|---|---|
Priority: | normal | Milestone: | Numerical methods |
Component: | Math | Version: | 2.0.1 |
Severity: | normal | Keywords: | BSR, CHOLMOD, BLAS, ATLAS |
Cc: |
Description
CHOLMOD efficiency decays for dense and block-dense matrices, and it's almost sure due to it is using a non optimiced BLAS library.
It's neccessary to force CHOLMOD to use ATLAS or other optimiced libraries.
Change History (6)
comment:1 Changed 16 years ago by
Owner: | changed from Víctor de Buen Remiro to Jorge |
---|---|
Status: | new → assigned |
comment:2 Changed 16 years ago by
Keywords: | BSR added |
---|
comment:3 Changed 16 years ago by
Priority: | high → normal |
---|---|
Severity: | critical → normal |
comment:4 Changed 16 years ago by
Milestone: | → BSR Numerical methods |
---|---|
Version: | → 2.0.1 |
comment:5 Changed 16 years ago by
In windows is posible to force CHOLMOD to use ATLAS libraries instead of non optimized ones (libblas32.dll and liblapack32.dll), by mean of this environment variables
set BLAPI_BLAS_LIB=libcblas.dll set BLAPI_LAPACK_LIB=libclapack.dll set BLAPI_VERBOSE=1
However, it causes some problems in release mode that cannot be reproduce in debug mode compiling trunk/tol_tests/tol/LinearAlgebra/vmatrix/test_0012/test.tol
comment:6 Changed 13 years ago by
Resolution: | → remind |
---|---|
Status: | assigned → closed |
A new analysis of Cholesky performance for semi-sparse matrices has been codified in TOL file
trunk/tol_tests/tol/LinearAlgebra/vmatrix/test_0014/performance.tol
Results are stored in trunk/tol_tests/tol/LinearAlgebra/vmatrix/test_0014/performance.ods
Main conclusion is that efficience decays specially when CholeskiFactor of X'X is calculated directly over a non square matrix X that have a density of non null zeros not very low. Then it's a better option to calculate explicitly X'X and then to apply CholeskiFactor. Obtained simulation results show that density frontier is above 13%.
Using BLAS optimized libraries can help in generic terms but it doesn't seem very probable that could solve the problem by itself.