visualdynamics 0.1.0a41__tar.gz → 0.1.0a42__tar.gz

{visualdynamics-0.1.0a41/src/visualdynamics.egg-info → visualdynamics-0.1.0a42}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: visualdynamics
-Version: 0.1.0a41
+Version: 0.1.0a42
 Summary: Structural Dynamics Visualization: a units-aware structural dynamics analysis toolset
 Author-email: Brandon Zwink <bzwink@gmail.com>
 License-Expression: MIT
@@ -73,11 +73,12 @@ file.
 
 **This is an alpha.** Every launch opens with a disclaimer to
 acknowledge: the software is offered as is, and its numbers are the
-user's to check against a tool they already trust. The project file,
-`.vdyn`, is provisional: the next alpha will save projects in the
-Engineering Sciences Common Data Format (`.h5`), a public standard,
-and will still open `.vdyn` files; the first release that is not an
-alpha will not. Open each `.vdyn` once in that alpha and save it again.
+user's to check against a tool they already trust. Projects are
+saved as `.vdyn` files, and Engineering Sciences Common Data Format
+(`.h5`) files are imported and exported for exchanging data with
+other tools. While Visual Dynamics is an alpha the `.vdyn` format may
+still change; the promise that a later version opens every earlier
+file starts with the first release that is not an alpha.
 
 **[Core principles](PRINCIPLES.md)** — thirteen principles
 that shape every part of it: data you can see and handle, a workflow

{visualdynamics-0.1.0a41 → visualdynamics-0.1.0a42}/README.md

@@ -16,11 +16,12 @@ file.
 
 **This is an alpha.** Every launch opens with a disclaimer to
 acknowledge: the software is offered as is, and its numbers are the
-user's to check against a tool they already trust. The project file,
-`.vdyn`, is provisional: the next alpha will save projects in the
-Engineering Sciences Common Data Format (`.h5`), a public standard,
-and will still open `.vdyn` files; the first release that is not an
-alpha will not. Open each `.vdyn` once in that alpha and save it again.
+user's to check against a tool they already trust. Projects are
+saved as `.vdyn` files, and Engineering Sciences Common Data Format
+(`.h5`) files are imported and exported for exchanging data with
+other tools. While Visual Dynamics is an alpha the `.vdyn` format may
+still change; the promise that a later version opens every earlier
+file starts with the first release that is not an alpha.
 
 **[Core principles](PRINCIPLES.md)** — thirteen principles
 that shape every part of it: data you can see and handle, a workflow

{visualdynamics-0.1.0a41 → visualdynamics-0.1.0a42}/src/visualdynamics/__init__.py

@@ -90,7 +90,7 @@ from .units import (
 #: alpha, said in the version itself so the dialog, the About text, the
 #: installers and latest.json all say the same thing, and the first
 #: non-alpha release is offered as an upgrade (Brandon, 2026-09-11)
-__version__ = '0.1.0a41'
+__version__ = '0.1.0a42'
 
 
 def launch_gui(*paths: str, in_process: bool | None = None) -> Any:

{visualdynamics-0.1.0a41 → visualdynamics-0.1.0a42}/src/visualdynamics/core/data.py

@@ -229,6 +229,18 @@ class _LogAbscissa:
         return cls._log_abscissa_default
 
 
+def _recast(obj: Any, cls: type) -> Any:
+    """`obj` as an instance of `cls`, every field copied — what reading
+    a PSD as a specification, or back, changes is the kind and nothing
+    else, so nothing is rebuilt through a constructor that might not
+    know a field (a declared-later unit, a dimension hint, a band)."""
+    import copy
+
+    out = cls.__new__(cls)
+    out.__dict__.update(copy.deepcopy(obj.__dict__))
+    return out
+
+
 class DataArray:
     """Base class; use a concrete subclass (TimeHistory, Spectrum, Frf, Psd).
 
@@ -3000,6 +3012,30 @@ class Psd(_Bands, DataArray):
             self.ordinate: np.ndarray = np.ascontiguousarray(self.ordinate.real)
         self._take_bandwidth(bandwidth)
 
+    def as_specification(self) -> Specification:
+        """This PSD read as a specification: the same numbers, records,
+        units and bands, held as the requirement a measurement is judged
+        against, with no limits.
+
+        Nothing in a file says a PSD is a requirement when the format
+        has no word for one — ESCDF's data types do not, and a target
+        written by another program arrives as a PSD with only its name
+        to tell — so the person who knows says so (Brandon,
+        2026-10-09). A banded PSD steps, as a banded specification does;
+        otherwise its frequencies decide whether it reads as a
+        controller's lines or as breakpoints (`Specification.reading_of`).
+
+        Returns
+        -------
+        Specification
+            A copy; this PSD is unchanged.
+        """
+        spec = _recast(self, Specification)
+        spec.limits = {}
+        spec.interpolation = ('bin' if self.bandwidth is not None
+                              else Specification.reading_of(self.abscissa))
+        return spec
+
     def to_octave(self, per_octave: int | None = None,
                   low: float | None = None,
                   high: float | None = None) -> Psd:
@@ -3433,6 +3469,31 @@ class Specification(Bounded, Psd):
         with np.errstate(invalid='ignore'):
             return (np.isfinite(magnitude) & (magnitude > 0.0)).any(axis=0)
 
+    def as_psd(self) -> Psd:
+        """This specification read as a plain PSD — `Psd.as_specification`
+        undone, for a PSD read as a requirement by mistake. Its reading
+        of its own points is kept.
+
+        Returns
+        -------
+        Psd
+            A copy; this specification is unchanged.
+
+        Raises
+        ------
+        ValueError
+            When it has limits: a PSD has nowhere to keep them, and
+            dropping them unasked would lose what the file said.
+            `limit(name)` hands one out as a PSD of its own.
+        """
+        if self.limits:
+            raise ValueError(
+                'it carries warning or abort limits, which a PSD cannot '
+                'hold; take one out with limit(name) instead')
+        plain = _recast(self, Psd)
+        plain.__dict__.pop('limits', None)
+        return plain
+
     def to_octave(self, per_octave: int | None = None,
                   low: float | None = None,
                   high: float | None = None) -> Specification:

{visualdynamics-0.1.0a41 → visualdynamics-0.1.0a42}/src/visualdynamics/core/filters.py

@@ -54,6 +54,7 @@ double integration invents are set out.
 from __future__ import annotations
 
 from dataclasses import dataclass, replace
+from functools import cache
 from typing import TYPE_CHECKING, Any
 
 import numpy as np
@@ -94,11 +95,27 @@ def threaded_sosfiltfilt(sos: Any, rows: Any, axis: int = -1) -> np.ndarray:
     if rows.ndim < 2 or axis not in (-1, rows.ndim - 1) or rows.shape[0] < 2 \
             or rows.size < 1_000_000:
         return sosfiltfilt(sos, rows, axis=axis)
+    return np.stack(list(_filter_pool().map(
+        lambda row: sosfiltfilt(sos, row), rows)))
+
+
+@cache
+def _filter_pool() -> Any:
+    """The threads `threaded_sosfiltfilt` shares, made once and kept.
+
+    A pool per call started up to eight threads every time the filter
+    preview redrew. One gate run stalled for an hour inside exactly
+    that: the main thread in `Thread.start`, waiting on a new worker
+    that never said it had begun, from a render of the filter view
+    (2026-10-09; CPython 3.13.2, not reproduced in 1 200 calls). A kept
+    pool starts its threads once a session, so no redraw waits on a
+    thread being born — and skips the start-up besides.
+    """
     import os
     from concurrent.futures import ThreadPoolExecutor
 
-    with ThreadPoolExecutor(min(8, os.cpu_count() or 1)) as pool:
-        return np.stack(list(pool.map(lambda row: sosfiltfilt(sos, row), rows)))
+    return ThreadPoolExecutor(min(8, os.cpu_count() or 1),
+                              thread_name_prefix='visualdynamics-filter')
 
 
 @dataclass(frozen=True, kw_only=True)

{visualdynamics-0.1.0a41 → visualdynamics-0.1.0a42}/src/visualdynamics/gui/disclaimer.py

@@ -35,12 +35,10 @@ TEXT = (
     'set of tests. Numbers it produces — modal parameters, spectra, '
     'levels, comparisons against a specification — are to be verified by '
     'independent means before anything depends on them.</p>'
-    '<p>The project file, <code>.vdyn</code>, is provisional. The next '
-    'alpha will save projects in the Engineering Sciences Common Data '
-    'Format (<code>.h5</code>), a public standard, and will still '
-    'open <code>.vdyn</code> files; the first release that is not an '
-    'alpha will not. Open each <code>.vdyn</code> once in that alpha and '
-    'save it again.</p>'
+    '<p>Projects are saved as <code>.vdyn</code> files. While this is '
+    'an alpha the <code>.vdyn</code> format may still change; the '
+    'promise that a later version opens every earlier file starts with '
+    'the first release that is not an alpha.</p>'
     '<p>Anything that looks wrong is worth reporting: '
     '<a href="mailto:contact@visualdynamics.org">contact@visualdynamics.org</a>.</p>'
 )

{visualdynamics-0.1.0a41 → visualdynamics-0.1.0a42}/src/visualdynamics/gui/icons.py

@@ -745,6 +745,30 @@ def _draw_curves(painter, color):
         for i in range(45)]))
 
 
+def _spectrum(dy=0.0):
+    """A PSD's shape across the glyph: a rise to a plateau and a fall."""
+    return _curve([(8 + i, 40 + dy - 14 * math.sin(min(i, 40) / 40 * math.pi))
+                   for i in range(49)])
+
+
+def _draw_as_specification(painter, color):
+    """The spectrum between two dashed limits: read this PSD as the
+    requirement it is."""
+    painter.setPen(_pen(color, 5))
+    painter.drawPath(_spectrum())
+    dashed = QPen(QColor(color), 3)
+    dashed.setStyle(Qt.PenStyle.DashLine)
+    painter.setPen(dashed)
+    painter.drawPath(_spectrum(-9))
+    painter.drawPath(_spectrum(9))
+
+
+def _draw_as_psd(painter, color):
+    """The spectrum alone: read this specification as a plain PSD."""
+    painter.setPen(_pen(color, 5))
+    painter.drawPath(_spectrum())
+
+
 def _draw_map(painter, color):
     """Banded rows: every channel at once, coherence as color."""
     painter.setPen(Qt.PenStyle.NoPen)
@@ -1723,7 +1747,9 @@ def control_icon(name: str) -> QIcon:
      'set_view': _draw_set_view, 'reset_view': _draw_reset_view,
      'sine': _draw_sine,
      'refresh': _draw_refresh,
-     'copy': _draw_copy, 'movie': _draw_movie}.get(
+     'copy': _draw_copy, 'movie': _draw_movie,
+     'as_specification': _draw_as_specification,
+     'as_psd': _draw_as_psd}.get(
          name, _draw_default)(painter, CONTROL_COLOR)
     painter.end()
     return QIcon(pixmap)

{visualdynamics-0.1.0a41 → visualdynamics-0.1.0a42}/src/visualdynamics/gui/main_window.py

@@ -3076,6 +3076,9 @@ class MainWindow(QMainWindow):
         ('add_plane', 'Add Plane', 'plane', 'add_plane_act'),
         ('add_block', 'Add Block', 'block', 'add_block_act'),
         ('set_view', 'Set Default View', 'set_view', 'set_view_act'),
+        ('read_as_specification', 'Read as Specification',
+         'as_specification', 'read_as_specification_act'),
+        ('read_as_psd', 'Read as PSD', 'as_psd', 'read_as_psd_act'),
     )
 
     def acts_for(self, names=None):
@@ -11807,6 +11810,38 @@ class MainWindow(QMainWindow):
             f'{added}: {shapes.num_shapes} modes to {top:g} Hz, {rigid} '
             f'rigid, {percent:g}% damping — linked to {name}')
 
+    def read_as_specification_act(self) -> None:
+        """The bar's Read as Specification: the selected PSD becomes the
+        requirement it is, in place (`Project.read_as_specification`)."""
+        self._read_as(Psd, 'Select a PSD to read as a specification',
+                      self.project.read_as_specification,
+                      'read as a specification — compare a measurement '
+                      'against it by selecting both')
+
+    def read_as_psd_act(self) -> None:
+        """The bar's Read as PSD: the selected specification is plain
+        data again, in place (`Project.read_as_psd`)."""
+        self._read_as(Specification, 'Select a specification to read as a PSD',
+                      self.project.read_as_psd, 'read as a PSD')
+
+    def _read_as(self, kind, refusal, verb, said) -> None:
+        """Swap the selected object for its other reading, under the
+        same name: the tree row, its grid and the view follow, the way a
+        recompute's in-place swap does."""
+        acted = self._act_on(kind, refusal, verb)
+        if acted is None:
+            return
+        name = acted[0]
+        obj = self.objects[name]
+        item = self._item_for_object(name)
+        if item is not None:
+            self._refresh_item(item, obj)
+            item.setToolTip(0, self._object_tooltip(name, obj))
+            self._build_children(item, obj, name)
+        self._report_content_changed()
+        self.render_current()
+        self._show_status(f'{name} {said}')
+
     def integrate_history(self) -> None:
         """One integration of the selected time history — acceleration
         to velocity, velocity to displacement — added beside it."""

{visualdynamics-0.1.0a41 → visualdynamics-0.1.0a42}/src/visualdynamics/project.py

@@ -2501,6 +2501,53 @@ class Project(dict):
         return self._derive(source, shapes, name or f'{source} Modes',
                             recipe=('solve_modes', params))
 
+    def read_as_specification(self, name: Any) -> str:
+        """Read this PSD as a specification — the requirement it is
+        judged against — under its own name, its object group and every
+        reference to it kept.
+
+        A file whose format has no word for a requirement (ESCDF's data
+        types do not) can only say so by its name, and a name is no
+        evidence; the person who knows says so instead.
+
+        Parameters
+        ----------
+        name : str or Psd
+            The PSD.
+
+        Returns
+        -------
+        str
+            Its name.
+        """
+        name = self.name_of(name)
+        obj = self[name]
+        if not isinstance(obj, Psd) or isinstance(obj, Specification):
+            raise TypeError(f'{name} is not a PSD read as data')
+        self[name] = obj.as_specification()
+        return name
+
+    def read_as_psd(self, name: Any) -> str:
+        """Read this specification as a plain PSD again, under its own
+        name — `read_as_specification` undone.
+
+        Parameters
+        ----------
+        name : str or Specification
+            The specification.
+
+        Returns
+        -------
+        str
+            Its name.
+        """
+        name = self.name_of(name)
+        obj = self[name]
+        if not isinstance(obj, Specification):
+            raise TypeError(f'{name} is not a specification')
+        self[name] = obj.as_psd()
+        return name
+
     def scale_family(self, measured: Any) -> list[str]:
         """The measured PSDs that share one comparison scale: the object
         and every non-specification PSD in its object group — its
@@ -4039,6 +4086,12 @@ _VERB_APPLIES: tuple = (
         for other in p.values()))),
     # a specification bands too, limits and all (Brandon, 2026-09-18)
     ('compute_octave', lambda p, o: isinstance(o, Psd)),
+    # what a file could not say: a PSD that is a requirement, and back
+    # — but not a specification with limits, which a PSD cannot hold
+    ('read_as_specification', lambda p, o: (
+        isinstance(o, Psd) and not isinstance(o, Specification))),
+    ('read_as_psd', lambda p, o: (isinstance(o, Specification)
+                                  and not o.has_limits)),
     ('fit_modes', lambda p, o: isinstance(o, Frf)),
     ('transform', lambda p, o: _reads_as(p, o, 'physical')),
     ('expand', lambda p, o: _reads_as(p, o, 'modal')),
@@ -4142,6 +4195,7 @@ _JOURNALED_VERBS = (
     'merge_coincident_nodes', 'new_geometry', 'add_plane', 'add_block',
     'tie_elements',
     'merge_groups', 'set_view', 'set_comparison_scale',
+    'read_as_specification', 'read_as_psd',
     'author_specification',
     'transform', 'expand', 'project_onto_basis', 'match_modes',
     'extract_sine', 'refresh', 'refresh_stale', 'work_up',

{visualdynamics-0.1.0a41 → visualdynamics-0.1.0a42}/src/visualdynamics/viz/paired.py

@@ -29,6 +29,7 @@ from ._quiet import one_render
 from .waterfall import (
     POINT_BUDGET,
     STAGE,
+    drawn_span,
     finish_scene,
     place_camera,
     waterfall_arrays,
@@ -143,10 +144,13 @@ def paired_stage_curves(loud: Any, quiet: Any,
 
     stations = arrays['stations']
     n = len(stations)
-    xs = [float(np.nanmin(c[which][0])) for c in stations
-          for which in ('loud', 'quiet') if len(c[which][0])]
-    xe = [float(np.nanmax(c[which][0])) for c in stations
-          for which in ('loud', 'quiet') if len(c[which][0])]
+    # where something is drawn, not every line a curve is carried on
+    # (`waterfall.drawn_span`)
+    spans = [drawn_span([c[which]]) for c in stations
+             for which in ('loud', 'quiet')
+             if len(c[which][0]) and np.isfinite(c[which][1]).any()]
+    xs = [low for low, _high in spans]
+    xe = [high for _low, high in spans]
     zs = np.concatenate(
         [c[which][1][np.isfinite(c[which][1])] for c in stations
          for which in ('loud', 'quiet')]) if stations else np.array([0.0, 1.0])

{visualdynamics-0.1.0a41 → visualdynamics-0.1.0a42}/src/visualdynamics/viz/waterfall.py

@@ -412,6 +412,29 @@ def nice_axis(low: float, high: float,
     return float(wide_low), float(wide_high), intervals + 1
 
 
+def drawn_span(curves: Any, default: tuple[float, float] = (0.0, 1.0)
+               ) -> tuple[float, float]:
+    """(low, high) of the abscissa where something is drawn: points
+    whose level is finite, never the whole abscissa a curve is carried
+    on. A PSD written as a controller writes a specification — every
+    line to Nyquist, zero outside its band — draws only over the band
+    (zero is a gap on a log level), and a stage that spanned all its
+    lines drew 20 to 400 Hz as a sliver of a 0 to 2048 Hz floor
+    (Brandon, 2026-10-09)."""
+    lows, highs = [], []
+    for cx, cz in curves:
+        cx, cz = np.asarray(cx, dtype=float), np.asarray(cz, dtype=float)
+        if cx.size != cz.size:
+            cx = cx[:cz.size]
+        shown = cx[np.isfinite(cz) & np.isfinite(cx)] if cz.size else cx[:0]
+        if shown.size:
+            lows.append(float(shown.min()))
+            highs.append(float(shown.max()))
+    if not lows:
+        return default
+    return min(lows), max(highs)
+
+
 def stage_curves(arrays: dict[str, Any],
                  ordinate_limits: tuple[float, float] | None = None
                  ) -> dict[str, Any]:
@@ -432,10 +455,7 @@ def stage_curves(arrays: dict[str, Any],
     curves = arrays['curves']
     n = len(curves)
     sx, sy, sz = STAGE
-    x0 = min((float(np.nanmin(cx)) for cx, _cz in curves if len(cx)),
-             default=0.0)
-    x1 = max((float(np.nanmax(cx)) for cx, _cz in curves if len(cx)),
-             default=1.0)
+    x0, x1 = drawn_span(curves)
     # widened to round tick values where that is nearly free, before
     # anything is normalized against the range — the axis, the curves,
     # the marks and the report's canvas all read these extents, so the

{visualdynamics-0.1.0a41 → visualdynamics-0.1.0a42/src/visualdynamics.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: visualdynamics
-Version: 0.1.0a41
+Version: 0.1.0a42
 Summary: Structural Dynamics Visualization: a units-aware structural dynamics analysis toolset
 Author-email: Brandon Zwink <bzwink@gmail.com>
 License-Expression: MIT
@@ -73,11 +73,12 @@ file.
 
 **This is an alpha.** Every launch opens with a disclaimer to
 acknowledge: the software is offered as is, and its numbers are the
-user's to check against a tool they already trust. The project file,
-`.vdyn`, is provisional: the next alpha will save projects in the
-Engineering Sciences Common Data Format (`.h5`), a public standard,
-and will still open `.vdyn` files; the first release that is not an
-alpha will not. Open each `.vdyn` once in that alpha and save it again.
+user's to check against a tool they already trust. Projects are
+saved as `.vdyn` files, and Engineering Sciences Common Data Format
+(`.h5`) files are imported and exported for exchanging data with
+other tools. While Visual Dynamics is an alpha the `.vdyn` format may
+still change; the promise that a later version opens every earlier
+file starts with the first release that is not an alpha.
 
 **[Core principles](PRINCIPLES.md)** — thirteen principles
 that shape every part of it: data you can see and handle, a workflow