pyfemtet 0.4.25__py3-none-any.whl → 0.5.1__py3-none-any.whl

pyfemtet-0.5.1.dist-info/RECORD

@@ -0,0 +1,115 @@
+pyfemtet/__init__.py,sha256=PQu9BZ4UZrI--6Rc-s0RhFPPEWVTqddrgsWkwYl87dI,21
+pyfemtet/core.py,sha256=3lqfBGJ5IuKz2Nqj5pRo7YQqKwx_0ZDL72u95Ur_1p0,1386
+pyfemtet/dispatch_extensions.py,sha256=XVZajbjh7mb6NG4Hq8qff2TJWab75r4Hd59cIvCRsVg,16213
+pyfemtet/logger.py,sha256=JYD0FvzijMS2NvZN7VT7vZA5hqtHEkvS93AHlIMDePw,2507
+pyfemtet/message/1. make_pot.bat,sha256=oS38xYsaUnQAuKwUR8hZJFgt3AKBU993fWFDSg2ROz4,570
+pyfemtet/message/2. make_mo.bat,sha256=nqUi3Cze7JGKkYItlch8ZG2gSbRNZiS2ltuCS7DbmG8,154
+pyfemtet/message/__init__.py,sha256=gE1-XX_PzHj9BbhqPaK5VcIHuv6_Tec5qlPMC3IRiBg,100
+pyfemtet/message/babel.cfg,sha256=rlvsm_EjDU2ki-OLVPG3toWspVuVA0JFAFvFp2U-By4,72
+pyfemtet/message/locales/ja/LC_MESSAGES/messages.mo,sha256=hbW3KO524_lFYIhU-QGMhHMW4LHnxyecnptFD55ECMo,18332
+pyfemtet/message/locales/ja/LC_MESSAGES/messages.po,sha256=f8rimJEqmLdlrE7JBwETm-Gq2-USOZcSEpnzMqBqk88,23617
+pyfemtet/message/locales/messages.pot,sha256=_gVMy6qaSU7Qu6zIY8KQizb1XpTZ0ijLKbnnk4lAAdg,13635
+pyfemtet/message/messages.py,sha256=i0cOOV8NnBeT1gG-tqgnu_3dIym1hZG8Hzg38sh_Qeo,13187
+pyfemtet/opt/__init__.py,sha256=MPrUWeLZLrJ-ApVckn8dsn3QmRH13aPzit5JgaoshG8,696
+pyfemtet/opt/_femopt.py,sha256=V1D0n8WgzDRE6NbvRzJdhxXiqHh7GGRAmvS7jyQDGfA,26611
+pyfemtet/opt/_femopt_core.py,sha256=ocn7qOMCDRXuS-3DRXgGIo1hGnewUYSUxE7BuvqCXMY,27858
+pyfemtet/opt/_test_utils/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+pyfemtet/opt/_test_utils/control_femtet.py,sha256=Oy2MmNS-LhUXF9rKLa8AXAfJhppIQI8Nha8LmEZflmk,1169
+pyfemtet/opt/_test_utils/hyper_sphere.py,sha256=nQhw8EIY0DwvcTqrbKhkxiITLZifr4-nG77E-_6ggmA,700
+pyfemtet/opt/_test_utils/record_history.py,sha256=qg5kNURiFHr40YN4f4s9Q6JNed-1y5ysQFXwIbYpKZw,2512
+pyfemtet/opt/interface/__init__.py,sha256=qz5BszPuU3jZIoDnPjkPDAgvgHLlx1sYhuqh5ID798k,480
+pyfemtet/opt/interface/_base.py,sha256=I4pJttLeRW-6WWMuCNynwxWPwriiGZk20vHLVcfixZY,2332
+pyfemtet/opt/interface/_femtet.py,sha256=PjmwRpkY99i4Tstva6ahlG8beqQytkJsAbT9nU_d64A,35956
+pyfemtet/opt/interface/_femtet_parametric.py,sha256=KDG8SB43AgwuhpCStjvx10G0RzyHhga6k4dfvp0gvYU,2175
+pyfemtet/opt/interface/_femtet_with_nx/__init__.py,sha256=-6W2g2FDEcKzGHmI5KAKQe-4U5jDpMj0CXuma-GZca0,83
+pyfemtet/opt/interface/_femtet_with_nx/_interface.py,sha256=cibPOhsUAvTDJa1xdEhVClZjwPkEx2kBmi61vyTJHUs,5514
+pyfemtet/opt/interface/_femtet_with_nx/update_model.py,sha256=P7VH0i_o-X9OUe6AGaLF1fACPeHNrMjcrOBCA3MMrI4,3092
+pyfemtet/opt/interface/_femtet_with_sldworks.py,sha256=Ldr8Esa5xZ-D_E5uIXBTF2DHslXVMFDYOsqTd8FhY1M,6242
+pyfemtet/opt/optimizer/__init__.py,sha256=wdz7PXkcSJ9Z2OHoegfGvYrnj_OAyU39BvDyiy3QnVw,407
+pyfemtet/opt/optimizer/_base.py,sha256=lKCDGoAtOLwqNaP7qfqYJ69AHxq9YZVivTP14VHYxZU,12136
+pyfemtet/opt/optimizer/_optuna.py,sha256=MojfB0oQCjSzS-hJNiD_Wn-vPDJbsMUydMq2ht6Ww9g,11570
+pyfemtet/opt/optimizer/_optuna_botorchsampler_parameter_constraint_helper.py,sha256=6um5tdhpMirzwAPmp_K_YluguP0W9dWPwHJIqWICeOM,11908
+pyfemtet/opt/optimizer/_scipy.py,sha256=M3q25VuRazDzXUWC25z7hGbNq2Qg9158NjTSNBvi2pk,4347
+pyfemtet/opt/optimizer/_scipy_scalar.py,sha256=YmErX4s_ywVIn861-D06d-TvDWi0rp843Jit3x80qtU,3202
+pyfemtet/opt/parameter.py,sha256=YLE9lmYRaZA8isnTPJnbYXpUn6zsJFW4xg03QaSWey8,3950
+pyfemtet/opt/prediction/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+pyfemtet/opt/prediction/base.py,sha256=q4lDqrOtAkiWG-OblQEddnEVj29Q-EJE3-O5QTOb3Q4,1761
+pyfemtet/opt/prediction/single_task_gp.py,sha256=VbsVllzXzCGqkM1fC61Ouqwuv3ddS5INbFwDG8v-d6g,3303
+pyfemtet/opt/samples/femprj_sample/.gitignore,sha256=hx-5Hhaf7kpHe1wvWWfJqjPfObg-zf9CTI4joNh2Hk4,28
+pyfemtet/opt/samples/femprj_sample/ParametricIF.femprj,sha256=9BtDHmc3cdom0Zq33DTdZ0mDAsIUY6i8SRkkg-n7GO0,442090
+pyfemtet/opt/samples/femprj_sample/ParametricIF.py,sha256=oXzchBZEbH69xacDht5HDnbZzKwapXsn6bp9qihY17Y,707
+pyfemtet/opt/samples/femprj_sample/ParametricIF_test_result.reccsv,sha256=TiOAqEDMub6SCGYClBv1JvQxphDOY3iIdr_pMmGgJ9M,2859
+pyfemtet/opt/samples/femprj_sample/cad_ex01_NX.femprj,sha256=KC8JlHqHzqgyKriK911QSnQByQpRlw-SX5OSQ_GNe5M,149193
+pyfemtet/opt/samples/femprj_sample/cad_ex01_NX.prt,sha256=3okHLeMdslrRA_wkhppZtxIe-2-ZPMfNqWCdQwUV31o,226626
+pyfemtet/opt/samples/femprj_sample/cad_ex01_NX.py,sha256=_wrpHQ6NlNJ9js0PwR3tfuxXt5fxORTa9h8BvjiWOcg,4621
+pyfemtet/opt/samples/femprj_sample/cad_ex01_NX_test_result.reccsv,sha256=hFxjat32G71F3XE_UtC9zzi5HZxu0Q0MXjU9e1xVIDM,4710
+pyfemtet/opt/samples/femprj_sample/cad_ex01_SW.SLDPRT,sha256=jjBi4aRRwZPK-4-YRKDC4eO_Ch2jwd7M7chvJlnBbZU,97158
+pyfemtet/opt/samples/femprj_sample/cad_ex01_SW.femprj,sha256=knN0bBTHm5CqExLdmxdJvPldJ6ahnQesKt974qRjWh4,126837
+pyfemtet/opt/samples/femprj_sample/cad_ex01_SW.py,sha256=JoGQSwH3yJnABxyd-WJfrwMkhd1UV0yYF2L2RvMFXmc,4559
+pyfemtet/opt/samples/femprj_sample/cad_ex01_SW_test_result.reccsv,sha256=NS0Zik2c1mbMdGa0hGJaRQdCD08Bltx84n9QzP5CjPo,4736
+pyfemtet/opt/samples/femprj_sample/constrained_pipe.femprj,sha256=MAl-VQfethwYvl49RkuW7FQlFCQ9_mYvc03SsqBCad0,57414
+pyfemtet/opt/samples/femprj_sample/constrained_pipe.py,sha256=-9LYdLQ4MEsdt9FAzaqhtRPwBQN5mnsXcOMyNqsv5Iw,3217
+pyfemtet/opt/samples/femprj_sample/constrained_pipe_test_result.reccsv,sha256=AbAMXLtxrEBHLjdM0HnGQ7j3uJtVHZpxOZxOqNmO1CA,1796
+pyfemtet/opt/samples/femprj_sample/gal_ex58_parametric.femprj,sha256=dbanN3W2eX_ciZ0wZGqK60mc4edSVh5G2OqbbMKyFng,74898
+pyfemtet/opt/samples/femprj_sample/gal_ex58_parametric.py,sha256=RxLSXFs0SqUjlug_JZAKlkJhqJdQCY3Y3F-DtSQRtVM,2458
+pyfemtet/opt/samples/femprj_sample/gal_ex58_parametric_test_result.reccsv,sha256=NOaYmpmrhn9WMbIaLWFlR0IvRheGqcqu0J9Nf3JQnfo,1131
+pyfemtet/opt/samples/femprj_sample/gau_ex08_parametric.femprj,sha256=Yb9ILeTEKx5xfJGk8IZH_DVlgkpGB33Vy9-LGIEQboY,279251
+pyfemtet/opt/samples/femprj_sample/gau_ex08_parametric.py,sha256=5KazqJ5wRbGs0dBMJslZ1eRCUWq8j3k1mqlhyB8M0g8,1929
+pyfemtet/opt/samples/femprj_sample/gau_ex08_parametric_test_result.reccsv,sha256=yZ9aHthiKIBY_NMOz94Jl2dyHIH-GWMvukgHk4ZeT_o,3474
+pyfemtet/opt/samples/femprj_sample/her_ex40_parametric.femprj,sha256=LLAUDlUo1dIpRzlKPs1lvACzJQxjnWW3xAGAodYEqRM,117221
+pyfemtet/opt/samples/femprj_sample/her_ex40_parametric.py,sha256=47Lqul9LDlbMB2Bvcn_B6wzrdNvQp1ksYSDRUZfR76Q,4846
+pyfemtet/opt/samples/femprj_sample/her_ex40_parametric_test_result.reccsv,sha256=0yWqTpmpAtFvYRRyk2zneAVnl_5qJDeVwG4aeIWxXv8,3679
+pyfemtet/opt/samples/femprj_sample/paswat_ex1_parametric.femprj,sha256=z2NHFJWiuiLv_zhxjFpLpmRbYVvQ43bAMj_NLioQGsk,262283
+pyfemtet/opt/samples/femprj_sample/paswat_ex1_parametric.py,sha256=2oUBg2MvKuewO5Acs5FO3dkXX0QIZGRR5RRvKAb4CMk,2361
+pyfemtet/opt/samples/femprj_sample/paswat_ex1_parametric_parallel.py,sha256=V6XYS3RblruTi8T43lEnNl_3pgdsobknspsD_K2iEpU,2472
+pyfemtet/opt/samples/femprj_sample/paswat_ex1_parametric_test_result.reccsv,sha256=6PjbmLfwDx01Z0pywDL9eT5lfRKgZDU3r_H2SfZYaT4,3885
+pyfemtet/opt/samples/femprj_sample/wat_ex14_parametric.femprj,sha256=F-yu2dGrsbrIA1Lhizu2aHTjQFTohyBmOuJv-Iyl8jk,179596
+pyfemtet/opt/samples/femprj_sample/wat_ex14_parametric.py,sha256=K9ZlIkmD05NOF2dXnv4AmYsvM7klhm8Ss7irAr74MVQ,2213
+pyfemtet/opt/samples/femprj_sample/wat_ex14_parametric_parallel.py,sha256=UjPMpoZY99rJY9mqudsK7MOGURhzwOulwDkrRiclZzI,2289
+pyfemtet/opt/samples/femprj_sample/wat_ex14_parametric_test_result.reccsv,sha256=bvHgcuFZ1iS8sb_FT5W7MuquODxmv3IcFkiNOdSMiK8,3450
+pyfemtet/opt/samples/femprj_sample_jp/ParametricIF_jp.femprj,sha256=9BtDHmc3cdom0Zq33DTdZ0mDAsIUY6i8SRkkg-n7GO0,442090
+pyfemtet/opt/samples/femprj_sample_jp/ParametricIF_jp.py,sha256=oXzchBZEbH69xacDht5HDnbZzKwapXsn6bp9qihY17Y,707
+pyfemtet/opt/samples/femprj_sample_jp/cad_ex01_NX_jp.femprj,sha256=0RBhOGhtiFAp0QSCTBYEaDY9EZymn9hJYchAOJ6PaBA,143533
+pyfemtet/opt/samples/femprj_sample_jp/cad_ex01_NX_jp.py,sha256=B7wVemkiVzh0NJXDzG3fdigWBOb58ZxrJZUT4NRvW9Q,4899
+pyfemtet/opt/samples/femprj_sample_jp/cad_ex01_SW_jp.femprj,sha256=ZZhT9XjB9Xu9YwHWv4gbvKBiUWlOFKEoHjAcGWb3vvQ,128026
+pyfemtet/opt/samples/femprj_sample_jp/cad_ex01_SW_jp.py,sha256=e-XaHBZ8syWlG9pObcDDkHGLC5t338zAt_NyRXty338,4837
+pyfemtet/opt/samples/femprj_sample_jp/constrained_pipe_jp.py,sha256=QhAwJOFEknf6Yk3mMgr1gdqB_Db8akjJGZNBepsrdgA,3654
+pyfemtet/opt/samples/femprj_sample_jp/gal_ex58_parametric_jp.femprj,sha256=PzqtNVde6VnMHFpedRBsOq9JVhCY0ymQPVs54EKsNLw,75668
+pyfemtet/opt/samples/femprj_sample_jp/gal_ex58_parametric_jp.py,sha256=97np4uH-UQqpv4UDwJS0doFYA7TOkXnbhLdkZExdNek,2461
+pyfemtet/opt/samples/femprj_sample_jp/gau_ex08_parametric_jp.femprj,sha256=TTXw_8YT8pzHQlu4ufGzTq1IFYSwcWWt4GA6sIY1YPM,295600
+pyfemtet/opt/samples/femprj_sample_jp/gau_ex08_parametric_jp.py,sha256=8Op_zwz9SD0NfGg4TFlcNvs-ZlU0bxgs5oaaI9UtlRU,2087
+pyfemtet/opt/samples/femprj_sample_jp/her_ex40_parametric_jp.femprj,sha256=OJ7f8iw0z1BZqanuNn71uEaoM2Kgb93ptUU8iYwYON0,129783
+pyfemtet/opt/samples/femprj_sample_jp/her_ex40_parametric_jp.py,sha256=MlFk6KRCQeCX1J0DNOjph75qjCUHg5UQPNTcHxIEnoo,5279
+pyfemtet/opt/samples/femprj_sample_jp/paswat_ex1_parametric_jp.femprj,sha256=y7eURFBdqh6PmD4zbelGuB458HmfihVht0K4wVI-mik,265368
+pyfemtet/opt/samples/femprj_sample_jp/paswat_ex1_parametric_jp.py,sha256=ibXVhya5JV7CfCD2srQceXOj1sYjH3kWtaIEfyu98u4,2489
+pyfemtet/opt/samples/femprj_sample_jp/paswat_ex1_parametric_parallel_jp.py,sha256=cY9RhoAFEnVUkfGhbXpn3LavT_NRp_OwVxGbL7jSbBQ,2605
+pyfemtet/opt/samples/femprj_sample_jp/wat_ex14_parametric_jp.femprj,sha256=dMwQMt6yok_PbZLyxPYdmg5wJQwgQDZ4RhS76zdGLGk,177944
+pyfemtet/opt/samples/femprj_sample_jp/wat_ex14_parametric_jp.py,sha256=vMy-KUP1wEMV9Rt6yXjkE40Fcs1t1cpQK-nQJK8hHao,2284
+pyfemtet/opt/samples/femprj_sample_jp/wat_ex14_parametric_parallel_jp.py,sha256=4X0cl3YWpYarcNBCH79mrlyFuKUYSqvnGzokEbv9ILk,2335
+pyfemtet/opt/visualization/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+pyfemtet/opt/visualization/_create_wrapped_components.py,sha256=XpOv_iEEB3lFT3vYuCnCCbwUwRFbCYgvsZDpaiSzmZY,2511
+pyfemtet/opt/visualization/base.py,sha256=XhZodLEr5PPHmtNai6uZcPJobKurENepuGCUaqDOI7I,7588
+pyfemtet/opt/visualization/complex_components/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+pyfemtet/opt/visualization/complex_components/alert_region.py,sha256=IIsfw4Dwlag7ked94X7MF2UxEODJCBeYruDOUQsKzGQ,2104
+pyfemtet/opt/visualization/complex_components/control_femtet.py,sha256=b4gn_pQMXe3NL3aqfEgsPTvhmdWmg-tRBvaHG659lZg,6240
+pyfemtet/opt/visualization/complex_components/main_figure_creator.py,sha256=ueVyDFnY_KPIXuLM585bRToAGDCUXmd-cy_ZNMSFWs4,6949
+pyfemtet/opt/visualization/complex_components/main_graph.py,sha256=iPxH7XdE0vWByzZfVka_93mLM70svv20OTpFMwAB-UE,21528
+pyfemtet/opt/visualization/complex_components/pm_graph.py,sha256=n1p-yruqo9NW93UOUIdxPF2QdO4qDSuH37h6l8C68BQ,24921
+pyfemtet/opt/visualization/complex_components/pm_graph_creator.py,sha256=ukv9XqW_Mzdsw9Oxc3JXLwKJTyXPSrt9DiB3dvKv4uk,7302
+pyfemtet/opt/visualization/process_monitor/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+pyfemtet/opt/visualization/process_monitor/application.py,sha256=ir0xVTsbqFWq2FlPw75XNcMpFEUQNUQr5OcAW0r-FuY,7922
+pyfemtet/opt/visualization/process_monitor/pages.py,sha256=HSjBeZLnVEgm-VP_NJIpjRUsGhphNuV4eMqK7xC77es,15110
+pyfemtet/opt/visualization/result_viewer/.gitignore,sha256=ryvb4aqbbsHireHWlPQfxxqDHQJo6YkVYhE9imKt0b8,6
+pyfemtet/opt/visualization/result_viewer/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+pyfemtet/opt/visualization/result_viewer/application.py,sha256=99Zn-i2mJLIN0DeKjewCBW9UCa78YpyRoMLXzbMd4Qg,1687
+pyfemtet/opt/visualization/result_viewer/pages.py,sha256=wzTbd-3XHDZ0rK98sE55tXQ90nb1IbUhav3sIVFfjOY,32147
+pyfemtet/opt/visualization/wrapped_components/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+pyfemtet/opt/visualization/wrapped_components/dbc.py,sha256=wzR1ZMOb4uwPNTMFn5upLJ9tXeuJ22mLGynyJlA2cho,42161
+pyfemtet/opt/visualization/wrapped_components/dcc.py,sha256=hcW7SR6VIMn4S4-JMyohvOzdc0Aw8A4chIeHqQEUbFU,17499
+pyfemtet/opt/visualization/wrapped_components/html.py,sha256=sE2XHTDY1GvA1NW7y6SKWf-WglVXFKKvXhU9h3z53_g,95652
+pyfemtet/opt/visualization/wrapped_components/str_enum.py,sha256=NZqbh2jNEAckvJyZv__MWeRs2F2Q-dkJCWo30rU2rrM,1383
+pyfemtet-0.5.1.dist-info/LICENSE,sha256=sVQBhyoglGJUu65-BP3iR6ujORI6YgEU2Qm-V4fGlOA,1485
+pyfemtet-0.5.1.dist-info/METADATA,sha256=PLu9EtBr3EFxUh2vMonVk6klrI_tyaS_h9hfXO9qfRE,3287
+pyfemtet-0.5.1.dist-info/WHEEL,sha256=sP946D7jFCHeNz5Iq4fL4Lu-PrWrFsgfLXbbkciIZwg,88
+pyfemtet-0.5.1.dist-info/entry_points.txt,sha256=ZfYqRaoiPtuWqFi2_msccyrVF0LurMn-IHlYamAegZo,104
+pyfemtet-0.5.1.dist-info/RECORD,,

pyfemtet/_test_util.py

@@ -1,135 +0,0 @@
-import os
-import csv
-from shutil import copy
-from time import sleep
-from subprocess import run
-from multiprocessing import Process
-from glob import glob
-
-import numpy as np
-import pandas as pd
-from tqdm import tqdm
-from win32com.client import Dispatch
-from femtetutils import util
-
-from pyfemtet.opt import FEMOpt
-
-
-class SuperSphere:
-    def __init__(self, n):
-        self.n = n
-
-    def x(self, radius, *angles):
-        assert len(angles) == self.n - 1, 'invalid angles length'
-
-        out = []
-
-        for i in range(self.n):
-            if i == 0:
-                out.append(radius * np.cos(angles[0]))
-            elif i < self.n - 1:
-                product = radius
-                for j in range(i):
-                    product *= np.sin(angles[j])
-                product *= np.cos(angles[i])
-                out.append(product)
-            else:
-                product = radius
-                for j in range(i):
-                    product *= np.sin(angles[j])
-                out.append(product)
-        return out
-
-
-def _open_femprj(femprj_path):
-    Femtet = Dispatch('FemtetMacro.Femtet')
-    for _ in tqdm(range(5), 'wait for dispatch Femtet'):
-        sleep(1)
-    Femtet.LoadProject(femprj_path, True)
-
-
-def launch_femtet(femprj_path):
-    # launch Femtet externally
-    print('Launching Femtet...')
-    util.execute_femtet()
-    pid = util.get_last_executed_femtet_process_id()
-    for _ in tqdm(range(8), 'Wait for launch Femtet.'):
-        sleep(1)
-
-    # open femprj in a different process
-    # to release Femtet for sample program
-    print('Opening femprj...')
-    if femprj_path:
-        p = Process(
-            target=_open_femprj,
-            args=(femprj_path,),
-        )
-        p.start()
-        p.join()
-
-
-def taskkill_femtet():
-    for _ in tqdm(range(3), 'wait before taskkill Femtet'):
-        sleep(1)
-    run(['taskkill', '/f', '/im', 'Femtet.exe'])
-    for _ in tqdm(range(3), 'wait after taskkill Femtet'):
-        sleep(1)
-
-
-def find_latest_csv(dir_path=None):
-    if dir_path is None: dir_path = '.'
-    target = os.path.join(dir_path, '*.csv')
-    files = [(f, os.path.getmtime(f)) for f in glob(target)]
-    out = sorted(files, key=lambda files: files[1])[-1]
-    return os.path.abspath(out[0])
-
-
-def py_to_reccsv(py_path, suffix=''):
-    dst_csv_path = py_path + suffix
-    dst_csv_path = dst_csv_path.replace(f'.py{suffix}', f'{suffix}.reccsv')
-    return dst_csv_path
-
-
-def record_result(src: FEMOpt or str, py_path, suffix=''):
-    """Record the result csv for `is_equal_result`."""
-
-    if isinstance(src, FEMOpt):  # get df directory
-        src_csv_path = src.history_path
-    else:
-        src_csv_path = os.path.abspath(src)
-
-    dst_csv_path = py_to_reccsv(py_path, suffix)
-    copy(src_csv_path, dst_csv_path)
-
-
-def _get_obj_from_csv(csv_path, encoding='cp932'):
-    df = pd.read_csv(csv_path, encoding=encoding, header=2)
-    columns = df.columns
-    with open(csv_path, mode='r', encoding=encoding, newline='\n') as f:
-        reader = csv.reader(f, delimiter=',')
-        meta = reader.__next__()
-    obj_indices = np.where(np.array(meta) == 'obj')[0]
-    out = df.iloc[:, obj_indices]
-    return out, columns
-
-
-def is_equal_result(ref_path, dif_path, log_path):
-    """Check the equality of two result csv files."""
-    ref_df, ref_columns = _get_obj_from_csv(ref_path)
-    dif_df, dif_columns = _get_obj_from_csv(dif_path)
-
-    with open(log_path, 'a', newline='\n') as f:
-        f.write('\n\n===== 結果の分析 =====\n\n')
-        f.write(f'        \tref\tdif\n')
-        f.write(f'---------------------\n')
-        f.write(f'len(col)\t{len(ref_columns)}\t{len(dif_columns)}\n')
-        f.write(f'len(df) \t{len(ref_df)}\t{len(dif_df)}\n')
-        try:
-            difference = (np.abs(ref_df.values - dif_df.values) / np.abs(dif_df.values)).mean()
-            f.write(f'diff    \t{int(difference*100)}%\n')
-        except Exception:
-            f.write(f'diff    \tcannot calc\n')
-
-    assert len(ref_columns) == len(dif_columns), '結果 csv の column 数が異なります。'
-    assert len(ref_df) == len(dif_df), '結果 csv の row 数が異なります。'
-    assert difference <= 0.05, '前回の結果との平均差異が 5% を超えています。'

pyfemtet/opt/femprj_sample/her_ex40_parametric_test_result.reccsv

@@ -1,18 +0,0 @@
-"{""femprj_path"": null, ""model_name"": ""Harm""}",prm,prm_lb,prm_ub,prm,prm_lb,prm_ub,prm,prm_lb,prm_ub,obj,obj_direction,,cns,cns_lb,cns_ub,cns,cns_lb,cns_ub,,,,
-,,,,,,,,,,,,,,,,,,,,,,
-trial,antenna_radius,antenna_radius_lower_bound,antenna_radius_upper_bound,substrate_w,substrate_w_lower_bound,substrate_w_upper_bound,port_x,port_x_lower_bound,port_x_upper_bound,first resonant frequency(Hz),first resonant frequency(Hz)_direction,non_domi,antenna and substrate clearance,antenna and substrate clearance_lower_bound,antenna and substrate clearance_upper_bound,antenna and port clearance,antenna and port clearance_lower_bound,antenna and port clearance_upper_bound,feasible,hypervolume,message,time
-1,10.0,5.0,20.0,50.0,40.0,60.0,5.0,1.0,20.0,4289000000.0,3000000000.0,False,15.0,1,,5.0,1,,True,-1.0,initial,2024-07-24 10:30:06.939869
-2,13.979877262955549,5.0,20.0,43.12037280884873,40.0,60.0,3.96389588638785,1.0,20.0,3050000000.0,3000000000.0,False,7.580309099999999,1,,10.01598141,1,,True,-1.0,,2024-07-24 10:30:19.275568
-3,17.48663961200633,5.0,20.0,44.246782213565524,40.0,60.0,4.4546743769349115,1.0,20.0,2460000000.0,3000000000.0,False,4.636751499999999,1,,13.03196522,1,,True,-1.0,,2024-07-24 10:30:35.798531
-4,11.841049763255539,5.0,20.0,55.70351922786027,40.0,60.0,4.793801861008835,1.0,20.0,3581000000.0,3000000000.0,False,16.0107098,1,,7.04724794,1,,True,-1.0,,2024-07-24 10:30:51.243382
-5,12.713516576204174,5.0,20.0,51.84829137724085,40.0,60.0,1.8825578416799567,1.0,20.0,2932000000.0,3000000000.0,False,13.2106291,1,,10.83095876,1,,True,-1.0,,2024-07-24 10:31:01.335664
-6,14.113172778521575,5.0,20.0,43.41048247374583,40.0,60.0,2.235980266720311,1.0,20.0,2991000000.0,3000000000.0,True,7.592068450000001,1,,11.877192529999999,1,,True,-1.0,,2024-07-24 10:31:14.045279
-7,19.233283058799998,5.0,20.0,59.31264066149119,40.0,60.0,16.35954961421276,1.0,20.0,2283000000.0,3000000000.0,False,10.42303725,1,,2.8737335,1,,True,-1.0,,2024-07-24 10:31:40.560920
-8,11.60228740609402,5.0,20.0,42.44076469689558,40.0,60.0,10.408361292114133,1.0,20.0,3758000000.0,3000000000.0,False,9.618094950000001,1,,1.1939261000000005,1,,True,-1.0,,2024-07-24 10:31:57.496815
-9,14.93783426530973,5.0,20.0,46.23422152178822,40.0,60.0,10.881292402378406,1.0,20.0,2932000000.0,3000000000.0,False,8.179276449999998,1,,4.056541900000001,1,,True,-1.0,,2024-07-24 10:32:18.322964
-10,13.968499682166277,5.0,20.0,58.437484700462335,40.0,60.0,2.6813575389864703,1.0,20.0,2991000000.0,3000000000.0,True,15.250242649999999,1,,11.28714216,1,,True,-1.0,,2024-07-24 10:32:32.049325
-11,14.058429837987106,5.0,20.0,51.01029232478569,40.0,60.0,1.738898798219626,1.0,20.0,2814000000.0,3000000000.0,False,11.44671635,1,,12.319531000000001,1,,True,-1.0,,2024-07-24 10:32:44.208324
-12,14.260807151712159,5.0,20.0,57.66392990209854,40.0,60.0,8.792627280104398,1.0,20.0,3050000000.0,3000000000.0,False,14.57115775,1,,5.468179920000001,1,,True,-1.0,,2024-07-24 10:33:18.982843
-13,13.508818923150182,5.0,20.0,59.99999999999999,40.0,60.0,4.887635502556277,1.0,20.0,3109000000.0,3000000000.0,False,16.4911811,1,,8.6211834,1,,True,-1.0,,2024-07-24 10:33:35.259275
-14,14.885966644225132,5.0,20.0,40.0,40.0,60.0,4.993357418116345,1.0,20.0,2814000000.0,3000000000.0,False,5.1140334,1,,9.892609180000001,1,,True,-1.0,,2024-07-24 10:34:06.248951
-15,19.999999999999996,5.0,20.0,59.99999999999999,40.0,60.0,1.0,1.0,20.0,4525000000.0,3000000000.0,False,10.0,1,,19.0,1,,True,-1.0,,2024-07-24 10:34:37.774576

pyfemtet/opt/opt/__init__.py

@@ -1,12 +0,0 @@
-from pyfemtet.opt.opt._base import AbstractOptimizer, logger, OptimizationMethodChecker
-from pyfemtet.opt.opt._optuna import OptunaOptimizer
-from pyfemtet.opt.opt._scipy import ScipyOptimizer
-from pyfemtet.opt.opt._scipy_scalar import ScipyScalarOptimizer
-
-__all__ = [
-    'ScipyScalarOptimizer',
-    'ScipyOptimizer',
-    'OptunaOptimizer',
-    'AbstractOptimizer',
-    'logger',
-]

pyfemtet/opt/opt/_optuna_botorch_helper.py

@@ -1,209 +0,0 @@
-from typing import Optional, List, Tuple, Callable
-from functools import partial
-import inspect
-
-import numpy as np
-import optuna.study
-import torch
-from torch import Tensor
-from botorch.optim.initializers import gen_batch_initial_conditions
-from botorch.utils.transforms import unnormalize
-from optuna.study import Study
-from botorch.acquisition import AcquisitionFunction
-
-from pyfemtet.opt.opt import AbstractOptimizer
-from pyfemtet.opt.parameter import ExpressionEvaluator
-
-# module to monkey patch
-import optuna_integration
-
-
-# モンキーパッチを実行するため、optimize_acqf の引数を MonkyPatch クラスで定義し optuna に上書きされないようにするためのクラス
-class NonOverwritablePartial(partial):
-    def __call__(self, /, *args, **keywords):
-        stored_kwargs = self.keywords
-        keywords.update(stored_kwargs)
-        return self.func(*self.args, *args, **keywords)
-
-
-# prm_name を引数に取る関数を optimize_acqf の nonlinear_inequality_constraints に入れられる形に変換する関数
-class ConvertedConstraintFunction:
-    def __init__(self, fun, prm_args, kwargs, variables: ExpressionEvaluator, study: optuna.study.Study):
-        self.fun = fun
-        self.prm_args = prm_args
-        self.kwargs = kwargs
-        self.variables = variables
-        self.study = study
-
-        self.bounds = None
-        self.prm_name_seq = None
-
-        # fun の prm として使う引数が指定されていなければ fun の引数を取得
-        if self.prm_args is None:
-            signature = inspect.signature(fun)
-            prm_inputs = set([a.name for a in signature.parameters.values()])
-        else:
-            prm_inputs = set(self.prm_args)
-
-        # 引数の set から kwargs の key を削除
-        self.prm_arg_names = prm_inputs - set(kwargs.keys())
-
-        # 変な引数が残っていないか確認
-        assert all([(arg in variables.get_parameter_names()) for arg in self.prm_arg_names])
-
-    def __call__(self, x: Tensor or np.ndarray):
-        # x: all of normalized parameters whose sequence is sorted by optuna
-
-        if not isinstance(x, Tensor):
-            x = torch.tensor(np.array(x)).double()
-
-        x = unnormalize(x, self.bounds)
-
-        # fun で使うパラメータのみ value を取得
-        kwargs = self.kwargs
-        kwargs.update(
-            {k: v for k, v in zip(self.prm_name_seq, x) if k in self.prm_arg_names}
-        )
-
-        return self.fun(**kwargs)
-
-
-# 与えられた獲得関数に拘束を満たさない場合 0 を返すよう加工された獲得関数
-class AcqWithConstraint(AcquisitionFunction):
-
-    # noinspection PyAttributeOutsideInit
-    def set(self, _org_acq_function: AcquisitionFunction, nonlinear_constraints):
-        self._org_acq_function = _org_acq_function
-        self._nonlinear_constraints = nonlinear_constraints
-
-    def forward(self, X: Tensor) -> Tensor:
-        base = self._org_acq_function.forward(X)
-
-        is_feasible = all([cons(X[0][0]) > 0 for cons, _ in self._nonlinear_constraints])
-        if is_feasible:
-            return base
-        else:
-            # penalty = torch.Tensor(size=base.shape)
-            # penalty = torch.fill(penalty, -1e10)
-            # return base * penalty
-            return base * 0.
-
-
-def remove_infeasible(_ic_batch, nonlinear_constraints):
-    # infeasible なものを削除
-    remove_indices = []
-    for i, ic in enumerate(_ic_batch):  # ic: 1 x len(params) tensor
-        # cons: Callable[["Tensor"], "Tensor"]
-        is_feasible = all([cons(ic[0]) > 0 for cons, _ in nonlinear_constraints])
-        if not is_feasible:
-            # ic_batch[i] = torch.nan  # これで無視にならない
-            remove_indices.append(i)
-    for i in remove_indices[::-1]:
-        _ic_batch = torch.cat((_ic_batch[:i], _ic_batch[i + 1:]))
-    return _ic_batch
-
-
-class OptunaBotorchWithParameterConstraintMonkeyPatch:
-
-    def __init__(self, study: Study, opt: AbstractOptimizer):
-        self.num_restarts: int = 20
-        self.raw_samples_additional: int = 512
-        self.eta: float = 2.0
-        self.study = study
-        self.opt = opt
-        self.nonlinear_inequality_constraints = []
-        self.additional_kwargs = dict()
-        self.bounds = None
-        self.prm_name_seq = None
-
-    def add_nonlinear_constraint(self, fun, prm_args, kwargs):
-        f = ConvertedConstraintFunction(
-            fun,
-            prm_args,
-            kwargs,
-            self.opt.variables,
-            self.study,
-        )
-
-        # 初期化
-        self.nonlinear_inequality_constraints = self.nonlinear_inequality_constraints or []
-
-        # 自身に追加
-        self.nonlinear_inequality_constraints.append((f, True))
-
-        # optimize_acqf() に渡す引数に追加
-        self.additional_kwargs.update(
-            nonlinear_inequality_constraints=self.nonlinear_inequality_constraints
-        )
-
-    def _detect_prm_seq_if_needed(self):
-        # study から distribution の情報を復元する。
-        if self.bounds is None or self.prm_name_seq is None:
-            from optuna._transform import _transform_search_space
-            # sample_relative の後に呼ばれているから最後の trial は search_space を持つはず
-            search_space: dict = self.study.sampler.infer_relative_search_space(self.study, self.study.trials[-1])
-            self.bounds = _transform_search_space(search_space, False, False)[0].T
-            self.prm_name_seq = list(search_space.keys())
-
-            for cns in self.nonlinear_inequality_constraints:
-                cns[0].bounds = torch.tensor(self.bounds)
-                cns[0].prm_name_seq = self.prm_name_seq
-
-    def generate_initial_conditions(self, *args, **kwargs):
-        self._detect_prm_seq_if_needed()
-
-        # acqf_function を 上書きし、拘束を満たさないならば 0 を返すようにする
-        org_acq_function = kwargs['acq_function']
-        new_acqf = AcqWithConstraint(None)
-        new_acqf.set(org_acq_function, self.nonlinear_inequality_constraints)
-        kwargs['acq_function'] = new_acqf
-
-        # initial condition の提案 batch を作成
-        # ic: `num_restarts x q x d` tensor of initial conditions.
-        # q = 1, d = len(params)
-        ic_batch = gen_batch_initial_conditions(*args, **kwargs)
-
-        # 拘束を満たさないものを削除
-        ic_batch = remove_infeasible(ic_batch, self.nonlinear_inequality_constraints)
-
-        # 全部なくなっているならばランダムに生成
-        if len(ic_batch) == 0:
-            print('拘束を満たす組み合わせがなかったのでランダムサンプリングします')
-        while len(ic_batch) == 0:
-            size = ic_batch.shape
-            ic_batch = torch.rand(size=[100, *size[1:]])  # 正規化された変数の組合せ
-            ic_batch = remove_infeasible(ic_batch, self.nonlinear_inequality_constraints)
-
-        return ic_batch
-
-    def do_monkey_patch(self):
-        """optuna_integration.botorch には optimize_acqf に constraints を渡す方法が用意されていないので、モンキーパッチして渡す
-
-        モンキーパッチ自体は最適化実行前のどの時点で呼んでも機能するが、additional_kwargs の更新後に
-        モンキーパッチを呼ぶ必要があるのでコンストラクタにこの処理は入れない。
-        各 add_constraint に入れるのはいいかも。
-
-        """
-
-        # === reconstruct argument ``options`` for optimize_acqf ===
-        options = dict()  # initialize
-
-        # for nonlinear-constraint
-        options.update(dict(batch_limit=1))
-
-        # for gen_candidates_scipy()
-        # use COBYLA or SLSQP only.
-        options.update(dict(method='SLSQP'))
-
-        # make partial of optimize_acqf used in optuna_integration.botorch and replace to it.
-        original_fun = optuna_integration.botorch.optimize_acqf
-        overwritten_fun = NonOverwritablePartial(
-            original_fun,
-            q=1,  # for nonlinear constraints
-            options=options,
-            num_restarts=20,  # gen_batch_initial_conditions に渡すべきで、self.generate_initial_conditions に渡される変数。
-            raw_samples=512,  # gen_batch_initial_conditions に渡すべきで、self.generate_initial_conditions に渡される変数。
-            nonlinear_inequality_constraints=self.nonlinear_inequality_constraints,
-            ic_generator=self.generate_initial_conditions,
-        )
-        optuna_integration.botorch.optimize_acqf = overwritten_fun