https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Sparse-least-squares/m-p/850927#M6565 <BR />Hi!<BR /><BR />You can use MKL LAPACK to solve a least squares problem, but in this case you should convert your data from sparse format into dense one. This approach leads to extensive usage of memory, so it is limited by the size of physical memory and size of solved task. <BR />To solve this problem in sparse format, you can use SpaseBLAS functions to compute A^T A and A^T b. For example you can use mkl_?csrmm function to compute matrix-matrix product in the CSR format and mkl_?csrgemv to compute matrix-vector product (also in CSR format, see MKL manual). Then MKL PARDISO (or DSS simplified interface to PARDISO) could be used to solve obtained system of equations. Note that PARDISO demonstrated good scalability on shared memory systems.<BR /><BR />- Best regards, Sergey<BR /> Wed, 11 Feb 2009 10:32:01 GMT Sergey_P_Intel2 2009-02-11T10:32:01Z https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Sparse-least-squares/m-p/850926#M6564 So, I was wondering what the best way to solve a sparse least squares problem (Ax=b) with the MKL is.<BR />The two ways that I can see are:<BR />1 - use LAPACK (which I don't believe is sparse, so probably not usable) <BR />2 - compute ATAx=ATb sparsely (perhaps using sparse BLAS?) and then use one of the sparse solvers, like PARDISO or dss.<BR /><BR />I don't have a lot of applied math experience, so if there's some better way, I would love to hear it. Right now I'm using CSparse's qr solver, looking to see if there's a simple way to get better performance.<BR /><BR />Thanks Mon, 02 Feb 2009 18:49:36 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Sparse-least-squares/m-p/850926#M6564 tfinniga 2009-02-02T18:49:36Z https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Sparse-least-squares/m-p/850927#M6565 <BR />Hi!<BR /><BR />You can use MKL LAPACK to solve a least squares problem, but in this case you should convert your data from sparse format into dense one. This approach leads to extensive usage of memory, so it is limited by the size of physical memory and size of solved task. <BR />To solve this problem in sparse format, you can use SpaseBLAS functions to compute A^T A and A^T b. For example you can use mkl_?csrmm function to compute matrix-matrix product in the CSR format and mkl_?csrgemv to compute matrix-vector product (also in CSR format, see MKL manual). Then MKL PARDISO (or DSS simplified interface to PARDISO) could be used to solve obtained system of equations. Note that PARDISO demonstrated good scalability on shared memory systems.<BR /><BR />- Best regards, Sergey<BR /> Wed, 11 Feb 2009 10:32:01 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Sparse-least-squares/m-p/850927#M6565 Sergey_P_Intel2 2009-02-11T10:32:01Z