KratosShapeOptimizationApplication 10.4.4__py3-none-any.whl

KratosMultiphysics/ShapeOptimizationApplication/loggers/design_logger_base.py

@@ -0,0 +1,28 @@
+# ==============================================================================
+#  KratosShapeOptimizationApplication
+#
+#  License:         BSD License
+#                   license: ShapeOptimizationApplication/license.txt
+#
+#  Main authors:    Baumgaertner Daniel, https://github.com/dbaumgaertner
+#                   Geiser Armin, https://github.com/armingeiser
+#
+# ==============================================================================
+
+
+# ==============================================================================
+class DesignLogger:
+
+    # --------------------------------------------------------------------------
+    def InitializeLogging( self ):
+        raise RuntimeError("Design logger base class is called. Please check your implementation of the function >> InitializeLogging << .")
+
+    # --------------------------------------------------------------------------
+    def LogCurrentDesign( self, optimizationIteration ):
+        raise RuntimeError("Design logger base class is called. Please check your implementation of the function >> LogCurrentDesign << .")
+
+    #---------------------------------------------------------------------------
+    def FinalizeLogging( self ):
+        raise RuntimeError("Design logger base class is called. Please check your implementation of the function >> FinalizeLogging << .")
+
+# ==============================================================================

KratosMultiphysics/ShapeOptimizationApplication/loggers/design_logger_gid.py

@@ -0,0 +1,112 @@
+# ==============================================================================
+#  KratosShapeOptimizationApplication
+#
+#  License:         BSD License
+#                   license: ShapeOptimizationApplication/license.txt
+#
+#  Main authors:    Baumgaertner Daniel, https://github.com/dbaumgaertner
+#                   Geiser Armin, https://github.com/armingeiser
+#
+# ==============================================================================
+
+
+# Kratos Core and Apps
+import KratosMultiphysics as KM
+
+# For GID output
+from KratosMultiphysics.gid_output_process import GiDOutputProcess
+
+# Import logger base classes
+from KratosMultiphysics.ShapeOptimizationApplication.loggers.design_logger_base import DesignLogger
+
+# ==============================================================================
+class DesignLoggerGID( DesignLogger ):
+
+    # --------------------------------------------------------------------------
+    def __init__( self, model_part_controller, optimization_settings ):
+        self.output_settings = optimization_settings["output"]
+        minimal_gid_parameters = KM.Parameters("""
+        {
+            "name"           : "gid",
+            "gid_parameters" : { }
+        }""")
+
+        self.output_settings["output_format"].ValidateAndAssignDefaults(minimal_gid_parameters)
+        self.output_settings["output_format"]["gid_parameters"].RecursivelyValidateAndAssignDefaults(GiDOutputProcess.defaults)
+
+        self.optimization_model_part = model_part_controller.GetOptimizationModelPart()
+        self.design_surface = model_part_controller.GetDesignSurface()
+
+        self.write_design_surface = False
+        self.write_optimization_model_part = False
+        self.design_history_filename = None
+
+        self.__DetermineOutputMode()
+        self.__ModifySettingsToMatchDefaultGiDOutputProcess()
+        self.__CreateGiDIO()
+
+    # --------------------------------------------------------------------------
+    def __DetermineOutputMode( self ):
+        output_mode = self.output_settings["design_output_mode"].GetString()
+
+        if output_mode == "write_design_surface":
+            self.write_design_surface = True
+            self.design_history_filename = self.design_surface.Name
+        elif output_mode == "write_optimization_model_part":
+            if self.optimization_model_part.NumberOfElements() == 0:
+                raise NameError("Output of optimization model part in Gid-format requires definition of elements. No elements are given in current mdpa! You may change the design output mode.")
+            self.write_optimization_model_part = True
+            self.design_history_filename = self.optimization_model_part.Name
+        else:
+            raise NameError("The following design output mode is not defined within a GiD output (name may be misspelled): " + output_mode)
+
+    # --------------------------------------------------------------------------
+    def __ModifySettingsToMatchDefaultGiDOutputProcess( self ):
+        gid_parameters = self.output_settings["output_format"]["gid_parameters"]
+
+        # Add nodal results
+        self.output_settings["output_format"]["gid_parameters"]["result_file_configuration"]["nodal_results"] = self.output_settings["nodal_results"]
+
+        # Set condition flag
+        if not gid_parameters["result_file_configuration"]["gidpost_flags"].Has("WriteConditionsFlag"):
+            gid_parameters["result_file_configuration"]["gidpost_flags"].AddEmptyValue("WriteConditionsFlag")
+
+        if self.write_design_surface:
+            gid_parameters["result_file_configuration"]["gidpost_flags"]["WriteConditionsFlag"].SetString("WriteConditions")
+        elif self.write_optimization_model_part:
+            gid_parameters["result_file_configuration"]["gidpost_flags"]["WriteConditionsFlag"].SetString("WriteElementsOnly")
+
+    # --------------------------------------------------------------------------
+    def __CreateGiDIO( self ):
+
+        gid_config = self.output_settings["output_format"]["gid_parameters"]
+        results_directory = self.output_settings["output_directory"].GetString()
+        design_history_file_path =  results_directory+"/"+self.design_history_filename
+
+        if self.write_design_surface:
+            self.gid_io = GiDOutputProcess(self.design_surface, design_history_file_path, gid_config)
+        elif self.write_optimization_model_part:
+            self.gid_io = GiDOutputProcess(self.optimization_model_part, design_history_file_path, gid_config)
+
+    # --------------------------------------------------------------------------
+    def InitializeLogging( self ):
+        self.gid_io.ExecuteInitialize()
+        self.gid_io.ExecuteBeforeSolutionLoop()
+
+    # --------------------------------------------------------------------------
+    def LogCurrentDesign( self, optimizationIteration ):
+        OriginalTime = self.optimization_model_part.ProcessInfo[KM.TIME]
+        self.optimization_model_part.ProcessInfo[KM.TIME] = optimizationIteration
+
+        self.gid_io.ExecuteInitializeSolutionStep()
+        if(self.gid_io.IsOutputStep()):
+            self.gid_io.PrintOutput()
+        self.gid_io.ExecuteFinalizeSolutionStep()
+
+        self.optimization_model_part.ProcessInfo[KM.TIME] = OriginalTime
+
+    # --------------------------------------------------------------------------
+    def FinalizeLogging( self ):
+        self.gid_io.ExecuteFinalize()
+
+# ==============================================================================

KratosMultiphysics/ShapeOptimizationApplication/loggers/design_logger_unv.py

@@ -0,0 +1,73 @@
+# ==============================================================================
+#  KratosShapeOptimizationApplication
+#
+#  License:         BSD License
+#                   license: ShapeOptimizationApplication/license.txt
+#
+#  Main authors:    Baumgaertner Daniel, https://github.com/dbaumgaertner
+#                   Geiser Armin, https://github.com/armingeiser
+#
+# ==============================================================================
+
+
+# Kratos Core and Apps
+import KratosMultiphysics.ShapeOptimizationApplication as KSO
+
+# Import logger base classes
+from KratosMultiphysics.ShapeOptimizationApplication.loggers.design_logger_base import DesignLogger
+
+# ==============================================================================
+class DesignLoggerUNV( DesignLogger ):
+
+    # --------------------------------------------------------------------------
+    def __init__( self, model_part_controller, optimization_settings ):
+        self.optimization_model_part = model_part_controller.GetOptimizationModelPart()
+        self.design_surface = model_part_controller.GetDesignSurface()
+        self.output_settings = optimization_settings["output"]
+
+        self.write_design_surface = False
+        self.write_optimization_model_part = False
+        self.design_history_filename = None
+
+        self.__DetermineOutputMode()
+        self.__CreateUNVIO()
+
+    # --------------------------------------------------------------------------
+    def __DetermineOutputMode( self ):
+        output_mode = self.output_settings["design_output_mode"].GetString()
+
+        if output_mode == "write_design_surface":
+            self.write_design_surface = True
+            self.design_history_filename = self.design_surface.Name
+        elif output_mode == "write_optimization_model_part":
+            if self.optimization_model_part.NumberOfElements() == 0:
+                raise NameError("Output of optimization model part in UNV-format requires definition of elements. No elements are given in current mdpa! You may change the design output mode.")
+            self.write_optimization_model_part = True
+            self.design_history_filename = self.optimization_model_part.Name
+        else:
+            raise NameError("The following design output mode is not defined within a UNV output (name may be misspelled): " + OutputMode)
+
+    # --------------------------------------------------------------------------
+    def __CreateUNVIO( self ):
+        results_directory = self.output_settings["output_directory"].GetString()
+        design_history_file_path =  results_directory+"/"+self.design_history_filename
+        nodal_results = self.output_settings["nodal_results"]
+
+        if self.write_design_surface:
+            self.UNVIO = KSO.UniversalFileIO( self.design_surface, design_history_file_path, "WriteConditionsOnly", nodal_results )
+        elif self.write_optimization_model_part:
+            self.UNVIO = KSO.UniversalFileIO( self.optimization_model_part, design_history_file_path, "WriteElementsOnly", nodal_results )
+
+    # --------------------------------------------------------------------------
+    def InitializeLogging( self ):
+        self.UNVIO.InitializeLogging()
+
+    # --------------------------------------------------------------------------
+    def LogCurrentDesign( self, OptimizationIteration ):
+        self.UNVIO.LogNodalResults( OptimizationIteration )
+
+    # --------------------------------------------------------------------------
+    def FinalizeLogging( self ):
+        pass
+
+# ==============================================================================

KratosMultiphysics/ShapeOptimizationApplication/loggers/design_logger_vtk.py

@@ -0,0 +1,96 @@
+# ==============================================================================
+#  KratosShapeOptimizationApplication
+#
+#  License:         BSD License
+#                   license: ShapeOptimizationApplication/license.txt
+#
+#  Main authors:    Baumgaertner Daniel, https://github.com/dbaumgaertner
+#                   Geiser Armin, https://github.com/armingeiser
+#
+# ==============================================================================
+
+
+# Kratos Core and Apps
+import KratosMultiphysics as KM
+
+# Additional imports
+from KratosMultiphysics.ShapeOptimizationApplication.loggers.design_logger_base import DesignLogger
+from KratosMultiphysics.vtk_output_process import VtkOutputProcess
+
+# ==============================================================================
+class DesignLoggerVTK( DesignLogger ):
+
+    # --------------------------------------------------------------------------
+    def __init__( self, model_part_controller, optimization_settings ):
+        self.output_settings = optimization_settings["output"]
+        minimal_vtk_settings = KM.Parameters("""
+        {
+            "name"       : "vtk",
+            "vtk_parameters" : {
+                "model_part_name"        : "PLEASE_SPECIFY_MODEL_PART_NAME",
+                "write_ids"              : false,
+                "file_format"            : "binary",
+                "output_sub_model_parts" : false,
+                "output_path"            : "Optimization_Results"
+            }
+        }""")
+
+        output_format = self.output_settings["output_format"]
+        if not output_format.Has("vtk_parameters"):
+            output_format.AddValue("vtk_parameters", minimal_vtk_settings["vtk_parameters"])
+        else:
+            if output_format["vtk_parameters"].Has("model_part_name"):
+                KM.Logger.PrintWarning("ShapeOpt::DesignLoggerVTK", "vtk output parameter `model_part_name` will be overwritten!")
+            if output_format["vtk_parameters"].Has("output_path"):
+                KM.Logger.PrintWarning("ShapeOpt::DesignLoggerVTK", "vtk output parameter `output_path` will be overwritten!")
+
+        output_format["vtk_parameters"].ValidateAndAssignDefaults(minimal_vtk_settings["vtk_parameters"])
+
+        self.model = model_part_controller.GetModel()
+        self.optimization_model_part = model_part_controller.GetOptimizationModelPart()
+        self.design_surface = model_part_controller.GetDesignSurface()
+
+        self.vtk_io = self.__CreateVTKIO()
+
+    # --------------------------------------------------------------------------
+    def __CreateVTKIO( self ):
+
+        vtk_parameters = self.output_settings["output_format"]["vtk_parameters"]
+        output_mode = self.output_settings["design_output_mode"].GetString()
+
+        nodal_results = self.output_settings["nodal_results"]
+        vtk_parameters.AddValue("nodal_solution_step_data_variables", nodal_results)
+
+        if output_mode == "write_design_surface":
+            vtk_parameters["model_part_name"].SetString(self.design_surface.FullName())
+        elif output_mode == "write_optimization_model_part":
+            vtk_parameters["model_part_name"].SetString(self.optimization_model_part.FullName())
+        else:
+            raise NameError("The following design output mode is not defined within a VTK output (name may be misspelled): " + output_mode)
+
+        vtk_parameters["output_path"].SetString(self.output_settings["output_directory"].GetString())
+
+        return VtkOutputProcess(self.model, vtk_parameters)
+
+    # --------------------------------------------------------------------------
+    def InitializeLogging( self ):
+        self.vtk_io.ExecuteInitialize()
+        self.vtk_io.ExecuteBeforeSolutionLoop()
+
+    # --------------------------------------------------------------------------
+    def LogCurrentDesign( self, optimizationIteration ):
+        OriginalTime = self.optimization_model_part.ProcessInfo[KM.TIME]
+        self.optimization_model_part.ProcessInfo[KM.TIME] = optimizationIteration
+
+        self.vtk_io.ExecuteInitializeSolutionStep()
+        if(self.vtk_io.IsOutputStep()):
+            self.vtk_io.PrintOutput()
+        self.vtk_io.ExecuteFinalizeSolutionStep()
+
+        self.optimization_model_part.ProcessInfo[KM.TIME] = OriginalTime
+
+    # --------------------------------------------------------------------------
+    def FinalizeLogging( self ):
+        self.vtk_io.ExecuteFinalize()
+
+# ==============================================================================

KratosMultiphysics/ShapeOptimizationApplication/loggers/sensitivity_heatmap_logger.py

@@ -0,0 +1,91 @@
+# ==============================================================================
+#  KratosShapeOptimizationApplication
+#
+#  License:         BSD License
+#                   license: ShapeOptimizationApplication/license.txt
+#
+#  Main authors:    Schmölz David, https://github.com/dschmoelz
+#
+# ==============================================================================
+
+# importing the Kratos Library
+import KratosMultiphysics as Kratos
+from KratosMultiphysics.ShapeOptimizationApplication.utilities.custom_sensitivity_heatmap import SensitivityHeatmapCalculator
+
+# ==============================================================================
+class SensitivityHeatmapLoggerBase:
+
+    # --------------------------------------------------------------------------
+    def __init__( self, model_part_controller, optimization_settings ):
+
+        default_heatmap_settings = Kratos.Parameters("""
+        {
+            "norm_type": "l2",
+            "sensitivity_weighting": true,
+            "mapping" : true,
+            "relaxation_method": "reciprocal",
+            "relaxation_coefficient": 0.5
+        }""")
+
+        self.heatmap_settings = optimization_settings["output"]["sensitivity_heatmap_settings"]
+        if self.heatmap_settings.Has("relaxation_method") and self.heatmap_settings["relaxation_method"].GetString() == "constant":
+            if not self.heatmap_settings.Has("relaxation_coefficient"):
+                raise RuntimeError("'relaxation_coefficient' is missing for relaxation method 'constant'!")
+        self.heatmap_settings.ValidateAndAssignDefaults(default_heatmap_settings)
+
+        self.optimization_model_part = model_part_controller.GetOptimizationModelPart()
+        self.design_surface = model_part_controller.GetDesignSurface()
+
+        self.objectives = optimization_settings["objectives"]
+        self.constraints = optimization_settings["constraints"]
+
+        self.sensitivity_heatmap_calculator = None
+
+    # --------------------------------------------------------------------------
+    def InitializeLogging ( self ):
+        raise RuntimeError("Sensitivity Heatmap logger base class is called. Please check your implementation of the function >> InitializeLogging << .")
+
+    # --------------------------------------------------------------------------
+    def LogSensitivityHeatmap(self, optimization_iteration, mapper):
+        self.sensitivity_heatmap_calculator.ComputeHeatmaps(optimization_iteration, mapper)
+
+# ==============================================================================
+class SensitivityHeatmapLoggerSteepestDescent(SensitivityHeatmapLoggerBase):
+
+    def InitializeLogging(self):
+        self.heatmap_settings.AddString("design_variable_name", "X")
+        self.heatmap_settings.AddInt("design_variable_dimension", 3)
+        self.sensitivity_heatmap_calculator = SensitivityHeatmapCalculator(self.design_surface, self.objectives, self.constraints, self.heatmap_settings)
+
+
+# ==============================================================================
+class SensitivityHeatmapLoggerPenalizedProjection(SensitivityHeatmapLoggerBase):
+
+    def InitializeLogging(self):
+        self.heatmap_settings.AddString("design_variable_name", "X")
+        self.heatmap_settings.AddInt("design_variable_dimension", 3)
+        self.sensitivity_heatmap_calculator = SensitivityHeatmapCalculator(self.design_surface, self.objectives, self.constraints, self.heatmap_settings)
+
+# ==============================================================================
+class SensitivityHeatmapLoggerGradientProjection(SensitivityHeatmapLoggerBase):
+
+    def InitializeLogging(self):
+        self.heatmap_settings.AddString("design_variable_name", "X")
+        self.heatmap_settings.AddInt("design_variable_dimension", 3)
+        self.sensitivity_heatmap_calculator = SensitivityHeatmapCalculator(self.design_surface, self.objectives, self.constraints, self.heatmap_settings)
+
+# ==============================================================================
+class SensitivityHeatmapLoggerTrustRegion(SensitivityHeatmapLoggerBase):
+
+    def InitializeLogging(self):
+        self.heatmap_settings.AddString("design_variable_name", "X")
+        self.heatmap_settings.AddInt("design_variable_dimension", 3)
+        self.sensitivity_heatmap_calculator = SensitivityHeatmapCalculator(self.design_surface, self.objectives, self.constraints, self.heatmap_settings)
+
+# ==============================================================================
+class SensitivityHeatmapLoggerBeadOptimization(SensitivityHeatmapLoggerBase):
+
+    def InitializeLogging(self):
+        self.heatmap_settings.AddString("design_variable_name", "ALPHA")
+        self.heatmap_settings.AddInt("design_variable_dimension", 1)
+        self.sensitivity_heatmap_calculator = SensitivityHeatmapCalculator(self.design_surface, self.objectives, self.constraints, self.heatmap_settings)

KratosMultiphysics/ShapeOptimizationApplication/loggers/value_logger_base.py

@@ -0,0 +1,128 @@
+# ==============================================================================
+#  KratosShapeOptimizationApplication
+#
+#  License:         BSD License
+#                   license: ShapeOptimizationApplication/license.txt
+#
+#  Main authors:    Baumgaertner Daniel, https://github.com/dbaumgaertner
+#                   Geiser Armin, https://github.com/armingeiser
+#
+# ==============================================================================
+
+# importing the Kratos Library
+import KratosMultiphysics as KM
+
+# Import additional libraries
+import os
+
+# ==============================================================================
+class ValueLogger():
+
+    # --------------------------------------------------------------------------
+    def __init__( self, communicator, optimization_settings ):
+        self.communicator = communicator
+        self.optimization_settings = optimization_settings
+
+        self.objectives = optimization_settings["objectives"]
+        self.constraints = optimization_settings["constraints"]
+
+        self.complete_log_file_name = self.__CreateCompleteLogFileName( optimization_settings )
+
+        self.obj_reference_value = 0
+        self.current_index = 0
+        self.previos_index = 0
+
+        self.history = { "response_value"              : {},
+                         "standardized_response_value" : {},
+                         "abs_change_objective"        : {},
+                         "rel_change_objective"        : {} }
+
+        for itr in range(self.objectives.size()):
+            objective_id = self.objectives[itr]["identifier"].GetString()
+            self.history["response_value"][objective_id] = {}
+            self.history["standardized_response_value"][objective_id] = {}
+
+        for itr in range(self.constraints.size()):
+            constraint_id = self.constraints[itr]["identifier"].GetString()
+            self.history["response_value"][constraint_id] = {}
+            self.history["standardized_response_value"][constraint_id] = {}
+
+        self.predefined_keys = list(self.history.keys())
+
+    # --------------------------------------------------------------------------
+    def InitializeLogging( self ):
+        pass
+
+    # --------------------------------------------------------------------------
+    def LogCurrentValues( self, index, additional_values_dictionary ):
+        self.current_index = index
+
+        self.__LogObjectiveValuesToHistory()
+        self.__LogConstraintValuesToHistory()
+        self.__LogAdditionalValuesToHistory( additional_values_dictionary )
+
+        self._WriteCurrentValuesToConsole()
+        self._WriteCurrentValuesToFile()
+
+        self.previos_index = index
+
+    # --------------------------------------------------------------------------
+    def FinalizeLogging( self ):
+        pass
+
+    # --------------------------------------------------------------------------
+    def GetValues( self, key ):
+        return self.history[key]
+
+    # --------------------------------------------------------------------------
+    @staticmethod
+    def __CreateCompleteLogFileName( optimization_settings ):
+        resultsDirectory = optimization_settings["output"]["output_directory"].GetString()
+        optimizationLogFilename = optimization_settings["output"]["optimization_log_filename"].GetString() + ".csv"
+        return os.path.join( resultsDirectory, optimizationLogFilename )
+
+    # --------------------------------------------------------------------------
+    def __LogObjectiveValuesToHistory( self ):
+        objective_id = self.objectives[0]["identifier"].GetString()
+        is_first_log = self.history["response_value"][objective_id] == {}
+
+        self.history["response_value"][objective_id][self.current_index] = self.communicator.getValue( objective_id )
+        self.history["standardized_response_value"][objective_id][self.current_index] = self.communicator.getStandardizedValue( objective_id )
+
+        if is_first_log:
+            self.obj_reference_value = self.communicator.getValue( objective_id )
+
+            if abs(self.obj_reference_value)<1e-12:
+                KM.Logger.PrintWarning("ShapeOpt::ValueLogger", "Objective reference value < 1e-12!! Therefore, standard reference value of 1 is assumed! ")
+                self.obj_reference_value = 1.0
+
+            self.history["abs_change_objective"] = {self.current_index: 0.0}
+            self.history["rel_change_objective"] = {self.current_index: 0.0}
+        else:
+            current_obj_value = self.history["response_value"][objective_id][self.current_index]
+            previos_obj_value = self.history["response_value"][objective_id][self.previos_index]
+            abs_change = 100*(current_obj_value-self.obj_reference_value) / abs(self.obj_reference_value)
+            rel_change = 100*(current_obj_value-previos_obj_value) / abs(self.obj_reference_value)
+
+            self.history["abs_change_objective"][self.current_index] = abs_change
+            self.history["rel_change_objective"][self.current_index] = rel_change
+
+    # --------------------------------------------------------------------------
+    def __LogConstraintValuesToHistory( self ):
+        for itr in range(self.constraints.size()):
+            constraint_id = self.constraints[itr]["identifier"].GetString()
+            self.history["response_value"][constraint_id][self.current_index] = self.communicator.getValue( constraint_id )
+            self.history["standardized_response_value"][constraint_id][self.current_index] = self.communicator.getStandardizedValue( constraint_id )
+
+    # --------------------------------------------------------------------------
+    def __LogAdditionalValuesToHistory( self, additional_values_dictionary ):
+        for key, value in additional_values_dictionary.items():
+            if key in self.predefined_keys:
+                raise NameError("ValueLogger: The key \""+key+"\" is already predefined and may not be overwritten! Predefined keys are: "+str(self.predefined_keys))
+
+            if key not in self.history.keys():
+                self.history[key] = {}
+
+            self.history[key][self.current_index] = value
+
+# ==============================================================================

KratosMultiphysics/ShapeOptimizationApplication/loggers/value_logger_bead_optimization.py

@@ -0,0 +1,79 @@
+# ==============================================================================
+#  KratosShapeOptimizationApplication
+#
+#  License:         BSD License
+#                   license: ShapeOptimizationApplication/license.txt
+#
+#  Main authors:    Baumgaertner Daniel, https://github.com/dbaumgaertner
+#                   Geiser Armin, https://github.com/armingeiser
+#
+# ==============================================================================
+
+# importing the Kratos Library
+import KratosMultiphysics as KM
+
+# Import logger base classes
+from KratosMultiphysics.ShapeOptimizationApplication.loggers.value_logger_base import ValueLogger
+
+# Import additional libraries
+import csv
+from KratosMultiphysics.ShapeOptimizationApplication.utilities.custom_timer import Timer
+
+# ==============================================================================
+class ValueLoggerBeadOptimization( ValueLogger ):
+    # --------------------------------------------------------------------------
+    def InitializeLogging( self ):
+        with open(self.complete_log_file_name, 'w') as csvfile:
+            historyWriter = csv.writer(csvfile, delimiter=',',quotechar='|',quoting=csv.QUOTE_MINIMAL)
+            row = []
+            row.append("{:>12s}".format("overall_itr"))
+            row.append("{:>10s}".format("outer_itr"))
+            row.append("{:>10s}".format("inner_itr"))
+            row.append("{:>13s}".format("l"))
+            row.append("{:>13s}".format("dl_rel[%]"))
+            row.append("{:>13s}".format("f"))
+            row.append("{:>13s}".format("df_abs[%]"))
+            row.append("{:>13s}".format("df_rel[%]"))
+            row.append("{:>14s}".format("max_norm_dfda"))
+            row.append("{:>13s}".format("lambda_p"))
+            row.append("{:>13s}".format("p"))
+            row.append("{:>13s}".format("p_scaling"))
+            row.append("{:>13s}".format("p_fac"))
+            row.append("{:>13s}".format("step_size"))
+            row.append("{:>25s}".format("time_stamp"))
+            historyWriter.writerow(row)
+
+    # --------------------------------------------------------------------------
+    def _WriteCurrentValuesToConsole( self ):
+        objective_id = self.objectives[0]["identifier"].GetString()
+        KM.Logger.Print("")
+        KM.Logger.PrintInfo("ShapeOpt", "Current value of objective = ","{:> .5E}".format(self.history["response_value"][objective_id][self.current_index]))
+
+        KM.Logger.PrintInfo("ShapeOpt", "Absolut change of objective = ","{:> .5E}".format(self.history["abs_change_objective"][self.current_index])," [%]")
+        KM.Logger.PrintInfo("ShapeOpt", "Relative change of objective = ","{:> .5E}".format(self.history["rel_change_objective"][self.current_index])," [%]\n")
+        KM.Logger.PrintInfo("ShapeOpt", "Value of penalty term = ","{:> .5E}".format(self.history["penalty_value"][self.current_index]),"\n")
+
+    # --------------------------------------------------------------------------
+    def _WriteCurrentValuesToFile( self ):
+        with open(self.complete_log_file_name, 'a') as csvfile:
+            historyWriter = csv.writer(csvfile, delimiter=',',quotechar='|',quoting=csv.QUOTE_MINIMAL)
+            row = []
+            row.append("{:>12d}".format(self.current_index))
+            row.append("{:>10d}".format(self.history["outer_iteration"][self.current_index]))
+            row.append("{:>10d}".format(self.history["inner_iteration"][self.current_index]))
+            objective_id = self.objectives[0]["identifier"].GetString()
+            row.append(" {:> .5E}".format(self.history["lagrange_value"][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["lagrange_value_relative_change"][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["response_value"][objective_id][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["abs_change_objective"][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["rel_change_objective"][self.current_index]))
+            row.append("  {:> .5E}".format(self.history["max_norm_objective_gradient"][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["penalty_lambda"][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["penalty_value"][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["penalty_scaling"][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["penalty_factor"][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["step_size"][self.current_index]))
+            row.append("{:>25}".format(Timer().GetTimeStamp()))
+            historyWriter.writerow(row)
+
+# ==============================================================================