At the moment the functor for scal takes as template parameter the type of the
scalar.
This means that when doing something like this:
cusp::array1d<float,cusp::host_memory> y(1);
y[0] = 2.3;
cusp::blas::scal(y,2);
cusp::print_matrix(y);
The elements of y are casted to int during the multiply resulting in unexpected
results.
Original issue reported on code.google.com by filipe.c...@gmail.com on 28 Feb 2011 at 8:40
Original issue reported on code.google.com by
filipe.c...@gmail.com
on 28 Feb 2011 at 8:40