Tobias Ribizel
Tobias Ribizel
Examples where bitwise equivalence is important: * MPI runs where all nodes need to do the exactly same calculations * Debugging numerical issues is harder if non-determinism is involved
I would suggest adding a section to each LinOp (if it concerns apply) and function (if it concerns something else), and maybe a table giving an overview?
Stealing a bit of MATLAB's notation, we could also add a (defaulted) parameter to compute_norm2 that for a nxm matrix specifies whether we want column norms (output 1xm), row norms...
> compute_norm2 can relay on the output size, so we do not need to introduce another parameter Hm... while that is true, I am not sure whether that makes for...
All of these operations are heavily motivated by the solver design, and most of the uses of Dense right now are solver-related as multivectors, so I think we need to...
label!
check-format!
@thoasm These are some good points. As @pratikvn already suggested, I would probably move the whole `factorization` namespace to `preconditioner` instead of building two separate hierarchies. As for the ISAI,...
We had a discussion on how to combine the ILU preconditioner interface with preconditioners for the triangular solves (like ISAI #489). The simplest way in my opinion is to use...
Is the implementation intended for single-executor or distributed use? The interface seems to imply the former