Sanjay

Hello Karl, viennacl::linalg::cuda::prod does work for int and long. I had to change the iteration bounds of the GPU-CPU result check to iterate over ublas_C, instead of ublas_A. -Sanjay

I'm looking at using viennacl for integer matrices, especially for matrix multiplication. I've been using the 1.7.1 release. Changing the ScalarType to int or long in examples/tutorials/blas3.cu results in both ublas::prod and viennacl::prod to compute the wrong answer. Surprisingly, the random number generator used to initialise the matrices generates numbers between 0 and 1, is static_cast to an int or long. This fills the matrices with zeros. I'm perplexed as to what's causing ublas's and viennacl's...

I'm looking at using viennacl for integer matrices, especially for matrix multiplication. I've been using the 1.7.1 release. Changing the ScalarType to int or long in examples/tutorials/blas3.cu results in both ublas::prod and viennacl::prod to compute the wrong answer. Surprisingly, the random number generator used to initialise the matrices generates numbers between 0 and 1, is static_cast to an int or long. This fills the matrices with zeros. I'm perplexed as to what's causing ublas's and viennacl's...

Hello Karl, Can I declare 2 viennacl::matrix<int>es and run linalg::cuda::prod_impl (matmul) on them and expect a correct answer, as of v1.7.1 ? I went through the source, and there seems to be a limitation on the datatypes only when prod is accessed through the BLAS like APIs, i.e.ViennaCLgemm and ViennaCLCUDAgemm_impl. I'm not 100% sure as I have not perused the code. Thanks !

Hello Karl, Can I declare 2 viennacl::matrix<int>es and run prod (matmul) on them as of v1.7.1 ?