Added OpenCL kernels for element-wise operations.

Kernels for x += element_op(y) and x -= element_op(y) are included as well,
even though the expression template parser does not take them into account yet.