Update MeshedRegion.grid coordinates when necessary and clear data before plotting.

.ci/build_doc.bat

@@ -3,10 +3,11 @@ call sphinx-apidoc -o ../docs/source/api ../ansys ../ansys/dpf/core/log.py ^
 ../ansys/dpf/core/help.py ../ansys/dpf/core/mapping_types.py ../ansys/dpf/core/ipconfig.py ^
 ../ansys/dpf/core/field_base.py ../ansys/dpf/core/cache.py ../ansys/dpf/core/misc.py ^
 ../ansys/dpf/core/check_version.py ../ansys/dpf/core/operators/build.py ../ansys/dpf/core/operators/specification.py ^
-../ansys/dpf/core/vtk_helper.py ../ansys/dpf/core/label_space.py ^
+../ansys/dpf/core/vtk_helper.py ../ansys/dpf/core/label_space.py ../ansys/dpf/core/examples/python_plugins/* ^
+../ansys/dpf/core/examples/examples.py ^
  -f --implicit-namespaces --separate  --no-headings
 pushd .
 cd ../docs/
 call make clean
-call make html
+call make html -v -v -v -P
 popd

.github/workflows/ci.yml

@@ -243,6 +243,10 @@ jobs:
       - name: "Setup Graphviz"
         uses: ts-graphviz/setup-graphviz@v1
 
+      - name: "Install OS packages"
+        run: |
+          choco install pandoc
+
       - name: "Install documentation packages for Python"
         run: |
           pip install -r requirements/requirements_docs.txt
@@ -297,12 +301,14 @@ jobs:
         with:
           files: docs/build/html
           dest: HTML-doc-${{env.PACKAGE_NAME}}.zip
+        if: always()
 
       - name: "Upload HTML Documentation"
         uses: actions/upload-artifact@v3
         with:
           name: HTML-doc-${{env.PACKAGE_NAME}}
           path: HTML-doc-${{env.PACKAGE_NAME}}.zip
+        if: always()
 
   run_examples:
     name: "Run Examples with/without bin"

ansys/dpf/core/meshed_region.py

@@ -18,6 +18,21 @@
 from ansys.dpf.gate import meshed_region_capi, meshed_region_grpcapi
 
 
+def update_grid(func):
+    # Decorate mesh setters to centralize the update logic of pyvista objects.
+    def wrapper(*args, **kwargs):
+        mesh = args[0]
+        if mesh._full_grid is not None:
+            # Treat each setter separately to improve performance by updating the minimum required.
+            if func.__name__ == 'set_coordinates_field':
+                # When setting node coordinates
+                from ansys.dpf.core.vtk_helper import vtk_update_coordinates
+                vtk_update_coordinates(vtk_grid=mesh._full_grid, coordinates_array=args[1].data)
+
+        return func(*args, **kwargs)
+    return wrapper
+
+
 @class_handling_cache
 class MeshedRegion:
     """
@@ -275,6 +290,7 @@ def set_property_field(self, property_name, value):
         else:
             self._api.meshed_region_set_property_field(self, property_name, value)
 
+    @update_grid
     @version_requires("3.0")
     def set_coordinates_field(self, coordinates_field):
         """

ansys/dpf/core/plotter.py

@@ -391,6 +391,8 @@ def add_mesh(self, meshed_region, deform_by=None, scale_factor=1.0, **kwargs):
         >>> pl.add_mesh(mesh)
 
         """
+        if meshed_region.grid is not None:
+            meshed_region.grid.clear_data()
         self._internal_plotter.add_mesh(meshed_region=meshed_region,
                                         deform_by=deform_by,
                                         scale_factor=scale_factor,
@@ -767,6 +769,7 @@ def plot_contour(
             self._internal_plotter.add_scale_factor_legend(scale_factor, **kwargs)
         else:
             grid = mesh.grid
+        grid.clear_data()
         self._internal_plotter._plotter.add_mesh(grid, scalars=overall_data, **kwargs_in)
 
         background = kwargs.pop("background", None)

ansys/dpf/core/vtk_helper.py

@@ -226,3 +226,8 @@ def compute_offset():
         offset = compute_offset()
 
     return pv.UnstructuredGrid(offset, cells, vtk_cell_type, nodes)
+
+
+def vtk_update_coordinates(vtk_grid, coordinates_array):
+    from copy import copy
+    vtk_grid.points = copy(coordinates_array)

docs/make.bat

@@ -11,6 +11,7 @@ set SOURCEDIR=source
 set BUILDDIR=build
 
 if "%1" == "" goto help
+if "%1" == "clean" goto clean
 
 %SPHINXBUILD% >NUL 2>NUL
 if errorlevel 9009 (
@@ -28,6 +29,30 @@ if errorlevel 9009 (
 %SPHINXBUILD% -M %1 %SOURCEDIR% %BUILDDIR% %SPHINXOPTS% %O%
 goto end
 
+:clean
+echo.Cleaning files form previous build...
+IF EXIST "build" (
+    rmdir "build" /s /q
+)
+IF EXIST "source\images\auto-generated" (
+    rmdir "source\images\auto-generated" /s /q
+)
+IF EXIST "source\examples\07-python-operators\plugins" (
+    robocopy "source\examples\07-python-operators\plugins" "source\_temp\plugins" /E >nul 2>&1
+)
+IF EXIST "source\examples" (
+    rmdir "source\examples" /s /q
+)
+IF EXIST "source\_temp\plugins" (
+    robocopy "source\_temp\plugins" "source\examples\07-python-operators\plugins" /E >nul 2>&1
+)
+IF EXIST "source\_temp" (
+    rmdir "source\_temp" /s /q
+)
+
+echo.Done.
+goto end
+
 :help
 %SPHINXBUILD% -M help %SOURCEDIR% %BUILDDIR% %SPHINXOPTS% %O%
 

docs/source/conf.py

@@ -67,6 +67,7 @@
 
 autosummary_generate = True
 
+autodoc_mock_imports = ["ansys.dpf.core.examples.python_plugins"]
 
 # Add any paths that contain templates here, relative to this directory.
 # templates_path = ['_templates']
@@ -85,7 +86,7 @@
 #
 # This is also used if you do content translation via gettext catalogs.
 # Usually you set "language" from the command line for these cases.
-language = None
+language = "en"
 
 # List of patterns, relative to source directory, that match files and
 # directories to ignore when looking for source files.
@@ -99,6 +100,28 @@
 # -- Sphinx Gallery Options
 from sphinx_gallery.sorting import FileNameSortKey
 
+
+def reset_servers(gallery_conf, fname, when):
+    import psutil
+    from ansys.dpf.core import server
+    import gc
+
+    gc.collect()
+    server.shutdown_all_session_servers()
+
+    proc_name = "Ans.Dpf.Grpc"
+    nb_procs = 0
+    for proc in psutil.process_iter():
+        try:
+            # check whether the process name matches
+            if proc_name in proc.name():
+                # proc.kill()
+                nb_procs += 1
+        except psutil.NoSuchProcess:
+            pass
+    print(f"Counted {nb_procs} {proc_name} processes {when} the example.")
+
+
 sphinx_gallery_conf = {
     # convert rst to md for ipynb
     "pypandoc": True,
@@ -118,18 +141,24 @@
     # 'first_notebook_cell': ("%matplotlib inline\n"
     #                         "from pyvista import set_plot_theme\n"
     #                         "set_plot_theme('document')"),
+    "reset_modules_order": 'both',
+    "reset_modules": (reset_servers,),
 }
 
 autodoc_member_order = "bysource"
 
 
 # -- Options for HTML output -------------------------------------------------
+html_short_title = html_title = "PyDPF-Core"
 html_theme = "ansys_sphinx_theme"
 html_logo = pyansys_logo_black
 html_theme_options = {
-    "github_url": "https://github.com/pyansys/DPF-Core",
+    "github_url": "https://github.com/pyansys/pydpf-core",
     "show_prev_next": False,
-    "logo_link": "https://dpfdocs.pyansys.com/"  # navigate to the main page
+    "show_breadcrumbs": True,
+    "additional_breadcrumbs": [
+        ("PyAnsys", "https://docs.pyansys.com/"),
+    ],
 }
 
 

docs/source/user_guide/operators.rst

@@ -126,7 +126,7 @@ You can create data sources in two ways:
 
 
 Because several other examples use the ``Model`` class, this example uses the
-``DataSources``class:
+``DataSources`` class:
 
 .. code-block:: python
 
@@ -168,7 +168,7 @@ can also be connected to work on a temporal subset:
 
 Evaluate operators
 ~~~~~~~~~~~~~~~~~~
-With all the required inputs assigned, you can output the :class:`ansys.dpf.core.fields_container`_
+With all the required inputs assigned, you can output the :class:`ansys.dpf.core.fields_container`
 class from the operator:
 
 .. code-block:: python

examples/03-advanced/04-extrapolation_stress_3d.py

@@ -69,7 +69,7 @@
 ###############################################################################
 # Nodal stress result of integration points
 ###############################################################################
-# The MAPLD command ``ERESX,NO``is used to copy directly the
+# The MAPLD command ``ERESX,NO`` is used to copy directly the
 # Gaussian (integration) points results to the nodes, instead of the
 # results at nodes or elements (which are interpolation of results at a
 # few gauss points).

examples/03-advanced/05-extrapolation_strain_2d.py

@@ -53,7 +53,7 @@
 ###############################################################################
 # Extrapolate from integration points for elastic strain result
 # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-# This example uses the ``gauss_to_node_fc``operator to compute nodal component
+# This example uses the ``gauss_to_node_fc`` operator to compute nodal component
 # elastic strain results from the elastic strain at the integration points.
 
 # Create elastic strain operator to get strain result of integration points

examples/03-advanced/10-asme_secviii_divtwo.py

@@ -41,7 +41,7 @@
 #
 # - ``alfasl`` = input("Introduce ``alfasl`` parameter from ASME\n")
 # - ``alfasl`` = float(alfasl)
-# -``m2`` = input("Introduce ``m2`` parameter from ASME\n")
+# - ``m2`` = input("Introduce ``m2`` parameter from ASME\n")
 # - ``m2`` = float(m2)
 #
 # For this exercise, ``alfasl`` = 2.2 and ``m2`` = .288, which is the same

examples/08-averaging/00-compute_and_average.py

@@ -146,8 +146,7 @@ def average_then_compute_von_mises(analysis):
 # Plot the results
 # ~~~~~~~~~~~~~~~~
 # Plot both von Mises stress fields side by side to compare them.
-# - The first plot displays the results when the equivalent stresses are calculated
-#   first.
+# - The first plot displays the results when the equivalent stresses are calculated first.
 # - The second plot shows the results when the averaging is done first.
 #
 

examples/08-averaging/README.txt

@@ -1,4 +1,4 @@
-.. _averaging_examples
+.. _averaging_examples:
 
 Averaging examples
 ==================

tests/test_meshregion.py

@@ -121,7 +121,15 @@ def test_set_coordinates_field_meshedregion(simple_bar_model):
     field_coordinates.data = new_data
     mesh.set_coordinates_field(field_coordinates)
     field_coordinates = mesh.nodes.coordinates_field
+    point_0 = list(mesh.grid.points[0])
     assert np.allclose(field_coordinates.data[0], [1.0, 1.0, 1.0])
+    assert np.allclose(point_0, [1.0, 1.0, 1.0])
+    field_coordinates.data = field_coordinates.data * 2.0
+    mesh.set_coordinates_field(field_coordinates)
+    field_coordinates = mesh.nodes.coordinates_field
+    point_1 = list(mesh.grid.points[0])
+    assert np.allclose(field_coordinates.data[0], [2.0, 2.0, 2.0])
+    assert np.allclose(point_1, [2.0, 2.0, 2.0])
 
 
 def test_get_element_types_field_meshedregion(simple_bar_model):