Lower triangular matrix of a symmetric sparse matrix

Hi,

Is there any way to get lower (or upper) triangular matrix of a symmetric sparse matrix from pardiso solver? If not, what is the best way to do that using MKL intel in FORTRAN? I know I can get lower triangular matrix of a symmetric matrix with pbtrf() function but I could not find an function for sparse matrices!

The size of symmetric sparse matrix is around 10000*10000. 

Thank you,

Vahid