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<NumericT> compressed_A;
  viennacl::matrix<NumericT, FactorLayoutT>  B1(std_A.size(), cols_rhs);
  viennacl::matrix_base<NumericT> B1_ref(B1);
  viennacl::matrix_base<NumericT> 
C2(viennacl::linalg::prod(compressed_A, B1_ref));

compiles cleanly. Could you please provide a code snippet demonstrating 
the problem you are encountering?

Thanks and best regards,
Karli



On 08/05/2016 09:04 PM, Andrew Palumbo wrote:
> Hi Karl,
>
>
> I've been trying to implement tests for:
>
>
>      matrix_base<double> C = compressed_matrix<double> A %*%
>
>      matrix_base<double,row_major> B.
>
>
> I cant find in the code or the documentation any constructor for
> matrix_base<T>(
>
> matrix_expression<const viennacl::compressed_matrix<T>, const
> viennacl::matrix_base<T>, viennacl::op_prod>)
>
> ie. a mixed expression of compressed_matrix and matrix_base
>
> and get a compilation error when I try to instantiate a:
>
>      matrix_base<double>(matrix_expression<const
>      viennacl::compressed_matrix<double>, const
> viennacl::matrix_base<double>,
>      viennacl::op_prod>)
>
> Is there a transformation that I need to do from this
>
>      matrix_expression<compressed_matrix<double>, matrix_base<double>,
> op_prod>
>
> to something else so that I may be able to initialize a matrix_base (or
> possibly even a compressed_matrix) from it?
>
> The compilation error that i get is below.
>
> Thanks,
>
> Andy
>


------------------------------------------------------------------------------
_______________________________________________
ViennaCL-devel mailing list
ViennaCL-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/viennacl-devel

Reply via email to