def CreateMetricProcess(main_model_part, adaptative_remesh_parameters):
if main_model_part.ProcessInfo[KratosMultiphysics.DOMAIN_SIZE] == 2:
metric_process = MeshingApplication.MetricErrorProcess2D(main_model_part, adaptative_remesh_parameters["metric_error_parameters"])
else:
metric_process = MeshingApplication.MetricErrorProcess3D(main_model_part, adaptative_remesh_parameters["metric_error_parameters"])
return metric_process
def _create_metric_process(self):
if (self.main_model_part.ProcessInfo[KratosMultiphysics.DOMAIN_SIZE] == 2):
metric_process = MeshingApplication.MetricErrorProcess2D(self.main_model_part, self.adaptative_remesh_parameters["metric_error_parameters"])
else:
metric_process = MeshingApplication.MetricErrorProcess3D(self.main_model_part, self.adaptative_remesh_parameters["metric_error_parameters"])
return metric_process
def _CreateMetricsProcess(self):
self.metric_processes = []
if self.strategy == "LevelSet":
level_set_parameters = KratosMultiphysics.Parameters("""{}""")
level_set_parameters.AddValue("minimal_size",self.settings["minimal_size"])
level_set_parameters.AddValue("maximal_size",self.settings["maximal_size"])
level_set_parameters.AddValue("sizing_parameters",self.settings["sizing_parameters"])
level_set_parameters.AddValue("enforce_current",self.settings["enforce_current"])
level_set_parameters.AddValue("anisotropy_remeshing",self.settings["anisotropy_remeshing"])
level_set_parameters.AddValue("anisotropy_parameters",self.settings["anisotropy_parameters"])
level_set_parameters["anisotropy_parameters"].RemoveValue("boundary_layer_min_size_ratio")
if self.domain_size == 2:
self.metric_processes.append(MeshingApplication.ComputeLevelSetSolMetricProcess2D(self.main_model_part, self.gradient_variable, level_set_parameters))
else:
self.metric_processes.append(MeshingApplication.ComputeLevelSetSolMetricProcess3D(self.main_model_part, self.gradient_variable, level_set_parameters))
elif self.strategy == "Hessian":
hessian_parameters = KratosMultiphysics.Parameters("""{}""")
hessian_parameters.AddValue("minimal_size",self.settings["minimal_size"])
hessian_parameters.AddValue("maximal_size",self.settings["maximal_size"])
hessian_parameters.AddValue("enforce_current",self.settings["enforce_current"])
hessian_parameters.AddValue("hessian_strategy_parameters",self.settings["hessian_strategy_parameters"])
hessian_parameters["hessian_strategy_parameters"].RemoveValue("metric_variable")
hessian_parameters.AddValue("anisotropy_remeshing",self.settings["anisotropy_remeshing"])
hessian_parameters.AddValue("enforce_anisotropy_relative_variable",self.settings["enforce_anisotropy_relative_variable"])
hessian_parameters.AddValue("enforced_anisotropy_parameters",self.settings["anisotropy_parameters"])
hessian_parameters["enforced_anisotropy_parameters"].RemoveValue("boundary_layer_min_size_ratio")
for current_metric_variable in self.metric_variable:
self.metric_processes.append(MeshingApplication.ComputeHessianSolMetricProcess(self.main_model_part, current_metric_variable, hessian_parameters))
elif self.strategy == "superconvergent_patch_recovery":
if not structural_dependencies:
raise Exception("You need to compile the StructuralMechanicsApplication in order to use this criteria")
# We compute the error
error_compute_parameters = KratosMultiphysics.Parameters("""{}""")
error_compute_parameters.AddValue("stress_vector_variable", self.settings["compute_error_extra_parameters"]["stress_vector_variable"])
error_compute_parameters.AddValue("echo_level", self.settings["echo_level"])
if self.domain_size == 2:
self.error_compute = StructuralMechanicsApplication.SPRErrorProcess2D(self.main_model_part, error_compute_parameters)
else:
self.error_compute = StructuralMechanicsApplication.SPRErrorProcess3D(self.main_model_part, error_compute_parameters)
# Now we compute the metric
error_metric_parameters = KratosMultiphysics.Parameters("""{}""")
error_metric_parameters.AddValue("minimal_size",self.settings["minimal_size"])
error_metric_parameters.AddValue("maximal_size",self.settings["maximal_size"])
error_metric_parameters.AddValue("target_error",self.settings["error_strategy_parameters"]["error_metric_parameters"]["interpolation_error"])
error_metric_parameters.AddValue("set_target_number_of_elements", self.settings["error_strategy_parameters"]["set_target_number_of_elements"])
error_metric_parameters.AddValue("target_number_of_elements", self.settings["error_strategy_parameters"]["target_number_of_elements"])
error_metric_parameters.AddValue("perform_nodal_h_averaging", self.settings["error_strategy_parameters"]["perform_nodal_h_averaging"])
error_metric_parameters.AddValue("echo_level", self.settings["echo_level"])
if self.domain_size == 2:
self.metric_process = MeshingApplication.MetricErrorProcess2D(self.main_model_part, error_metric_parameters)
else:
self.metric_process = MeshingApplication.MetricErrorProcess3D(self.main_model_part, error_metric_parameters)
def _create_metric_process(self):
if (self.main_model_part.ProcessInfo[KM.DOMAIN_SIZE] == 2):
metric_process = MA.MetricErrorProcess2D(
self.main_model_part,
self.adaptative_remesh_parameters["metric_error_parameters"])
else:
metric_process = MA.MetricErrorProcess3D(
self.main_model_part,
self.adaptative_remesh_parameters["metric_error_parameters"])
return metric_process
def _CreateMetricsProcess(self):
""" This method is responsible to create the metrics of the process
Keyword arguments:
self -- It signifies an instance of a class.
"""
self.metric_processes = []
if self.strategy == "LevelSet":
level_set_parameters = KratosMultiphysics.Parameters("""{}""")
level_set_parameters.AddValue("minimal_size",self.settings["minimal_size"])
level_set_parameters.AddValue("maximal_size",self.settings["maximal_size"])
level_set_parameters.AddValue("sizing_parameters",self.settings["sizing_parameters"])
level_set_parameters.AddValue("enforce_current",self.settings["enforce_current"])
level_set_parameters.AddValue("anisotropy_remeshing",self.settings["anisotropy_remeshing"])
level_set_parameters.AddValue("anisotropy_parameters",self.settings["anisotropy_parameters"])
level_set_parameters["anisotropy_parameters"].RemoveValue("boundary_layer_min_size_ratio")
if self.domain_size == 2:
self.metric_processes.append(MeshingApplication.ComputeLevelSetSolMetricProcess2D(self.main_model_part, self.gradient_variable, level_set_parameters))
else:
self.metric_processes.append(MeshingApplication.ComputeLevelSetSolMetricProcess3D(self.main_model_part, self.gradient_variable, level_set_parameters))
elif self.strategy == "Hessian":
hessian_parameters = KratosMultiphysics.Parameters("""{}""")
hessian_parameters.AddValue("minimal_size",self.settings["minimal_size"])
hessian_parameters.AddValue("maximal_size",self.settings["maximal_size"])
hessian_parameters.AddValue("enforce_current",self.settings["enforce_current"])
hessian_parameters.AddValue("hessian_strategy_parameters",self.settings["hessian_strategy_parameters"])
hessian_parameters["hessian_strategy_parameters"].RemoveValue("metric_variable")
hessian_parameters["hessian_strategy_parameters"].RemoveValue("non_historical_metric_variable")
hessian_parameters["hessian_strategy_parameters"].AddEmptyValue("non_historical_metric_variable")
hessian_parameters["hessian_strategy_parameters"].RemoveValue("normalization_factor")
hessian_parameters["hessian_strategy_parameters"].AddEmptyValue("normalization_factor")
hessian_parameters["hessian_strategy_parameters"].RemoveValue("normalization_alpha")
hessian_parameters["hessian_strategy_parameters"].AddEmptyValue("normalization_alpha")
hessian_parameters["hessian_strategy_parameters"].RemoveValue("normalization_method")
hessian_parameters["hessian_strategy_parameters"].AddEmptyValue("normalization_method")
hessian_parameters.AddValue("anisotropy_remeshing",self.settings["anisotropy_remeshing"])
hessian_parameters.AddValue("enforce_anisotropy_relative_variable",self.settings["enforce_anisotropy_relative_variable"])
hessian_parameters.AddValue("enforced_anisotropy_parameters",self.settings["anisotropy_parameters"])
hessian_parameters["enforced_anisotropy_parameters"].RemoveValue("boundary_layer_min_size_ratio")
for current_metric_variable, non_historical_metric_variable, normalization_factor, normalization_alpha, normalization_method in zip(self.metric_variables, self.non_historical_metric_variable, self.normalization_factor, self.normalization_alpha, self.normalization_method):
hessian_parameters["hessian_strategy_parameters"]["non_historical_metric_variable"].SetBool(non_historical_metric_variable)
hessian_parameters["hessian_strategy_parameters"]["normalization_factor"].SetDouble(normalization_factor)
hessian_parameters["hessian_strategy_parameters"]["normalization_alpha"].SetDouble(normalization_alpha)
hessian_parameters["hessian_strategy_parameters"]["normalization_method"].SetString(normalization_method)
self.metric_processes.append(MeshingApplication.ComputeHessianSolMetricProcess(self.main_model_part, current_metric_variable, hessian_parameters))
elif self.strategy == "superconvergent_patch_recovery" or self.strategy == "SPR":
# Generate SPR process
self.error_compute = self._GenerateErrorProcess()
# Now we compute the metric
error_metric_parameters = KratosMultiphysics.Parameters("""{"error_strategy_parameters":{}}""")
error_metric_parameters.AddValue("minimal_size",self.settings["minimal_size"])
error_metric_parameters.AddValue("maximal_size",self.settings["maximal_size"])
error_metric_parameters["error_strategy_parameters"].AddValue("target_error",self.settings["error_strategy_parameters"]["error_metric_parameters"]["interpolation_error"])
error_metric_parameters["error_strategy_parameters"].AddValue("set_target_number_of_elements", self.settings["error_strategy_parameters"]["set_target_number_of_elements"])
error_metric_parameters["error_strategy_parameters"].AddValue("target_number_of_elements", self.settings["error_strategy_parameters"]["target_number_of_elements"])
error_metric_parameters["error_strategy_parameters"].AddValue("perform_nodal_h_averaging", self.settings["error_strategy_parameters"]["perform_nodal_h_averaging"])
error_metric_parameters.AddValue("echo_level", self.settings["echo_level"])
if self.domain_size == 2:
self.metric_process = MeshingApplication.MetricErrorProcess2D(self.main_model_part, error_metric_parameters)
else:
self.metric_process = MeshingApplication.MetricErrorProcess3D(self.main_model_part, error_metric_parameters)
