I put that typecast back in and now i have things properly offset gain. I look forward to your fix. I need the SpMV to be asynchronous, and the typecast is a barrier.
-E On Wed, Jul 24, 2013 at 4:25 PM, Evan Bollig <bol...@gmail.com> wrote: > Karl, I have some issues with that new project operation. I want to > assign (in Matlab notation): > > b(5:10) = A(1:5,1:10) x(1:10) > > With: > > viennacl::range r4(115, 334); > project(out_deriv, r4) = viennacl::linalg::prod(DM_b, u_gpu); > > the assignment fills vector elements out_deriv(0:219) . The matrix in > is unprojected and the vector in is unprojected, so the output > projection is not offsetting the write properly. > > -E > > On Mon, Jul 22, 2013 at 2:59 PM, Evan Bollig <bol...@gmail.com> wrote: >> Thanks! Ill pull in the latest dev and test against that. >> >> Cheers, >> -E >> >> On Mon, Jul 22, 2013 at 1:47 PM, Karl Rupp <r...@iue.tuwien.ac.at> wrote: >>> Hi Evan, >>> >>> I just pushed support for sparse matrix-vector products when using >>> vector-ranges and vector-slices. >>> >>> https://github.com/viennacl/viennacl-dev/commit/7e68860886d6461e483a5965e0ede68146c74189 >>> >>> You cannot take ranges and slices of the matrix at the moment. >>> >>> Best regards, >>> Karli >>> >>> >>> >>> On 07/21/2013 04:23 PM, Evan Bollig wrote: >>>> >>>> Btw Karl: >>>> >>>> vector_range_M = viennacl::linalg::prod(A_MxN, vector_N); >>>> >>>> does not work. That would be equivalent to: >>>> >>>> vector_range_M = project( viennacl::linalg::prod(A_MxN, vector_N) , 0,M ); >>>> >>>> -E >>>> >>>> >>>> >>>> >>>> On Sun, Jul 21, 2013 at 4:13 PM, Evan Bollig <bol...@gmail.com> wrote: >>>>> >>>>> Ah I see. Thanks for the clarification on projecting the LHS. I was >>>>> trying to assign the result to an unprojected vector assuming (vector >>>>> - vector_range) would return a vector and (vector_range - vector) >>>>> would return a vector_range (i.e., operation returns operand type). >>>>> >>>>> -E >>>>> >>>>> On Sun, Jul 21, 2013 at 4:09 PM, Karl Rupp <r...@iue.tuwien.ac.at> wrote: >>>>>> >>>>>> Hi Evan, >>>>>> >>>>>> >>>>>>> Can you refresh my memory on how ViennaCL behaves with this operation: >>>>>>> >>>>>>> >>>>>>> vector = vector - vector_range >>>>>>> >>>>>>> I would like to subtract an M < N sized vector_range from an N sized >>>>>>> vector (by operating on the first M elements). Is this supported in >>>>>>> 1.4.2? >>>>>> >>>>>> >>>>>> >>>>>> Yes, this is supported. You can just use >>>>>> >>>>>> project(x, range(a, b)) = project(y, range(c, d)) - z; >>>>>> >>>>>> where z is some existing vector_range. Mixing of vectors with ranges and >>>>>> slices is perfectly fine and supported. >>>>>> >>>>>> Best regards, >>>>>> Karli >>>>>> >>>>> >>>>> >>>>> >>>>> -- >>>>> -Evan Bollig >>>>> bol...@gmail.com >>>>> bol...@scs.fsu.edu >>>> >>>> >>>> >>>> >>> >> >> >> >> -- >> -Evan Bollig >> bol...@gmail.com >> bol...@scs.fsu.edu > > > > -- > -Evan Bollig > bol...@gmail.com > bol...@scs.fsu.edu -- -Evan Bollig bol...@gmail.com bol...@scs.fsu.edu ------------------------------------------------------------------------------ See everything from the browser to the database with AppDynamics Get end-to-end visibility with application monitoring from AppDynamics Isolate bottlenecks and diagnose root cause in seconds. Start your free trial of AppDynamics Pro today! http://pubads.g.doubleclick.net/gampad/clk?id=48808831&iu=/4140/ostg.clktrk _______________________________________________ ViennaCL-devel mailing list ViennaCL-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/viennacl-devel
