Add support for f-order (column-major) arrays to linalg

This PR adds support for fortran-order arrays to `scikits.cuda.linalg.dot` and `scikits.cuda.linalg.cho_solve`. 

In the latter case, it also makes it possible to have a full B-matrix instead of just a single vector (i.e., compute inv(A)*B  when B is a matrix). However this will only work if both matrices are column-major.