I stumbeled across a different method to make the linear system less i'll conditioned. The approach in eq. (30) [1] seems to be very simple and easy to implement. It is based in the ideas discussed in [2]. Do you think there will be time on your side to investigate this? Removing the iterative refinement seems to increase the performance since only a single linear system has to be solve in each iteration.
NB: Thanks for this great library.
I stumbeled across a different method to make the linear system less i'll conditioned. The approach in eq. (30) [1] seems to be very simple and easy to implement. It is based in the ideas discussed in [2]. Do you think there will be time on your side to investigate this? Removing the iterative refinement seems to increase the performance since only a single linear system has to be solve in each iteration.
NB: Thanks for this great library.