Robust gradual incrementing of boundary conditions.

[tupek2]

Add a recursive call to advanceStep() to handle cases where the solver fails due to the boundary conditions changing rapidly. Also, add utilities to check if a matrix is symmetric and if to print matrix norms to diagnose matrix asymmetries.

[tupek2]

/style