haasad / PyPardiso

Python interface to the Intel MKL Pardiso library to solve large sparse linear systems of equations
BSD 3-Clause "New" or "Revised" License
135 stars 20 forks source link

Always convert b to fortran array if b is not fortran contiguous #40

Closed haasad closed 2 years ago

haasad commented 2 years ago

This fixes a bug where pypardiso produces wrong results if b is a one-dimensional slice of c-contiguous array

Fixes #39