Denis Demidov

Another step to better dense matrix support is 3e452b7. It implements a tensordot() operation, borrowed from python's [numpy.tensordot](http://docs.scipy.org/doc/numpy/reference/generated/numpy.tensordot.html). The dense matrix-matrix product implemented above with reshape/reduce combination becomes ``` .cpp...

I am sorry for the shameless plug, but you could use [vexcl](https://github.com/ddemidov/vexcl) with the boost.compute backend (see http://vexcl.readthedocs.io/en/latest/interop.html). This will solve the problem with scalar kernel arguments, since vexcl by...

[Sure](https://github.com/ddemidov/vexcl/releases/tag/1.4.2), thanks for reminding me :).

The main problem here is that algorithms that use `vex::symbolic` are not allowed to use any data-dependent branching. This is because C++ does not allow overloading conditional operators. So I...

> Is there any reason why dt cannot be vex::symbolic? I don't see any right away, but with the logic you showed the update to `dt` will be wasted anyway...

> Why is that? Wouldn't the updated dt persist between kernels calls? I would then do something like: Ok, if you make dt a vector parameter, then sure. > I...

> Would you implement the loop in the kernel by manually updating the output stream That was what I meant, yes. > Based on your experience do you think I...

It is not possible to send a `vector_view` to a custom kernel as a pointer, because a `vector_view` is an expression that is evaluated on the fly by vexcl. You...

> There does not seem to be a way to define a VEX_FUNCTION from a string at runtime. There is `VEX_FUNCTION_S`: https://github.com/ddemidov/vexcl/blob/192137a2c1f74de97dc2859d8af76b822a9a4dd6/vexcl/function.hpp#L205-L210 > Would it be possible to expose this...

> not sure if this is meant by "static string expression" This should be an expression that always returns the same string. Also, VEX_FUNCTION macro unpacks into something like ```cpp...