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

float32 support #42

Closed xper0418 closed 1 year ago

xper0418 commented 2 years ago

Hi. As far as I know, PyPardiso only supports float64. I hope to use float32 to save memory. Any update plans?

haasad commented 2 years ago

Hi there, no plans to support float32 as of now. I assume the accuracy of the result will suffer a lot when using float32. Have you tried to solve your system with scipy's spsolve to check if the result is acceptable in terms of accuracy and computation time?

xper0418 commented 2 years ago

Thanks for the response. I tested it with sicpy, but the calculation was not finished, so I stopped in the middle. Since I already store the results as a float32, I don't think there will be any accuracy issues.