Something which should be discussed is where to handle all the numerical issues. At the moment I agree with @h-g-s that all results should be handled with "care". But a lot of corrections are already applied during calculation in the generation process, e.g. values smaller than eps are handled as 0. I don't know if this is necessary, isn't that something the core solver must maintain? I guess some tests should be added to see how to make inputs and results numerical save.
Something which should be discussed is where to handle all the numerical issues. At the moment I agree with @h-g-s that all results should be handled with "care". But a lot of corrections are already applied during calculation in the generation process, e.g. values smaller than eps are handled as 0. I don't know if this is necessary, isn't that something the core solver must maintain? I guess some tests should be added to see how to make inputs and results numerical save.