jmrp.io / PyOctaveBand

0

[READ-ONLY] Mirror of https://github.com/jmrplens/PyOctaveBand. [Python3] Octave-Band and Fractional Octave-Band filter. For signal in time domain. jmrplens.github.io/PyOctaveBand/

acoustics audio filter frequency frequency-analysis frequency-domain octave python3 signal time-domain

0

PyOctaveBand / docs /

at fix/audit-p1 29 files

José M. Requena Plens fix: first pass of audit findings (docs, figures, plotting) 15d ago

fix: first pass of audit findings (docs, figures, plotting) Documentation: - fix the Spanish filter-figure regression: _showfilter gains a close flag and the generator saves through themed_path after drawing, so the translation pass runs on the finished figure; restore the two _ES_EXACT keys removed in #82 (they are emitted by filter_design, not by the script) - repair the broken building-acoustics.md link in the enclosed-space guide (the repo split lives in room-acoustics.md) - make the multiple-shock and intensity snippets self-contained: define az/fs with a synthetic shock train (verified outputs 20.94 m/s2, R=0.46, probability 0.03) and import numpy with stated outputs (F2=3.41, Ld=8.0), in the repo doc and both site languages - decimal points in the English Formula texts (1.07, 55.3) that had kept the Spanish comma Figures: - stop drawing data series in the grid colour (invisible on light backgrounds): time-weighting input burst, hearing-threshold 20 yr and NIPTS 10 yr curves, and the unweighted/partial bars now use a visible neutral grey; regenerate the 15 affected figure sets Plotting layer: - plot_facade_insulation forwards user kwargs to the primary D2m,nT curve only, so label=/color= no longer raise or merge the four curves; the kwargs-forwarding test now covers the facade plot - harden three edge cases: all-masked SII (zero contribution no longer divides 0/0), empty per-impulse prominence, and nearest-band lookup for the NC governing marker instead of float equality - sharpness normalization sanity guards raise RuntimeError instead of AssertionError

2 weeks ago

EN 12354-6: sound absorption in enclosed spaces (#107) * feat: EN 12354-6 sound absorption in enclosed spaces Add the EN 12354-6:2003 prediction of a room's total equivalent sound absorption area and reverberation time from the absorption of its surfaces and objects (the normative Clause 4 model). The absorption area sums the surface contributions, the object areas and the air absorption (Formula 1), with the air term Aair = 4*m*V*(1 - psi) (Formula 2, Table 1), the object fraction psi = sum(Vobj)/V (Formula 3) and the empirical hard-object area Aobj = Vobj^(2/3) (Formula 4). The reverberation time follows as T = 55.3/c0 * V*(1 - psi)/A (Formula 5, the 0.16 factor at c0 = 345.6 m/s). The informative Annex D method for irregular spaces is out of scope. New module phonometry.en12354_6 with equivalent_absorption_area, object_fraction, hard_object_absorption, air_absorption_area, reverberation_time and the per-band enclosed_space_reverberation returning a ReverberationResult with .plot(). Conformance report gains an enclosed-space absorption domain (2 checks, 90 total). Adds the enclosed-space-absorption guide (repo + EN/ES site) with a figure and a flow diagram. Validated digit-exact against the three worked cases of Annex E (A = 2.26/5.03/10.21 m2, T = 2.1/0.9/0.5 s, and the air-absorption note Aair = 0.12 m2 -> T = 2.0 s). * refactor: address review feedback on EN 12354-6 - validation helpers _require_positive/_require_fraction: single error literal and no negated comparison, still rejecting NaN (SonarCloud S1192/S1940) - object_fraction: reject negative object volumes and a fraction >= 1 (Gemini) - air_absorption_area: validate the object fraction range and a non-negative attenuation coefficient (Gemini) - equivalent_absorption_area: reject negative air_area, absorption coefficients and object areas (Gemini) - reverberation_time: validate the object fraction range (Gemini) - enclosed_space_reverberation: clear error when a built-in air_condition is combined with non-standard frequencies (Gemini) - diagram + alt text: the speed-of-sound symbol is c0, not co (CodeRabbit); the diagram now renders it as a proper subscript - tests: cover the speed_of_sound guard (CodeRabbit), NaN volume, the physical-range validations and the standard-bands guard * docs: label the object-array sum in the Formula 1 prose The third summand of Formula 1 is the standard's object arrays k (groups of identical objects treated as an absorbing surface), not a duplicated surface term; name all three indices in the surrounding prose so the equation reads unambiguously (CodeRabbit).

2 weeks ago

api-reference.md

feat: plottable impulse-response results and ISO 3382 measurement setup (#98) impulse_response / mls_impulse_response now return an ImpulseResponseResult with .plot() (waveform + log-magnitude envelope + Schroeder decay), staying a drop-in for the raw IR array via __array__ so existing callers are unchanged. Adds plot_excitation() for the sweep and MLS excitation signals, and an ISO 3382 measurement-setup diagram (source/microphone positions and ISO 3382-1 placement rules, plus the ISO 3382-2 minimum-position grades). Room-acoustics guide §1 (repo + EN/ES site) shows the new signal/IR plots and the diagram.

2 weeks ago

block-processing.md

Docs usability: measurement-role snippet clarity, one-line result plots, responsive figures (#86) * fix: responsive figure sizing — full-width on mobile, 92% plots on desktop (plan 16 T4) Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: snippet clarity wave — recording/excitation convention (plan 16 T1) Apply the recording/excitation naming convention from the P16 audit across the calibration, intensity, block-processing, time-weighting, levels, weighting and room-acoustics guides (docs/ + EN guides + ES guides). - calibration: ref_signal->calibrator_recording, signal->recording (J1 cascade) - intensity: probe-recording synth + swap-in disclosure on p1/p2 - block-processing / time-weighting: audio_blocks provenance comments - levels / time-weighting / weighting: signal->recording, provenance upgraded to convention - room-acoustics: play/record framing on sweep/recorded; synthetic-IR stand-in note EN docs/guides kept code-identical; ES mirrored with UNE wording. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * feat: one-line canonical plots on result objects (plan 16 T2) Add a thin plot(ax=None, **kwargs) method to each result dataclass that renders the essence of its guide's canonical figure in one line, following the sklearn/statsmodels/librosa pattern. matplotlib stays a soft dependency: importing phonometry and computing works without it; .plot() raises a clear ImportError with install guidance. Rendering lives in a new lazy _plotting module; dataclass methods are thin delegates. Covers ZwickerLoudness, STIResult, WeightedRatingResult/ImpactRatingResult, RoomAcousticsResult, SoundPowerResult/ReverberationSoundPowerResult/ SoundPowerIntensityResult, IntensityResult and a new frozen DecayCurve (returned by decay_curve, tuple-unpackable for backward compat). Rating results gain band_centers/measured/shifted_reference fields for the figure. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: plot kwargs forwarding, deterministic invalid-band test, octave-impact annotation - _plot_rating now accepts **kwargs and forwards to the measured curve, so result.plot(linewidth=2) no longer raises TypeError; parametrized test covers every public .plot() forwarding a benign kwarg to its primary artist. - Replace skip-if-all-valid room-acoustics invalid-band test with a directly constructed RoomAcousticsResult (500 Hz band flagged invalid); assert exactly the flagged bars are hatched and greyed. - plot_impact_rating keeps the shifted-reference curve normatively honest (ref - shift) and marks/annotates the 500 Hz read value with the ISO 717-2 Clause 4.3.2 -5 dB octave rule instead of distorting the curve. - Move _bar_width next to its only consumer, plot_intensity. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: .plot() one-liners and figure-code collapsibles (plan 16 T3) Add one-line res.plot() mentions to every guide section whose result type gained .plot() in T2, and figure-code collapsibles (both the one-liner and the manual matplotlib path) under the five .plot()-backed data figures. Document .plot() availability and the DecayCurve return type in the API reference and CHANGELOG. EN docs + EN/ES site; llms-full.txt regenerated. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: collapsibles for STI-vs-T60 and open-plan figures Controller follow-up to plan 16 T3: the two figures are reproducible with public API only. STI-vs-T60 sweeps sti_from_impulse_response over synthetic exponential IRs at a 10-point T60 grid; the open-plan collapsible rebuilds the D2,S regression and rD/rP crossings from the OpenPlanResult fields. EN docs + EN/ES site; executed on Agg; llms-full.txt regenerated. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: final-review polish — result-field docs, comment hygiene, naming Document band_centers/measured/shifted_reference on WeightedRatingResult and ImpactRatingResult (CHANGELOG + EN/ES API tables), strip internal plan tags from theme-images.css and test_result_plots.py, note that [plot] also powers result .plot() methods, tighten the plot_zwicker_loudness docstring to say the two-panel layout only occurs when ax is None, sentence-case the intensity provenance comments, and rename the shadowed impact rating var to res_imp on the room-acoustics pages. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: PR 86 review round 1 — log-axis bars, style kwargs overrides, snippet physics and notation Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf

2 weeks ago

Docs usability: measurement-role snippet clarity, one-line result plots, responsive figures (#86) * fix: responsive figure sizing — full-width on mobile, 92% plots on desktop (plan 16 T4) Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: snippet clarity wave — recording/excitation convention (plan 16 T1) Apply the recording/excitation naming convention from the P16 audit across the calibration, intensity, block-processing, time-weighting, levels, weighting and room-acoustics guides (docs/ + EN guides + ES guides). - calibration: ref_signal->calibrator_recording, signal->recording (J1 cascade) - intensity: probe-recording synth + swap-in disclosure on p1/p2 - block-processing / time-weighting: audio_blocks provenance comments - levels / time-weighting / weighting: signal->recording, provenance upgraded to convention - room-acoustics: play/record framing on sweep/recorded; synthetic-IR stand-in note EN docs/guides kept code-identical; ES mirrored with UNE wording. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * feat: one-line canonical plots on result objects (plan 16 T2) Add a thin plot(ax=None, **kwargs) method to each result dataclass that renders the essence of its guide's canonical figure in one line, following the sklearn/statsmodels/librosa pattern. matplotlib stays a soft dependency: importing phonometry and computing works without it; .plot() raises a clear ImportError with install guidance. Rendering lives in a new lazy _plotting module; dataclass methods are thin delegates. Covers ZwickerLoudness, STIResult, WeightedRatingResult/ImpactRatingResult, RoomAcousticsResult, SoundPowerResult/ReverberationSoundPowerResult/ SoundPowerIntensityResult, IntensityResult and a new frozen DecayCurve (returned by decay_curve, tuple-unpackable for backward compat). Rating results gain band_centers/measured/shifted_reference fields for the figure. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: plot kwargs forwarding, deterministic invalid-band test, octave-impact annotation - _plot_rating now accepts **kwargs and forwards to the measured curve, so result.plot(linewidth=2) no longer raises TypeError; parametrized test covers every public .plot() forwarding a benign kwarg to its primary artist. - Replace skip-if-all-valid room-acoustics invalid-band test with a directly constructed RoomAcousticsResult (500 Hz band flagged invalid); assert exactly the flagged bars are hatched and greyed. - plot_impact_rating keeps the shifted-reference curve normatively honest (ref - shift) and marks/annotates the 500 Hz read value with the ISO 717-2 Clause 4.3.2 -5 dB octave rule instead of distorting the curve. - Move _bar_width next to its only consumer, plot_intensity. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: .plot() one-liners and figure-code collapsibles (plan 16 T3) Add one-line res.plot() mentions to every guide section whose result type gained .plot() in T2, and figure-code collapsibles (both the one-liner and the manual matplotlib path) under the five .plot()-backed data figures. Document .plot() availability and the DecayCurve return type in the API reference and CHANGELOG. EN docs + EN/ES site; llms-full.txt regenerated. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: collapsibles for STI-vs-T60 and open-plan figures Controller follow-up to plan 16 T3: the two figures are reproducible with public API only. STI-vs-T60 sweeps sti_from_impulse_response over synthetic exponential IRs at a 10-point T60 grid; the open-plan collapsible rebuilds the D2,S regression and rD/rP crossings from the OpenPlanResult fields. EN docs + EN/ES site; executed on Agg; llms-full.txt regenerated. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: final-review polish — result-field docs, comment hygiene, naming Document band_centers/measured/shifted_reference on WeightedRatingResult and ImpactRatingResult (CHANGELOG + EN/ES API tables), strip internal plan tags from theme-images.css and test_result_plots.py, note that [plot] also powers result .plot() methods, tighten the plot_zwicker_loudness docstring to say the two-panel layout only occurs when ax is None, sentence-case the intensity provenance comments, and rename the shadowed impact rating var to res_imp on the room-acoustics pages. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: PR 86 review round 1 — log-axis bars, style kwargs overrides, snippet physics and notation Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf

2 weeks ago

enclosed-space-absorption.md

fix: first pass of audit findings (docs, figures, plotting) Documentation: - fix the Spanish filter-figure regression: _showfilter gains a close flag and the generator saves through themed_path after drawing, so the translation pass runs on the finished figure; restore the two _ES_EXACT keys removed in #82 (they are emitted by filter_design, not by the script) - repair the broken building-acoustics.md link in the enclosed-space guide (the repo split lives in room-acoustics.md) - make the multiple-shock and intensity snippets self-contained: define az/fs with a synthetic shock train (verified outputs 20.94 m/s2, R=0.46, probability 0.03) and import numpy with stated outputs (F2=3.41, Ld=8.0), in the repo doc and both site languages - decimal points in the English Formula texts (1.07, 55.3) that had kept the Spanish comma Figures: - stop drawing data series in the grid colour (invisible on light backgrounds): time-weighting input burst, hearing-threshold 20 yr and NIPTS 10 yr curves, and the unweighted/partial bars now use a visible neutral grey; regenerate the 15 affected figure sets Plotting layer: - plot_facade_insulation forwards user kwargs to the primary D2m,nT curve only, so label=/color= no longer raise or merge the four curves; the kwargs-forwarding test now covers the facade plot - harden three edge cases: all-masked SII (zero contribution no longer divides 0/0), empty per-impulse prominence, and nearest-band lookup for the NC governing marker instead of float equality - sharpness normalization sanity guards raise RuntimeError instead of AssertionError

2 weeks ago

filter-banks.md

Audit-driven accuracy, robustness and didactic improvements across the library (#82) * fix: audit wave A — inter-sample LCpeak, cheby2 class-1 default, 144 kHz weighting target Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: audit wave A — decay validity thresholds, variant sharpness anchors, N5 phase, STIPA warning, Farina guard Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: audit wave B — self-contained snippets, tables, sharpness/open-plan figures, orphans Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: audit wave C — remaining minors and optimizations across the library Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: document the audit-wave parameters across guides and API tables Cover the public parameters added in the audit fix waves across all three documentation trees (docs/, site EN, site ES): - lc_peak(..., oversample=8): inter-sample peak recovery (levels guide + API) - calculate_sensitivity(..., narrowband=False): coherent tone estimator (calibration guide + API) - sound_intensity(..., bias_correct=False): finite-difference bias correction (intensity guide + API) - room_parameters/decay_curve(..., zero_phase=False): 125 Hz short-T bias (room-acoustics guide + API) - OctaveFilterBank/octavefilter attenuation default 60 -> 72 and cheby2 class-1 note (filter-banks guide + API) - STIPA < 15 s UserWarning (psychoacoustics guide + API) Also refresh dependent behaviour notes: ln_levels attack skip 2*tau -> 5*tau, T20/T30 validity thresholds 35/45 -> 46/54 dB, and the zero-phase broadband ~0.2-0.3 dB band-level caveat. Regenerated llms-full.txt. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: truthfulness fixes from the final audit-branch review Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: 2000 Hz, not 2 ms, in the N5/N10 docstring Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: PR 82 review round 1 — bias-correction cutoff, parabolic-peak guard, complexity, snippet self-containment Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf

2 weeks ago

getting-started.md

Docs usability: measurement-role snippet clarity, one-line result plots, responsive figures (#86) * fix: responsive figure sizing — full-width on mobile, 92% plots on desktop (plan 16 T4) Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: snippet clarity wave — recording/excitation convention (plan 16 T1) Apply the recording/excitation naming convention from the P16 audit across the calibration, intensity, block-processing, time-weighting, levels, weighting and room-acoustics guides (docs/ + EN guides + ES guides). - calibration: ref_signal->calibrator_recording, signal->recording (J1 cascade) - intensity: probe-recording synth + swap-in disclosure on p1/p2 - block-processing / time-weighting: audio_blocks provenance comments - levels / time-weighting / weighting: signal->recording, provenance upgraded to convention - room-acoustics: play/record framing on sweep/recorded; synthetic-IR stand-in note EN docs/guides kept code-identical; ES mirrored with UNE wording. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * feat: one-line canonical plots on result objects (plan 16 T2) Add a thin plot(ax=None, **kwargs) method to each result dataclass that renders the essence of its guide's canonical figure in one line, following the sklearn/statsmodels/librosa pattern. matplotlib stays a soft dependency: importing phonometry and computing works without it; .plot() raises a clear ImportError with install guidance. Rendering lives in a new lazy _plotting module; dataclass methods are thin delegates. Covers ZwickerLoudness, STIResult, WeightedRatingResult/ImpactRatingResult, RoomAcousticsResult, SoundPowerResult/ReverberationSoundPowerResult/ SoundPowerIntensityResult, IntensityResult and a new frozen DecayCurve (returned by decay_curve, tuple-unpackable for backward compat). Rating results gain band_centers/measured/shifted_reference fields for the figure. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: plot kwargs forwarding, deterministic invalid-band test, octave-impact annotation - _plot_rating now accepts **kwargs and forwards to the measured curve, so result.plot(linewidth=2) no longer raises TypeError; parametrized test covers every public .plot() forwarding a benign kwarg to its primary artist. - Replace skip-if-all-valid room-acoustics invalid-band test with a directly constructed RoomAcousticsResult (500 Hz band flagged invalid); assert exactly the flagged bars are hatched and greyed. - plot_impact_rating keeps the shifted-reference curve normatively honest (ref - shift) and marks/annotates the 500 Hz read value with the ISO 717-2 Clause 4.3.2 -5 dB octave rule instead of distorting the curve. - Move _bar_width next to its only consumer, plot_intensity. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: .plot() one-liners and figure-code collapsibles (plan 16 T3) Add one-line res.plot() mentions to every guide section whose result type gained .plot() in T2, and figure-code collapsibles (both the one-liner and the manual matplotlib path) under the five .plot()-backed data figures. Document .plot() availability and the DecayCurve return type in the API reference and CHANGELOG. EN docs + EN/ES site; llms-full.txt regenerated. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: collapsibles for STI-vs-T60 and open-plan figures Controller follow-up to plan 16 T3: the two figures are reproducible with public API only. STI-vs-T60 sweeps sti_from_impulse_response over synthetic exponential IRs at a 10-point T60 grid; the open-plan collapsible rebuilds the D2,S regression and rD/rP crossings from the OpenPlanResult fields. EN docs + EN/ES site; executed on Agg; llms-full.txt regenerated. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: final-review polish — result-field docs, comment hygiene, naming Document band_centers/measured/shifted_reference on WeightedRatingResult and ImpactRatingResult (CHANGELOG + EN/ES API tables), strip internal plan tags from theme-images.css and test_result_plots.py, note that [plot] also powers result .plot() methods, tighten the plot_zwicker_loudness docstring to say the two-panel layout only occurs when ax is None, sentence-case the intensity provenance comments, and rename the shadowed impact rating var to res_imp on the room-acoustics pages. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: PR 86 review round 1 — log-axis bars, style kwargs overrides, snippet physics and notation Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf

2 weeks ago

gum-uncertainty.md

feat: measurement uncertainty (GUM ISO/IEC Guide 98-3:2008 + Monte Carlo Supplement 1) (#102) Adds a `phonometry.uncertainty` module implementing the two propagation methods of the GUM: * Law of propagation of uncertainty (ISO/IEC Guide 98-3:2008, clause 5) — `combine_uncertainty` builds the combined standard uncertainty from central-difference sensitivity coefficients, with optional input correlations (validated for shape, symmetry, unit diagonal and positive-semidefiniteness), Welch-Satterthwaite effective degrees of freedom (Annex G.4) and the expanded uncertainty U = k*uc (clause 6). * Monte Carlo method (ISO/IEC Guide 98-3-1:2008, Supplement 1, clause 7) — `monte_carlo` propagates the input PDFs and returns the estimate, its standard uncertainty and the probabilistically symmetric coverage interval (clause 7.7). Type B quantities from a half-width: `rectangular` (a/sqrt3), `triangular` (a/sqrt6), `u_shaped` (a/sqrt2). Results expose `.expanded()` and `.plot()` (uncertainty budget). Validated against the Guides' worked examples: additive uc=2.0 (Suppl. 1, 9.2), coverage factor k=2.92 at p=0.99/v=16 (Table G.2), Welch-Satterthwaite v_eff=40 (Annex G.4), Monte Carlo interval [-3.88, 3.88] (Suppl. 1, 9.2.3). Conformance 79/79. Includes the budget/Monte-Carlo figure, the two-lane GUM-vs-Monte-Carlo diagram (EN/ES, light/dark), the repo doc and the EN/ES site guides under a new "Metrology & uncertainty" section.

2 weeks ago

hearing-threshold.md

feat: hearing threshold (ISO 7029:2017 age model, ISO 389-7:2006 reference) (#101) Add phonometry.hearing: age_threshold (ISO 7029:2017 statistical age distribution — median, spreads, population fractile) and reference_threshold (ISO 389-7:2006 free/diffuse-field audiometric zero) over 125 Hz - 8000 Hz. Coefficients parsed digit-exact from the standard. AgeThresholdResult with .plot(); 11 tests; 3 conformance checks (76/76); figure + diagram (EN/ES, light/dark); repo doc + site guides EN/ES.

2 weeks ago

human-vibration.md

Human vibration exposure (ISO 8041-1, ISO 2631, ISO 5349, Directive 2002/44/EC) (#94) New `human_vibration` module implementing the ISO 8041-1:2017 master weighting cascade (all nine weightings Wb/Wc/Wd/We/Wf/Wh/Wj/Wk/Wm), ISO 2631-1/-2/-4 whole-body metrics (running RMS, MTVV, VDV, MSDV, crest factor, vibration total value, energy-equivalent acceleration), the ISO 5349-1/-2 hand-arm A(8) arithmetic and vibration-white-finger latency, and the Directive 2002/44/EC EAV/ELV assessment. Every result object exposes .plot(); validated to <0.05% against the ISO 8041-1 Annex B tables and the ISO 5349-2 Annex E worked examples. Includes docs page, figures, setup diagram, and 7 conformance checks (100%).

2 weeks ago

impulse-prominence.md

NT ACOU 112: impulsive-sound prominence and LAeq adjustment (#105) * feat: NT ACOU 112 impulsive-sound prominence and LAeq adjustment Add the Nordtest NT ACOU 112:2002 method for the prominence of impulsive sounds: the predicted prominence P = 3 lg(OR) + 2 lg(LD) from each impulse's onset rate and level difference (clause 7, Formula 1), the graduated adjustment KI = 1.8 (P - 5) dB for P > 5 (clause 8, Formula 2), and the rating level over a reference time interval (clause 8, Note 1). New module phonometry.ntacou112 with predicted_prominence, impulse_adjustment, impulse_prominence (ImpulseProminence result with .plot()) and rating_level. Conformance report gains an impulsive-sound domain (2 checks). Adds the impulse-prominence guide (repo + EN/ES site) with figure and flow diagram. * refactor: address review feedback on NT ACOU 112 - impulse_prominence: ravel inputs so multi-dimensional arrays flatten to a 1-D impulse sequence (Gemini) - figure: use a distinct neutral for the threshold line instead of the grid color, matching plot_impulse_prominence (Gemini) - tests: use the shared NTACOU_PROMINENCE / NTACOU_ADJUSTMENT_P10 reference constants instead of hardcoded literals (CodeRabbit) * refactor: second-pass review polish for NT ACOU 112 - figure: give the impulse markers a distinct light blue (#aec7e8), matching plot_impulse_prominence, so they read apart from the neutral threshold line (CodeRabbit) - tests: split the chained magnitude/approx assertion into two clearer assertions (CodeRabbit)

2 weeks ago

fix: first pass of audit findings (docs, figures, plotting) Documentation: - fix the Spanish filter-figure regression: _showfilter gains a close flag and the generator saves through themed_path after drawing, so the translation pass runs on the finished figure; restore the two _ES_EXACT keys removed in #82 (they are emitted by filter_design, not by the script) - repair the broken building-acoustics.md link in the enclosed-space guide (the repo split lives in room-acoustics.md) - make the multiple-shock and intensity snippets self-contained: define az/fs with a synthetic shock train (verified outputs 20.94 m/s2, R=0.46, probability 0.03) and import numpy with stated outputs (F2=3.41, Ld=8.0), in the repo doc and both site languages - decimal points in the English Formula texts (1.07, 55.3) that had kept the Spanish comma Figures: - stop drawing data series in the grid colour (invisible on light backgrounds): time-weighting input burst, hearing-threshold 20 yr and NIPTS 10 yr curves, and the unweighted/partial bars now use a visible neutral grey; regenerate the 15 affected figure sets Plotting layer: - plot_facade_insulation forwards user kwargs to the primary D2m,nT curve only, so label=/color= no longer raise or merge the four curves; the kwargs-forwarding test now covers the facade plot - harden three edge cases: all-masked SII (zero contribution no longer divides 0/0), empty per-impulse prominence, and nearest-band lookup for the NC governing marker instead of float equality - sharpness normalization sanity guards raise RuntimeError instead of AssertionError

2 weeks ago

Outdoor sound propagation and occupational noise exposure (ISO 9613-1, ISO 9613-2, ISO 9612) (#89) * feat: ISO 9613-1 atmospheric sound absorption Full alpha(f, T, RH, p) model per ISO 9613-1:1993 clauses 6.2-6.4: oxygen/nitrogen relaxation frequencies, water-vapour molar concentration from relative humidity, classical + rotational + vibrational terms. Validated against a 13-point Table 1 grid (worst deviation 0.39 %, rounding-limited, vs the standard's stated +/-10 % accuracy) with the exact midband convention of Eq. (6) Note 5. air_attenuation_m() feeds the ISO 354 m parameter directly, closing the user-supplied gap. Two conformance checks added (34/34). Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * refactor: C1 review nits — test name and comment direction Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * test: extend ISO 9613-1 Table 1 oracle against the full standard text The full 37-page ISO 9613-1:1993 text (now in hand) confirms every implemented constant of Eqs. (3)-(6), the Annex B psychrometric conversion and the clause 4.2 reference conditions exactly as coded. Oracle extended from 13 to 23 Table 1 points, adding the 20-50 degC sub-tables 1(i)-1(p) unavailable in the interim source. Worst deviation unchanged at 0.385 % (rounding-limited). Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * feat: ISO 9613-2 outdoor sound propagation Implement the ISO 9613-2:1996 general method for octave-band outdoor sound attenuation: geometrical divergence (Eq. 7), atmospheric absorption (Eq. 8, via ISO 9613-1 air_absorption), ground effect (general per-region method of 7.3.1 with the Table 3 a'/b'/c'/d' functions, Eq. 9, plus the alternative method of 7.3.2, Eq. 10, and the DOmega index Eq. 11), and screening (Eq. 12-18 with the C2/C3 factors, single/double diffraction, Kmet and the 20/25 dB caps). Adds the Cmet meteorological correction (Eq. 21/22) and the LfT(DW)=Lw+Dc-A receiver-level composition (Eq. 3/6). Frozen dataclass result exposes the per-term breakdown (Adiv/Aatm/Agr/Abar). New module, tests and exports. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * feat: ISO 9612 occupational noise exposure strategies and uncertainty Add phonometry.occupational_exposure implementing the three ISO 9612:2009 measurement strategies and the normative Annex C uncertainty budget: - task_based_exposure (Clause 9): per-task energy average (Eq 7), task contributions (Eq 8), daily level (Eq 9/10); uncertainty Eq C.3 with sampling u1a (Eq C.6), duration u1b (Eq C.7) and sensitivity coefficients c1a (Eq C.4) / c1b (Eq C.5). - job_based_exposure (Clause 10): effective-day energy average (Eq 11), daily level (Eq 12); Table 1 minimum-duration check; uncertainty Eq C.9. - full_day_exposure (Clause 11): whole-day averaging (Eq 11/13) with the Clause 11.3 three-measurement / 3 dB spread rule. - Table C.4 (job/full-day c1*u1, bilinear interpolation), Table C.5 (u2), u3 = 1.0 dB, coverage factor k = 1.65 (one-sided 95 %). Frozen dataclasses (Task, TaskContribution, ExposureResult), ExposureWarning advisories for the sampling rules, and exports appended to __init__. Tests reproduce the Annex D (task), E (job) and F (full-day) worked examples to the standard's precision, plus closed-form oracles (single 85 dB/8 h task = 85; equal-level tasks; -3.01 dB per halved duration; k = 1.65) and the Table C.4 / Table 1 lookups. make check green: 1167 passed, 12 skipped. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: ISO 9612 review corrections — Table C.4 digits and Eq C.7 denominator Two Table C.4 cells verified against the page image at 200 dpi and corrected (N=6/u1=3.5: 3.4 -> 3.3; N=20/u1=2.0: 0.6 -> 0.5), the test that echoed the wrong N=6 cell now pins the printed value plus the corrected N=20 cell. The duration_samples branch of Eq C.7 divides by J(J-1) like Eq C.6 (standard error of the mean), with a regression test pinning u1b = 1.0 h for samples (4, 6) h. Module docstring now states the Lp,Cpeak uncertainty scope cut (Table C.5, NOTE 1). Claude-Session: https://claude.ai/code/session_01HUHMx4dxDes2t5UPQ25djP * test: ISO 9613-2 and ISO 9612 conformance checks Add a new "Outdoor propagation & occupational exposure" conformance domain with seven closed-form / worked-example checks backed by shared constants in tests/reference_data.py: - ISO 9613-2:1996 Eq. (7) geometrical divergence Adiv = 51 dB at 100 m. - ISO 9613-2:1996 Table 3 ground limit b'(0) = 10,1 -> Agr(250 Hz) = 17,2 dB (porous ground, hs = hr = 0, dp -> inf). - ISO 9613-2:1996 clause 7.4 single/double diffraction 20/25 dB caps. - ISO 9612:2009 Annex D/E/F daily exposure LEX,8h and expanded uncertainty U (84,3/2,7, 88,1/3,8, 90,1/3,4 dB), the Task objects rebuilt from the shared input table. Registry grows 34 -> 41 checks across 8 domains; docs/CONFORMANCE.md regenerated via make conformance. Registration and reference-data consistency tests pin the constants to the published values. Claude-Session: https://claude.ai/code/session_01HUHMx4dxDes2t5UPQ25djP * docs: figures for outdoor propagation and occupational exposure Add four documentation figures (each EN/ES x light/dark) plus a setup diagram: - air_absorption_alpha: ISO 9613-1 alpha(f) over frequency at four temperature/humidity combinations (relaxation roll-off, f^2 growth). - outdoor_attenuation_breakdown: ISO 9613-2 per-term stacked bar (Adiv, Aatm, Agr, Abar) with the total A, for a 200 m path over porous ground with a 4 m barrier. - exposure_uncertainty: ISO 9612 Annex D task contributions with the daily LEX,8h and the one-sided 95 % upper limit LEX,8h + U. - diagram_outdoor_geometry: source-barrier-receiver geometry (dss/dsr/z, the blocked direct ray and the diffracted ray over the top edge). Spanish translations added to the generate_graphs / generate_diagrams dictionaries (decimal commas, UNE terminology). Every image visually audited. Claude-Session: https://claude.ai/code/session_01HUHMx4dxDes2t5UPQ25djP * docs: ISO 9613-1/2 outdoor propagation guide and ISO 9612 exposure (EN) Add the didactic Outdoor Sound Propagation guide (docs/outdoor-propagation.md) covering ISO 9613-1 atmospheric absorption and the ISO 9613-2 general method (divergence, atmospheric, ground and barrier terms, the OutdoorAttenuation breakdown, the alternative 7.3.2 ground method, Cmet), with executed snippets and honest scope notes (Amisc/reflections not implemented, 7.3.2 not auto-wired, per-band Cmet a convenience). Extend the Levels guide with the ISO 9612 occupational-exposure section (task/job/full-day strategies, the C.6 I(I-1) vs C.12 N-1 asymmetry, Table C.4, k = 1,65, LEX,8h + U upper limit, Lp,Cpeak uncertainty out of scope), reproducing the Annex D/E/F worked examples. Add the matching theory sections, API-reference rows for every new public name, README feature list + docs table, docs index and CHANGELOG entries. llms-full.txt / llms.txt regenerated (make llms) with the new page wired into gen_llms.py. GitHub-safe math throughout; every snippet executed with pasted outputs. Claude-Session: https://claude.ai/code/session_01HUHMx4dxDes2t5UPQ25djP * docs: outdoor propagation and exposure site twins (EN + ES) Add the English and Spanish Starlight twins of the new content: - guides/outdoor-propagation.md (EN + ES) — code-identical snippets to the docs page; Spanish prose uses UNE terminology ("ponderación temporal"), decimal commas outside math and code, and the _es figure variants. - The ISO 9612 section added to guides/levels.md (EN + ES). - reference/theory.md and reference/api.md rows (EN + ES). - Landing cards and standards list bumped 34 -> 37 (index.mdx EN + ES); guides/outdoor-propagation added to the sidebar. Claude-Session: https://claude.ai/code/session_01HUHMx4dxDes2t5UPQ25djP * fix: final branch review — API rows, precision phrasing and ES links Adds the five missing public-symbol rows (TaskContribution, table_c4_contribution, minimum_cumulative_duration_hours, COVERAGE_FACTOR, INSTRUMENT_U2) to the three API references, states the ISO 9613-2 default atmosphere (20 °C / 70 %, a Table 2 reference atmosphere) in the API rows and guides, and softens "digit-for-digit" to "the standard's printed precision" wherever it covered Annex E, whose final level differs only by the standard's own pre-rounding. ES pages now use the /es/ link prefix consistently. Code nits from the review: dead arg > 1 guard removed (z > 0 implies arg >= 3), _with_advisory uses dataclasses.replace, and the float-keyed _PRIME_BY_BAND gets a safety comment. llms files regenerated. Claude-Session: https://claude.ai/code/session_01HUHMx4dxDes2t5UPQ25djP * fix: validate propagation/exposure inputs and docs lint (PR review) Addresses the PR #89 bot review: - Barrier.__post_init__ rejects negative edge distances and a zero edge_separation (which previously divided by zero in the C3 double- diffraction factor); Task.__post_init__ rejects empty samples and a non-positive duration. - ground_attenuation and barrier_attenuation reject non-positive frequencies (a zero frequency gave an infinite wavelength / division by zero); ground_attenuation also checks a positive distance and non-negative heights. job_based_exposure rejects a non-positive sample_duration_hours. - Rewrite the ISO 9613-2 7.3.2 closed form in the guides and API tables so Markdownlint no longer parses `(...)[...]` as a reversed link (MD011). - Regression tests for every new guard; docstrings note the new raises. Claude-Session: https://claude.ai/code/session_01HUHMx4dxDes2t5UPQ25djP

2 weeks ago

docs: correct impedance-tube inputs in materials reference (#96) transfer_matrix_two_load / transfer_matrix_one_load take the four microphone transfer functions (H1, H2, H3, H4) of each load (ASTM E2611), not "wave amplitudes". Aligns docs/materials.md with the API docstring and the documentation site.

2 weeks ago

multichannel.md

Audit-driven accuracy, robustness and didactic improvements across the library (#82) * fix: audit wave A — inter-sample LCpeak, cheby2 class-1 default, 144 kHz weighting target Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: audit wave A — decay validity thresholds, variant sharpness anchors, N5 phase, STIPA warning, Farina guard Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: audit wave B — self-contained snippets, tables, sharpness/open-plan figures, orphans Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: audit wave C — remaining minors and optimizations across the library Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: document the audit-wave parameters across guides and API tables Cover the public parameters added in the audit fix waves across all three documentation trees (docs/, site EN, site ES): - lc_peak(..., oversample=8): inter-sample peak recovery (levels guide + API) - calculate_sensitivity(..., narrowband=False): coherent tone estimator (calibration guide + API) - sound_intensity(..., bias_correct=False): finite-difference bias correction (intensity guide + API) - room_parameters/decay_curve(..., zero_phase=False): 125 Hz short-T bias (room-acoustics guide + API) - OctaveFilterBank/octavefilter attenuation default 60 -> 72 and cheby2 class-1 note (filter-banks guide + API) - STIPA < 15 s UserWarning (psychoacoustics guide + API) Also refresh dependent behaviour notes: ln_levels attack skip 2*tau -> 5*tau, T20/T30 validity thresholds 35/45 -> 46/54 dB, and the zero-phase broadband ~0.2-0.3 dB band-level caveat. Regenerated llms-full.txt. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: truthfulness fixes from the final audit-branch review Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: 2000 Hz, not 2 ms, in the N5/N10 docstring Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: PR 82 review round 1 — bias-correction cutoff, parabolic-peak guard, complexity, snippet self-containment Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf

2 weeks ago

multiple-shock-vibration.md

fix: first pass of audit findings (docs, figures, plotting) Documentation: - fix the Spanish filter-figure regression: _showfilter gains a close flag and the generator saves through themed_path after drawing, so the translation pass runs on the finished figure; restore the two _ES_EXACT keys removed in #82 (they are emitted by filter_design, not by the script) - repair the broken building-acoustics.md link in the enclosed-space guide (the repo split lives in room-acoustics.md) - make the multiple-shock and intensity snippets self-contained: define az/fs with a synthetic shock train (verified outputs 20.94 m/s2, R=0.46, probability 0.03) and import numpy with stated outputs (F2=3.41, Ld=8.0), in the repo doc and both site languages - decimal points in the English Formula texts (1.07, 55.3) that had kept the Spanish comma Figures: - stop drawing data series in the grid colour (invisible on light backgrounds): time-weighting input burst, hearing-threshold 20 yr and NIPTS 10 yr curves, and the unweighted/partial bars now use a visible neutral grey; regenerate the 15 affected figure sets Plotting layer: - plot_facade_insulation forwards user kwargs to the primary D2m,nT curve only, so label=/color= no longer raise or merge the four curves; the kwargs-forwarding test now covers the facade plot - harden three edge cases: all-masked SII (zero contribution no longer divides 0/0), empty per-impulse prominence, and nearest-band lookup for the NC governing marker instead of float equality - sharpness normalization sanity guards raise RuntimeError instead of AssertionError

2 weeks ago

noise-induced-hearing-loss.md

feat: noise-induced hearing loss (ISO 1999:2013 NIPTS and HTLAN) (#103) Adds a `phonometry.iso1999` module estimating occupational noise-induced hearing loss, extending the ISO 7029 age-related threshold with a noise component. * `nipts(l_ex, years, fractile)` — noise-induced permanent threshold shift over 500-6000 Hz: median N50 (clause 6.3.1, Formula 2/3, Table 1) and population fractiles from the half-Gaussian spreads du/dl (clause 6.3.2, Formulae 4-7, Tables 2-3), clamped at zero. * `htlan(age, sex, l_ex, years, fractile)` — hearing threshold level associated with age and noise, combining the age component (HTLA = ISO 7029) with the NIPTS by H' = H + N - H*N/120 (clause 6.1, Formula 1). Validated against ISO 1999 Annex D (Tables D.1-D.4, 85-100 dB, 10-40 years): reproduces every tabulated dB value. Conformance 82/82; full suite 1686 passed. Figure + two-lane NIHL diagram (EN/ES, light/dark), repo doc and EN/ES site guides.

2 weeks ago

outdoor-propagation.md

Outdoor sound propagation and occupational noise exposure (ISO 9613-1, ISO 9613-2, ISO 9612) (#89) * feat: ISO 9613-1 atmospheric sound absorption Full alpha(f, T, RH, p) model per ISO 9613-1:1993 clauses 6.2-6.4: oxygen/nitrogen relaxation frequencies, water-vapour molar concentration from relative humidity, classical + rotational + vibrational terms. Validated against a 13-point Table 1 grid (worst deviation 0.39 %, rounding-limited, vs the standard's stated +/-10 % accuracy) with the exact midband convention of Eq. (6) Note 5. air_attenuation_m() feeds the ISO 354 m parameter directly, closing the user-supplied gap. Two conformance checks added (34/34). Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * refactor: C1 review nits — test name and comment direction Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * test: extend ISO 9613-1 Table 1 oracle against the full standard text The full 37-page ISO 9613-1:1993 text (now in hand) confirms every implemented constant of Eqs. (3)-(6), the Annex B psychrometric conversion and the clause 4.2 reference conditions exactly as coded. Oracle extended from 13 to 23 Table 1 points, adding the 20-50 degC sub-tables 1(i)-1(p) unavailable in the interim source. Worst deviation unchanged at 0.385 % (rounding-limited). Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * feat: ISO 9613-2 outdoor sound propagation Implement the ISO 9613-2:1996 general method for octave-band outdoor sound attenuation: geometrical divergence (Eq. 7), atmospheric absorption (Eq. 8, via ISO 9613-1 air_absorption), ground effect (general per-region method of 7.3.1 with the Table 3 a'/b'/c'/d' functions, Eq. 9, plus the alternative method of 7.3.2, Eq. 10, and the DOmega index Eq. 11), and screening (Eq. 12-18 with the C2/C3 factors, single/double diffraction, Kmet and the 20/25 dB caps). Adds the Cmet meteorological correction (Eq. 21/22) and the LfT(DW)=Lw+Dc-A receiver-level composition (Eq. 3/6). Frozen dataclass result exposes the per-term breakdown (Adiv/Aatm/Agr/Abar). New module, tests and exports. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * feat: ISO 9612 occupational noise exposure strategies and uncertainty Add phonometry.occupational_exposure implementing the three ISO 9612:2009 measurement strategies and the normative Annex C uncertainty budget: - task_based_exposure (Clause 9): per-task energy average (Eq 7), task contributions (Eq 8), daily level (Eq 9/10); uncertainty Eq C.3 with sampling u1a (Eq C.6), duration u1b (Eq C.7) and sensitivity coefficients c1a (Eq C.4) / c1b (Eq C.5). - job_based_exposure (Clause 10): effective-day energy average (Eq 11), daily level (Eq 12); Table 1 minimum-duration check; uncertainty Eq C.9. - full_day_exposure (Clause 11): whole-day averaging (Eq 11/13) with the Clause 11.3 three-measurement / 3 dB spread rule. - Table C.4 (job/full-day c1*u1, bilinear interpolation), Table C.5 (u2), u3 = 1.0 dB, coverage factor k = 1.65 (one-sided 95 %). Frozen dataclasses (Task, TaskContribution, ExposureResult), ExposureWarning advisories for the sampling rules, and exports appended to __init__. Tests reproduce the Annex D (task), E (job) and F (full-day) worked examples to the standard's precision, plus closed-form oracles (single 85 dB/8 h task = 85; equal-level tasks; -3.01 dB per halved duration; k = 1.65) and the Table C.4 / Table 1 lookups. make check green: 1167 passed, 12 skipped. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: ISO 9612 review corrections — Table C.4 digits and Eq C.7 denominator Two Table C.4 cells verified against the page image at 200 dpi and corrected (N=6/u1=3.5: 3.4 -> 3.3; N=20/u1=2.0: 0.6 -> 0.5), the test that echoed the wrong N=6 cell now pins the printed value plus the corrected N=20 cell. The duration_samples branch of Eq C.7 divides by J(J-1) like Eq C.6 (standard error of the mean), with a regression test pinning u1b = 1.0 h for samples (4, 6) h. Module docstring now states the Lp,Cpeak uncertainty scope cut (Table C.5, NOTE 1). Claude-Session: https://claude.ai/code/session_01HUHMx4dxDes2t5UPQ25djP * test: ISO 9613-2 and ISO 9612 conformance checks Add a new "Outdoor propagation & occupational exposure" conformance domain with seven closed-form / worked-example checks backed by shared constants in tests/reference_data.py: - ISO 9613-2:1996 Eq. (7) geometrical divergence Adiv = 51 dB at 100 m. - ISO 9613-2:1996 Table 3 ground limit b'(0) = 10,1 -> Agr(250 Hz) = 17,2 dB (porous ground, hs = hr = 0, dp -> inf). - ISO 9613-2:1996 clause 7.4 single/double diffraction 20/25 dB caps. - ISO 9612:2009 Annex D/E/F daily exposure LEX,8h and expanded uncertainty U (84,3/2,7, 88,1/3,8, 90,1/3,4 dB), the Task objects rebuilt from the shared input table. Registry grows 34 -> 41 checks across 8 domains; docs/CONFORMANCE.md regenerated via make conformance. Registration and reference-data consistency tests pin the constants to the published values. Claude-Session: https://claude.ai/code/session_01HUHMx4dxDes2t5UPQ25djP * docs: figures for outdoor propagation and occupational exposure Add four documentation figures (each EN/ES x light/dark) plus a setup diagram: - air_absorption_alpha: ISO 9613-1 alpha(f) over frequency at four temperature/humidity combinations (relaxation roll-off, f^2 growth). - outdoor_attenuation_breakdown: ISO 9613-2 per-term stacked bar (Adiv, Aatm, Agr, Abar) with the total A, for a 200 m path over porous ground with a 4 m barrier. - exposure_uncertainty: ISO 9612 Annex D task contributions with the daily LEX,8h and the one-sided 95 % upper limit LEX,8h + U. - diagram_outdoor_geometry: source-barrier-receiver geometry (dss/dsr/z, the blocked direct ray and the diffracted ray over the top edge). Spanish translations added to the generate_graphs / generate_diagrams dictionaries (decimal commas, UNE terminology). Every image visually audited. Claude-Session: https://claude.ai/code/session_01HUHMx4dxDes2t5UPQ25djP * docs: ISO 9613-1/2 outdoor propagation guide and ISO 9612 exposure (EN) Add the didactic Outdoor Sound Propagation guide (docs/outdoor-propagation.md) covering ISO 9613-1 atmospheric absorption and the ISO 9613-2 general method (divergence, atmospheric, ground and barrier terms, the OutdoorAttenuation breakdown, the alternative 7.3.2 ground method, Cmet), with executed snippets and honest scope notes (Amisc/reflections not implemented, 7.3.2 not auto-wired, per-band Cmet a convenience). Extend the Levels guide with the ISO 9612 occupational-exposure section (task/job/full-day strategies, the C.6 I(I-1) vs C.12 N-1 asymmetry, Table C.4, k = 1,65, LEX,8h + U upper limit, Lp,Cpeak uncertainty out of scope), reproducing the Annex D/E/F worked examples. Add the matching theory sections, API-reference rows for every new public name, README feature list + docs table, docs index and CHANGELOG entries. llms-full.txt / llms.txt regenerated (make llms) with the new page wired into gen_llms.py. GitHub-safe math throughout; every snippet executed with pasted outputs. Claude-Session: https://claude.ai/code/session_01HUHMx4dxDes2t5UPQ25djP * docs: outdoor propagation and exposure site twins (EN + ES) Add the English and Spanish Starlight twins of the new content: - guides/outdoor-propagation.md (EN + ES) — code-identical snippets to the docs page; Spanish prose uses UNE terminology ("ponderación temporal"), decimal commas outside math and code, and the _es figure variants. - The ISO 9612 section added to guides/levels.md (EN + ES). - reference/theory.md and reference/api.md rows (EN + ES). - Landing cards and standards list bumped 34 -> 37 (index.mdx EN + ES); guides/outdoor-propagation added to the sidebar. Claude-Session: https://claude.ai/code/session_01HUHMx4dxDes2t5UPQ25djP * fix: final branch review — API rows, precision phrasing and ES links Adds the five missing public-symbol rows (TaskContribution, table_c4_contribution, minimum_cumulative_duration_hours, COVERAGE_FACTOR, INSTRUMENT_U2) to the three API references, states the ISO 9613-2 default atmosphere (20 °C / 70 %, a Table 2 reference atmosphere) in the API rows and guides, and softens "digit-for-digit" to "the standard's printed precision" wherever it covered Annex E, whose final level differs only by the standard's own pre-rounding. ES pages now use the /es/ link prefix consistently. Code nits from the review: dead arg > 1 guard removed (z > 0 implies arg >= 3), _with_advisory uses dataclasses.replace, and the float-keyed _PRIME_BY_BAND gets a safety comment. llms files regenerated. Claude-Session: https://claude.ai/code/session_01HUHMx4dxDes2t5UPQ25djP * fix: validate propagation/exposure inputs and docs lint (PR review) Addresses the PR #89 bot review: - Barrier.__post_init__ rejects negative edge distances and a zero edge_separation (which previously divided by zero in the C3 double- diffraction factor); Task.__post_init__ rejects empty samples and a non-positive duration. - ground_attenuation and barrier_attenuation reject non-positive frequencies (a zero frequency gave an infinite wavelength / division by zero); ground_attenuation also checks a positive distance and non-negative heights. job_based_exposure rejects a non-positive sample_duration_hours. - Rewrite the ISO 9613-2 7.3.2 closed form in the guides and API tables so Markdownlint no longer parses `(...)[...]` as a reversed link (MD011). - Regression tests for every new guard; docstrings note the new raises. Claude-Session: https://claude.ai/code/session_01HUHMx4dxDes2t5UPQ25djP

2 weeks ago

psychoacoustics.md

Advanced psychoacoustics (ISO 532-2/3 Moore-Glasberg, ECMA-418-2 Sottek loudness/tonality/roughness) + numerical conformance report (#87) * feat: ECMA-418-2 Sottek Hearing Model loudness and auditory front-end Implement the ECMA-418-2:2025 (4th ed.) psychoacoustic loudness metric and its shared auditory front-end (Sottek Hearing Model): - Clause 5 front-end: outer/middle-ear cascade (5.1.3, Table 1), 53-band gammatone-like auditory filter bank (5.1.4), band-dependent segmentation (5.1.5, Table 4), rectification/RMS (5.1.6-5.1.7), compressive nonlinearity (5.1.8, Formula 23, Table 2) and threshold in quiet (5.1.9, Table 3). - Clause 6.2.2-6.2.7 ACF-based tonal/noise specific loudness (reusable by the later tonality/roughness metrics). - Clause 8 assembly: tonal/noise power average, average specific loudness, time-dependent and single representative loudness value. The front-end helpers are factored for reuse. Calibration verified: a 1 kHz 40 dB SPL sinusoid yields 0.996 sone_HMS (target 1.0). Adds a frozen EcmaLoudness result with lazy .plot(), exports, and 12 tests. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * feat: ECMA-418-2 tonality (Sottek Hearing Model) Implements the psychoacoustic tonality metric of ECMA-418-2:2025 (Clause 6.2.8-6.2.11) on top of the committed Sottek front-end and ACF tonal/noise decomposition. - tonality_ecma() -> EcmaTonality with the single value T (Formula 63), average specific tonality T'(z) (Formula 53), tonal frequencies f_ton,z(z) (Formula 55), time-dependent tonality T(l) (Formula 61) and its frequency f_ton(l) (Formula 62); optional user band [f_low,f_high]. - Full Clause 6.2.3 band averaging with cross-block-size-group ACF recomputation (loudness keeps its simplified, valid-for-loudness path; boundary marker added there). Loudness 0.996 sone guard unchanged. - _tonal_estimate now also returns the DFT-peak tonal frequency, with the argmax restricted to the lower DFT half to avoid Nyquist-mirror aliasing. - .plot() via _plotting.plot_ecma_tonality; exports in __init__. - 13 conformance tests: 1 kHz/40 dB -> 1 tu_HMS calibration, pure tone >> noise, f_ton tracks the tone, silence -> 0, band restriction. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * feat: ECMA-418-2 roughness (Sottek Hearing Model) Implement the ECMA-418-2:2025 (4th ed.) psychoacoustic roughness metric (Clause 7), the third Sottek Hearing Model output and new library functionality. Reuses the Clause 5 auditory front-end from loudness_ecma (ear filter, 53-band gammatone bank, specific basis loudness) and adds the roughness-specific chain: start-only zero-padding and fixed 16384/4096 segmentation (5.1.2.2, 5.1.5.2, 7.1.1), the Hilbert envelope with x32 downsampling to 1500 Hz (7.1.2), the scaled envelope power spectrum (7.1.3), two-step noise reduction (7.1.4), the four-stage spectral weighting with quadratic-fit modulation-rate refinement, bias correction, high/low modulation-rate weighting and fundamental-rate estimation (7.1.5), interpolation to 50 Hz with the distribution-dependent nonlinear transform, tabulated calibration c_R and asymmetric smoothing (7.1.7), and the R'(z) / R(l50) / 90th-percentile aggregation (7.1.8-7.1.10). Uses the standard's tabulated c_R = 0.0180685 (Formula 104), not reverse-fit: the 1 kHz / 70 Hz / m=1 / 60 dB calibration signal yields R = 1.0735 asper (near-but-not-exactly 1). OCR-ambiguous formulae (66, 71, 78-81) cross-checked against the SQAT MATLAB reference; the optional entropy weighting (7.1.6, needs RPM telemetry) is documented but not implemented. Adds roughness_ecma()/EcmaRoughness, a two-panel .plot(), exports, a 10-test TDD suite and plan/notes-ecma418-2-roughness.md. Full suite 752 passed / 12 skipped; loudness (0.996) and tonality (0.99993) oracles unchanged (both modules untouched). Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * test: pin ECMA-418-2 roughness calibration value and document clean-room variance Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * feat: ISO 532-2 Moore-Glasberg loudness Add the stationary Moore-Glasberg loudness model of ISO 532-2:2017 as a third loudness method alongside the Zwicker (ISO 532-1) and Sottek (ECMA-418-2) implementations. The full Clause 7 chain is implemented: fixed outer/middle-ear transfer (Table 1), the level-dependent roex auditory filter bank and excitation pattern on the ERB-number scale (Formulae 1-6), the compressive specific-loudness transform (Formulae 7-9, Tables 2-4, C=0.0617 sone/Cam), binaural inhibition (Formulae 10-13) and integration to total loudness with the phon mapping of Table 5. Public API: loudness_moore_glasberg_from_spectrum (exact 5.2/5.4), loudness_moore_glasberg_from_third_octave (29 bands, 5.5) and loudness_moore_glasberg (signal wrapper), returning a frozen MooreGlasbergLoudness with .plot() over the ERB-number scale. Validated against the Annex B reference signals (tones, white/pink noise, multi-tone complexes, tone-in-noise) to ~1-4%, well inside the standard's 2.8 phon expanded uncertainty; the definitional anchor (1 kHz/40 dB/free/binaural = 1.000 sone) holds end-to-end from the tabulated constant. 50 new tests; full suite 802 passed / 12 skipped. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * feat: ISO 532-3 Moore-Glasberg-Schlittenlacher time-varying loudness Add loudness_moore_glasberg_time() implementing the ISO 532-3:2023 method: a six-FFT multi-resolution running short-term spectrum (clause 7.3) feeding the ISO 532-2 excitation and specific-loudness model on the 0.25 Cam grid with the 532-3 constants (C=0.063, Tables 2/4, E_THRQ/E0=2.307), integrated by asymmetric attack/release smoothing into short-term (clauses 7.6-7.8) and long-term (clause 7.9) loudness traces, with peak long-term loudness, percentiles and a lazy STL/LTL .plot(). Validated against the Annex C.1 tone table (1 kHz 10-80 dB, 3 kHz, 4 kHz, 100 Hz and monaural earphone) to within 0.4 phon, well inside the 2.8 phon expanded uncertainty; the 1 kHz/40 dB anchor is exact at 1.000 sone. The ISO 532-2 and ECMA-418-2 oracles are unchanged (831 passed / 12 skipped). Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: ISO 532-3 per-ear LTL and Formula-14 kernel width (dichotic path) Two transcription errors in the ISO 532-3:2023 dichotic (per-ear) binaural path of the time-varying Moore-Glasberg-Schlittenlacher module: 1. Clause 7.9: the long-term loudness must be computed per ear by the attack/release averager (Formulae 18-21) and then summed. The code ran a single AGC on the binaural sum S'_L + S'_R; because the attack/release branch is nonlinear, AGC(S'_L) + AGC(S'_R) != AGC(S'_L + S'_R) for dichotic input. Keep two LTL states and sum the two long-term loudnesses. For diotic/monaural the averager is linear in the (scaled) input so the result is bit-for-bit identical to the old sum-then-AGC. 2. Formula (14)/(15): the smoothing weight is exp(-(B*Di)^2) with B = 0.08 and Di in Cam. The code evaluated exp(-((0.08 * Di/0.1)^2)) = exp(-(0.8*Di)^2), a Gaussian ~10x too narrow, and truncated the taps at +/-1.75 Cam. Remove the spurious /0.1 and widen the kernel to the standard's +/-18 Cam (145 taps on the 0.25-Cam grid). For diotic/monaural the L/R excitation ratio is 1 regardless of the kernel, so all oracles are unaffected. No current validation signal exercises the dichotic path (Annex C signals are diotic/monaural), so the 532-3 tone oracle and every other oracle remain byte-identical; added focused dichotic tests plus a byte-identical anchor pin. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * feat: numerical conformance report in CI PR comment Reframe the PR comment into a persistent numerical conformance report: each implemented standard -> normative expected value/range -> computed result -> pass/fail with deviation, grouped by domain. - scripts/conformance_report.py: registry of 21 checks across 17 standards and 6 domains (filters/weightings, levels/dosimetry, psychoacoustics, speech, intensity/power, room/building). Emits a headline summary, a "Numerical validation - filters & weightings" section (per-architecture IEC 61260-1 class margins + A/C/G weighting deviation vs the normative curve) and one conformance table per domain. Deterministic, ~2.7 s. - tests/reference_data.py: single source of truth for the shared normative tables (IEC 61672-1 Table 3, ISO 7196 Table 2, ISO 717-1 Annex C R), imported by both the tests and the report so they cannot drift. - tests/test_conformance_report.py: smoke test asserting every registered check passes and the Markdown is well-formed. - comment_pr.py / workflow: conformance report leads the comment; the test/coverage table is collapsed below. Retire benchmark_filters.py (numerical content absorbed into the report). make check green (740 passed, 12 skipped); ruff, mypy --strict, bandit clean. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * refactor: block-A review minors — docstring notes and tonality band-range Four Minor findings from the plan-17 block-A psychoacoustics reviews. No change to any validated numeric output; all five oracles stay byte-identical. 1. ISO 532-2 _source_levels: docstring note that clause 7.4 specifies a roex weighting for X where the code uses a rectangular +/-ERB_n/2 window — identical for isolated tones, immaterial for the Annex B broadband cases. 2. ISO 532-2 Formula-10 binaural smoothing kernel: the PDF (clause 8.1) states "Di is changed in steps of 0,1" over -18..+18 Cam; the code confused the 0.1-Cam index offset with Di-in-Cam (the same /0.1 width bug just fixed in 532-3 Formula-14), making the Gaussian 10x too narrow and truncating to +/-1.8 Cam. Widen to +/-18 Cam (361 taps) with Di = tap*_I_STEP. Provably inert: all Annex B cases are diotic (ratio 1) or monaural (unit inhibition). 3. ECMA tonality _band_range: implement the exact Formulae 56-57 edge-midpoint boundary conditions instead of a centre-frequency threshold; only affects the optional user-band feature. Added a user-band edge test. 4. ECMA tonality reach: comment that the _CBF-1-band upper-edge term is a dead defensive guard (N_B>0 only for low bands 0..24, so it never limits reach). Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * chore: move conformance design doc out of repo root * feat: clearer conformance report (value-vs-range) + block-A checks + audit fixes Table 1 (filters): show, per architecture, the measured relative attenuation at the binding band vs the class-1 limit it must clear (value + range, not just a margin). Replace the "none" verdict (read as a failure) with a "By design (ripple/soft rolloff)" label + footnote; Butterworth (default) and Chebyshev-II are the only mask-compliant architectures. The measured value and limit are re-derived with the public class_limits on the same designed SOS; a smoke-test guard asserts the re-derived class-1 margin matches verify_filter_class (single source of truth). Table 2 (weightings): separate the informational max-deviation-from-nominal (at a frequency extreme with wide, asymmetric tolerance) from the compliance margin at the binding frequency, where the deviation, the +/- tolerance band and the headroom are co-located so "value vs range" is unambiguous and in-spec. Audit: verified all 21 checks (expected value/range vs the standard, real library call, honest tolerance). One fix - ISO 226 anchored at 60 phon @ 2 kHz -> 60 dB, a genuine Table B.1 value but coincidentally equal to the phon number (reads as the trivial 1 kHz identity); re-anchored to the non-coincidental 60 phon @ 100 Hz -> 78.5 dB, sourced from reference_data. No library bug found. Block-A psychoacoustics (+5, single-source expected in reference_data): ECMA-418-2 loudness (1 kHz/40 dB -> 1 sone_HMS, c_N), tonality (-> 1 tu_HMS, c_T), roughness (1 kHz/70 Hz/m=1/60 dB -> clean-room 1.0735 asper with a note that the standard target is 1.0); ISO 532-2 (1 kHz/40 dB -> 1 sone, C=0.0617); ISO 532-3 (steady 1 kHz/40 dB -> peak LTL 1 sone). 26/26 checks, 20 standards. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * feat: versioned docs/CONFORMANCE.md with make target, CI staleness check and README link Persist the numerical conformance report as a committed, versioned file so every reader sees the library's per-standard validation (expected vs computed) without running CI. - docs/CONFORMANCE.md: committed report (26/26 checks, 6 domains, 20 standards) with an auto-generated "do not hand-edit" header linking the standards docs. - scripts/conformance_report.py: add a --file-header flag that emits that header for the committed file (PR-comment body stays header-free); coarsen the informational delta column to 3 decimals so BLAS/FFT-dependent sub-milli residuals do not churn the file (status/computed columns keep regressions visible). - Makefile: `make conformance` regenerates the file; `make install-hooks` installs the optional pre-commit hook. - CI: new fast `conformance` job (runtime deps only) fails a PR if the committed file is stale, telling the author to run `make conformance` and commit. - hooks/pre-commit: opt-in convenience hook, scoped to src/scripts/reference changes; CI is the enforcement. - README + CONTRIBUTING: link and document the generated report. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: advanced psychoacoustics guide, theory, API (EN) Document the plan-17 block-A advanced psychoacoustics models: Moore-Glasberg loudness (ISO 532-2/532-3), Sottek Hearing Model loudness, tonality and roughness (ECMA-418-2:2025). - Psychoacoustics guide: new "Advanced loudness & sound-quality models" section with a model-comparison table, executed snippets, parameter tables, res.plot() one-liners and figure-code collapsibles (docs + site twin, code blocks byte-identical). - Theory: excitation-pattern/specific-loudness chain, Sottek front-end, ACF tonality and envelope-modulation roughness (GitHub-safe math). - API reference: rows for the 5 functions (+ 2 MG spectral entry points) and 5 result dataclasses with their .plot(). - README, EN landing (standards 26 -> 29: ISO 532-2, ISO 532-3, ECMA-418-2), CHANGELOG Unreleased/Added, llms regen. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: advanced-psychoacoustics figures (loudness models, Sottek, tonality/roughness, time-varying) Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: gemelas ES de psicoacústica avanzada (Moore-Glasberg, Sottek, tonalidad, aspereza) Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: ISO 532-2 signal-input path builds a narrowband spectrum (pure-tone loudness) The loudness_moore_glasberg(x, fs) wrapper routed the signal through the one-third-octave method, which smeared a pure tone across a whole 1/3-octave band and used a coherent-gain window normalisation that over-counted power by 1.76 dB. A calibrated 1 kHz tone at 40 dB SPL therefore returned 1.327 sone instead of the definitional anchor of 1.000 sone. ISO 532-2 is spectrum-based and its exact input (clauses 5.2/5.4) is a set of discrete sinusoidal components; the excitation pattern (Formula 5) is computed per component, so a tone must enter as a single line. The wrapper now forms the signal's narrowband (FFT) line spectrum with a power-preserving (Parseval) window normalisation and feeds it to loudness_moore_glasberg_from_spectrum. A 1 kHz/40 dB tone now yields 1.0001 sone / 40 phon and multi-tone signals match the exact spectrum path. The validated spectrum and third-octave paths are byte-identical. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: correct MG signal-path mechanism text and presentation default Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: PR 87 review round 1 — float-equality guards, 532-3 FFT window length, tonality short-signal, report/figure fixes - roughness_ecma: replace three float == 0.0 checks with abs>eps guards (SonarCloud gate); roughness oracle byte-identical (1.0735). - loudness_moore_glasberg_time: per-window FFT length (>= window length) so the 64 ms window is not truncated at 44.1/48 kHz (was ~0.74 dB low in the 20-80 Hz band). 32 kHz anchor byte-identical (n_max 1.0000044713237626). - tonality_ecma: short signals fall back to averaging all blocks (matches loudness_ecma), not the final block only; validated anchor unaffected. - conformance_report: document deliberate 273.0 (speed of sound) vs 273.15 (C1/C2) split matching the library; complete bind_side docstring. - generate_graphs: roughness annotation uses computed peak fmod; Sottek metrics use tu_HMS/sone_HMS units; ES translations + 8 figures regenerated. - hooks/pre-commit: drop set -e so a benign regen failure does not block the commit (CI staleness check is the gate). - test_iec_weighting_table3 docstring, Makefile .PHONY completion. All 5 psychoacoustics oracles byte-identical; make check 868 passed / 12 skipped; ruff + mypy + bandit clean; conformance 26/26. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: localize decimal commas in Sottek-unit figure labels (underscore guard) The `_` guard in `_translate_figure` skipped decimal-comma localization for any label containing an underscore. The new Sottek/HMS unit tokens (`sone_HMS`, `tu_HMS`, `sonios_HMS`) carry underscores, so Spanish figures regressed to dots (`8.0`, `0.42`) instead of commas (`8,0`, `0,42`). Drop the underscore guard: the decimal-comma regex only rewrites a bare `digit.digit` not adjacent to more digits/dots, so unit identifiers and version numbers (5.3.3) stay intact while genuine decimals get commas. Mathtext is still skipped via the `$` guard. Regenerated the affected ES figures (sottek_specific_loudness, tonality_roughness_demo). Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf

2 weeks ago

room-acoustics.md

feat: plottable impulse-response results and ISO 3382 measurement setup (#98) impulse_response / mls_impulse_response now return an ImpulseResponseResult with .plot() (waveform + log-magnitude envelope + Schroeder decay), staying a drop-in for the raw IR array via __array__ so existing callers are unchanged. Adds plot_excitation() for the sweep and MLS excitation signals, and an ISO 3382 measurement-setup diagram (source/microphone positions and ISO 3382-1 placement rules, plus the ISO 3382-2 minimum-position grades). Room-acoustics guide §1 (repo + EN/ES site) shows the new signal/IR plots and the diagram.

2 weeks ago

fix: SonarCloud smells + CodeRabbit review on room-noise (#100) Hoist the repeated "upper right" legend location into a Final constant (S1192); build the RC spectral tag with explicit concatenation (S5797); show the interpolated NC rating with :g so the figure title and prose agree (NC-42.5), regenerating the figure and alt text.

2 weeks ago

Scattering/diffusion, in-situ road absorption & precision sound power (ISO 17497-1/-2, ISO 13472-1/-2, ISO 3745, ISO 9614-3) (#92) Clean-room ISO 17497-1/-2 (scattering & diffusion), ISO 13472-1/-2 (in-situ road absorption) and ISO 3745 / 9614-3 (precision sound power). Plottable result objects with .plot(), single-concept documentation figures with dual-snippet <details>, 12 new conformance checks (60/60), and 6 experimental-setup diagrams. All subagent-reviewed and visually validated. https://claude.ai/code/session_01JJLhQzF5hTdnEi3bfmqnqm

2 weeks ago

speech-intelligibility.md

feat: SII standard speech spectra for raised, loud and shout vocal efforts (#104) Completes the ANSI S3.5-1997 Table 3 standard speech spectra in the SII module (#97), which shipped with only the normal-effort spectrum. Adds raised, loud and shout so `standard_speech_spectrum` and `speech_intelligibility_index` accept all four vocal-effort names. The three 18-band Table 3 arrays were cross-verified digit-exact against independent reference implementations (Google speech_intelligibility_index, R CRAN SII) and their reconstructed overall levels match the known ANSI vocal-effort overalls (normal 62.35, raised 68.3, loud 74.86, shout 82.36 dB SPL), strictly increasing. Speaking louder raises the index in a fixed noise (0.12 -> 0.79). Conformance 83/83; full suite 1691 passed. Adds a vocal-effort figure (EN/ES, light/dark) and a "vocal effort" section in the repo and site SII guides.

2 weeks ago

surface-scattering.md

Scattering/diffusion, in-situ road absorption & precision sound power (ISO 17497-1/-2, ISO 13472-1/-2, ISO 3745, ISO 9614-3) (#92) Clean-room ISO 17497-1/-2 (scattering & diffusion), ISO 13472-1/-2 (in-situ road absorption) and ISO 3745 / 9614-3 (precision sound power). Plottable result objects with .plot(), single-concept documentation figures with dual-snippet <details>, 12 new conformance checks (60/60), and 6 experimental-setup diagrams. All subagent-reviewed and visually validated. https://claude.ai/code/session_01JJLhQzF5hTdnEi3bfmqnqm

2 weeks ago

Outdoor sound propagation and occupational noise exposure (ISO 9613-1, ISO 9613-2, ISO 9612) (#89) * feat: ISO 9613-1 atmospheric sound absorption Full alpha(f, T, RH, p) model per ISO 9613-1:1993 clauses 6.2-6.4: oxygen/nitrogen relaxation frequencies, water-vapour molar concentration from relative humidity, classical + rotational + vibrational terms. Validated against a 13-point Table 1 grid (worst deviation 0.39 %, rounding-limited, vs the standard's stated +/-10 % accuracy) with the exact midband convention of Eq. (6) Note 5. air_attenuation_m() feeds the ISO 354 m parameter directly, closing the user-supplied gap. Two conformance checks added (34/34). Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * refactor: C1 review nits — test name and comment direction Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * test: extend ISO 9613-1 Table 1 oracle against the full standard text The full 37-page ISO 9613-1:1993 text (now in hand) confirms every implemented constant of Eqs. (3)-(6), the Annex B psychrometric conversion and the clause 4.2 reference conditions exactly as coded. Oracle extended from 13 to 23 Table 1 points, adding the 20-50 degC sub-tables 1(i)-1(p) unavailable in the interim source. Worst deviation unchanged at 0.385 % (rounding-limited). Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * feat: ISO 9613-2 outdoor sound propagation Implement the ISO 9613-2:1996 general method for octave-band outdoor sound attenuation: geometrical divergence (Eq. 7), atmospheric absorption (Eq. 8, via ISO 9613-1 air_absorption), ground effect (general per-region method of 7.3.1 with the Table 3 a'/b'/c'/d' functions, Eq. 9, plus the alternative method of 7.3.2, Eq. 10, and the DOmega index Eq. 11), and screening (Eq. 12-18 with the C2/C3 factors, single/double diffraction, Kmet and the 20/25 dB caps). Adds the Cmet meteorological correction (Eq. 21/22) and the LfT(DW)=Lw+Dc-A receiver-level composition (Eq. 3/6). Frozen dataclass result exposes the per-term breakdown (Adiv/Aatm/Agr/Abar). New module, tests and exports. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * feat: ISO 9612 occupational noise exposure strategies and uncertainty Add phonometry.occupational_exposure implementing the three ISO 9612:2009 measurement strategies and the normative Annex C uncertainty budget: - task_based_exposure (Clause 9): per-task energy average (Eq 7), task contributions (Eq 8), daily level (Eq 9/10); uncertainty Eq C.3 with sampling u1a (Eq C.6), duration u1b (Eq C.7) and sensitivity coefficients c1a (Eq C.4) / c1b (Eq C.5). - job_based_exposure (Clause 10): effective-day energy average (Eq 11), daily level (Eq 12); Table 1 minimum-duration check; uncertainty Eq C.9. - full_day_exposure (Clause 11): whole-day averaging (Eq 11/13) with the Clause 11.3 three-measurement / 3 dB spread rule. - Table C.4 (job/full-day c1*u1, bilinear interpolation), Table C.5 (u2), u3 = 1.0 dB, coverage factor k = 1.65 (one-sided 95 %). Frozen dataclasses (Task, TaskContribution, ExposureResult), ExposureWarning advisories for the sampling rules, and exports appended to __init__. Tests reproduce the Annex D (task), E (job) and F (full-day) worked examples to the standard's precision, plus closed-form oracles (single 85 dB/8 h task = 85; equal-level tasks; -3.01 dB per halved duration; k = 1.65) and the Table C.4 / Table 1 lookups. make check green: 1167 passed, 12 skipped. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: ISO 9612 review corrections — Table C.4 digits and Eq C.7 denominator Two Table C.4 cells verified against the page image at 200 dpi and corrected (N=6/u1=3.5: 3.4 -> 3.3; N=20/u1=2.0: 0.6 -> 0.5), the test that echoed the wrong N=6 cell now pins the printed value plus the corrected N=20 cell. The duration_samples branch of Eq C.7 divides by J(J-1) like Eq C.6 (standard error of the mean), with a regression test pinning u1b = 1.0 h for samples (4, 6) h. Module docstring now states the Lp,Cpeak uncertainty scope cut (Table C.5, NOTE 1). Claude-Session: https://claude.ai/code/session_01HUHMx4dxDes2t5UPQ25djP * test: ISO 9613-2 and ISO 9612 conformance checks Add a new "Outdoor propagation & occupational exposure" conformance domain with seven closed-form / worked-example checks backed by shared constants in tests/reference_data.py: - ISO 9613-2:1996 Eq. (7) geometrical divergence Adiv = 51 dB at 100 m. - ISO 9613-2:1996 Table 3 ground limit b'(0) = 10,1 -> Agr(250 Hz) = 17,2 dB (porous ground, hs = hr = 0, dp -> inf). - ISO 9613-2:1996 clause 7.4 single/double diffraction 20/25 dB caps. - ISO 9612:2009 Annex D/E/F daily exposure LEX,8h and expanded uncertainty U (84,3/2,7, 88,1/3,8, 90,1/3,4 dB), the Task objects rebuilt from the shared input table. Registry grows 34 -> 41 checks across 8 domains; docs/CONFORMANCE.md regenerated via make conformance. Registration and reference-data consistency tests pin the constants to the published values. Claude-Session: https://claude.ai/code/session_01HUHMx4dxDes2t5UPQ25djP * docs: figures for outdoor propagation and occupational exposure Add four documentation figures (each EN/ES x light/dark) plus a setup diagram: - air_absorption_alpha: ISO 9613-1 alpha(f) over frequency at four temperature/humidity combinations (relaxation roll-off, f^2 growth). - outdoor_attenuation_breakdown: ISO 9613-2 per-term stacked bar (Adiv, Aatm, Agr, Abar) with the total A, for a 200 m path over porous ground with a 4 m barrier. - exposure_uncertainty: ISO 9612 Annex D task contributions with the daily LEX,8h and the one-sided 95 % upper limit LEX,8h + U. - diagram_outdoor_geometry: source-barrier-receiver geometry (dss/dsr/z, the blocked direct ray and the diffracted ray over the top edge). Spanish translations added to the generate_graphs / generate_diagrams dictionaries (decimal commas, UNE terminology). Every image visually audited. Claude-Session: https://claude.ai/code/session_01HUHMx4dxDes2t5UPQ25djP * docs: ISO 9613-1/2 outdoor propagation guide and ISO 9612 exposure (EN) Add the didactic Outdoor Sound Propagation guide (docs/outdoor-propagation.md) covering ISO 9613-1 atmospheric absorption and the ISO 9613-2 general method (divergence, atmospheric, ground and barrier terms, the OutdoorAttenuation breakdown, the alternative 7.3.2 ground method, Cmet), with executed snippets and honest scope notes (Amisc/reflections not implemented, 7.3.2 not auto-wired, per-band Cmet a convenience). Extend the Levels guide with the ISO 9612 occupational-exposure section (task/job/full-day strategies, the C.6 I(I-1) vs C.12 N-1 asymmetry, Table C.4, k = 1,65, LEX,8h + U upper limit, Lp,Cpeak uncertainty out of scope), reproducing the Annex D/E/F worked examples. Add the matching theory sections, API-reference rows for every new public name, README feature list + docs table, docs index and CHANGELOG entries. llms-full.txt / llms.txt regenerated (make llms) with the new page wired into gen_llms.py. GitHub-safe math throughout; every snippet executed with pasted outputs. Claude-Session: https://claude.ai/code/session_01HUHMx4dxDes2t5UPQ25djP * docs: outdoor propagation and exposure site twins (EN + ES) Add the English and Spanish Starlight twins of the new content: - guides/outdoor-propagation.md (EN + ES) — code-identical snippets to the docs page; Spanish prose uses UNE terminology ("ponderación temporal"), decimal commas outside math and code, and the _es figure variants. - The ISO 9612 section added to guides/levels.md (EN + ES). - reference/theory.md and reference/api.md rows (EN + ES). - Landing cards and standards list bumped 34 -> 37 (index.mdx EN + ES); guides/outdoor-propagation added to the sidebar. Claude-Session: https://claude.ai/code/session_01HUHMx4dxDes2t5UPQ25djP * fix: final branch review — API rows, precision phrasing and ES links Adds the five missing public-symbol rows (TaskContribution, table_c4_contribution, minimum_cumulative_duration_hours, COVERAGE_FACTOR, INSTRUMENT_U2) to the three API references, states the ISO 9613-2 default atmosphere (20 °C / 70 %, a Table 2 reference atmosphere) in the API rows and guides, and softens "digit-for-digit" to "the standard's printed precision" wherever it covered Annex E, whose final level differs only by the standard's own pre-rounding. ES pages now use the /es/ link prefix consistently. Code nits from the review: dead arg > 1 guard removed (z > 0 implies arg >= 3), _with_advisory uses dataclasses.replace, and the float-keyed _PRIME_BY_BAND gets a safety comment. llms files regenerated. Claude-Session: https://claude.ai/code/session_01HUHMx4dxDes2t5UPQ25djP * fix: validate propagation/exposure inputs and docs lint (PR review) Addresses the PR #89 bot review: - Barrier.__post_init__ rejects negative edge distances and a zero edge_separation (which previously divided by zero in the C3 double- diffraction factor); Task.__post_init__ rejects empty samples and a non-positive duration. - ground_attenuation and barrier_attenuation reject non-positive frequencies (a zero frequency gave an infinite wavelength / division by zero); ground_attenuation also checks a positive distance and non-negative heights. job_based_exposure rejects a non-positive sample_duration_hours. - Rewrite the ISO 9613-2 7.3.2 closed form in the guides and API tables so Markdownlint no longer parses `(...)[...]` as a reversed link (MD011). - Regression tests for every new guard; docstrings note the new raises. Claude-Session: https://claude.ai/code/session_01HUHMx4dxDes2t5UPQ25djP

2 weeks ago

time-weighting.md

Docs usability: measurement-role snippet clarity, one-line result plots, responsive figures (#86) * fix: responsive figure sizing — full-width on mobile, 92% plots on desktop (plan 16 T4) Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: snippet clarity wave — recording/excitation convention (plan 16 T1) Apply the recording/excitation naming convention from the P16 audit across the calibration, intensity, block-processing, time-weighting, levels, weighting and room-acoustics guides (docs/ + EN guides + ES guides). - calibration: ref_signal->calibrator_recording, signal->recording (J1 cascade) - intensity: probe-recording synth + swap-in disclosure on p1/p2 - block-processing / time-weighting: audio_blocks provenance comments - levels / time-weighting / weighting: signal->recording, provenance upgraded to convention - room-acoustics: play/record framing on sweep/recorded; synthetic-IR stand-in note EN docs/guides kept code-identical; ES mirrored with UNE wording. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * feat: one-line canonical plots on result objects (plan 16 T2) Add a thin plot(ax=None, **kwargs) method to each result dataclass that renders the essence of its guide's canonical figure in one line, following the sklearn/statsmodels/librosa pattern. matplotlib stays a soft dependency: importing phonometry and computing works without it; .plot() raises a clear ImportError with install guidance. Rendering lives in a new lazy _plotting module; dataclass methods are thin delegates. Covers ZwickerLoudness, STIResult, WeightedRatingResult/ImpactRatingResult, RoomAcousticsResult, SoundPowerResult/ReverberationSoundPowerResult/ SoundPowerIntensityResult, IntensityResult and a new frozen DecayCurve (returned by decay_curve, tuple-unpackable for backward compat). Rating results gain band_centers/measured/shifted_reference fields for the figure. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: plot kwargs forwarding, deterministic invalid-band test, octave-impact annotation - _plot_rating now accepts **kwargs and forwards to the measured curve, so result.plot(linewidth=2) no longer raises TypeError; parametrized test covers every public .plot() forwarding a benign kwarg to its primary artist. - Replace skip-if-all-valid room-acoustics invalid-band test with a directly constructed RoomAcousticsResult (500 Hz band flagged invalid); assert exactly the flagged bars are hatched and greyed. - plot_impact_rating keeps the shifted-reference curve normatively honest (ref - shift) and marks/annotates the 500 Hz read value with the ISO 717-2 Clause 4.3.2 -5 dB octave rule instead of distorting the curve. - Move _bar_width next to its only consumer, plot_intensity. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: .plot() one-liners and figure-code collapsibles (plan 16 T3) Add one-line res.plot() mentions to every guide section whose result type gained .plot() in T2, and figure-code collapsibles (both the one-liner and the manual matplotlib path) under the five .plot()-backed data figures. Document .plot() availability and the DecayCurve return type in the API reference and CHANGELOG. EN docs + EN/ES site; llms-full.txt regenerated. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: collapsibles for STI-vs-T60 and open-plan figures Controller follow-up to plan 16 T3: the two figures are reproducible with public API only. STI-vs-T60 sweeps sti_from_impulse_response over synthetic exponential IRs at a 10-point T60 grid; the open-plan collapsible rebuilds the D2,S regression and rD/rP crossings from the OpenPlanResult fields. EN docs + EN/ES site; executed on Agg; llms-full.txt regenerated. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: final-review polish — result-field docs, comment hygiene, naming Document band_centers/measured/shifted_reference on WeightedRatingResult and ImpactRatingResult (CHANGELOG + EN/ES API tables), strip internal plan tags from theme-images.css and test_result_plots.py, note that [plot] also powers result .plot() methods, tighten the plot_zwicker_loudness docstring to say the two-panel layout only occurs when ax is None, sentence-case the intensity provenance comments, and rename the shadowed impact rating var to res_imp on the room-acoustics pages. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: PR 86 review round 1 — log-axis bars, style kwargs overrides, snippet physics and notation Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf

2 weeks ago

Docs usability: measurement-role snippet clarity, one-line result plots, responsive figures (#86) * fix: responsive figure sizing — full-width on mobile, 92% plots on desktop (plan 16 T4) Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: snippet clarity wave — recording/excitation convention (plan 16 T1) Apply the recording/excitation naming convention from the P16 audit across the calibration, intensity, block-processing, time-weighting, levels, weighting and room-acoustics guides (docs/ + EN guides + ES guides). - calibration: ref_signal->calibrator_recording, signal->recording (J1 cascade) - intensity: probe-recording synth + swap-in disclosure on p1/p2 - block-processing / time-weighting: audio_blocks provenance comments - levels / time-weighting / weighting: signal->recording, provenance upgraded to convention - room-acoustics: play/record framing on sweep/recorded; synthetic-IR stand-in note EN docs/guides kept code-identical; ES mirrored with UNE wording. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * feat: one-line canonical plots on result objects (plan 16 T2) Add a thin plot(ax=None, **kwargs) method to each result dataclass that renders the essence of its guide's canonical figure in one line, following the sklearn/statsmodels/librosa pattern. matplotlib stays a soft dependency: importing phonometry and computing works without it; .plot() raises a clear ImportError with install guidance. Rendering lives in a new lazy _plotting module; dataclass methods are thin delegates. Covers ZwickerLoudness, STIResult, WeightedRatingResult/ImpactRatingResult, RoomAcousticsResult, SoundPowerResult/ReverberationSoundPowerResult/ SoundPowerIntensityResult, IntensityResult and a new frozen DecayCurve (returned by decay_curve, tuple-unpackable for backward compat). Rating results gain band_centers/measured/shifted_reference fields for the figure. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: plot kwargs forwarding, deterministic invalid-band test, octave-impact annotation - _plot_rating now accepts **kwargs and forwards to the measured curve, so result.plot(linewidth=2) no longer raises TypeError; parametrized test covers every public .plot() forwarding a benign kwarg to its primary artist. - Replace skip-if-all-valid room-acoustics invalid-band test with a directly constructed RoomAcousticsResult (500 Hz band flagged invalid); assert exactly the flagged bars are hatched and greyed. - plot_impact_rating keeps the shifted-reference curve normatively honest (ref - shift) and marks/annotates the 500 Hz read value with the ISO 717-2 Clause 4.3.2 -5 dB octave rule instead of distorting the curve. - Move _bar_width next to its only consumer, plot_intensity. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: .plot() one-liners and figure-code collapsibles (plan 16 T3) Add one-line res.plot() mentions to every guide section whose result type gained .plot() in T2, and figure-code collapsibles (both the one-liner and the manual matplotlib path) under the five .plot()-backed data figures. Document .plot() availability and the DecayCurve return type in the API reference and CHANGELOG. EN docs + EN/ES site; llms-full.txt regenerated. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: collapsibles for STI-vs-T60 and open-plan figures Controller follow-up to plan 16 T3: the two figures are reproducible with public API only. STI-vs-T60 sweeps sti_from_impulse_response over synthetic exponential IRs at a 10-point T60 grid; the open-plan collapsible rebuilds the D2,S regression and rD/rP crossings from the OpenPlanResult fields. EN docs + EN/ES site; executed on Agg; llms-full.txt regenerated. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: final-review polish — result-field docs, comment hygiene, naming Document band_centers/measured/shifted_reference on WeightedRatingResult and ImpactRatingResult (CHANGELOG + EN/ES API tables), strip internal plan tags from theme-images.css and test_result_plots.py, note that [plot] also powers result .plot() methods, tighten the plot_zwicker_loudness docstring to say the two-panel layout only occurs when ax is None, sentence-case the intensity provenance comments, and rename the shadowed impact rating var to res_imp on the room-acoustics pages. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: PR 86 review round 1 — log-axis bars, style kwargs overrides, snippet physics and notation Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf

2 weeks ago

why-phonometry.md

Didactic documentation: theory everywhere, bilingual figures, setup diagrams, parameter tables (#77) * docs: harden diagram engine per visual audit - dim() draws real witness lines from the measured points to the offset dimension line (the offset=0 path emitted zero-length lines); vertical labels sit left of the line, clear of masts and figures - env diagram: the 4 m dimension now measures the mic axis itself and the flush-mounted label moved to the clear zone between masts - tonality top view: the 1.00 m dimension connects face and microphone Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: bilingual diagram variants (EN/ES x light/dark) Every setup diagram is now emitted in four variants; Spanish strings (including decimal commas) come from a translation table applied by the SVG engine, and layouts were adjusted where the longer Spanish labels collided (env titles, processing-chain box width, calibration output label, tonality dimension label side). Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: bilingual figure variants for the ES site (EN/ES x light/dark) generate_graphs.py now emits every figure in four variants: a Spanish translation table (~75 exact strings + parametrized patterns) is applied to the figure's Text artists at savefig time, so the generator functions stay single-language internally. Spanish output uses UNE terminology (ponderación temporal, líneas isofónicas, fonios, ráfaga) and a uniform decimal comma (annotations and tick labels included); code identifiers and normative mode names (Fast/Slow/Impulse) stay untranslated. The ES site pages now reference the _es/_es_dark variants (56 URLs). The ES set passed a 10-figure visual audit (10/10 approved, no blockers). Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: didactic calibration page - theory, setup diagram, stability figure, parameter table Adds the sensitivity-factor derivation with the calibration-chain SVG, a new autogenerated calibration_stability figure (stable vs 3% AM tone against the IEC 60942:2017 class 1 deviation band; 4 language/theme variants; the ES decimal-comma pass now preserves clause numbers like 5.3.3), and a full parameter table for calculate_sensitivity - EN and ES in parity. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: SEL concept, Lden 24h profile and annotated TNR spectrum figures Three new didactic figures for the levels guide, each in four language/theme variants: sel_concept (pass-by event vs its Leq and the equal-energy 1 s SEL block, dBFS mode), lden_profile (synthetic urban 24 h LAeq with colored day/evening/night bands and the +5/+10 weighted period levels clearly separated from the composite Lden), and tonality_spectrum (averaged Hann spectrum with the critical band, the detected tone and the live TNR value vs its criterion, computed by the library itself). Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: didactic levels guide, theory extensions and refreshed why-phonometry - levels guide (EN docs+site): energy-mean rationale for Leq, percentile semantics for LN, SEL and noise-dose theory with the new figures embedded, TNR critical-band explanation with the annotated spectrum, Lden profile figure, and parameter tables (type/units/range/default) for every level function including the spectrogram - theory page: five new sections (G-weighting pole/zero derivation, ISO 226 formulae, ECMA-418-1 TNR/PR math, event/dose metrics, ISO 1996-1 descriptors) with cross-references to the guides - why-phonometry: rewritten around the conformance-testing story with the per-standard verification table and honest ecosystem positioning Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: weighting/time-weighting theory and parameter tables; fix math-block delimiters The physiological origin of A/C (inverted 40/100-phon contours, the four IEC corner frequencies) and the exponential-detector derivation with the tau semantics, each with full parameter tables. Display math inserted by the docs expansion now always uses block-style $$ delimiters on their own lines - mid-line multi-line $$ blocks silently derailed remark parsing and truncated every heading after them (caught because the built pages lost their anchor ids); also replaced a raw < in a table. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: processing-chain and multirate diagrams, filter-bank math, streaming patterns - getting-started embeds the processing-chain diagram as a guided map of the guides - filter-banks gains the IEC 61260-1 base-10 band math (G, fm, band edges), a new multirate-decimation diagram (d5, four language/theme variants) and the full octavefilter/OctaveFilterBank parameter table - block-processing gains the canonical real-time level-meter loop (validated against the actual API - TimeWeighting is inherently stateful) and a stateful-constraints table - multichannel gains the accepted-shapes table Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: Spanish parity for the didactic expansion Every enriched page now has its Spanish twin: levels (theory, three new figures via _es variants, six parameter tables), weighting and time-weighting (curve-origin and exponential-detector sections), filter-banks (fractional-band math, multirate diagram, full parameter table - terminology unified on 'decimación'), getting-started (processing-chain map), block-processing (real-time meter pattern and stateful constraints), multichannel (shapes table), theory (five new sections) and a full rewrite of why-phonometry around the conformance story. API reference now covers all 27 public exports in EN and ES (added ToneAssessment and CalibrationWarning rows). UNE terminology and Spanish decimal commas throughout; site builds with all links valid. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs: audited fixes - Spanish tick commas and multirate note Final visual-audit round: the d5 diagram now shows the no-decimation note on the 16 kHz row in both languages, and Spanish figures render every numeric tick with a decimal comma (FixedFormatter sequences are rewritten in place; set_xticklabels-style FuncFormatters and plain linear ScalarFormatters are wrapped - the first naive wrapper broke ScalarFormatter ticks entirely and was caught visually). Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * docs(site): language selector in the mobile header, scrollable display math Custom Header component (default markup plus a compact icon-only LanguageSelect visible below the md breakpoint, where Starlight hides the whole right group behind the hamburger menu) and overflow-x on katex-display so wide formulas scroll inside their box on phones. Verified with Playwright at 390 px in EN and ES. Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf * fix: address review feedback on the didactic docs - levels tables: spectrogram calibration_factor/dbfs marked constructor- only, overlap range as 0 <= overlap < 1, composite hours constraints, ln_levels shared-semantics row filled; ES translation leftovers and the 'delante' typo fixed - block-processing: y[..., -1] so the meter loop is multichannel-safe; high_accuracy wording matches the real behavior (default resolves to False, explicit True raises) - API: CalibrationWarning row states the fs/validate preconditions (3 files) - weighting: high_accuracy default column and oversampling claim made precise (8x cap; silently ignored for G); ES row in parity - filter-banks: 25 Hz relative bandwidth corrected to 0.024% of Nyquist (EN/ES); verify_filter_class classes documented as 1/2/None - generate_graphs: dropped the global plt.savefig monkey-patch - the translation now runs inside themed_path(), which every generator already calls at save time (verified output identical) Claude-Session: https://claude.ai/code/session_013kkVt3nxi9svp1an28uHxf

2 weeks ago

README.md

phonometry Documentation#

Full documentation for phonometry. Also available as a website: https://jmrplens.github.io/phonometry/

Guides#

Getting Started — installation and first analysis
Filter Banks — architectures, responses, band decomposition
Frequency Weighting — A, C, Z curves
Time Weighting — Fast, Slow, Impulse ballistics
Integrated & Statistical Levels — Leq, LAeq, L10/L50/L90, noise dose and occupational exposure (ISO 9612), octave spectrogram
Psychoacoustics and Speech Intelligibility — Zwicker loudness, sharpness, STI/STIPA
Speech Intelligibility Index — the ANSI S3.5-1997 one-third-octave-band SII: band-importance weighting (Table 3), self-speech and upward spread of masking, band audibility, and the index in noise and hearing loss
Room-noise criteria — the ANSI/ASA S12.2-2019 room-noise ratings: the NC tangency method (Table 1) and the RC Mark II rating with its rumble/hiss/neutral spectral tag (Annex D)
Hearing threshold — the age-related hearing threshold distribution (ISO 7029:2017) and the free-field/diffuse-field reference threshold of hearing (ISO 389-7:2006)
Noise-induced hearing loss — the ISO 1999:2013 noise-induced permanent threshold shift (NIPTS) and its population distribution, and the combination with age into the hearing threshold level associated with age and noise (HTLAN)
Impulsive-sound prominence — the NT ACOU 112:2002 predicted prominence of impulsive sounds (onset rate and level difference) and the graduated adjustment KI added to LAeq
Measurement uncertainty — the GUM law of propagation of uncertainty and the Monte Carlo method (ISO/IEC Guide 98-3:2008 and Supplement 1): combined and expanded uncertainty, Welch–Satterthwaite effective degrees of freedom, and probabilistically symmetric coverage intervals
Sound Intensity (p-p) — two-microphone intensity and field indicators
Room and Building Acoustics — impulse-response acquisition, reverberation and room parameters, open-plan speech metrics, field airborne/impact/façade sound insulation (ISO 16283-1/2/3), laboratory characterisation (ISO 10140), flanking-transmission prediction (EN 12354-1/2), measurement uncertainty (ISO 12999-1), sound absorption (ISO 354)
Sound absorption in enclosed spaces — the EN 12354-6:2003 prediction of a room's total equivalent absorption area and reverberation time from its surfaces and objects (Clause 4)
Outdoor Sound Propagation — atmospheric absorption α(f) (ISO 9613-1) and the ISO 9613-2 general method: divergence, atmospheric absorption, ground effect and barrier screening (occupational exposure ISO 9612 lives in Levels)
Sound Power — sound power level by enveloping surface (ISO 3744/3746), reverberation room (ISO 3741), intensity scanning (ISO 9614-2), and the precision grades in an anechoic room (ISO 3745) and by precision intensity scanning (ISO 9614-3)
Acoustic Materials — sound-absorption rating α_w and classes (ISO 11654), airflow resistance static and alternating methods (ISO 9053-1/-2), and impedance-tube measurement of absorption, surface impedance and transmission loss (ISO 10534-1/-2, ASTM E2611)
Surface Scattering, Diffusion and In-situ Absorption — random-incidence scattering (ISO 17497-1), free-field diffusion coefficient (ISO 17497-2), and in-situ road-surface absorption by the extended-surface subtraction technique (ISO 13472-1) and the spot method (ISO 13472-2)
Human Vibration — whole-body and hand-arm frequency weightings (ISO 8041-1), weighted r.m.s. acceleration, running r.m.s./MTVV/VDV and crest factor (ISO 2631-1), vibration in buildings (ISO 2631-2), vibration total value and daily exposure A(8) (ISO 5349-1/-2), and the exposure action/limit values of Directive 2002/44/EC
Multiple-shock whole-body vibration — the ISO 2631-5:2018 spinal-response model: the seat-to-spine transfer function, the acceleration and daily dose from the response peaks (Clause 5), and the compressive stress, stress variable R and Weibull probability of lumbar injury (Annex C)
Calibration and dBFS — physical SPL and digital analysis
Block Processing — stateful real-time workflows
Multichannel — vectorized multichannel analysis

Reference#

API Reference — every public function and class
Theory — standards, math and design decisions
Why phonometry — IEC compliance vs other libraries

Development#

Run the test suite with pytest tests/, the full quality gate with make check (ruff + mypy + bandit + tests), and regenerate the documentation images with make graphs. See CONTRIBUTING.md.