Sparse Matrix Inversion / Thomas method

I have differential equation that I am solving with FV method. From FV I have tri-diagonal sparse matrix and I use

SprD=sparse(DD); %DD is my FV matrix u=SprD\RHST % RHST is my Right hand side of the differential equation and u is the solution.

What I realized; Matlab does not work in parallel to solve this system. If I did not put sparse(DD) it would be parallel. I am trying to find if Matlab has an option to solve sparse Matrix system in parallel, or a parallel Thomas method.

Thank you.