Update to 0.21.0

Author: dprohe

docs/source/example_problems/substructuring.ipynb

@@ -1335,7 +1335,7 @@
     "tags": []
    },
    "source": [
-    "In the previous subsection, we saw the typical substructuring workflow: concatenate the :py:class:`System<sdynpy.core.sdynpy_system.System>` objects together, select degrees of freedom by closest position, apply constraints, and visualize results.  However, because this type of substructuring calculation, where degrees of freedom already exist co-located spatially in the two systems, is common, SDynPy has a shortcut method: \\subspos.  This method automatically searches out common degrees of freedom in the :py:class:`System<sdynpy.core.sdynpy_system.System>` objects and applies constraints between them, meaning the entirety of Case 1a can be performed in a single function call."
+    "In the previous subsection, we saw the typical substructuring workflow: concatenate the :py:class:`System<sdynpy.core.sdynpy_system.System>` objects together, select degrees of freedom by closest position, apply constraints, and visualize results.  However, because this type of substructuring calculation, where degrees of freedom already exist co-located spatially in the two systems, is common, SDynPy has a shortcut method: `substructure_by_position`.  This method automatically searches out common degrees of freedom in the :py:class:`System<sdynpy.core.sdynpy_system.System>` objects and applies constraints between them, meaning the entirety of Case 1a can be performed in a single function call."
    ]
   },
   {

pyproject.toml

@@ -39,7 +39,7 @@ dependencies = [
     ]
 
 requires-python = '>=3.10'
-version = '0.19.0'
+version = '0.21.0'
 
 classifiers = [
     'Natural Language :: English',
@@ -78,8 +78,8 @@ all = [
 ]
 
 [project.urls]
-Documentation = 'http://structmechtools.cee-gitlab.lan/structural-dynamics-python-libraries/'
-Repository = 'https://cee-gitlab.sandia.gov/structMechTools/structural-dynamics-python-libraries'
+Documentation = 'https://sandialabs.github.io/sdynpy/'
+Repository = 'https://github.com/sandialabs/sdynpy'
 
 [tool.pytest.ini_options]
 junit_family = 'xunit2'

src/sdynpy/__init__.py

@@ -17,7 +17,7 @@
 You should have received a copy of the GNU General Public License
 along with this program.  If not, see <https://www.gnu.org/licenses/>.
 """
-from .core import coordinate, colors, array, geometry, shape, data, system, matrix_mod
+from .core import coordinate, colors, array, geometry, shape, data, system, matrix_mod, system_nl
 from .fileio import unv, uff, rattlesnake, vic, tshaker, pdf3D, escdf
 from .fem.sdynpy_exodus import Exodus, ExodusInMemory, read_sierra_matlab_map_file, read_sierra_matlab_matrix_file
 from .fem import sdynpy_beam as beam
@@ -59,6 +59,8 @@
 power_spectral_density_array = data.power_spectral_density_array
 PowerSpectralDensityArray = data.PowerSpectralDensityArray
 spectrum_array = data.spectrum_array
+ShockResponseSpectrumArray = data.ShockResponseSpectrumArray
+shock_response_spectrum_array = data.shock_response_spectrum_array
 SpectrumArray = data.SpectrumArray
 GUIPlot = data.GUIPlot
 CPSDPlot = data.CPSDPlot
@@ -67,3 +69,4 @@
 matrix_plot = correlation.matrix_plot
 Matrix = matrix_mod.Matrix
 matrix = matrix_mod.matrix
+SystemNL = system_nl.SystemNL

src/sdynpy/core/__init__.py

@@ -24,3 +24,4 @@
 from . import sdynpy_array as array
 from . import sdynpy_system as system
 from . import sdynpy_matrix as matrix_mod
+from . import sdynpy_nonlinear_system as system_nl

src/sdynpy/core/sdynpy_coordinate.py

@@ -246,6 +246,40 @@ def offset_node_ids(self, offset_value):
         output = self.copy()
         output.node += offset_value
         return output
+    
+    def find_indices(self,coordinate):
+        """Finds the indices and signs of the specified coordinates
+
+        Parameters
+        ----------
+        coordinate : CoordinateArray
+            A CoordinateArray containing the coordinates to find.
+
+        Returns
+        -------
+        indices : tuple of ndarray
+            An tuple of arrays of indices into this CoordinateArray for each entry in
+            `coordinate`, which can be used directly to index this CoordinateArray.
+        signs : ndarray
+            An array of sign flips between this CoordinateArray at the indices in
+            `indices` for each entry in `coordinate`
+
+        Raises
+        ------
+        ValueError
+            If more than one coordinate matches, or no coordinate matches.
+        """
+        indices = np.zeros((self.ndim,)+coordinate.shape,dtype=int)
+        signs = np.zeros(coordinate.shape,dtype=int)
+        for check_index,coord in coordinate.ndenumerate():
+            index = np.array(np.where(abs(self)==abs(coord)))
+            if index.shape[-1] == 0:
+                raise ValueError(f'Coordinate {coordinate} not found in {self}')
+            if index.shape[-1] > 1:
+                raise ValueError(f'Coordinate {coordinate} found {index.shape[-1]} times in {self}')
+            indices[(slice(None),)+check_index] = index[...,0]
+            signs[check_index] = self[tuple(index)].sign()*coord.sign()
+        return tuple(indices),signs
 
     @classmethod
     def from_matlab_cellstr(cls, cellstr_data):