direct_solve: opencl doesnt work, bug in prod(), benchmark extended on matrix-vector