Try to make log output more readable

[paveltomin]

Trying to make log output a bit more readable.

• Try to make more consistent indentation: time level messages - no indent, nonlinear iteration level - two or four spaces, and 8 spaces for residuals, deltas, etc and linear solver.
• Make residual outputs more consistent between different solvers.
• Add main residual output in line search.
• Remove some messages that appear debugish.
• Replace cout with GEOS_LOG/GEOS_FMT in some places.
• For aperture table: add message about adding aperture transition and print the modified table.
• Add fracture propagation iteration output for HydrofractureSolver

Example:

Time: 3.00e+00 s, dt: 1 s, Cycle: 3
  Fracture propagation iteration 0
    Attempt:  0, ConfigurationIter:  0, NewtonIter:  0
        ( R ) = ( 3.02e-01 )
        hydrofracture: Global solution scaling factor = 1
    Attempt:  0, ConfigurationIter:  0, NewtonIter:  1
        ( R ) = ( 5.39e-08 )
        Last LinSolve(iter,res) = (   1, 9.19e-10 )
++ Fracture propagation. Re-entering Newton Solve.
  Fracture propagation iteration 1
    Attempt:  0, ConfigurationIter:  0, NewtonIter:  0
        ( R ) = ( 3.92e+00 )
        hydrofracture: Global solution scaling factor = 1
    Attempt:  0, ConfigurationIter:  0, NewtonIter:  1
        ( R ) = ( 6.26e+00 )
        Last LinSolve(iter,res) = (   1, 1.91e-16 )
        Line search @ 0.500:              ( R ) = ( 2.48e+00 )
        hydrofracture: Global solution scaling factor = 1
    Attempt:  0, ConfigurationIter:  0, NewtonIter:  2
        ( R ) = ( 6.38e+00 )
        Last LinSolve(iter,res) = (   1, 1.08e-16 )
        Line search @ 0.500:              ( R ) = ( 2.48e+00 )
        hydrofracture: Global solution scaling factor = 1
    Attempt:  0, ConfigurationIter:  0, NewtonIter:  3
        ( R ) = ( 2.22e+01 )
        Last LinSolve(iter,res) = (   1, 2.79e-15 )
        Line search @ 0.500:              ( R ) = ( 1.47e+01 )
        Line search @ 0.250:              ( R ) = ( 5.28e+01 )
        Line search @ 0.125:              ( R ) = ( 5.94e+01 )
        Line search @ 0.062:              ( R ) = ( 5.56e+00 )
        Line search failed to produce reduced residual. Accepting iteration.
        hydrofracture: Global solution scaling factor = 1
    Attempt:  0, ConfigurationIter:  0, NewtonIter:  4
        ( R ) = ( 1.63e+00 )
        Last LinSolve(iter,res) = (   1, 1.95e-16 )
        hydrofracture: Global solution scaling factor = 1
    Attempt:  0, ConfigurationIter:  0, NewtonIter:  5
        ( R ) = ( 3.80e-02 )
        Last LinSolve(iter,res) = (   1, 9.69e-11 )
        hydrofracture: Global solution scaling factor = 1
    Attempt:  0, ConfigurationIter:  0, NewtonIter:  6
        ( R ) = ( 1.14e-03 )
        Last LinSolve(iter,res) = (   1, 4.44e-11 )
        hydrofracture: Global solution scaling factor = 1
    Attempt:  0, ConfigurationIter:  0, NewtonIter:  7
        ( R ) = ( 3.38e-05 )
        Last LinSolve(iter,res) = (   1, 5.82e-11 )
Time: 4.00e+00 s, dt: 1 s, Cycle: 4