Improve Examples and Testing (ansys#301)

* use rcParams default

* use rcParams for default plotting

* update for additional kwargs and dependencies

* fixed load_table

* using shallow copy for plotting

Author: akaszynski

docs/conf.py

@@ -99,6 +99,7 @@
 from sphinx_gallery.sorting import FileNameSortKey
 
 sphinx_gallery_conf = {
+    'pypandoc': True,  # convert rst to md for ipynb
     # path to your examples scripts
     "examples_dirs": [
         "../examples/",

examples/02-mapdl-examples/2d_plate_with_a_hole.py

@@ -17,7 +17,6 @@
 
 import pyansys
 
-# os.environ['I_MPI_SHM_LMT'] = 'shm'  # necessary on Ubuntu without "smp"
 mapdl = pyansys.launch_mapdl(override=True, additional_switches='-smp',
                              loglevel='ERROR')
 
@@ -126,7 +125,7 @@
 mapdl.nsel('S', 'LOC', 'X', length)
 
 # Verify that only the nodes at length have been selected:
-assert np.unique(mapdl.mesh.nodes[:, 0]) == length
+assert np.allclose(mapdl.mesh.nodes[:, 0], length)
 
 # Next, couple the DOF for these nodes.  This lets us provide a force
 # to one node that will be spread throughout all nodes in this coupled
@@ -148,7 +147,8 @@
 # Solve the static analysis
 mapdl.run('/SOLU')
 mapdl.antype('STATIC')
-mapdl.solve()
+output = mapdl.solve()
+print(output)
 
 ###############################################################################
 # Post-Processing
@@ -370,3 +370,10 @@ def compute_stress_con(ratio):
 plt.xlabel('Ratio of Hole Diameter to Width of Plate')
 plt.ylabel('Stress Concentration')
 plt.show()
+
+
+###############################################################################
+# Cleanup
+# ~~~~~~~
+# Close mapdl when complete
+mapdl.exit()

examples/02-mapdl-examples/3d_notch.py

@@ -16,7 +16,6 @@
 import numpy as np
 import pyansys
 
-# os.environ['I_MPI_SHM_LMT'] = 'shm'  # necessary on Ubuntu without "smp"
 mapdl = pyansys.launch_mapdl(override=True, additional_switches='-smp',
                              loglevel='ERROR')
 
@@ -305,3 +304,10 @@
 # plt.plot(ratios, k_t_exp, label=r'$K_t$ ANSYS')
 # plt.legend()
 # plt.show()
+
+
+###############################################################################
+# Cleanup
+# ~~~~~~~
+# Close mapdl when complete
+mapdl.exit()

examples/02-mapdl-examples/3d_plate_thermal.py

@@ -15,7 +15,6 @@
 import os
 import pyansys
 
-os.environ['I_MPI_SHM_LMT'] = 'shm'  # necessary on Ubuntu without "smp"
 mapdl = pyansys.launch_mapdl(loglevel='ERROR')
 
 ###############################################################################
@@ -67,6 +66,15 @@
 # Alternatively you could also use the result object that reads in the
 # result file using pyansys
 
-nnum, temp = mapdl.result.nodal_temperature(0)
+result = mapdl.result
+nnum, temp = result.nodal_temperature(0)
 # this is the same as pyansys.read_binary(mapdl._result_file)
+print(nnum, temp)
+
+
+###############################################################################
+# Cleanup
+# ~~~~~~~
+# Close mapdl when complete
+mapdl.exit()
 

examples/02-mapdl-examples/mapdl_3d_beam.py

@@ -10,12 +10,11 @@
 
 """
 
-import os
 from pyansys import examples
 import pyansys
 
-os.environ['I_MPI_SHM_LMT'] = 'shm'  # necessary on ubuntu
-mapdl = pyansys.launch_mapdl(override=True, additional_switches='-smp')
+
+mapdl = pyansys.launch_mapdl(loglevel='ERROR')
 
 mapdl.cdread('db', examples.hexarchivefile)
 mapdl.esel('s', 'ELEM', vmin=5, vmax=20)
@@ -40,9 +39,37 @@
 
 mapdl.mxpand(elcalc='YES')
 mapdl.modal_analysis(nmode=6)
-mapdl.exit()
 
-# view the results using pyansys's result viewer
+
+###############################################################################
+# View the results using the pyansys result object
 result = mapdl.result
-result.animate_nodal_solution(0, show_edges=True, loop=False, displacement_factor=10,
-                              movie_filename='demo.gif')
+print(result)
+
+
+###############################################################################
+# Access nodal displacement values
+nnum, disp = result.nodal_displacement(0)
+
+# print the nodes 50 - 59
+for i in range(49, 59):
+    print(nnum[i], disp[i])
+
+
+###############################################################################
+# Plot a modal result
+result.plot_nodal_displacement(0, show_edges=True)
+
+
+###############################################################################
+# Animate a modal result
+# result.animate_nodal_solution(0, show_edges=True, loop=False, displacement_factor=10,
+                              # movie_filename='demo.gif')
+
+
+###############################################################################
+# Cleanup
+# ~~~~~~~
+# Close mapdl when complete
+mapdl.exit()
+

examples/02-mapdl-examples/mapdl_beam.py

@@ -11,8 +11,7 @@
 ###############################################################################
 # Launch MAPDL with interactive plotting
 import pyansys
-mapdl = pyansys.launch_mapdl(override=True, loglevel='WARNING',
-                             interactive_plotting=True)
+mapdl = pyansys.launch_mapdl(interactive_plotting=True, loglevel='ERROR')
 
 # reduce pixel count for larger font documentation plots
 mapdl.gfile(1200)

examples/02-mapdl-examples/path_operations.py

@@ -1,5 +1,5 @@
 """
-.. _ref_plane_stress_concentration:
+.. _ref_path_operation:
 
 Path Operations within pyansys and MAPDL
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
@@ -19,7 +19,6 @@
 
 import pyansys
 
-os.environ['I_MPI_SHM_LMT'] = 'shm'  # necessary on Ubuntu
 mapdl = pyansys.launch_mapdl(loglevel='ERROR')
 
 ###############################################################################

examples/02-mapdl-examples/pyvista_mesh.py

@@ -14,8 +14,7 @@
 import pyansys
 
 # launch MAPDL and run a modal analysis
-os.environ['I_MPI_SHM_LMT'] = 'shm'  # necessary for Ubuntu
-mapdl = pyansys.launch_mapdl(loglevel='WARNING', override=True)
+mapdl = pyansys.launch_mapdl(loglevel='WARNING')
 
 # Create a simple plane mesh centered at (0, 0, 0) on the XY plane
 mesh = pv.Plane(i_resolution=100, j_resolution=100)

examples/02-mapdl-examples/transient_thermal.py

@@ -16,14 +16,12 @@
 """
 # sphinx_gallery_thumbnail_number = 4
 
-import os
 import numpy as np
 
 import matplotlib.pyplot as plt
 import pyansys
 
-os.environ['I_MPI_SHM_LMT'] = 'shm'  # necessary on Ubuntu
-mapdl = pyansys.launch_mapdl()
+mapdl = pyansys.launch_mapdl(loglevel='ERROR')
 
 mapdl.clear()
 mapdl.prep7()
@@ -234,3 +232,4 @@
 plt.xlabel('Time (seconds)')
 plt.ylabel('Temperature ($^\circ$F)')
 plt.show()
+

pyansys/__init__.py

@@ -14,6 +14,13 @@
 from pyansys.misc import Report, _configure_pyvista, _check_has_ansys
 from pyansys.examples.downloads import *
 
+# optional www.pyansys.com feature
+if 'ANSJUPHUB_VER' in os.environ:
+    try:
+        from ansys.mapdl import launch_mapdl
+    except ImportError:
+        raise RuntimeError("you're not supposed to be here")
+
 
 _HAS_ANSYS = _check_has_ansys()