Add an iterative scalar solver

An iterative diagonal solver was already implemented, this adds a scalar solver. It appears to be accurate and fast so far :-).