I have a huge sparse function L (10 million*10 million). It is not symmetric. Most off-diagonals are zeros. I want to solve for sol=L\R with matlab. It is very slow. So I tried Cholosky factorization using "klu" function from SuiteSparse. However, it reports "klu mexFunction not found". Does anyone know what's the problem, or have any better suggestions?