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

2016-08-16 Thread Dmitriy Lyubimov
Ok, so the problem here seems to narrow down to what include files are included and in what order. In particular, sparse operations are broken iff both compressed_matrix.hpp and prod.hpp are included (and maybe forward.hpp). Strange that it doesn't break dense template operations but it is confir

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

2016-08-16 Thread Dmitriy Lyubimov
Karl, i can independently confirm the problem with prod_impl instantiation over expression of compressed times base_matrix into matrix type. I understand there are tests examples but something goes wrong with the straightforward code. We are compiling for open cl and open mp at the same time.

Re: [ViennaCL-devel] Interfacing with clMAGMA?

2016-08-16 Thread Charles Determan
Karl, I have tried now with the regular .size() calls and setting the leading dimensions with .internal_size(). Same error. The error is reported on line 274 of the clBLAS/src/libary/blas/xgemm.cc file. It appears this loop is where it is causing the problem. for (unsigned int i = 0; i < numKe

Re: [ViennaCL-devel] Interfacing with clMAGMA?

2016-08-16 Thread Karl Rupp
Hi, On 08/15/2016 08:56 PM, Charles Determan wrote: > Karl, > > I have the OpenCL backend enabled and I have tried: > > cl_mem bufA = A.handle().opencl_handle().get() > cl_mem bufB = B.handle().opencl_handle().get() > cl_mem bufC = C.handle().opencl_handle().get() > > cl_command_queue queue = >