Iterative solver for turbine model #10
Description
Currently, the turbine model solves the steady shallow water equations using a full LU decomposition at each nonlinear iteration. This is fine at low resolution, but gets expensive because the LU factorisation is O(n^3). It would be good if we could devise an alternative iterative solver strategy that is more scalable.