KratosShapeOptimizationApplication 10.4.4__py3-none-any.whl

KratosMultiphysics/ShapeOptimizationApplication/loggers/value_logger_gradient_projection.py

@@ -0,0 +1,80 @@
+# ==============================================================================
+#  KratosShapeOptimizationApplication
+#
+#  License:         BSD License
+#                   license: ShapeOptimizationApplication/license.txt
+#
+#  Main authors:    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 ValueLoggerGradientProjection( 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("{:>4s}".format("itr"))
+            row.append("{:>13s}".format("f"))
+            row.append("{:>13s}".format("df_abs[%]"))
+            row.append("{:>13s}".format("df_rel[%]"))
+
+            for itr in range(self.constraints.size()):
+                con_type = self.constraints[itr]["type"].GetString()
+                row.append("{:>13s}".format("c"+str(itr+1)+": "+con_type))
+                row.append("{:>13s}".format("c"+str(itr+1)+"_ref"))
+
+            row.append("{:>13s}".format("step_size"))
+            row.append("{:>13s}".format("inf_norm_s"))
+            row.append("{:>13s}".format("inf_norm_c"))
+            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", "Absolute 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")
+
+        for itr in range(self.constraints.size()):
+            constraint_id = self.constraints[itr]["identifier"].GetString()
+            KM.Logger.PrintInfo("ShapeOpt", "Value of C"+str(itr+1)+" = ", "{:> .5E}".format(self.history["response_value"][constraint_id][self.current_index]))
+
+    # --------------------------------------------------------------------------
+    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("{:>4d}".format(self.current_index))
+
+            objective_id = self.objectives[0]["identifier"].GetString()
+            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]))
+
+            for itr in range(self.constraints.size()):
+                constraint_id = self.constraints[itr]["identifier"].GetString()
+                row.append(" {:> .5E}".format(self.history["response_value"][constraint_id][self.current_index]))
+                row.append(" {:> .5E}".format(self.communicator.getReferenceValue(constraint_id)))
+
+            row.append(" {:> .5E}".format(self.history["step_size"][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["inf_norm_s"][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["inf_norm_c"][self.current_index]))
+            row.append("{:>25}".format(Timer().GetTimeStamp()))
+            historyWriter.writerow(row)
+
+# ==============================================================================

KratosMultiphysics/ShapeOptimizationApplication/loggers/value_logger_penalized_projection.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 ValueLoggerPenalizedProjection( 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("{:>4s}".format("itr"))
+            row.append("{:>13s}".format("f"))
+            row.append("{:>13s}".format("df_abs[%]"))
+            row.append("{:>13s}".format("df_rel[%]"))
+
+            for itr in range(self.constraints.size()):
+                con_type = self.constraints[itr]["type"].GetString()
+                row.append("{:>13s}".format("c"+str(itr+1)+": "+con_type))
+                row.append("{:>13s}".format("c"+str(itr+1)+"_ref"))
+
+            row.append("{:>13s}".format("c_scaling"))
+            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")
+
+        for itr in range(self.constraints.size()):
+            constraint_id = self.constraints[itr]["identifier"].GetString()
+            KM.Logger.PrintInfo("ShapeOpt", "Value of C"+str(itr+1)+" = ", "{:> .5E}".format(self.history["response_value"][constraint_id][self.current_index]))
+
+    # --------------------------------------------------------------------------
+    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("{:>4d}".format(self.current_index))
+
+            objective_id = self.objectives[0]["identifier"].GetString()
+            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]))
+
+            for itr in range(self.constraints.size()):
+                constraint_id = self.constraints[itr]["identifier"].GetString()
+                row.append(" {:> .5E}".format(self.history["response_value"][constraint_id][self.current_index]))
+                row.append(" {:> .5E}".format(self.communicator.getReferenceValue(constraint_id)))
+
+            row.append(" {:> .5E}".format(self.history["correction_scaling"][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["step_size"][self.current_index]))
+            row.append("{:>25}".format(Timer().GetTimeStamp()))
+            historyWriter.writerow(row)
+
+# ==============================================================================

KratosMultiphysics/ShapeOptimizationApplication/loggers/value_logger_relaxed_gradient_projection.py

@@ -0,0 +1,108 @@
+# ==============================================================================
+#  KratosShapeOptimizationApplication
+#
+#  License:         BSD License
+#                   license: ShapeOptimizationApplication/license.txt
+#
+#  Main authors:    Ihar Antonau
+#
+# ==============================================================================
+
+# importing the Kratos Library
+import KratosMultiphysics as Kratos
+
+# Import logger base classes
+from .value_logger_base import ValueLogger
+
+# Import additional libraries
+import csv
+from KratosMultiphysics.ShapeOptimizationApplication.utilities.custom_timer import Timer
+
+# ==============================================================================
+class ValueLoggerRelaxedGradientProjection( 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("{:>4s}".format("itr"))
+            row.append("{:>13s}".format("f"))
+            row.append("{:>13s}".format("df_abs[%]"))
+            row.append("{:>13s}".format("df_rel[%]"))
+
+            for itr in range(self.constraints.size()):
+                con_type = self.constraints[itr]["type"].GetString()
+                row.append("{:>13s}".format("c"+str(itr+1)+": "+con_type))
+                row.append("{:>13s}".format("c"+str(itr+1)+"_ref"))
+
+            row.append("{:>13s}".format("step_size"))
+            row.append("{:>13s}".format("inf_norm_p"))
+            row.append("{:>13s}".format("inf_norm_c"))
+            row.append("{:>13s}".format("projection_norm"))
+
+            for itr in range(self.constraints.size()):
+                row.append("{:>13s}".format("c"+str(itr+1)+"_buffer_value"))
+                row.append("{:>13s}".format("c"+str(itr+1)+"_buffer_size"))
+                row.append("{:>13s}".format("c"+str(itr+1)+"_buffer_size_factor"))
+                row.append("{:>13s}".format("c"+str(itr+1)+"_central_buffer_value"))
+                row.append("{:>13s}".format("c"+str(itr+1)+"_lower_buffer_value"))
+                row.append("{:>13s}".format("c"+str(itr+1)+"_upper_buffer_value"))
+
+            row.append("{:>25s}".format("time_stamp"))
+            historyWriter.writerow(row)
+
+    # --------------------------------------------------------------------------
+    def _WriteCurrentValuesToConsole( self ):
+        objective_id = self.objectives[0]["identifier"].GetString()
+        Kratos.Logger.Print("")
+        Kratos.Logger.PrintInfo("ShapeOpt", "Current value of objective = ", "{:> .5E}".format(self.history["response_value"]
+                                                                                               [objective_id][self.current_index]))
+
+        Kratos.Logger.PrintInfo("ShapeOpt", "Absolute change of objective = ","{:> .5E}".format(self.history["abs_change_objective"][self.current_index]),
+                                                                                                " [%]")
+        Kratos.Logger.PrintInfo("ShapeOpt", "Relative change of objective = ","{:> .5E}".format(self.history["rel_change_objective"][self.current_index]),
+                                                                                                " [%]\n")
+
+        for itr in range(self.constraints.size()):
+            constraint_id = self.constraints[itr]["identifier"].GetString()
+            Kratos.Logger.PrintInfo("ShapeOpt", f"Value of C{str(itr+1)} = ","{:> .5E}".format(
+                                    self.history["response_value"][constraint_id][self.current_index]))
+
+    # --------------------------------------------------------------------------
+    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("{:>4d}".format(self.current_index))
+
+            objective_id = self.objectives[0]["identifier"].GetString()
+            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]))
+
+            for itr in range(self.constraints.size()):
+                constraint_id = self.constraints[itr]["identifier"].GetString()
+                row.append(" {:> .5E}".format(self.history["response_value"][constraint_id][self.current_index]))
+                row.append(" {:> .5E}".format(self.communicator.getReferenceValue(constraint_id)))
+
+            row.append(" {:> .5E}".format(self.history["step_size"][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["inf_norm_p"][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["inf_norm_c"][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["projection_norm"][self.current_index]))
+
+            for itr in range(self.constraints.size()):
+                constraint_id = self.constraints[itr]["identifier"].GetString()
+                row.append(" {:> .5E}".format(self.history["c"+str(itr+1)+"_buffer_value"][self.current_index]))
+                row.append(" {:> .5E}".format(self.history["c"+str(itr+1)+"_buffer_size"][self.current_index]))
+                row.append(" {:> .5E}".format(self.history["c"+str(itr+1)+"_buffer_size_factor"][self.current_index]))
+                row.append(" {:> .5E}".format(self.history["c"+str(itr+1)+"_central_buffer_value"][self.current_index] +
+                                              self.communicator.getReferenceValue(constraint_id)))
+                row.append(" {:> .5E}".format(self.history["c"+str(itr+1)+"_lower_buffer_value"][self.current_index]
+                                              + self.communicator.getReferenceValue(constraint_id)))
+                row.append(" {:> .5E}".format(self.history["c"+str(itr+1)+"_upper_buffer_value"][self.current_index]
+                                              + self.communicator.getReferenceValue(constraint_id)))
+
+            row.append("{:>25}".format(Timer().GetTimeStamp()))
+            historyWriter.writerow(row)
+
+# ==============================================================================

KratosMultiphysics/ShapeOptimizationApplication/loggers/value_logger_shape_fraction_optimization.py

@@ -0,0 +1,92 @@
+# ==============================================================================
+#  KratosShapeOptimizationApplication
+#
+#  License:         BSD License
+#                   license: ShapeOptimizationApplication/license.txt
+#
+#  Main authors:    Schmölz David, https://github.com/dschmoelz
+#
+# ==============================================================================
+
+# 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 ValueLoggerShapeFractionOptimization( 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("{:>4s}".format("itr"))
+            row.append("{:>13s}".format("f"))
+            row.append("{:>13s}".format("df_abs[%]"))
+            row.append("{:>13s}".format("df_rel[%]"))
+            row.append("{:>13s}".format("penalty_value"))
+            row.append("{:>13s}".format("penalty_factor"))
+            row.append("{:>13s}".format("f_p"))
+            row.append("{:>13s}".format("df_rel_p"))
+            row.append("{:>13s}".format("shape_fraction"))
+
+            for itr in range(self.constraints.size()):
+                con_type = self.constraints[itr]["type"].GetString()
+                row.append("{:>13s}".format("c"+str(itr+1)+": "+con_type))
+                row.append("{:>13s}".format("c"+str(itr+1)+"_ref"))
+
+            row.append("{:>13s}".format("step_size"))
+            row.append("{:>13s}".format("inf_norm_s"))
+            row.append("{:>13s}".format("inf_norm_c"))
+            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", "Absolute 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")
+
+        for itr in range(self.constraints.size()):
+            constraint_id = self.constraints[itr]["identifier"].GetString()
+            KM.Logger.PrintInfo("ShapeOpt", "Value of C"+str(itr+1)+" = ", "{:> .5E}".format(self.history["response_value"][constraint_id][self.current_index]))
+
+    # --------------------------------------------------------------------------
+    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("{:>4d}".format(self.current_index))
+
+            objective_id = self.objectives[0]["identifier"].GetString()
+            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["penalty_value"][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["penalty_factor"][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["f_p"][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["df_rel_p"][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["shape_fraction"][self.current_index]))
+
+
+            for itr in range(self.constraints.size()):
+                constraint_id = self.constraints[itr]["identifier"].GetString()
+                row.append(" {:> .5E}".format(self.history["response_value"][constraint_id][self.current_index]))
+                row.append(" {:> .5E}".format(self.communicator.getReferenceValue(constraint_id)))
+
+            row.append(" {:> .5E}".format(self.history["step_size"][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["inf_norm_s"][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["inf_norm_c"][self.current_index]))
+            row.append("{:>25}".format(Timer().GetTimeStamp()))
+            historyWriter.writerow(row)
+
+# ==============================================================================

KratosMultiphysics/ShapeOptimizationApplication/loggers/value_logger_steepest_descent.py

@@ -0,0 +1,64 @@
+# ==============================================================================
+#  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 ValueLoggerSteepestDescent( 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("{:>4s}".format("itr"))
+            row.append("{:>13s}".format("f"))
+            row.append("{:>13s}".format("df_abs[%]"))
+            row.append("{:>13s}".format("df_rel[%]"))
+            row.append("{:>13s}".format("norm_df"))
+            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")
+
+    # --------------------------------------------------------------------------
+    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("{:>4d}".format(self.current_index))
+
+            objective_id = self.objectives[0]["identifier"].GetString()
+            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["norm_objective_gradient"][self.current_index]))
+            row.append(" {:> .5E}".format(self.history["step_size"][self.current_index]))
+            row.append("{:>25}".format(Timer().GetTimeStamp()))
+            historyWriter.writerow(row)
+
+# ==============================================================================

KratosMultiphysics/ShapeOptimizationApplication/loggers/value_logger_trust_region.py

@@ -0,0 +1,100 @@
+# ==============================================================================
+#  KratosShapeOptimizationApplication
+#
+#  License:         BSD License
+#                   license: ShapeOptimizationApplication/license.txt
+#
+#  Main authors:    Baumgaertner Daniel, https://github.com/dbaumgaertner
+#
+# ==============================================================================
+
+# 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 ValueLoggerTrustRegion( 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("{:>4s}".format("itr"))
+            row.append("{:>13}".format("f"))
+            row.append("{:>12s}".format("df_abs[%]"))
+            row.append("{:>12s}".format("df_rel[%]"))
+
+            for itr in range(self.constraints.size()):
+                con_type = self.constraints[itr]["type"].GetString()
+                row.append("{:>13}".format("c"+str(itr+1)+": "+con_type))
+                row.append("{:>13}".format("c"+str(itr+1)+"_ref"))
+                row.append("{:>12s}".format("len_C"+str(itr+1)))
+                row.append("{:>12s}".format("adj_len_C"+str(itr+1)))
+
+            row.append("{:>12s}".format("bi_itrs"))
+            row.append("{:>12s}".format("bi_err"))
+            row.append("{:>17}".format("test_norm_dX_bar"))
+            row.append("{:>12s}".format("norm_dX"))
+            row.append("{:>12s}".format("step_length"))
+            row.append("{:>25s}".format("time_stamp"))
+            historyWriter.writerow(row)
+
+    # --------------------------------------------------------------------------
+    def _WriteCurrentValuesToConsole( self ):
+        KM.Logger.Print("")
+        KM.Logger.Print("-------------------------------------------------------\n")
+
+        objective_id = self.objectives[0]["identifier"].GetString()
+        KM.Logger.PrintInfo("ShapeOpt", "Current value of objective = ", round(self.history["response_value"][objective_id][self.current_index],12))
+
+        KM.Logger.PrintInfo("ShapeOpt", "Absolut change of objective = ",round(self.history["abs_change_objective"][self.current_index],4)," [%]")
+        KM.Logger.PrintInfo("ShapeOpt", "Relative change of objective = ",round(self.history["rel_change_objective"][self.current_index],4)," [%]\n")
+
+        for itr in range(self.constraints.size()):
+            constraint_id = self.constraints[itr]["identifier"].GetString()
+            KM.Logger.PrintInfo("ShapeOpt", "Value of C"+str(itr+1)+" = ", round(self.history["response_value"][constraint_id][self.current_index],12), "\n")
+
+        KM.Logger.PrintInfo("ShapeOpt", "NormInf3D of dX = ", round(self.history["norm_dX"][self.current_index],6), "\n")
+
+        KM.Logger.PrintInfo("ShapeOpt", "len_bar_obj = ", round(self.history["len_bar_obj"][self.current_index],6))
+        KM.Logger.PrintInfo("ShapeOpt", "adj_len_bar_obj = ", round(self.history["adj_len_bar_obj"][self.current_index],6), "\n")
+
+        KM.Logger.PrintInfo("ShapeOpt", "len_bar_cons = ", [round(entry, 6) for entry in self.history["len_bar_cons"][self.current_index]])
+        KM.Logger.PrintInfo("ShapeOpt", "adj_len_bar_cons = ", [round(entry, 6) for entry in self.history["adj_len_bar_cons"][self.current_index]], "\n")
+
+        KM.Logger.Print("-------------------------------------------------------")
+
+    # --------------------------------------------------------------------------
+    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("{:>4d}".format(self.current_index))
+
+            objective_id = self.objectives[0]["identifier"].GetString()
+            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]))
+
+            for itr in range(self.constraints.size()):
+                constraint_id = self.constraints[itr]["identifier"].GetString()
+                row.append(" {:> .5E}".format(self.history["response_value"][constraint_id][self.current_index]))
+                row.append(" {:> .5E}".format(self.communicator.getReferenceValue(constraint_id)))
+                row.append("{:>12f}".format(self.history["len_bar_cons"][self.current_index][itr]))
+                row.append("{:>12f}".format(self.history["adj_len_bar_cons"][self.current_index][itr]))
+
+            row.append("{:>12d}".format(self.history["bi_itrs"][self.current_index]))
+            row.append("{:>12f}".format(self.history["bi_err"][self.current_index]))
+            row.append("{:>17f}".format(self.history["test_norm_dX_bar"][self.current_index]))
+            row.append("{:>12f}".format(self.history["norm_dX"][self.current_index]))
+            row.append(" {:>.5E}".format(self.history["step_length"][self.current_index]))
+            row.append("{:>25}".format(Timer().GetTimeStamp()))
+            historyWriter.writerow(row)
+
+# ==============================================================================

KratosMultiphysics/ShapeOptimizationApplication/mapper_factory.py

@@ -0,0 +1,100 @@
+# ==============================================================================
+#  KratosShapeOptimizationApplication
+#
+#  License:         BSD License
+#                   license: ShapeOptimizationApplication/license.txt
+#
+#  Main authors:    Baumgaertner Daniel, https://github.com/dbaumgaertner
+#
+# ==============================================================================
+
+
+# Kratos Core and Apps
+import KratosMultiphysics as KM
+import KratosMultiphysics.ShapeOptimizationApplication as KSO
+from .mapping import in_plane_vertex_morphing_mapper, sliding_vertex_morphing_mapper
+
+# ==============================================================================
+def CreateMapper(origin_model_part, destination_model_part, mapper_settings):
+    default_settings = KM.Parameters("""
+    {
+        "filter_function_type"       : "linear",
+        "filter_radius"              : 0.000000000001,
+        "max_nodes_in_filter_radius" : 10000,
+        "matrix_free_filtering"      : false,
+        "consistent_mapping"         : false,
+        "improved_integration"       : false,
+        "integration_method"         : "gauss_integration",
+        "number_of_gauss_points"     : 5,
+        "in_plane_morphing"          : false,
+        "in_plane_morphing_settings" : {},
+        "sliding_morphing"           : false,
+        "sliding_morphing_settings"  : {},
+        "plane_symmetry"             : false,
+        "plane_symmetry_settings"    : {
+            "point" : [0.0, 0.0, 0.0],
+            "normal": [0.0, 0.0, 0.0]
+        },
+        "revolution"                 : false,
+        "revolution_settings"        : {
+            "point" : [0.0, 0.0, 0.0],
+            "normal": [0.0, 0.0, 0.0]
+        },
+        "adaptive_filter_settings"   : {
+            "adaptive_filter_method": "curvature_based",
+            "radius_function": "analytic",
+            "radius_function_parameter": 1,
+            "minimum_filter_radius": 1e-3,
+            "curvature_limit": 1e-3,
+            "filter_radius_smoothing_iterations": 10
+        }
+    }""")
+
+    mapper_vertex_morphing_matrix_free = KSO.MapperVertexMorphingMatrixFree
+    mapper_vertex_morphing_improved_integration = KSO.MapperVertexMorphingImprovedIntegration
+    mapper_vertex_morphing_symmetric = KSO.MapperVertexMorphingSymmetric
+    mapper_vertex_morphing = KSO.MapperVertexMorphing
+    if mapper_settings.Has("filter_radius") and mapper_settings["filter_radius"].IsString():
+        if mapper_settings["filter_radius"].GetString() == "adaptive":
+            if mapper_settings.Has("adaptive_filter_method") and mapper_settings["adaptive_filter_method"].GetString() != "curvature_based":
+                raise Exception("For now, only \"curvature_based\" is available for \"adaptive_filter_method\".")
+            else:
+                mapper_vertex_morphing_matrix_free = KSO.MapperVertexMorphingMatrixFreeAdaptiveRadius
+                mapper_vertex_morphing_improved_integration = KSO.MapperVertexMorphingImprovedIntegrationAdaptiveRadius
+                mapper_vertex_morphing_symmetric = KSO.MapperVertexMorphingSymmetricAdaptiveRadius
+                mapper_vertex_morphing = KSO.MapperVertexMorphingAdaptiveRadius
+
+                if mapper_settings.Has("in_plane_morphing"):
+                    if mapper_settings["in_plane_morphing"].GetBool():
+                        raise Exception("\"in_plane_morphing\" is not yet supported with \"adaptive\" filter radius.")
+                mapper_settings["filter_radius"].SetDouble(-1.0)
+        else:
+            raise Exception("\"filter_radius\" either should be double value or \"adaptive\".")
+
+    mapper_settings.ValidateAndAssignDefaults(default_settings)
+
+    if mapper_settings["in_plane_morphing"].GetBool():
+        return in_plane_vertex_morphing_mapper.InPlaneVertexMorphingMapper(origin_model_part, destination_model_part, mapper_settings)
+    elif mapper_settings["sliding_morphing"].GetBool():
+        return sliding_vertex_morphing_mapper.SlidingVertexMorphingMapper(origin_model_part, destination_model_part, mapper_settings)
+    elif mapper_settings["matrix_free_filtering"].GetBool():
+        if mapper_settings["consistent_mapping"].GetBool():
+             raise ValueError ("Matrix free mapper has no option to map consistently yet!")
+        if mapper_settings["improved_integration"].GetBool():
+             raise ValueError ("Matrix free mapper does not yet allow for an improved integration!")
+        else:
+            return mapper_vertex_morphing_matrix_free(origin_model_part, destination_model_part, mapper_settings)
+    else:
+        if mapper_settings["revolution"].GetBool() and mapper_settings["plane_symmetry"].GetBool():
+            raise RuntimeError("revolution and plane_symmetry can not be combined!")
+
+        if mapper_settings["improved_integration"].GetBool():
+            return mapper_vertex_morphing_improved_integration(origin_model_part, destination_model_part, mapper_settings)
+        elif mapper_settings["plane_symmetry"].GetBool():
+            return mapper_vertex_morphing_symmetric(origin_model_part, destination_model_part, mapper_settings)
+        elif mapper_settings["revolution"].GetBool():
+            return mapper_vertex_morphing_symmetric(origin_model_part, destination_model_part, mapper_settings)
+        else:
+            return mapper_vertex_morphing(origin_model_part, destination_model_part, mapper_settings)
+
+# ==============================================================================