Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-17 Thread Karl Rupp
Hi Charles, > Just adding my opinion here as I have been following this thread. Would > it be possible to have both the .so library and header only options > available or is it a strictly 'this-or-that' scenario? a header-only version should be possible, but exposes the user to all kinds of com

Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-17 Thread Charles Determan
Just adding my opinion here as I have been following this thread. Would it be possible to have both the .so library and header only options available or is it a strictly 'this-or-that' scenario? Regards, Charles On Wed, Aug 17, 2016 at 2:53 PM, Karl Rupp wrote: > Hi Dmitriy, > > > We could

Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-17 Thread Karl Rupp
Hi Dmitriy, > We could (and probably should?) add such a convenience header file > at the expense of increased compilation times (and reduced > encapsulation of source code against compiler issues). > > > +1 on single header! :) thanks for the feedback: https://github.com/viennacl/vie

Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-17 Thread Dmitriy Lyubimov
gt;> // set the ptrs for A >>> A_compressed_matrix.set(A_row_ptr, A_col_ptr, A_values_ptr, m, >>> s, NNz_A); >>> >>> // B is dense s so we only need s x n values. >>> double B_values[s * n] = {0}; >

Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-17 Thread Dmitriy Lyubimov
major dense_matrix, B can set the double values >> in the construcor >> // this is currently the constructor that we're using through >> scala/javacpp. >> const viennacl::matrix >>B_dense_matrix(B_v

Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-17 Thread Karl Rupp
B_dense_matrix(B_values_ptr, > viennacl::MAIN_MEMORY, s, n); > > > // perform multiplication and inside of a compressed_matrix > constructor > viennacl::matrix > C_dense_matrix(viennacl::linalg::prod(A_compressed_matrix ,

Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-16 Thread Dmitriy Lyubimov
mber; >> for (int i = 0; i< s * n; i++) { >> B_values[i] = randomNumber(); >> } >> >> double* B_values_ptr = B_values; >> >> >> // for our row_major dense_matrix, B can set the double values in the >> construcor >> // this is currently t

Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-16 Thread Dmitriy Lyubimov
const viennacl::matrix >B_dense_matrix(B_values_ptr, viennacl::MAIN_MEMORY, > s, n); > > > // perform multiplication and inside of a compressed_matrix constructor > viennacl::matrix C_dense_matrix(viennacl:: > linalg::prod(A_compressed_matrix

Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-08 Thread Andrew Palumbo
dl; // just exit with success for now if there are no runtime errors. return EXIT_SUCCESS; } ____ From: Karl Rupp Sent: Sunday, August 7, 2016 2:20:26 PM To: Andrew Palumbo; viennacl-devel@lists.sourceforge.net Subject: Re: [ViennaCL-devel] compressed_matrix %*%

Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-07 Thread Karl Rupp
Hi Andy, the relevant tests for sparse matrices times dense matrices are in tests/spmdm.cpp. In particular, I recreated a test case based on your description and couldn't find any issues: viennacl::compressed_matrix compressed_A; viennacl::matrix B1(std_A.size(), cols_rhs); viennacl::mat

Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-05 Thread Andrew Palumbo
Thanks Karli, much appreciated as always. Original message From: Karl Rupp Date: 08/05/2016 5:13 PM (GMT-05:00) To: Andrew Palumbo , viennacl-devel@lists.sourceforge.net Subject: Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base Hi Andy, hmm, a first look doesn't r

Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-05 Thread Karl Rupp
------- > *From:* Andrew Palumbo > *Sent:* Thursday, August 4, 2016 11:06:33 AM > *To:* Karl Rupp; viennacl-devel@lists.sourceforge.net > *Subject:* Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base > > Thanks Karl- Appreci

Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-05 Thread Andrew Palumbo
unsigned int; DistanceT = long int]’ {...} From: Andrew Palumbo Sent: Thursday, August 4, 2016 11:06:33 AM To: Karl Rupp; viennacl-devel@lists.sourceforge.net Subject: Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base Thanks Karl- Appreciate it! From:

Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-04 Thread Andrew Palumbo
Thanks Karl- Appreciate it! From: Karl Rupp Sent: Thursday, August 4, 2016 3:36:53 AM To: Andrew Palumbo; viennacl-devel@lists.sourceforge.net Subject: Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base Hi Andrew, On 08/04/2016 01:33 AM, Andrew Palumbo

Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-04 Thread Karl Rupp
li > > *From:* Andrew Palumbo > *Sent:* Wednesday, August 3, 2016 6:44:10 PM > *To:* Karl Rupp; viennacl-devel@lists.sourceforge.net > *Subject:* Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base >

Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-03 Thread Andrew Palumbo
first of the matrix before converting it to CSR, etc? Curious, as we may be able to shave a small amount of time if so. Thanks! Andy From: Karl Rupp Sent: Wednesday, August 3, 2016 5:28:58 PM To: Andrew Palumbo; viennacl-devel@lists.sourceforge.net Subject: Re:

Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-03 Thread Andrew Palumbo
e may be able to shave a small amount of time if so. Thanks! Andy From: Karl Rupp Sent: Wednesday, August 3, 2016 5:28:58 PM To: Andrew Palumbo; viennacl-devel@lists.sourceforge.net Subject: Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base Hi Andrew, > I'm having some tr

Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-03 Thread Karl Rupp
Hi Andrew, > I'm having some trouble with sparse `compressed_matrix` `matrix`(base) > matrix multiplication. This is supported, correct? Yes. Could you please let us know what you have tried already? It shouldn't be any more code to write than viennacl::compressed_matrix A(...); viennacl::

[ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-03 Thread Andrew Palumbo
Hello, I'm having some trouble with sparse `compressed_matrix` `matrix`(base) matrix multiplication. This is supported, correct? I've been trying to use the: template< typename SparseMatrixType, typename SCALARTYPE> typename viennacl::enable_if< viennacl::is_any_sparse_matrix::value