Add preconditioner using iterative solver + use in a prec
Add inner preconditioner: gmres+ilu
Move some wave solution code into their own functions

(So that we can call them from elsewhere)
Fix constness on exact_solution(..) in problem classes
Initial versions of error norms for wave solution