Added requested changes

Author: arnoldkk13

include/aspect/simulator.h

@@ -1340,10 +1340,16 @@ namespace aspect
       /**
        * Perform a Newton line search to determine the optimal step length
        * along the search direction. After the update, the current_linearization_point
-      * is set to the old current_linearlization_point plus a suitable update.
+       * is set to the old current_linearlization_point plus a suitable update.
+       *
+      * @param dcr The defect correction residuals associated with the current nonlinear
+      * iteration.
+      * @param use_picard Whether a Picard iteration was used to update the nonlinear
+      * iteration (true) or a Newton update (false).
+      * @param search_direction The proposed update direction for the solution vector.
       *
-      * @return This function returns the updated residual after the line
-      * search is performed.
+       * @return This function returns the updated residual after the line
+       * search is performed.
        *
        * This function is implemented in
        * <code>source/simulator/helper_functions.cc</code>

source/simulator/helper_functions.cc

@@ -2629,9 +2629,8 @@ namespace aspect
     return initial_newton_residual;
   }
 
-  /**
-   * Performs a Newton line search for the do_one_defect_correction_Stokes_step function
-   */
+
+
   template <int dim>
   double Simulator<dim>::perform_line_search(const DefectCorrectionResiduals &dcr,
                                              const bool use_picard,
@@ -2640,6 +2639,7 @@ namespace aspect
     LinearAlgebra::BlockVector backup_linearization_point(introspection.index_sets.stokes_partitioning, mpi_communicator);
     double step_length_factor = 1.0;
     unsigned int line_search_iteration = 0;
+
     // Many parts of the solver depend on the block layout (velocity = 0,
     // pressure = 1). For example the linearized_stokes_initial_guess vector or the StokesBlock matrix
     // wrapper. Let us make sure that this holds (and shorten their names):
@@ -2650,9 +2650,9 @@ namespace aspect
     Assert(velocity_block_index == 0, ExcNotImplemented());
     Assert(pressure_block_index == 1, ExcNotImplemented());
     (void) pressure_block_index;
+
     // Do the loop for the line search at least once with the full step length.
     // If line search is disabled we will exit the loop in the first iteration.
-
     do
       {
         if (line_search_iteration == 0)
@@ -2717,6 +2717,7 @@ namespace aspect
   }
 
 
+
   template <int dim>
   double
   Simulator<dim>::compute_Eisenstat_Walker_linear_tolerance(const bool EisenstatWalkerChoiceOne,

source/simulator/solver_schemes.cc

@@ -709,10 +709,7 @@ namespace aspect
                                                             dcr.stokes_residuals.second,
                                                             dcr.residual,
                                                             dcr.residual_old);
-
-                pcout << "   The linear solver tolerance is set to "
-                      << parameters.linear_stokes_solver_tolerance
-                      << ". ";
+                pcout << "   The linear solver tolerance is set to " << parameters.linear_stokes_solver_tolerance << std::endl;
               }
           }