|
Hi I have a C program which does process on a large symmetric matrix with CBLAS and I need to allocate least possible memory in it. So I wanna pass triangular matrix which needs n(n+1)/2 memory, rather than full matrix with n*n. As far as I concened, CBLAS routines such as cblas_dtrsv takes full matrix as input argument, which just process wither lower or upper triangular part. Is it possible to use REAL triangular matrix in CBLAS? Thanks in advance -- Arya |
find it! packed format routines xTPSV (e.g. http://www.netlib.org/blas/dtpsv.f) should be used.