I noticed that the (iterated advection and) defect correction and newton solver have diverged, both in output and code, which should not have happened. In this pull request I combine them and let the defect correction version call the newton one with a parameter to state that it is not using the newton solver. This should prevent these kind of issue in the future.
I noticed that the (iterated advection and) defect correction and newton solver have diverged, both in output and code, which should not have happened. In this pull request I combine them and let the defect correction version call the newton one with a parameter to state that it is not using the newton solver. This should prevent these kind of issue in the future.