Added Cholesky factorisation and equation-solving to linalg. Modified linalg.diag to extract entries of main-diagonal if input array is 2D
Fix bugs in new Cholesky routines, make examples run as doctests, add documentation references. Fix eye doctest.
Add missing CULA exceptions.