Improved direct solver interface, now also accepts matrix_range and...

Improved direct solver interface, now also accepts matrix_range and matrix_proxy. Single-threaded module works well, OpenCL part needs further debugging.