> http://arxiv.org/pdf/math/0307321v2.pdf > Can this be implemented? I am not aware of Group theory :) > Apparently, this has a complexity of O(n^2) !
Well, O(n^2) does not have any practical value without specifying the constants. For example, 10000000000*n^2 won't be faster than 2*n^3 for any practical matrix sizes ;-) Either way, we are happy about any contributions. :-) Best regards, Karli > ------------------------------------------------------------------------ > *From:* Karl Rupp <r...@iue.tuwien.ac.at> > *To:* Sumit Kumar <dost_4_e...@yahoo.com> > *Cc:* "viennacl-devel@lists.sourceforge.net" > <viennacl-devel@lists.sourceforge.net> > *Sent:* Monday, August 3, 2015 1:23 AM > *Subject:* Re: [ViennaCL-devel] ViennaCL reductions > > Hi Sumit, > > > I was trying to run sparse matrix multiplication, but one of my explicit > > template typedefs had Int. After some digging, I found out that CSR only > > supported Float or double. Is there any reason for this? Can we also > > have support for other templates? (like int ?) > > It is technically possible, but we haven't implemented it yet. After 5 > years you are the first to even ask for it ;-) > > > > Another thing would be the alignment order. Suppose I have a Row-major > > Sparse Eigen matrix, then I can copy it to a (Row-Major ?) VCL > > compressed matrix. What about a column-major sparse Eigen matrix? > > If you can point me to a fast, massively parallel column-major > matrix-vector multiplication routine, I look into it. However, as far as > I know, there is no such routine for general sparse matrices, hence it > does not make sense for us to support it. > > Best regards, > Karli > > > > > > > > ------------------------------------------------------------------------ > > *From:* Karl Rupp <r...@iue.tuwien.ac.at <mailto:r...@iue.tuwien.ac.at>> > > *To:* Sumit Kumar <dost_4_e...@yahoo.com <mailto:dost_4_e...@yahoo.com>> > > *Cc:* "viennacl-devel@lists.sourceforge.net > <mailto:viennacl-devel@lists.sourceforge.net>" > > <viennacl-devel@lists.sourceforge.net > <mailto:viennacl-devel@lists.sourceforge.net>> > > *Sent:* Friday, July 31, 2015 9:04 PM > > *Subject:* Re: [ViennaCL-devel] ViennaCL reductions > > > > Hi Sumit, > > > > > I am aware that Eigen can do it for its matrices and I am also aware > > > that VCL cannot do it natively. My question was this: > > > In your example of interfacing with Eigen, you have shown a VCL dense > > > matrix interfacing with an Eigen dense matrix. Do you have any example > > > of interfacing an Eigen Sparse matrix with a VCL dense matrix? > > > > No, we don't have this. > > > > > > > > > > Best regards, > > Karli > > > > > > > > > > > ------------------------------------------------------------------------------ _______________________________________________ ViennaCL-devel mailing list ViennaCL-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/viennacl-devel