pyfemtet 0.4.9__py3-none-any.whl → 0.4.11__py3-none-any.whl

pyfemtet/__init__.py

@@ -1 +1 @@
-__version__ = "0.4.9"
+__version__ = "0.4.11"

pyfemtet/opt/_femopt.py

@@ -13,7 +13,7 @@ from dask.distributed import LocalCluster, Client
 # pyfemtet relative
 from pyfemtet.opt.interface import FEMInterface, FemtetInterface
 from pyfemtet.opt.opt import AbstractOptimizer, OptunaOptimizer
-from pyfemtet.opt.visualization._monitor import ProcessMonitorApp
+from pyfemtet.opt.visualization.process_monitor.application import main as process_monitor_main
 from pyfemtet.opt._femopt_core import (
     _check_bound,
     _is_access_gogh,
@@ -574,6 +574,12 @@ def _start_monitor_server(
         host=None,
         port=None,
 ):
-    monitor = ProcessMonitorApp(history, status, worker_addresses, worker_status_list)
-    monitor.start_server(host, port)
+    process_monitor_main(
+        history,
+        status,
+        worker_addresses,
+        worker_status_list,
+        host,
+        port,
+    )
     return 'Exit monitor server process gracefully'

pyfemtet/opt/femprj_sample/wat_ex14_parametric_parallel.py

@@ -0,0 +1,66 @@
+"""Multi-objective optimization: heating element on substrate.
+
+Using Femtet's heat conduction analysis solver, we will design
+to reduce the chip temperature and shrink the board size.
+
+Corresponding project: wat_ex14_parametric.femprj
+"""
+from pyfemtet.opt import FEMOpt
+
+
+def chip_temp(Femtet, chip_name):
+    """Obtain the maximum temperature of the chip.
+
+    Note:
+        The objective or constraint function should take Femtet
+        as its first argument and return a float as the output.
+
+    Params:
+        Femtet: An instance for manipulating Femtet with macros. For detailed information, please refer to "Femtet Macro Help".
+        chip_name (str): The body attribute name defined in femprj. Valid values are 'MAINCHIP' or 'SUBCHIP'.
+
+    Returns:
+        float: The maximum temperature of the body with the specified body attribute name.
+    """
+    Gogh = Femtet.Gogh
+
+    max_temperature, min_temperature, mean_temperature = Gogh.Watt.GetTemp(chip_name)
+
+    return max_temperature  # unit: degree
+
+
+def substrate_size(Femtet):
+    """Calculate the occupied area on the XY plane of the substrate."""
+    substrate_w = Femtet.GetVariableValue('substrate_w')
+    substrate_d = Femtet.GetVariableValue('substrate_d')
+    return substrate_w * substrate_d  # unit: mm2
+
+
+if __name__ == '__main__':
+
+    # Initialize the FEMOpt object.
+    # (establish connection between the optimization problem and Femtet)
+    femopt = FEMOpt()
+
+    # Add design variables to the optimization problem.
+    # (Specify the variables registered in the femprj file.)
+    femopt.add_parameter("substrate_w", 40, lower_bound=22, upper_bound=40)
+    femopt.add_parameter("substrate_d", 60, lower_bound=33, upper_bound=60)
+
+    # Add the objective function to the optimization problem.
+    # The target bending angle is 90 degrees.
+    femopt.add_objective(chip_temp, name='max temp. of<br>MAINCHIP (degree)', direction='minimize', args=('MAINCHIP',))
+    femopt.add_objective(chip_temp, name='max temp. of<br>SUBCHIP (degree)', direction='minimize', args=('SUBCHIP',))
+    femopt.add_objective(substrate_size, name='substrate size')
+
+    # Run optimization.
+    femopt.set_random_seed(42)
+    femopt.optimize(n_trials=15, n_parallel=3)  # This line is the only difference with no parallel pattern.
+
+    # Stop script to keep process alive
+    # while you check the result in process monitor.
+    print('================================')
+    print('Finished. Press Enter to quit...')
+    print('================================')
+    input()
+    femopt.terminate_all()

pyfemtet/opt/femprj_sample_jp/wat_ex14_parametric_parallel_jp.py

@@ -0,0 +1,64 @@
+"""多目的最適化: プリント基板上ICの発熱
+
+Femtetの熱伝導解析ソルバを使用して、ICチップの発熱を抑えつつ
+基板サイズを小さくする設計を行います。
+
+対応プロジェクト: wat_ex14_parametric_jp.femprj
+"""
+from pyfemtet.opt import FEMOpt
+
+
+def chip_temp(Femtet, chip_name):
+    """チップの最高温度を取得します。
+
+    Note:
+        目的関数または制約関数は、
+        第一引数としてFemtetを受け取り、
+        戻り値としてfloat型を返す必要があります。
+
+    Params:
+        Femtet: Femtet をマクロで操作するためのインスタンスです。詳細な情報については、「Femtet マクロヘルプ」をご覧ください。
+        chip_name (str): femprj 内で定義されているボディ属性名です。有効な値は 'MAINCHIP' 又は 'SUBCHIP' です。
+
+    Returns:
+        float: 指定されたボディ属性名のボディの最高温度です。
+    """
+    Gogh = Femtet.Gogh
+
+    max_temperature, min_temperature, mean_temperature = Gogh.Watt.GetTemp(chip_name)
+
+    return max_temperature  # 単位: 度
+
+
+def substrate_size(Femtet):
+    """基板のXY平面上での専有面積を計算します。"""
+    substrate_w = Femtet.GetVariableValue('substrate_w')
+    substrate_d = Femtet.GetVariableValue('substrate_d')
+    return substrate_w * substrate_d  # 単位: mm2
+
+
+if __name__ == '__main__':
+
+    # FEMOpt オブジェクトの初期化 (最適化問題とFemtetとの接続を行います)
+    femopt = FEMOpt()
+
+    # 設計変数を最適化問題に追加 (femprj ファイルに登録されている変数を指定してください)
+    femopt.add_parameter("substrate_w", 40, lower_bound=22, upper_bound=60)
+    femopt.add_parameter("substrate_d", 60, lower_bound=34, upper_bound=60)
+
+    # 目的関数を最適化問題に追加
+    femopt.add_objective(chip_temp, name='MAINCHIP<br>最高温度(度)', direction='minimize', args=('MAINCHIP',))
+    femopt.add_objective(chip_temp, name='SUBCHIP<br>最高温度(度)', direction='minimize', args=('SUBCHIP',))
+    femopt.add_objective(substrate_size, name='基板サイズ(mm2)')
+
+    # 最適化を実行
+    femopt.set_random_seed(42)
+    femopt.optimize(n_trials=15, n_parallel=3)  # この部分のみ変更します
+
+    # プロセスモニタで結果を確認するために
+    # Enter キーが押されるまで処理を停止します。
+    print('================================')
+    print('Finished. Press Enter to quit...')
+    print('================================')
+    input()
+    femopt.terminate_all()

pyfemtet/opt/interface/_femtet.py

@@ -3,15 +3,19 @@ from typing import Optional, List
 import os
 import sys
 from time import sleep, time
-import signal
 
 import pandas as pd
 import psutil
 from dask.distributed import get_worker
 
-from pywintypes import com_error
+# noinspection PyUnresolvedReferences
+from pywintypes import com_error, error
+# noinspection PyUnresolvedReferences
 from pythoncom import CoInitialize, CoUninitialize
+# noinspection PyUnresolvedReferences
 from win32com.client import constants
+import win32con
+import win32gui
 from femtetutils import util
 
 from pyfemtet.core import (
@@ -31,6 +35,10 @@ from pyfemtet.dispatch_extensions import (
 from pyfemtet.opt.interface import FEMInterface, logger
 
 
+def post_activate_message(hwnd):
+    win32gui.PostMessage(hwnd, win32con.WM_ACTIVATE, win32con.WA_ACTIVE, 0)
+
+
 class FemtetInterface(FEMInterface):
     """Concrete class for the interface with Femtet.
 
@@ -113,7 +121,7 @@ class FemtetInterface(FEMInterface):
         if self.original_femprj_path is None:
             # dask worker でなければ original のはず
             try:
-                worker = get_worker()
+                _ = get_worker()
             except ValueError:
                 self.original_femprj_path = self.femprj_path
 
@@ -252,7 +260,7 @@ class FemtetInterface(FEMInterface):
 
         """
 
-        # FIXME: Gaudi へのアクセスなど、self.Femtet.Gaudi.somefunc() のような場合、この関数を呼び出す前に Gaudi へのアクセスの時点で com_error が起こる
+        # FIXME: Gaudi へのアクセスなど、self.Femtet.Gaudi.SomeFunc() のような場合、この関数を呼び出す前に Gaudi へのアクセスの時点で com_error が起こる
         # FIXME: => 文字列で渡して eval() すればよい。
 
         if args is None:
@@ -282,9 +290,16 @@ class FemtetInterface(FEMInterface):
 
         # API を実行
         try:
-            returns = fun(*args, **kwargs)
-        except com_error:
+            # 解析結果を開いた状態で Gaudi.Activate して ReExecute する場合、ReExecute の前後にアクティブ化イベントが必要
+            if fun.__name__ == 'ReExecute':
+                post_activate_message(self.Femtet.hWnd)  # can raise pywintypes.error
+                returns = fun(*args, **kwargs)
+                post_activate_message(self.Femtet.hWnd)
+            else:
+                returns = fun(*args, **kwargs)
+        except (com_error, error):
             # パターン 2 エラーが生じたことは確定なのでエラーが起こるよう returns を作る
+            # com_error ではなく error の場合はおそらく Femtet が落ちている
             if ret_for_check_idx is None:
                 returns = return_value_if_failed
             else:
@@ -665,6 +680,13 @@ class FemtetInterface(FEMInterface):
         )
         return out
 
+    @staticmethod
+    def create_pdt_path(femprj_path, model_name, trial):
+        result_dir = femprj_path.replace('.femprj', '.Results')
+        pdt_path = os.path.join(result_dir, model_name + f'_trial{trial}.pdt')
+        return pdt_path
+
+    # noinspection PyMethodOverriding
     @staticmethod
     def postprocess_func(
             trial: int,
@@ -676,7 +698,7 @@ class FemtetInterface(FEMInterface):
     ):
         result_dir = original_femprj_path.replace('.femprj', '.Results')
         if pdt_file_content is not None:
-            pdt_path = os.path.join(result_dir, model_name + f'_trial{trial}.pdt')
+            pdt_path = FemtetInterface.create_pdt_path(original_femprj_path, model_name, trial)
             with open(pdt_path, 'wb') as f:
                 f.write(pdt_file_content)
 
@@ -710,13 +732,12 @@ class FemtetInterface(FEMInterface):
         jpg_path = os.path.join(result_dir, self.model_name + '.jpg')
 
         # モデル表示画面の設定
-        self.Femtet.SetWindowSize(200, 200)
+        self.Femtet.SetWindowSize(600, 600)
         self.Femtet.Fit()
-        self.Femtet.ViewNumeric.SetCoord(1, 1, 1)
 
         # ---モデルの画面を保存---
         self.Femtet.Redraw()  # 再描画
-        succeed = self.Femtet.SavePicture(jpg_path, 200, 200, 80)
+        succeed = self.Femtet.SavePicture(jpg_path, 600, 600, 80)
 
         self.Femtet.RedrawMode = True  # 逐一の描画をオン
 
@@ -745,6 +766,7 @@ class _UnPicklableNoFEM(FemtetInterface):
     Femtet = None
     quit_when_destruct = False
 
+    # noinspection PyMissingConstructor
     def __init__(self):
         CoInitialize()
         self.unpicklable_member = Dispatch('FemtetMacro.Femtet')

pyfemtet/opt/interface/_femtet_parametric.py

@@ -1,36 +1,15 @@
+from femtetutils import util
 from pyfemtet.dispatch_extensions import _get_pid
 
-import winreg
 import ctypes
 
 
-def _get_dll(Femtet):
-    femtet_major_version = _get_femtet_major_version(Femtet)
-    dll_path = _get_parametric_dll_path(femtet_major_version)
+def _get_dll():
+    femtet_exe_path = util.get_femtet_exe_path()
+    dll_path = femtet_exe_path.replace('Femtet.exe', 'ParametricIF.dll')
     return ctypes.cdll.LoadLibrary(dll_path)
 
 
-def _get_femtet_major_version(Femtet):
-    from pyfemtet.core import _version
-    version_int8 = _version(Femtet=Femtet)
-    return str(version_int8)[0:4]
-
-
-def _get_parametric_dll_path(femtet_major_version) -> str:
-    key = winreg.OpenKey(winreg.HKEY_LOCAL_MACHINE, r'SOFTWARE\FemtetInfo\InstallVersion\x64')
-    _, nValues, _ = winreg.QueryInfoKey(key)
-    for i in range(nValues):
-        name, data, _ = winreg.EnumValue(key, i)
-        if name == str(femtet_major_version):
-            winreg.CloseKey(key)
-            import os
-            dll_path = os.path.join(data, 'Program', 'ParametricIF.dll')
-            return dll_path
-    # ここまで来ていたら失敗
-    winreg.CloseKey(key)
-    raise RuntimeError('パラメトリック解析機能へのアクセスに失敗しました')
-
-
 def _get_dll_with_set_femtet(Femtet):
     dll = _get_dll(Femtet)
     pid = _get_pid(Femtet.hWnd)
@@ -40,6 +19,7 @@ def _get_dll_with_set_femtet(Femtet):
 
 
 def _get_prm_result_names(Femtet):
+    """Used by pyfemtet-opt-gui"""
     out = []
 
     # load dll and set target femtet

pyfemtet/opt/opt/_base.py

@@ -65,10 +65,15 @@ class AbstractOptimizer(ABC):
 
         # x の更新
         self.parameters['value'] = x
+        logger.info(f'Start calculation with input: {x}')
 
         # FEM の更新
         logger.debug('fem.update() start')
-        self.fem.update(self.parameters)
+        try:
+            self.fem.update(self.parameters)
+        except Exception as e:
+            logger.warning('An exception has occurred during FEM update.')
+            raise e
 
         # y, _y, c の更新
         logger.debug('calculate y start')
@@ -93,6 +98,9 @@ class AbstractOptimizer(ABC):
         )
 
         logger.debug('history.record end')
+
+        logger.info(f'End calculation with output: {_y}')
+
         return np.array(y), np.array(_y), np.array(c)
 
     def _reconstruct_fem(self, skip_reconstruct=False):

pyfemtet/opt/opt/_optuna.py

@@ -18,6 +18,10 @@ from pyfemtet.core import MeshError, ModelError, SolveError
 # filter warnings
 import warnings
 from optuna.exceptions import ExperimentalWarning
+
+
+optuna.logging.set_verbosity(optuna.logging.ERROR)
+
 warnings.filterwarnings('ignore', category=ExperimentalWarning)
 
 

pyfemtet/opt/visualization/__init__.py

@@ -1,7 +0,0 @@
-from pyfemtet.opt.visualization._monitor import ResultViewerApp
-from pyfemtet.opt.visualization._monitor import ProcessMonitorApp
-
-__all__ = [
-    'ResultViewerApp',
-    'ProcessMonitorApp',
-]

pyfemtet/opt/visualization/_create_wrapped_components.py

@@ -0,0 +1,93 @@
+"""Create autocompletable components"""
+import os
+import inspect
+
+from dash.development.base_component import Component
+
+
+# noinspection PyUnresolvedReferences
+from dash import html, dcc, dash_table
+# noinspection PyUnresolvedReferences
+import dash_bootstrap_components as dbc
+
+
+here, me = os.path.split(__file__)
+COMPONENT_FILE_DIR = os.path.join(here, 'wrapped_components')
+indent = '    '
+
+
+def create(module_name: str) -> str:
+    header = '''# auto created module
+from pyfemtet.opt.visualization.wrapped_components.str_enum import StrEnum
+# from enum import StrEnum
+import dash
+import dash_bootstrap_components
+
+
+'''
+    path = os.path.join(COMPONENT_FILE_DIR, module_name.replace('.py', '') + '.py')
+    with open(path, 'w', newline='\n') as f:
+        f.write(header)
+    return path
+
+
+def append(html_class, module_path: str):
+    print('==========')
+    print(html_class)
+
+    # get property names
+    init_signature = inspect.signature(html_class.__init__)
+    props = [param.name for param in init_signature.parameters.values() if (param.name != 'self') and (param.name != 'args') and (param.name != 'kwargs')]
+
+    # create class definition
+    class_definition = f'class {html_class.__name__}({html_class.__module__}):\n'  # library specific
+
+    # create id property
+    class_definition += indent + 'def _dummy(self):\n'
+    class_definition += indent*2 + '# noinspection PyAttributeOutsideInit\n'
+    class_definition += indent*2 + 'self.id = None\n\n'
+
+    # create Prop attribute
+    class_definition += indent + 'class Prop(StrEnum):\n'
+    for available_property in props:
+        property_definition = f'{available_property} = "{available_property}"'
+        try:
+            exec(property_definition)
+        except (SyntaxError, TypeError):
+            continue
+        class_definition += indent*2 + f'{property_definition}\n'
+
+    if class_definition.endswith(':\n'):
+        return
+
+    class_definition += '\n\n'
+
+    with open(module_path, 'a', newline='\n') as f:
+        f.write(class_definition)
+
+
+def sub(module_name):
+    # glob component classes
+    module_path = create(module_name)
+    class_names = dir(eval(module_name))
+    for class_name in class_names:
+        cls = eval(f'{module_name}.{class_name}')
+        if not inspect.isclass(cls):
+            continue
+        if issubclass(cls, Component):
+            append(cls, module_path)
+
+
+def main():
+    # html
+    sub('html')
+
+    # dcc
+    sub('dcc')
+
+    # dbc
+    sub('dbc')
+
+
+if __name__ == '__main__':
+    main()