Solves unsymmetric, structurally symmetric or symmetric systems, real or complex, positive definite or indefinite, hermitian.
LU decomposition with complete pivoting and extremely fast selected inversion strategies.
Parallel on SMPs and Cluster of SMPs.
Automatic combination of iterative and direct solver algorithms.
The PARDISO Version 6.0 was released in May 2018. It contains full support of multi-threaded Schur-complement computations and full support for parallel selected inversion. The new features are described in the following papers:
D. Kourounis, A. Fuchs, and O. Schenk, Towards the next generation of multiperiod optimal power flow solvers,
IEEE Transactions on Power Systems, vol. PP, no. 99, pp. 1–10, 2018.
F. Verbosio, A. D. Coninck, D. Kourounis, and O. Schenk, Enhancing the scalability of selected inversion factorization algorithms in genomic predictions,
Journal of Computational Science, vol. 22, no. Supplement C, pp. 99-108, 2017.
A. De Coninck, B. De Baets, D. Kourounis, F. Verbosio, O. Schenk, S. Maenhout, and J. Fostier, Needles: Toward large-scale genomic prediction with marker-by-environment interaction,
Genetics, vol. 203, no. 1, pp. 543–555, 2016.
Additional publications related to PARDISO are available on the webpage of the research group of Prof. Olaf Schenk, Advanced Computing Laboratory, Institute of Computational Science, USI Lugano, Switzerland.
Important: Please note that the Intel MKL version of PARDISO is based on our version from 2006 and that a lot of new features and improvements of PARDISO are not available in the Intel MKL library.
MKL PARDISO vs PARDISO 6 Performance Comparison
In the performance plot above we benchmarked PARDISO 6.0 against Intel MKL PARDISO version 2018 (Update 4) on up to 16 cores on a Intel-Xeon-CPU-E7email@example.comGHz (1TB, 16 cores). The matrices are from an industrial analog electronic circuit simulator based on nonlinear differential algebraic equations which are solved using implicit integration methods. The plots shows the performance acceleration for the factorization and solution phase against MKL PARDISO using one core.
How To Use PARDISO
- symmetric linear systems
- unsymmetric linear systems
- complex unsymmetric linear systems
- symmetric linear systems showing an example based on a discretization of the 2D Laplacian equation (laplace.c, laplace.h)
- scalable MPI C example for symmetric linear systems showing a parallel MPI example based on a discretization of the 2D Laplacian equation (laplace_mpi.tgz)
C++, Schur Complement
Julia interface for all types of linear systems in PARDISO
(by Kristoffer Carlsson , Chalmers University of Technology, Göteborg).
Binaries / Libraries
Current Available Libraries
Version 6.0 (Architecture x86-64, 64-bit)
|Compiler||Operating System||PARDISO Libraries|
|gcc/gfortran 7.2.0||Linux||libpardiso600-GNU720-X86-64.so||gcc/gfortran 8.0||Linux||libpardiso600-GNU800-X86-64.so|
|gcc/gfortran 8.1.0||MAC OSX 10.13.4 High Sierra||libpardiso600-MACOS-X86-64.dylib|
|Intel compiler, Windows (including optimzed BLAS)||Windows||libpardiso600-WIN-X86-64.dll, libpardiso600-WIN-X86-64.lib, libpardiso600-WIN-X86-64.exp|
|GNU compiler (MinGW)||Windows||libpardiso600-WIN-X86-64-MINGW.dll|
Other libraries can be compiled upon request. Please let us know in case that some of these libraries are not working for you. You can only download the files above if you are a registered user and agreeing to the license conditions.
Ipopt 3.12.9 (extended version), PARDISO 6.0 (threaded version)
You need to set the numbers of threads with, for example, if you would like to use 8 cores you should use
export OMP_NUM_THREADS=8 before running ipopt.mexa64 in Matlab.
In addition, you need a PARDISO license file pardiso.lic which can be obtained here.
|Architecture x86-64, 64-bit, Matlab R2016b||Linux||ipopt.mexa64, ipopt_auxdata.m, ipopt.opt, examplehs071.m|
|Architecture x86-64, 64-bit, Matlab R2016b||Mac OS X 10.13.4||ipopt.mexmaci64, ipopt_auxdata.m, ipopt.opt, examplehs071.m|
Academic Licenses / Corporate Licensing Program
Here you can dowload a cost-free academic license or a corporate license for the current release of PARDISO. The option corporate is only available for selected corporate or commercial entities. Please choose the desired license type:
Release of Version 6.0 for the R-INLA project including supernodal data compression and selected inversion software (release notes)
Release of new webpage, including user map and list of successful recent consulting projects.).
Release of Version 5.0.0, including new Schur-complement solver and selected inversion software (release notes).
Release of Version 4.1.2, including the distributed memory solver (release notes).
Release of Version 4.1.0 (release notes).
Map of pardiso users
Release of Version 4.0.0 (release notes).
The solver is now available for both academic and commercial use.
The solver PARDISO is now available under www.pardiso-project.org.
In case that you are using the new version 6.0 please cite:
- D. Kourounis, A. Fuchs, and O. Schenk, Towards the next generation of multiperiod optimal power flow solvers, IEEE Transactions on Power Systems, vol. PP, no. 99, pp. 1–10, 2018.
- F. Verbosio, A. D. Coninck, D. Kourounis, and O. Schenk, Enhancing the scalability of selected inversion factorization algorithms in genomic predictions, Journal of Computational Science, vol. 22, no. Supplement C, pp. 99 – 108, 2017.
- A. De Coninck, B. De Baets, D. Kourounis, F. Verbosio, O. Schenk, S. Maenhout, and J. Fostier, Needles: Toward large-scale genomic prediction with marker-by- environment interaction, Genetics, vol. 203, no. 1, pp. 543–555, 2016.
For support and consulting questions related to PARDISO 6.0, please send an e-mail to
We provide consulting services in the field of computational mathematics, scientific computing and high-performance computing including:
Benchmarking services and optimal use of the solver software
Effective use of high performance computational systems
Re-optimization of successfully implemented systems
Completed projects and a few selected consulting examples (since 2013)
Training and general solver recommendation for precision measuring technology product simulations.
Marathon Oil Company
Training and general solver recommendation for the Marathon Oil Reservoir Simulation product.
Personnel training in the area of OpenCL and sparse solvers.
Integration of PARDISO 6.0 into CST's electromagnetic design software.
Integration of PARDISO 5.0.0 into MOSEK's optimization software for solving mathematical optimization problems.
AutoForm Engineering AG
Integration of PARDISO 6.0 into AutoForm's automobile sheet metal forming tools.
Endress+Hauser Flowtec AG, Reinach, Switzerland
Integration of PARDISO 5.0.0 into finite-element tools at Endress+Hauser Flowtec AG.
U.S. Naval Observatory
Integration of PARDISO 5.0.0 into simulation tools at the U.S. Naval Observatory.
SMS Siemag AG
Integration of PARDISO 5.0.0 into simulation tools for metallurgical plant and rolling mill technology.
Integration of PARDISO 5.0.0 into Frustum's Distributed Optimization Service.
Integration of PARDISO 5.0.0 into Silvaco simulation tools.
Integration of PARDISO 5.0.0 into JCMsuite, a complete Finite Element software package for the computation of electromagnetic waves, elasticity and heat conduction.
Integration of PARDISO 4.1.2 into reservoir simulation tools at Saudi Aramco.
Map of PARDISO Users
Authors & Contact
Prof. Olaf Schenk
Institute of Computational Science, USI Lugano, Switzerland
USI Lugano, Switzerland, contributed with a block compression technique.
USI Lugano, Switzerland, contributed with a memory layer interface.
University of Minnessota, contributed with the METIS package
Technical University of Berlin, contributed with a matrix input checking tool
Dept. of Human Genetics, University of Chicago, contributed with an excellent MATLAB interface
Chalmers University of Technology, Göteborg, contributed with a Julia PARDISO interface