PyFemtet 1.8.3__tar.gz → 1.9.0__tar.gz

{pyfemtet-1.8.3 → pyfemtet-1.9.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pyfemtet
-Version: 1.8.3
+Version: 1.9.0
 Summary: Design parameter optimization using Femtet.
 License-Expression: MIT
 License-File: LICENSE

{pyfemtet-1.8.3 → pyfemtet-1.9.0}/pyfemtet/opt/femopt.py

@@ -49,8 +49,8 @@ class FEMOpt:
 
     def __init__(
             self,
-            fem: AbstractFEMInterface = None,
-            opt: AbstractOptimizer = None,
+            fem: AbstractFEMInterface | None = None,
+            opt: AbstractOptimizer | None = None,
     ):
         self.opt: AbstractOptimizer = opt or OptunaOptimizer()
         
@@ -160,7 +160,7 @@ class FEMOpt:
             names: str | list[str],
             fun: Callable[..., Sequence[float]],
             n_return: int,
-            directions: DIRECTION | Sequence[DIRECTION] = None,
+            directions: DIRECTION | Sequence[DIRECTION] | None = None,
             args: tuple | None = None,
             kwargs: dict | None = None,
     ):
@@ -208,7 +208,7 @@ class FEMOpt:
     ):
         self.opt.set_termination_condition(func)
 
-    def set_monitor_host(self, host: str = None, port: int = None):
+    def set_monitor_host(self, host: str | None = None, port: int | None = None):
         """Sets the host IP address and the port of the process monitor.
 
         Args:
@@ -240,14 +240,14 @@ class FEMOpt:
 
     def optimize(
             self,
-            n_trials: int = None,
+            n_trials: int | None = None,
             n_parallel: int = 1,
-            timeout: float = None,
+            timeout: float | None = None,
             wait_setup: bool = True,
             confirm_before_exit: bool = True,
-            history_path: str = None,
+            history_path: str | None = None,
             with_monitor: bool = True,
-            scheduler_address: str = None,
+            scheduler_address: str | None = None,
             seed: int | None = None,
     ):
 
@@ -258,7 +258,7 @@ class FEMOpt:
                 ver=pyfemtet.__version__,
             )
         )
-        client: Client
+        client: Client | DummyClient
 
         # set arguments
         self.opt.n_trials = n_trials or self.opt.n_trials
@@ -274,7 +274,6 @@ class FEMOpt:
             worker_name_base = 'Sub'
             if n_parallel == 1:
                 cluster = nullcontext()
-                # noinspection PyTypeChecker
                 client = DummyClient()
 
             else:
@@ -534,137 +533,3 @@ class FEMOpt:
         logger.info(_('All processes are terminated.'))
 
         return df
-
-
-def debug_1():
-    # noinspection PyUnresolvedReferences
-    from time import sleep
-    # from pyfemtet.opt.optimizer import InterruptOptimization
-    import optuna
-    from pyfemtet.opt.interface import AbstractFEMInterface, NoFEM
-
-    def _parabola(_fem: AbstractFEMInterface, _opt: AbstractOptimizer):
-        x = _opt.get_variables('values')
-        # print(os.getpid())
-        # raise RuntimeError
-        # raise Interrupt
-        # if get_worker() is None:
-        #     raise RuntimeError
-        return (x ** 2).sum()
-
-    def _cns(_fem: AbstractFEMInterface, _opt: AbstractOptimizer):
-        x = _opt.get_variables('values')
-        return x[0]
-
-    _opt = OptunaOptimizer()
-    _opt.sampler = optuna.samplers.TPESampler(seed=42)
-    _opt.n_trials = 10
-
-    _fem = NoFEM()
-    _opt.fem = _fem
-
-    _args = (_opt,)
-
-    _opt.add_parameter('x1', 1, -1, 1, step=0.2)
-    _opt.add_parameter('x2', 1, -1, 1, step=0.2)
-
-    _opt.add_constraint('cns', _cns, lower_bound=-0.5, args=_args)
-
-    _opt.add_objective('obj', _parabola, args=_args)
-
-    _femopt = FEMOpt(fem=_fem, opt=_opt)
-    _femopt.opt = _opt
-    # _femopt.opt.history.path = 'v1test/femopt-restart-test.csv'
-    _femopt.optimize(n_parallel=2)
-
-    print(os.path.abspath(_femopt.opt.history.path))
-
-
-def substrate_size(Femtet):
-    """基板のXY平面上での専有面積を計算します。"""
-    substrate_w = Femtet.GetVariableValue('substrate_w')
-    substrate_d = Femtet.GetVariableValue('substrate_d')
-
-    # assert get_worker() is not None
-
-    return substrate_w * substrate_d  # 単位: mm2
-
-
-def debug_2():
-    from pyfemtet.opt.interface import FemtetInterface
-    from pyfemtet.opt.optimizer import OptunaOptimizer
-
-    fem = FemtetInterface(
-        femprj_path=os.path.join(os.path.dirname(__file__), 'wat_ex14_parametric_jp.femprj'),
-    )
-
-    opt = OptunaOptimizer()
-
-    opt.fem = fem
-
-    opt.add_parameter(name="substrate_w", initial_value=40, lower_bound=22, upper_bound=60)
-    opt.add_parameter(name="substrate_d", initial_value=60, lower_bound=34, upper_bound=60)
-    opt.add_objective(name='基板サイズ(mm2)', fun=substrate_size)
-    opt.add_objective(name='obj2', fun=substrate_size)
-    opt.add_objective(name='obj3', fun=substrate_size)
-
-    opt.n_trials = 10
-    # opt.history.path = os.path.join(os.path.dirname(__file__), 'femtet-test.csv')
-
-    femopt = FEMOpt()
-
-    femopt.opt = opt
-
-    femopt.optimize(n_parallel=1)
-
-
-def debug_3():
-    # noinspection PyUnresolvedReferences
-    from time import sleep
-    # from pyfemtet.opt.optimizer import InterruptOptimization
-    import optuna
-    from pyfemtet.opt.interface import AbstractFEMInterface, NoFEM
-
-    def _parabola(_fem: AbstractFEMInterface, _opt: AbstractOptimizer):
-        x = _opt.get_variables('values')
-        # print(os.getpid())
-        # raise RuntimeError
-        # raise Interrupt
-        # if get_worker() is None:
-        #     raise RuntimeError
-        return (x ** 2).sum()
-
-    def _cns(_fem: AbstractFEMInterface, _opt: AbstractOptimizer):
-        x = _opt.get_variables('values')
-        return x[0]
-
-    _opt = OptunaOptimizer()
-    _opt.sampler = optuna.samplers.TPESampler(seed=42)
-
-    _fem = NoFEM()
-    _opt.fem = _fem
-
-    _args = (_opt,)
-
-    _opt.add_parameter('x1', 1, -1, 1, step=0.2)
-    _opt.add_parameter('x2', 1, -1, 1, step=0.2)
-
-    _opt.add_constraint('cns', _cns, lower_bound=-0.5, args=_args)
-
-    _opt.add_objective('obj', _parabola, args=_args)
-
-    _femopt = FEMOpt(fem=_fem, opt=_opt)
-    _femopt.optimize(
-        scheduler_address='<dask scheduler で起動したスケジューラの tcp をここに入力>',
-        n_trials=80,
-        n_parallel=6,
-        with_monitor=True,
-        confirm_before_exit=False,
-    )
-
-
-if __name__ == '__main__':
-    # for i in range(1):
-        debug_1()
-    # debug_2()
-    # debug_3()

{pyfemtet-1.8.3 → pyfemtet-1.9.0}/pyfemtet/opt/interface/_base_interface.py

@@ -244,6 +244,16 @@ class AbstractFEMInterface:
     def _get_additional_data(self) -> dict:
         return dict()
 
+    # name
+    @property
+    def name(self) -> str:
+        if hasattr(self, '_name'):
+            return self._name
+        return self.__class__.__name__
+
+    @name.setter
+    def name(self, value: str):
+        self._name = value
 
 class COMInterface(AbstractFEMInterface):
 

{pyfemtet-1.8.3 → pyfemtet-1.9.0}/pyfemtet/opt/interface/_excel_interface/excel_interface.py

@@ -662,7 +662,6 @@ class ExcelInterface(COMInterface):
                     upper_bound=ub,
                     step=step,
                     fix=not use,
-                    supress_duplicated_name_check=True,
                 )
 
             elif kind == 'Categorical':
@@ -675,7 +674,6 @@ class ExcelInterface(COMInterface):
                     initial_value=value,
                     choices=choices,
                     fix=not use,
-                    supress_duplicated_name_check=True,
                 )
 
             else:
@@ -722,7 +720,6 @@ class ExcelInterface(COMInterface):
                     direction=direction,
                     fun=_ScapeGoatObjective(),
                     kwargs=dict(name=name),
-                    supress_duplicated_name_check=True,
                 )
 
     def load_constraints(self, opt: OptimizationDataPerFEM, raise_if_no_keyword=False):
@@ -784,7 +781,6 @@ class ExcelInterface(COMInterface):
                     fun=_ScapeGoatObjective(),
                     kwargs=dict(name=name),
                     using_fem=not calc_before_solve,
-                    supress_duplicated_name_check=True,
                 )
 
     def objective_from_excel(self, _, name: str):

{pyfemtet-1.8.3 → pyfemtet-1.9.0}/pyfemtet/opt/interface/_femtet_interface/_femtet_parametric.py

@@ -223,7 +223,6 @@ def add_parametric_results_as_objectives(opt, Femtet, indexes, directions) -> bo
         # objective value function
         opt.add_objective(
             name=name, fun=_parametric_objective, direction=direction, args=(i,),
-            supress_duplicated_name_check=True,            
         )
     return True  # ここまで来たら成功
 

{pyfemtet-1.8.3 → pyfemtet-1.9.0}/pyfemtet/opt/interface/_femtet_interface/femtet_interface.py

@@ -185,6 +185,17 @@ class FemtetInterface(COMInterface):
 
     # ===== system =====
 
+    @property
+    def name(self) -> str:
+        name = ""
+        if self.femprj_path is not None:
+            name += os.path.basename(self.femprj_path)
+        if self.model_name is not None:
+            name += f' - {self.model_name}'
+        if name == "":
+            name = super().name
+        return name
+
     @property
     def object_pass_to_fun(self):
         """The object pass to the first argument of user-defined objective functions.

{pyfemtet-1.8.3 → pyfemtet-1.9.0}/pyfemtet/opt/interface/_surrogate_model_interface/base_surrogate_interface.py

@@ -122,7 +122,7 @@ class AbstractSurrogateModelInterfaceBase(AbstractFEMInterface):
             def dummy(*args, **kwargs):
                 assert False
 
-            ctx.add_objective(obj_name, dummy, direction, supress_duplicated_name_check=True)
+            ctx.add_objective(obj_name, dummy, direction)
         
         # ctx の目的関数のうち train に含まれるものを差し替え
         # directions の目的変数名を validation しているので漏れはないはず

{pyfemtet-1.8.3 → pyfemtet-1.9.0}/pyfemtet/opt/optimizer/_base_optimizer.py

@@ -84,6 +84,11 @@ class OptimizationData:
         self._initialize_objectives()
         self.constraints = Constraints()
         self.other_outputs = Functions()
+        self._suppress_duplicate_name_check = False
+
+    def _check_duplicated_name(self, name, names, suppress=False):
+        if not suppress and not self._suppress_duplicate_name_check:
+            _duplicated_name_check(name, names)
 
     def add_constant_value(
             self,
@@ -110,8 +115,7 @@ class OptimizationData:
         var.value = value
         var.pass_to_fem = pass_to_fem
         var.properties = properties if properties is not None else {}
-        if not supress_duplicated_name_check:
-            _duplicated_name_check(name, self.variable_manager.variables.keys())
+        self._check_duplicated_name(name, self.variable_manager.variables.keys(), supress_duplicated_name_check)
         self.variable_manager.set_variable(var)
         return var
 
@@ -141,8 +145,7 @@ class OptimizationData:
         prm.step = step
         prm.properties = properties
         prm.pass_to_fem = pass_to_fem
-        if not supress_duplicated_name_check:
-            _duplicated_name_check(name, self.variable_manager.variables.keys())
+        self._check_duplicated_name(name, self.variable_manager.variables.keys(), supress_duplicated_name_check)
         self.variable_manager.set_variable(prm)
         return prm
 
@@ -162,8 +165,7 @@ class OptimizationData:
         )
         var.properties = properties or dict()
         var.pass_to_fem = pass_to_fem
-        if not supress_duplicated_name_check:
-            _duplicated_name_check(name, self.variable_manager.variables.keys())
+        self._check_duplicated_name(name, self.variable_manager.variables.keys(), supress_duplicated_name_check)
         self.variable_manager.set_variable(var)
         return var
 
@@ -181,8 +183,7 @@ class OptimizationData:
         var._expr = ExpressionFromString.InternalClass(sympy_expr=sympy_expr)
         var.properties = properties or dict()
         var.pass_to_fem = pass_to_fem
-        if not supress_duplicated_name_check:
-            _duplicated_name_check(name, self.variable_manager.variables.keys())
+        self._check_duplicated_name(name, self.variable_manager.variables.keys(), supress_duplicated_name_check)
         self.variable_manager.set_variable(var)
         return var
 
@@ -204,8 +205,7 @@ class OptimizationData:
         var.kwargs = kwargs or dict()
         var.properties = properties or dict()
         var.pass_to_fem = pass_to_fem
-        if not supress_duplicated_name_check:
-            _duplicated_name_check(name, self.variable_manager.variables.keys())
+        self._check_duplicated_name(name, self.variable_manager.variables.keys(), supress_duplicated_name_check)
         self.variable_manager.set_variable(var)
         return var
 
@@ -231,8 +231,7 @@ class OptimizationData:
         prm.choices = choices
         prm.properties = properties
         prm.pass_to_fem = pass_to_fem
-        if not supress_duplicated_name_check:
-            _duplicated_name_check(name, self.variable_manager.variables.keys())
+        self._check_duplicated_name(name, self.variable_manager.variables.keys(), supress_duplicated_name_check)
         self.variable_manager.set_variable(prm)
         return prm
 
@@ -251,8 +250,7 @@ class OptimizationData:
         obj.kwargs = kwargs or {}
         obj.direction = direction
         obj.fem_ctx = None
-        if not supress_duplicated_name_check:
-            _duplicated_name_check(name, self.objectives.keys())
+        self._check_duplicated_name(name, self.objectives.keys(), supress_duplicated_name_check)
         self.objectives.update({name: obj})
         return obj  # Context で fem_ctx をセットするために返す
 
@@ -334,8 +332,7 @@ class OptimizationData:
         cns.hard = strict
         cns.fem_ctx = None
         cns.using_fem = using_fem
-        if not supress_duplicated_name_check:
-            _duplicated_name_check(name, self.constraints.keys())
+        self._check_duplicated_name(name, self.constraints.keys(), supress_duplicated_name_check)
         self.constraints.update({name: cns})
         return cns  # Context で fem_ctx をセットするために返す
 
@@ -352,8 +349,7 @@ class OptimizationData:
         other_func.args = args or ()
         other_func.kwargs = kwargs or {}
         other_func.fem_ctx = None
-        if not supress_duplicated_name_check:
-            _duplicated_name_check(name, self.other_outputs.keys())
+        self._check_duplicated_name(name, self.other_outputs.keys(), supress_duplicated_name_check)
         self.other_outputs.update({name: other_func})
         return other_func
 
@@ -473,8 +469,21 @@ class OptimizationDataPerFEM(OptimizationData):
         )
 
     def _y_common(self, fem: AbstractFEMInterface) -> TrialOutput:
+
+        if len(self.objectives) > 0:
+            logger.info(_(
+                en_message='Evaluating objective functions for {fem_name}...',
+                jp_message='{fem_name} の目的関数を評価しています...',
+                fem_name=fem.name,
+            ))
+
         out = TrialOutput()
         for name, obj in self.objectives.items():
+            logger.info(_(
+                en_message='  Evaluating objective function {obj_name}...',
+                jp_message='  目的関数 {obj_name} を評価しています...',
+                obj_name=name,
+            ))
             obj_result = ObjectiveResult(obj, fem)
             out.update({name: obj_result})
         return out
@@ -496,8 +505,28 @@ class OptimizationDataPerFEM(OptimizationData):
             fem: AbstractFEMInterface,
             hard: bool,
     ) -> TrialConstraintOutput:
+        if len([cns for cns in self.constraints.values() if cns.hard == hard]) > 0:
+            if hard:
+                logger.info(_(
+                    en_message='Evaluating hard constraint functions...',
+                    jp_message='{fem_name} の hard 拘束関数を評価しています...',
+                    fem_name=fem.name,
+                ))
+            else:
+                logger.info(_(
+                    en_message='Evaluating soft constraint functions...',
+                    jp_message='{fem_name} の soft 拘束関数を評価しています...',
+                    fem_name=fem.name,
+
+                ))
+
         for name, cns in self.constraints.items():
             if cns.hard == hard:
+                logger.info(_(
+                    en_message='  Evaluating constraint function {cns_name}...',
+                    jp_message='  制約関数 {cns_name} を評価しています...',
+                    cns_name=name,
+                ))
                 cns_result = ConstraintResult(cns, fem)
                 out.update({name: cns_result})
         return out
@@ -509,6 +538,14 @@ class OptimizationDataPerFEM(OptimizationData):
         return self._c_common(out, self.fem, hard=False)
 
     def _other_outputs_common(self, out: TrialFunctionOutput, fem: AbstractFEMInterface) -> TrialFunctionOutput:
+
+        if len(self.other_outputs) > 0:
+            logger.info(_(
+                en_message='Evaluating other outputs for {fem_name}...',
+                jp_message='{fem_name} のその他の出力を評価しています...',
+                fem_name=fem.name,
+            ))
+
         for name, other_func in self.other_outputs.items():
             other_func_result = FunctionResult(other_func, fem)
             out.update({name: other_func_result})
@@ -525,6 +562,16 @@ class OptimizationDataPerFEM(OptimizationData):
 
 class FEMListForGlobal(FEMListInterface):
 
+    @property
+    def name(self) -> str:
+        if len(self._fems) == 1:
+            return self._fems[0].name
+        else:
+            return _(
+                en_message='Entire',
+                jp_message='全体',
+            )
+
     # 単一 FEM しか使わない場合に list を意識させないため
     # 長さ 1 の時は要素を返す処理が方々で必要
     @property
@@ -989,34 +1036,47 @@ class AbstractOptimizer(OptimizationData):
 
                 from contextlib import nullcontext
 
+                def _get_pass_to_fem(ctx__):
+                    # 共通
+                    pass_to_fem__ = self.opt.fem_manager.global_data.variable_manager.get_variables(
+                        filter="pass_to_fem",
+                        format="raw",
+                    )
+                    # FEM ごとのものを追加
+                    pass_to_fem__.update(ctx__.variable_manager.get_variables(
+                        filter="pass_to_fem",
+                        format="raw",
+                    ))
+                    return pass_to_fem__
+
                 # start solve per FEM
                 if opt_.sub_fidelity_name != MAIN_FIDELITY_NAME:
                     logger.info('----------')
                     logger.info(_('fidelity: ({name})', name=opt_.sub_fidelity_name))
-                for ctx in opt_.fem_manager.contexts:
-                    logger.info(_('input variables:'))
-                    logger.info(parameters)
 
-                    # ===== update FEM parameter =====
-                    with nullcontext():
-                        # 共通
-                        pass_to_fem = self.opt.fem_manager.global_data.variable_manager.get_variables(
-                            filter="pass_to_fem",
-                            format="raw",
-                        )
-                        # FEM ごとのものを追加
-                        pass_to_fem.update(ctx.variable_manager.get_variables(
-                            filter="pass_to_fem",
-                            format="raw",
+                # ===== Phase 1: update_parameter + hard constraint for ALL contexts =====
+                for ctx in opt_.fem_manager.contexts:
+                    # ===== update FEM parameter (for hard constraint evaluation) =====
+                    if not isinstance(ctx.fem, FEMListForGlobal):
+                        pass_to_fem = _get_pass_to_fem(ctx)
+
+                        logger.info(_(
+                            en_message='Updating variables for {fem_name}...',
+                            jp_message='{fem_name} の変数を更新しています...',
+                            fem_name=ctx.fem.name,
                         ))
-                        logger.info(_('updating variables...'))
+                        logger.info(_(
+                            en_message='  Input: {variables}',
+                            jp_message='  入力: {variables}',
+                            variables=pass_to_fem,
+                        ))
+
                         ctx.fem.update_parameter(pass_to_fem)
                         opt_._check_and_raise_interruption()
 
                     # ===== evaluate hard constraint =====
                     with nullcontext():
                         ctx_hard_c = TrialConstraintOutput()
-                        logger.info(_('evaluating constraint functions...'))
 
                         try:
                             ctx._hard_c(ctx_hard_c)
@@ -1047,9 +1107,23 @@ class AbstractOptimizer(OptimizationData):
 
                             raise HardConstraintViolation
 
+                # ===== Phase 2: FEM execution + objectives + soft constraints + other outputs =====
+                for ctx in opt_.fem_manager.contexts:
+
+                    # ===== update FEM parameter (again, for Femtet parametric compatibility) =====
+                    with nullcontext():
+                        pass_to_fem = _get_pass_to_fem(ctx)
+                        ctx.fem.update_parameter(pass_to_fem)
+                        opt_._check_and_raise_interruption()
+
                     # ===== update FEM =====
                     with nullcontext():
-                        logger.info(_('Solving FEM...'))
+                        if not isinstance(ctx.fem, FEMListForGlobal):
+                            logger.info(_(
+                                en_message='Solving {fem_name}...',
+                                jp_message='{fem_name} の解析を実行しています...',
+                                fem_name=ctx.fem.name,
+                            ))
 
                         try:
                             ctx.fem.update()
@@ -1072,8 +1146,6 @@ class AbstractOptimizer(OptimizationData):
 
                     # ===== evaluate y =====
                     with nullcontext():
-                        logger.info(_('evaluating objective functions...'))
-
                         try:
                             ctx_y: TrialOutput = ctx._y()
                             record.y.update(ctx_y)
@@ -1097,7 +1169,6 @@ class AbstractOptimizer(OptimizationData):
                     # ===== evaluate soft constraint =====
                     with nullcontext():
                         ctx_soft_c = TrialConstraintOutput()
-                        logger.info(_('evaluating remaining constraints...'))
 
                         try:
                             ctx._soft_c(ctx_soft_c)
@@ -1117,8 +1188,6 @@ class AbstractOptimizer(OptimizationData):
 
                     # ===== evaluate other functions =====
                     with nullcontext():
-                        logger.info(_('evaluating other functions...'))
-
                         ctx_other_outputs = TrialFunctionOutput()
                         try:
                             ctx._other_outputs(ctx_other_outputs)
@@ -1141,8 +1210,11 @@ class AbstractOptimizer(OptimizationData):
 
                             raise e
 
-                logger.info(_('output:'))
-                logger.info(record.y)
+                logger.info(_(
+                    en_message='Trial finished. Objective values:',
+                    jp_message='トライアルが終了しました。目的関数の値:',
+                ))
+                logger.info(f'  出力: {record.y}')
                 record.state = TrialState.succeeded
 
                 opt_._check_and_raise_interruption()
@@ -1390,27 +1462,42 @@ class AbstractOptimizer(OptimizationData):
         return LoggingOutput()
 
     def _refresh_problem(self):
-        # ctx の内容が減っている場合でも検出できるように初期化
-        self._initialize_problem()
 
-        for ctx in self.fem_manager.contexts:
-            # femlist が _load すると各 FEMInterface が
-            # 想定するのと異なる object_pass_to_fun (=Sequence)が
-            # 渡されてしまうので、
-            # 各 context のみが _load_... を行う必要がある。
-            if not isinstance(ctx, GlobalOptimizationData):
-                ctx._load_problem_from_fem()
-
-            # optimizer に反映する
-            self.objectives.update(ctx.objectives)
-            self.constraints.update(ctx.constraints)
-            self.other_outputs.update(ctx.other_outputs)
-            self.variable_manager.variables.update(ctx.variable_manager.variables)
-
-        # 特殊処理が必要な場合は最後に fem の責任で行う
-        for ctx in self.fem_manager.contexts:
-            if not isinstance(ctx, GlobalOptimizationData):
-                ctx.fem.contact_to_optimizer(self, self.fem_manager.global_data, ctx)
+        class SuppressDuplicateNameCheck:
+
+            def __init__(self, opt: OptimizationData):
+                self._opt = opt
+
+            def __enter__(self_):
+                self_._opt._suppress_duplicate_name_check = True
+
+            def __exit__(self_, exc_type, exc_val, exc_tb):
+                self_._opt._suppress_duplicate_name_check = False
+
+        with SuppressDuplicateNameCheck(self):
+            # ctx の内容が減っている場合でも検出できるように初期化
+            self._initialize_problem()
+
+            for ctx in self.fem_manager.contexts:
+                with SuppressDuplicateNameCheck(ctx):
+                    # femlist が _load すると各 FEMInterface が
+                    # 想定するのと異なる object_pass_to_fun (=Sequence)が
+                    # 渡されてしまうので、
+                    # 各 context のみが _load_... を行う必要がある。
+                    if not isinstance(ctx, GlobalOptimizationData):
+                        ctx._load_problem_from_fem()
+
+                    # optimizer に反映する
+                    self.objectives.update(ctx.objectives)
+                    self.constraints.update(ctx.constraints)
+                    self.other_outputs.update(ctx.other_outputs)
+                    self.variable_manager.variables.update(ctx.variable_manager.variables)
+
+            # 特殊処理が必要な場合は最後に fem の責任で行う
+            for ctx in self.fem_manager.contexts:
+                if not isinstance(ctx, GlobalOptimizationData):
+                    with SuppressDuplicateNameCheck(ctx):
+                        ctx.fem.contact_to_optimizer(self, self.fem_manager.global_data, ctx)
 
     # noinspection PyMethodMayBeStatic
     def _get_additional_data(self) -> dict:

{pyfemtet-1.8.3 → pyfemtet-1.9.0}/pyproject.toml

@@ -34,7 +34,7 @@ dependencies = [
 
 #[project.optional-dependencies]
 #brep = ["brepmatching>=0.1.8,<0.2"]
-version = "1.8.3"
+version = "1.9.0"
 
 [project.urls]
 Repository = "https://github.com/pyfemtet/pyfemtet"