Bicgstabl from IterativeSolvers.jl gives different result every time

[ettersi]

@stevengj

@antoine-levitt

@ettersi of course it’s wrong, but it’s a useful simple mental model to have when looking at results, no?

Thinking of floating-point arithmetic as non-deterministic seems sufficiently wrong that it leads people to make weird programming choices.

This seems like an interesting topic to delve in deeper! I opened a new topic on this here.