Stack Overflow
Stack Overflow

Questions tagged [pardiso]

PARDISO (PARallel DIrect sparse SOlver) software packages from Intel MKL or Pardiso Project contain solution routines for solving sparse symmetric and unsymmetric linear systems.

PARDISO supplied via Intel MKL/Intel OneAPI and the Pardiso Project provides iterative and direct solution routines for solving large sparse symmetric and unsymmetric linear systems on shared memory multiprocessors. Both packages originate from a common predecessor.

38 questions
5
votes
1 answer

Replacement for (mkl) pardiso for arm64 (Apple Silicon)

For an academic numerical project we use as direct solvers mostly pardiso (from mkl) and cholmod. Now I switch to Apple Silicon (amazing fast, a quarter of the build time!) and cannot use Intel's mkl for native compiling any more. For BLAS openblas…
Fabian
  • 107
  • 1
  • 5
5
votes
1 answer

dask.delayed memory management when a single task can consume a lot of memory outside of python

I have some calculations calling the pardiso() solver from python. The solver allocates its own memory in a way that is opaque to python, but the pointers used to access that memory are stored in python. If I were to try and run these calculations…
cbf123
  • 51
  • 1
4
votes
1 answer

How to use (/install) the pardiso linear solver in ipopt using the pyomo modelling framework?

I am developing an optimization model using pyomo with python (python 3 - i manage packages with anaconda on windows). I need to use a non-linear solver, ipopt. It seems that the default linear solver used by ipopt (mumps or ma27) is relatively slow…
Jean Martin
  • 81
  • 3
4
votes
2 answers

Fast accessing elements of Compressed Sparse Row (CSR) sparse matrix

I want to test some of the newer sparse linear solvers and I want to know if there is a fast way of filling in the matrix. The format I'm interested is CSR (http://goo.gl/hLXYd). Let's say the matrix, in CSR format, is given by: values(num non-zero…
electrique
  • 431
  • 6
  • 18
3
votes
1 answer

'PyPardisoError: The Pardiso solver failed with error code -3. See Pardiso documentation for details'. How to solve it?

PyPardisoError: The Pardiso solver failed with error code -3. See Pardiso documentation for details. I was working on Jupyter Notebook and I found this error. What does it mean and how can I fix this problem?
Kevin
  • 35
  • 2
2
votes
0 answers

Julia exits when using Pardiso solver on Windows

I installed the Pardiso.jl with instructions on the package webpage. And when I run the test code, Julia just exited with the error message below. When I run from 568-574 it was not a problem. but when I try to use solve!() Julia just crushed. The…
CatieQQ
  • 21
  • 1
2
votes
0 answers

How to get intel MKL's cblas_dgemm and pardiso solvers in python?

I try to use Intel mkl from python. The only document I have found suggested the following: from ctypes import * # Load the share library mkl = cdll.LoadLibrary("./libmkl_rt.so") # For Intel MKL prior to version 10.3 us the created .so as below #…
zell
  • 9,830
  • 10
  • 62
  • 115
2
votes
0 answers

Problem building Pardiso in Julia on Windows

I am trying to use the Pardiso package in Julia together with the Pardiso library. However, I can’t build the package. When I look into the build log, I see this message: Looking in “C:\Users\Rodrigo.julia\packages\Pardiso\yZsYO\deps” for…
Rodrigo Pacheco
  • 65
  • 7
2
votes
0 answers

How to check the version of intel MKL pardiso in python?

I know that I can get the version of pardiso if I import MKL directly to python like this: import mkl mkl.get_version_string() But if I am using the following code to import intel MKL pardiso: import sys import ctypes def get_libmkl(): # linux…
f. c.
  • 1,095
  • 11
  • 26
2
votes
0 answers

Error in using getenv() to obtain number of threads

I am using the function getenv() to obtain the number of threads in the following manner: char* var; var = getenv("OMP_NUM_THREADS"); I get the following error: 'getenv': This function or variable may be unsafe. Consider using _dupenv_s instead. To…
Babaji
  • 398
  • 1
  • 4
  • 18
2
votes
1 answer

Solving sparse linear system with eigen/Intel MKL

I want to solve the Ax=b equation where A is a sparse matrix (1,964,568 x 1,964,568 nnz=75256446) and b is also sparse (1,964,568 x 1,964,568 nnz= 25354926) by using Eigen library in C++. At first I was trying to use the Eigen Sparse LU to solve my…
Vex994
  • 49
  • 1
  • 6
2
votes
2 answers

How to use .lib, .dll and .exp file in VC++?

I'm trying to solve linear systems of equations with pardiso. I have downloaded everything I might need from the website ( three files: .lib, .dll and .exp files). I know what's dll and lib, and know nothing about .exp. Moreover, I search the…
Alston
  • 33
  • 5
2
votes
2 answers

Compiling PARDISO linear solver test case with GCC

I am trying to compile a linear system solver using PARDISO. The test case (pardiso_sym.c) also downloaded from the same website above. I have the following files inside the directory: [gv@emerald my-pardiso]$ ls -lh total 1.3M -rw-r--r-- 1 gv…
neversaint
  • 60,904
  • 137
  • 310
  • 477
1
vote
1 answer

Using Eigen Pardiso with SparseMatrix which has 64 bit integer as StorageIndex

I recently got some problems when developing my code using Eigen Pardiso. Previously, I was using Eigen Pardiso to solve the 2D Poisson equation. The Matrix size is OK to be described by the C/C++ 32 bit int number. The code I was using is like…
Chen Cui
  • 11
  • 3
1
vote
0 answers

Accessing the pardiso solver from MKL 2020 in python

How can I access the pardiso solver from Intel's MKL in python? On anaconda3, I can import mkl. Here is my mkl's version. In [4]: mkl.get_version() Out[4]: {'MajorVersion': 2020, 'MinorVersion': 0, 'UpdateVersion': 2, 'ProductStatus':…
zell
  • 9,830
  • 10
  • 62
  • 115
1
2 3