Using MPI parMETIS with cluster_sparse_solver

Hello.

I am optimizing the `cluster_sparse_solver` runtime. In my case, the majority of the runtime is taken by phase `11`, symbolic factorization, with METIS. Additionally, only a single node is used in an otherwise `MPI`-enabled application.

I was wondering if there is a way to use `parMETIS` for fill-reducing ordering, in order to benefit from the cluster environment. One thing that would help tremendously is the source code for `cluster_sparse_solver`.

The version of MKL in question is mkl 11.2u3, which was bundled with composer_xe 2015 3.187.

Thanks!