morphopt 3.1.5__tar.gz → 3.2.0__tar.gz

morphopt-3.2.0/PKG-INFO

@@ -0,0 +1,151 @@
+Metadata-Version: 2.4
+Name: morphopt
+Version: 3.2.0
+Summary: A morphing optimization framework using PyTorch and Abaqus.
+Requires-Python: >=3.12
+Description-Content-Type: text/markdown
+Requires-Dist: torch>=2.9.1
+Requires-Dist: torchvision>=0.15.1
+Requires-Dist: numpy>=2.0.0
+Requires-Dist: pypardiso>=0.4.7
+Requires-Dist: pyside6>=6.8.0
+Requires-Dist: pyvista>=0.47.1
+Requires-Dist: pyvistaqt>=0.11.3
+Requires-Dist: vtk>=9.6.1
+Requires-Dist: gmsh>=4.15.0
+Requires-Dist: tabulate>=0.9.0
+Requires-Dist: imageio>=2.37.3
+Requires-Dist: bspmap>=1.0.3
+Requires-Dist: torchfea>=1.1.0
+Requires-Dist: trimesh>=4.11.5
+Requires-Dist: networkx>=3.6.1
+Requires-Dist: scipy>=1.17.1
+Requires-Dist: cpgeo>=1.0.12
+
+# MorphOpt — A Differentiable Optimization Framework for Soft Structure & Morphology Design
+
+**Developed by Zenan Song**
+
+A differentiable structural optimization framework that eliminates manual sensitivity derivation by fusing the adjoint method with automatic differentiation (AD). The sensitivity corresponds to the virtual work of residual force derivatives on the adjoint displacement field, computed via a single backpropagation through the residual graph. MorphOpt provides a unified interface for geometry, load, and material definition, calls torchfea for GPU-accelerated differentiable nonlinear FEA and adjoint-AD sensitivity, and employs a trust-region optimizer with L-BFGS for design updates. Validated on SIMP topology optimization with B-spline density fields and pneumatic soft robot shape optimization with deformation-dependent follower loads.
+
+## Features
+
+- Custom model definition covering geometry, material, and load parameters
+- Differentiable FEA via torchfea with automatic solving and sensitivity propagation
+- L-BFGS and trust-region methods for parameter updates within a subproblem framework
+- Fully extensible via inheritance — users customize every component in their definition scripts
+
+## Quick Start
+
+Install the package:
+
+```bash
+pip install morphopt
+```
+
+Requirements: Python ≥ 3.12; a PyTorch **float64** build is recommended.
+
+> **Optimizations are run as scripts.** Each job is an ordinary Python file that calls `morphopt.start_optimization(...)`. It runs headlessly in a terminal (local or remote), and a finished run can be stopped and resumed later from its result folder. Script execution is the recommended way to run jobs.
+
+1. **Run a minimal example** — execute it directly in a terminal:
+
+   ```bash
+   # Shape optimization (pneumatic bending actuator)
+   python examples/bendingactuator.py
+
+   # SIMP topology optimization (compliant gripper)
+   python examples/gripper.py
+   ```
+
+2. **Run your own job** — copy an example as a starting point, adapt the definition (geometry, material, load, updater, objective), and launch it as a script:
+
+   ```bash
+   python my_task.py
+   ```
+
+   Each run writes its results (logs, FEA files, restart script) into a timestamped result folder under the run root.
+
+3. **Optional — the GUI** is a visual aid for authoring a problem definition and inspecting finished runs in one window. It can also export a runnable script (`.py`), which you can either run directly in the UI or further customize before executing it via Step 2 in a terminal.
+
+   ```bash
+   morphopt-ui            # console command, available after `pip install morphopt`
+   # python -m morphopt.ui works too (handy when running from a source code directory)
+   ```
+
+The complete UI workflow (definition part + observer part) is described in [`docs/UI_usage.md`](docs/UI_usage.md).
+
+Further references:
+
+- Task definition guide: `docs/module_definition_guide.md`
+- Module API reference: `docs/module_reference.md`
+- Theory & papers: `docs/theory/`
+
+## Directory Overview
+
+- `src/morphopt/`
+  - `optcore` — Core package (controller, params, solver, objective, updaters)
+  - `shapeopt/` — B-spline shape optimization
+  - `simp/` — SIMP topology optimization
+  - `codesign/` — Shell co-design (shape + material)
+  - `ui/` — PySide6 graphical UI (definition + observer in one window)
+- `examples/` — Example tasks (bending actuator, gripper, etc.)
+- `docs/` — Usage guides, API reference, theory papers
+- `tests/` — Gradient checks, geometry tests, UI tests
+- `ui_runs/` — runs launched from the UI and their results (`<label>_T<timestamp>/`)
+
+## License
+
+Research use preferred. For production or commercial use, please evaluate and ensure numerical and engineering robustness first.
+
+## Citation
+
+If you use MorphOpt in your research, please cite:
+
+- TRO 2023
+
+```
+@ARTICLE{10296178,
+  author={Chen, Feifei and Song, Zenan and Chen, Shitong and Gu, Guoying and Zhu, Xiangyang},
+  journal={IEEE Transactions on Robotics}, 
+  title={Morphological Design for Pneumatic Soft Actuators and Robots With Desired Deformation Behavior}, 
+  year={2023},
+  volume={39},
+  number={6},
+  pages={4408-4428},
+  keywords={Shape;Surface morphology;Robots;Splines (mathematics);Pneumatic systems;Actuators;Soft robotics;Morphological design;shape optimization;pneumatic soft robots},
+  doi={10.1109/TRO.2023.3323825}}
+```
+
+- TRO 2026
+
+```
+@ARTICLE{11617347,
+  author={Song, Zenan and Chen, Feifei and Gu, Guoying and Zhu, Xiangyang},
+  journal={IEEE Transactions on Robotics}, 
+  title={Continuum Jacobian-Based Computational Morphogenesis for Soft Robotic Workspace Optimization}, 
+  year={2026},
+  volume={42},
+  number={},
+  pages={3048-3068},
+  keywords={Optimization;Design methodology;Surfaces;Modeling;End effectors;Shape;Soft robotics;Deformation;Vectors;Robots;Computational morphogenesis;Jacobian shape derivative;soft robotics;workspace optimization},
+  doi={10.1109/TRO.2026.3716051}}
+
+```
+
+- TMECH 2026
+
+```
+@ARTICLE{11536090,
+  author={Song, Zenan and Zhou, Guozhen and Chen, Feifei},
+  journal={IEEE/ASME Transactions on Mechatronics}, 
+  title={Contact-Aware Soft Robot Design via Differentiable Simulation and Morphological Optimization}, 
+  year={2026},
+  volume={31},
+  number={4},
+  pages={4576-4584},
+  keywords={Contacts;Optimization;Modeling;Design methodology;Soft robotics;Surfaces;Actuators;Materials;Force;Simulation;Differentiable contact modeling;morphological design;soft robotics},
+  doi={10.1109/TMECH.2026.3692920}}
+
+```
+
+

morphopt-3.2.0/README.md

@@ -0,0 +1,127 @@
+# MorphOpt — A Differentiable Optimization Framework for Soft Structure & Morphology Design
+
+**Developed by Zenan Song**
+
+A differentiable structural optimization framework that eliminates manual sensitivity derivation by fusing the adjoint method with automatic differentiation (AD). The sensitivity corresponds to the virtual work of residual force derivatives on the adjoint displacement field, computed via a single backpropagation through the residual graph. MorphOpt provides a unified interface for geometry, load, and material definition, calls torchfea for GPU-accelerated differentiable nonlinear FEA and adjoint-AD sensitivity, and employs a trust-region optimizer with L-BFGS for design updates. Validated on SIMP topology optimization with B-spline density fields and pneumatic soft robot shape optimization with deformation-dependent follower loads.
+
+## Features
+
+- Custom model definition covering geometry, material, and load parameters
+- Differentiable FEA via torchfea with automatic solving and sensitivity propagation
+- L-BFGS and trust-region methods for parameter updates within a subproblem framework
+- Fully extensible via inheritance — users customize every component in their definition scripts
+
+## Quick Start
+
+Install the package:
+
+```bash
+pip install morphopt
+```
+
+Requirements: Python ≥ 3.12; a PyTorch **float64** build is recommended.
+
+> **Optimizations are run as scripts.** Each job is an ordinary Python file that calls `morphopt.start_optimization(...)`. It runs headlessly in a terminal (local or remote), and a finished run can be stopped and resumed later from its result folder. Script execution is the recommended way to run jobs.
+
+1. **Run a minimal example** — execute it directly in a terminal:
+
+   ```bash
+   # Shape optimization (pneumatic bending actuator)
+   python examples/bendingactuator.py
+
+   # SIMP topology optimization (compliant gripper)
+   python examples/gripper.py
+   ```
+
+2. **Run your own job** — copy an example as a starting point, adapt the definition (geometry, material, load, updater, objective), and launch it as a script:
+
+   ```bash
+   python my_task.py
+   ```
+
+   Each run writes its results (logs, FEA files, restart script) into a timestamped result folder under the run root.
+
+3. **Optional — the GUI** is a visual aid for authoring a problem definition and inspecting finished runs in one window. It can also export a runnable script (`.py`), which you can either run directly in the UI or further customize before executing it via Step 2 in a terminal.
+
+   ```bash
+   morphopt-ui            # console command, available after `pip install morphopt`
+   # python -m morphopt.ui works too (handy when running from a source code directory)
+   ```
+
+The complete UI workflow (definition part + observer part) is described in [`docs/UI_usage.md`](docs/UI_usage.md).
+
+Further references:
+
+- Task definition guide: `docs/module_definition_guide.md`
+- Module API reference: `docs/module_reference.md`
+- Theory & papers: `docs/theory/`
+
+## Directory Overview
+
+- `src/morphopt/`
+  - `optcore` — Core package (controller, params, solver, objective, updaters)
+  - `shapeopt/` — B-spline shape optimization
+  - `simp/` — SIMP topology optimization
+  - `codesign/` — Shell co-design (shape + material)
+  - `ui/` — PySide6 graphical UI (definition + observer in one window)
+- `examples/` — Example tasks (bending actuator, gripper, etc.)
+- `docs/` — Usage guides, API reference, theory papers
+- `tests/` — Gradient checks, geometry tests, UI tests
+- `ui_runs/` — runs launched from the UI and their results (`<label>_T<timestamp>/`)
+
+## License
+
+Research use preferred. For production or commercial use, please evaluate and ensure numerical and engineering robustness first.
+
+## Citation
+
+If you use MorphOpt in your research, please cite:
+
+- TRO 2023
+
+```
+@ARTICLE{10296178,
+  author={Chen, Feifei and Song, Zenan and Chen, Shitong and Gu, Guoying and Zhu, Xiangyang},
+  journal={IEEE Transactions on Robotics}, 
+  title={Morphological Design for Pneumatic Soft Actuators and Robots With Desired Deformation Behavior}, 
+  year={2023},
+  volume={39},
+  number={6},
+  pages={4408-4428},
+  keywords={Shape;Surface morphology;Robots;Splines (mathematics);Pneumatic systems;Actuators;Soft robotics;Morphological design;shape optimization;pneumatic soft robots},
+  doi={10.1109/TRO.2023.3323825}}
+```
+
+- TRO 2026
+
+```
+@ARTICLE{11617347,
+  author={Song, Zenan and Chen, Feifei and Gu, Guoying and Zhu, Xiangyang},
+  journal={IEEE Transactions on Robotics}, 
+  title={Continuum Jacobian-Based Computational Morphogenesis for Soft Robotic Workspace Optimization}, 
+  year={2026},
+  volume={42},
+  number={},
+  pages={3048-3068},
+  keywords={Optimization;Design methodology;Surfaces;Modeling;End effectors;Shape;Soft robotics;Deformation;Vectors;Robots;Computational morphogenesis;Jacobian shape derivative;soft robotics;workspace optimization},
+  doi={10.1109/TRO.2026.3716051}}
+
+```
+
+- TMECH 2026
+
+```
+@ARTICLE{11536090,
+  author={Song, Zenan and Zhou, Guozhen and Chen, Feifei},
+  journal={IEEE/ASME Transactions on Mechatronics}, 
+  title={Contact-Aware Soft Robot Design via Differentiable Simulation and Morphological Optimization}, 
+  year={2026},
+  volume={31},
+  number={4},
+  pages={4576-4584},
+  keywords={Contacts;Optimization;Modeling;Design methodology;Soft robotics;Surfaces;Actuators;Materials;Force;Simulation;Differentiable contact modeling;morphological design;soft robotics},
+  doi={10.1109/TMECH.2026.3692920}}
+
+```
+
+

{morphopt-3.1.5 → morphopt-3.2.0}/pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name = "morphopt"
-version = "3.1.5"
+version = "3.2.0"
 description = "A morphing optimization framework using PyTorch and Abaqus."
 readme = "README.md"
 requires-python = ">=3.12"
@@ -9,7 +9,7 @@ dependencies = [
     "torchvision>=0.15.1",
     "numpy>=2.0.0",
     "pypardiso>=0.4.7",
-    "pyqt6>=6.11.0",
+    "pyside6>=6.8.0",
     "pyvista>=0.47.1",
     "pyvistaqt>=0.11.3",
     "vtk>=9.6.1",
@@ -17,9 +17,13 @@ dependencies = [
     "tabulate>=0.9.0",
     "imageio>=2.37.3",
     "bspmap>=1.0.3",
-    "torchfea>=1.0.17",
+    "torchfea>=1.1.0",
     "trimesh>=4.11.5",
     "networkx>=3.6.1",
     "scipy>=1.17.1",
     "cpgeo>=1.0.12",
 ]
+
+[project.scripts]
+# console entry for the graphical UI (after `pip install morphopt`)
+morphopt-ui = "morphopt.ui:run_app"

{morphopt-3.1.5 → morphopt-3.2.0}/src/morphopt/__init__.py

@@ -66,7 +66,7 @@ from .optcore.history import History
 from .optcore.baseobject import BaseObject
 # endregion
 
-from .opt_runner import start_optimization, debug_optimization, view_optimization_result
+from .opt_runner import start_optimization, debug_optimization
 
 # region modules
 from . import shapeopt

{morphopt-3.1.5 → morphopt-3.2.0}/src/morphopt/codesign/feaparams.py

@@ -40,7 +40,7 @@ class CodesignFEAParams(FEAParams):
             )
             surfcenter = surfnodes.mean(dim=1)
 
-            fe.assembly.get_instance('final_model').external_surface = 'surface_1_offset'
+            fe.assembly.get_instance('final_model').exterior_surface = 'surface_1_offset'
 
             mesh_1all = fe.assembly.get_instance('final_model').get_mesh(surf_name='surface_1_All')
             mesh_1offset = fe.assembly.get_instance('final_model').get_mesh(surf_name='surface_1_offset')
@@ -72,4 +72,4 @@ class CodesignFEAParams(FEAParams):
 
             raise ValueError("The Gaussian weights of the C3D6 elements are negative. Please check the mesh quality.")
         
-        return fe
+        return fe

morphopt-3.2.0/src/morphopt/opt_runner.py

@@ -0,0 +1,82 @@
+"""Run / debug optimization workers.
+
+Runs started from *code* never open a UI: the PySide6 observer lives inside the
+MorphOpt UI and simply polls the result folder on disk (per finished
+iteration).  The old ``no_gui`` flag and separate observer process were
+therefore removed; browse results from the observer page in the main UI.
+"""
+
+
+def start_optimization(device='cpu', restart_per_iteration: int = 20,
+                       path_result: str = None, target_iteration=None, **kwargs):
+    """
+    Start an optimization process (always headless).
+
+    Args:
+        device (str, optional): device to use for computation ('cpu'/'cuda:0').
+        restart_per_iteration (int): iterations between process restarts.
+        path_result (str, optional): existing result folder -> continue there,
+            otherwise a fresh timestamped run folder is created.
+        target_iteration (int, optional): when continuing, resume from this
+            step; None/0 -> resume from the last saved step.
+    """
+    import multiprocessing as mp
+
+    mp.set_start_method("spawn", force=True)
+
+    if path_result is None:
+        import __main__
+        main_filepath = __main__.__file__
+    else:
+        main_filepath = path_result + '/scripts/' + 'MAIN_SCRIPT_FOR_RESTART.py'
+
+    from .taskoptmization import TaskOptimization
+
+    process_optimization = mp.Process(
+        target=TaskOptimization.task_optimization,
+        kwargs={'path_result': path_result,
+                'main_filepath': main_filepath,
+                'device': device,
+                'target_iteration': target_iteration,
+                'restart_per_iteration': restart_per_iteration})
+    process_optimization.start()
+    process_optimization.join()
+
+
+def debug_optimization(device='cpu', restart_per_iteration: int = 20,
+                       path_result: str = None, target_iteration=None, **kwargs):
+    """
+    Run the optimization in the current process (debugger friendly).
+    """
+    import os
+    import sys
+    import __main__
+    import multiprocessing as mp
+
+    from .taskoptmization import TaskOptimization
+    import morphopt
+
+    mp.set_start_method("spawn", force=True)
+
+    if path_result is None:
+        main_filepath = __main__.__file__
+    else:
+        main_filepath = path_result + '/scripts/' + 'MAIN_SCRIPT_FOR_RESTART.py'
+
+    filename = os.path.splitext(os.path.basename(main_filepath))[0]
+    filepath = os.path.dirname(main_filepath)
+    os.chdir(filepath)
+    sys.path.append(os.getcwd())
+
+    Controller = getattr(__import__(filename), 'ThisController')
+
+    controller: morphopt.Controller = Controller()
+    controller.restart_per_iteration = restart_per_iteration
+    controller.optdevice = device
+    controller._debug_mode = True
+
+    TaskOptimization.optmain(device=device,
+                             path_result=path_result,
+                             target_iteration=target_iteration,
+                             restart_per_iteration=restart_per_iteration,
+                             main_filepath=main_filepath)

{morphopt-3.1.5 → morphopt-3.2.0}/src/morphopt/optcore/modelparams/geometry.py

@@ -85,6 +85,8 @@ class FixedGeometryINP(FixedGeometry):
         
         assembly = torchfea.Assembly()
 
-        assembly.add_part(part=fe_ext.assembly.get_part(self._part_name), name='final_model')
-        assembly.add_instance(instance=torchfea.Instance(part_name='final_model', external_surface='surface_0_All'), name='final_model')
-        return assembly
+        assembly.add_part(part=fe_ext.assembly.get_part(self._part_name), name='final_model')
+        instance = torchfea.Instance(part_name='final_model')
+        assembly.add_instance(instance=instance, name='final_model')
+        instance.exterior_surface = 'surface_0_All'
+        return assembly

{morphopt-3.1.5 → morphopt-3.2.0}/src/morphopt/shapeopt/geometryinterfaces/basesurfaceinterface.py

@@ -14,7 +14,7 @@ class BaseInterface():
     Class to handle the surface of the morphable model.
     """
     
-    def __init__(self, symmetric: list[int] = None, *args, **kwargs) -> None:
+    def __init__(self, *args, **kwargs) -> None:
         """
         Initialize the Surface class.
 
@@ -23,22 +23,6 @@ class BaseInterface():
             surf_type (int) : The type of the surface.
                 - 0: bspline surface
                 - 1: closed surface
-            symmetric (list[int]) : The symmetry of the surface.
-                - 0: no symmetry
-                - 1: axis symmetry
-                    0: x-axis symmetry
-                    1: y-axis symmetry
-                    2: z-axis symmetry
-        """
-
-        self.symmetric: list[int] = symmetric
-        """
-        the symmetry of the surface.
-        # 0: no symmetry
-        # 1: axis symmetry
-            ## 0: x-axis symmetry
-            ## 1: y-axis symmetry
-            ## 2: z-axis symmetry
         """
 
         self.surf_node_idx: np.ndarray
@@ -737,12 +721,6 @@ class CpBasedInterface(BaseInterface):
             surf_type (int) : The type of the surface.
                 - 0: bspline surface
                 - 1: closed surface
-            symmetric (list[int]) : The symmetry of the surface.
-                - 0: no symmetry
-                - 1: axis symmetry
-                    0: x-axis symmetry
-                    1: y-axis symmetry
-                    2: z-axis symmetry
         """
 
         super().__init__(*args, **kwargs)

{morphopt-3.1.5 → morphopt-3.2.0}/src/morphopt/shapeopt/geometryinterfaces/bspsurfaceinterface.py

@@ -36,7 +36,7 @@ class BspInterface(CpBasedInterface):
             
         def build(self):
             """
-            构建版本 - 改用 addBSplineSurface + Symmetric Wrapping + Trim
+            构建版本 - 改用 addBSplineSurface + Wrapping + Trim
             通过上下游双重延伸(Padding)确保接缝处C2连续性
             """
             self.volume_tag = None
@@ -392,8 +392,8 @@ class BspInterface(CpBasedInterface):
             gmsh.fltk.run()
 
 
-    def __init__(self, surface: bspmap.BSP, init_size: float, symmetric = [0], MaxR = 0.2, MaxFF = 0.1, MaxC = 1.0):
-        super().__init__(surface, symmetric)
+    def __init__(self, surface: bspmap.BSP, init_size: float, MaxR = 0.2, MaxFF = 0.1, MaxC = 1.0):
+        super().__init__(surface)
 
         self.model = surface
         """The B-spline surface model."""
@@ -691,7 +691,7 @@ class BspInterface(CpBasedInterface):
 
 
     @classmethod
-    def initialize_cylinder(cls, r0: float, length: float, seed_size: float, flip: bool, num_U_ratio: int = 1, num_V_ratio: int = 1, symmetric: list[int] = [0], degree = 3, init_location = [0.,0.,0.], maxR = 0.2, maxC = 1., maxFF = 0.2, perturbation_L = -1.):    
+    def initialize_cylinder(cls, r0: float, length: float, seed_size: float, flip: bool, num_U_ratio: int = 1, num_V_ratio: int = 1, degree = 3, init_location = [0.,0.,0.], maxR = 0.2, maxC = 1., maxFF = 0.2, perturbation_L = -1.):    
         """
         Initialize the B-spline surface for the optimization process.
 
@@ -701,7 +701,6 @@ class BspInterface(CpBasedInterface):
             seed_size (float): The size of the seed for the B-spline surface.
             num_U_ratio (int, optional): The ratio for the number of points in the U direction. Default is 1.
             num_V_ratio (int, optional): The ratio for the number of points in the V direction. Default is 1.
-            symmetric (list[int]): The symmetry of the surface.
             flip (bool): Whether to flip the surface or not.
             degree (int, optional): The degree of the B-spline surface. Default is 3.
             init_location (list[float], optional): The initial location of the surface. Default is [0., 0., 0.].
@@ -713,7 +712,6 @@ class BspInterface(CpBasedInterface):
         Returns:
             BSP (BSP_Surf): The initialized B-spline surface object.
             surf_type (int): The type of the surface (0 for B-spline surface).
-            symmetric (list[int]): The symmetry of the surface.
         """
 
         numU = round(r0 * 2 * np.pi / seed_size)
@@ -732,9 +730,6 @@ class BspInterface(CpBasedInterface):
         [y, x] = torch.meshgrid(y, x, indexing='ij')
         theta = 2 * torch.pi * x + (1 / numU) * torch.pi
 
-        if not flip:
-            theta = -theta
-
         P0[:, :, 0] = torch.cos(theta) * r0
         P0[:, :, 1] = torch.sin(theta) * r0
         P0[:, :, 2] = length * y
@@ -760,7 +755,7 @@ class BspInterface(CpBasedInterface):
                          size=[numV, numU],
                          control_points=P0.cpu().numpy().reshape([-1, 3]))
 
-        output = cls(bsp, init_size=seed_size, symmetric=symmetric, MaxR=maxR, MaxC=maxC, MaxFF=maxFF)
+        output = cls(bsp, init_size=seed_size, MaxR=maxR, MaxC=maxC, MaxFF=maxFF)
         output.flip = flip
 
         return output

{morphopt-3.1.5 → morphopt-3.2.0}/src/morphopt/shapeopt/geometryinterfaces/cpgeosurfaceinterface.py

@@ -20,8 +20,8 @@ class CPGEOInterface(CpBasedInterface):
     providing a flexible alternative to B-spline surfaces for complex geometries.
     """
 
-    def __init__(self, surface: cpgeo.CPGEO, init_size: float, symmetric = [0], MaxC = 1.0):
-        super().__init__(surface, symmetric)
+    def __init__(self, surface: cpgeo.CPGEO, init_size: float, MaxC = 1.0):
+        super().__init__(surface)
 
         self.model = surface
         """The CPGEO surface model."""
@@ -363,7 +363,7 @@ class CPGEOInterface(CpBasedInterface):
 
     @classmethod
     def initialize_cylinder(cls, r0: float, length: float, seed_size: float, flip: bool, 
-                          num_U_ratio: int = 1, num_V_ratio: int = 1, symmetric: list[int] = [0], 
+                          num_U_ratio: int = 1, num_V_ratio: int = 1, 
                           degree = 3, init_location = [0.,0.,0.], maxR = 0.2, maxC = 1., maxFF = 0.2, 
                           perturbation_L = -1.):
         """
@@ -376,7 +376,6 @@ class CPGEOInterface(CpBasedInterface):
             flip (bool): Whether to flip the surface normal or not.
             num_U_ratio (int, optional): The ratio for mesh density in circumferential direction. Default is 1.
             num_V_ratio (int, optional): The ratio for mesh density in axial direction. Default is 1.
-            symmetric (list[int]): The symmetry of the surface.
             degree (int, optional): Unused for CPGEO. Default is 3.
             init_location (list[float], optional): The initial location of the surface. Default is [0., 0., 0.].
             maxR (float, optional): The maximum curvature derivative constraint. Default is 0.2.
@@ -451,14 +450,14 @@ class CPGEOInterface(CpBasedInterface):
         )
         
         # Create interface
-        output = cls(cpgeo_model, init_size=seed_size, symmetric=symmetric, 
+        output = cls(cpgeo_model, init_size=seed_size, 
                     MaxR=maxR, MaxC=maxC, MaxFF=maxFF)
         output.flip = flip
 
         return output
     
     @classmethod
-    def initialize_Sphere(cls, r0: float, seed_size: float, flip: bool, symmetric: list[int] = [0], init_location = [0.,0.,0.], 
+    def initialize_Sphere(cls, r0: float, seed_size: float, flip: bool, init_location = [0.,0.,0.], 
                           MaxC = 1.0):
         """
         Initialize the CPGEO surface for sphere shape optimization.
@@ -467,7 +466,6 @@ class CPGEOInterface(CpBasedInterface):
             r0 (float): The radius of the sphere.
             seed_size (float): The size of the mesh seed.
             flip (bool): Whether to flip the surface normal or not.
-            symmetric (list[int]): The symmetry of the surface.
             init_location (list[float], optional): The initial location of the surface. Default is [0., 0., 0.].
             MaxC (float, optional): The maximum curvature constraint. Default is 1.0
         Returns:
@@ -505,7 +503,7 @@ class CPGEOInterface(CpBasedInterface):
         )
         
         # Create interface
-        output = cls(cpgeo_model, init_size=seed_size, symmetric=symmetric, MaxC=MaxC)
+        output = cls(cpgeo_model, init_size=seed_size, MaxC=MaxC)
         output.flip = flip
 
         return output

{morphopt-3.1.5 → morphopt-3.2.0}/src/morphopt/shapeopt/geometryparams.py

@@ -581,7 +581,7 @@ class GeometryParams(BaseGeometry):
 
     def apply_surface_constraints(self) -> None:
         """
-        Apply the constraints (e.g. the symmetric constraint) of the surfaces.
+        Apply the constraints (e.g. the surface constraint) of the surfaces.
         """
         pass
 
@@ -640,7 +640,9 @@ class GeometryParams(BaseGeometry):
         assembly = torchfea.Assembly()
 
         assembly.add_part(part=part, name='final_model')
-        assembly.add_instance(instance=torchfea.Instance(part_name='final_model', external_surface='surface_0_All'), name='final_model')
+        instance = torchfea.Instance(part_name='final_model')
+        assembly.add_instance(instance=instance, name='final_model')
+        instance.exterior_surface = 'surface_0_All'
         return assembly
 
     def _regenerate(self, path_result: str, pools = None):
@@ -776,4 +778,3 @@ class GeometryParams(BaseGeometry):
             nodes_new[part.mid_pt_idxmap_torch[:, 2]] = (nodes_new[part.mid_pt_idxmap_torch[:, 0]] + nodes_new[part.mid_pt_idxmap_torch[:, 1]]) / 2
 
         part.nodes = nodes_new
-