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$– Carl ChristianFeb 8, 2021 at 18:12
-
$\begingroup$ Do you need sparse LU or sparse Cholesky? $\endgroup$– Brian BorchersFeb 9, 2021 at 6:09
-
$\begingroup$ Cholesky is sufficient. $\endgroup$– Vítězslav ŠtemberaFeb 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$– Brian BorchersFeb 9, 2021 at 22:13