I need some (direct) sparse linear solver for fortran, which works with REAL16 data type. Any suggestions? Both Pardiso and MUMPS support only REAL8. (identical question: https://math.stackexchange.com/questions/4016856/sparse-linear-solver-in-fortran-working-with-real16?noredirect=1#comment8293130_4016856)

  • $\begingroup$ Possibly of use to you scicomp.stackexchange.com/q/8204/19136 $\endgroup$ Feb 8, 2021 at 18:12
  • $\begingroup$ Do you need sparse LU or sparse Cholesky? $\endgroup$ Feb 9, 2021 at 6:09
  • $\begingroup$ Cholesky is sufficient. $\endgroup$ Feb 9, 2021 at 12:01
  • 1
    $\begingroup$ It should be relatively easy to convert a simple older code like YSMP. I vaguely recall doing that back in the 90's. $\endgroup$ Feb 9, 2021 at 22:13


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.