correct for usage with map_flag implemented in 2021R1 (#18)

ansys/mapdl/reader/_version.py

@@ -1,7 +1,14 @@
-"""Version of ansys-mapdl-reader module."""
+"""Version of ansys-mapdl-reader module.
+
+On the ``master`` branch, use 'dev0' to denote a development version.
+For example:
+
+version_info = 0, 58, 'dev0'
+
+"""
 
 # major, minor, patch
-version_info = 0, 50, 6
+version_info = 0, 50, 'dev0'
 
 # Nice string for the version
 __version__ = '.'.join(map(str, version_info))

ansys/mapdl/reader/rst.py

@@ -411,10 +411,16 @@ def materials(self):
     def _load_section_data(self):
         """Loads the section data from the result file"""
         ptr_sec = self._geometry_header['ptrSEC']
-        sec_table = self.read_record(ptr_sec)
+        sec_table, sz = self.read_record(ptr_sec, True)
+
+        # Assemble section pointers relative to ptrSEC
+        if self._map_flag:  # required as of 2021R1
+            sec_pointers = self.read_record(ptr_sec + sz)
+        else:
+            sec_pointers = sec_table
 
         self._section_data = {}
-        for offset in sec_table:
+        for offset in sec_pointers:
             if offset:
                 table = self.read_record(ptr_sec + offset)
                 self._section_data[int(table[0])] = table[1:]

tests/test_rst.py

@@ -78,6 +78,7 @@
 temperature_rst = os.path.join(testfiles_path, 'temp_v13.rst')
 temperature_known_result = os.path.join(testfiles_path, 'temp_v13.npz')
 
+
 @pytest.fixture(scope='module')
 def hex_pipe_corner():
     filename = os.path.join(testfiles_path, 'rst', 'cyc_stress.rst')
@@ -96,6 +97,18 @@ def volume_rst():
     return pymapdl_reader.read_binary(rst_file)
 
 
+def test_map_flag_section_data():
+    # basic 4 element SHELL181 result files
+    shell181_2020r2 = os.path.join(testfiles_path, 'shell181_2020R2.rst')
+    shell181_2021r1 = os.path.join(testfiles_path, 'shell181_2021R1.rst')
+
+    rst_2020r2 = pymapdl_reader.read_binary(shell181_2020r2)
+    rst_2021r1 = pymapdl_reader.read_binary(shell181_2021r1)
+
+    for key in rst_2020r2.section_data:
+        assert np.allclose(rst_2020r2.section_data[key],
+                           rst_2021r1.section_data[key])
+
 def test_overwrite(tmpdir):
     tmp_path = str(tmpdir.mkdir("tmpdir"))
     rst = pymapdl_reader.read_binary(copy(examples.rstfile, tmp_path))