Update the PARDISO interface to match other sparse solvers.

- Add support for Upper or Lower inputs.
- Add supports for sparse RHS
- Remove transposed cases, remove ordering method interface
- Add full access to PARDISO parameters

1 files changed, 2 insertions, 1 deletions

diff --git a/test/sparse_solver.h b/test/sparse_solver.h
index 51bb33a92..1a5483050 100644
--- a/test/sparse_solver.h
+++ b/test/sparse_solver.h
@@ -101,6 +101,7 @@ template<typename Solver> void check_sparse_spd_solving(Solver& solver)
 {
   typedef typename Solver::MatrixType Mat;
   typedef typename Mat::Scalar Scalar;
+  typedef SparseMatrix<Scalar,ColMajor> SpMat;
   typedef Matrix<Scalar,Dynamic,Dynamic> DenseMatrix;
   typedef Matrix<Scalar,Dynamic,1> DenseVector;
 
@@ -112,7 +113,7 @@ template<typename Solver> void check_sparse_spd_solving(Solver& solver)
   // generate the right hand sides
   int rhsCols = internal::random<int>(1,16);
   double density = (std::max)(8./(size*rhsCols), 0.1);
-  Mat B(size,rhsCols);
+  SpMat B(size,rhsCols);
   DenseVector b = DenseVector::Random(size);
   DenseMatrix dB(size,rhsCols);
   initSparse<Scalar>(density, dB, B);