#  Copyright (c) 2020. Jose M. Requena-Plens
"""
Octave-Band and Fractional Octave-Band filter for signals in the time domain.
Implementation according to ANSI s1.11-2004 and IEC 61260-1-2014.
"""

from __future__ import annotations

from functools import lru_cache
from typing import Any, List, Tuple, overload, Literal

import numpy as np

from .calibration import CalibrationWarning, calculate_sensitivity, sensitivity
from .environmental import composite_rating_level, lden, ldn
from .environmental_measurement import (
    RepeatedMeasurementResult,
    ResidualCorrectionResult,
    TonalAssessmentResult,
    assess_tonal_audibility,
    combined_standard_uncertainty,
    critical_bandwidth,
    gaussian_residual_level,
    residual_correction_uncertainty,
    residual_sound_correction,
    tonal_adjustment,
    tonal_adjustment_from_mean_audibility,
    tonal_audibility,
    tonal_seeking_survey,
    uncertainty_from_repeated_measurements,
)
from .environmental_measurement import (
    expanded_uncertainty as environmental_expanded_uncertainty,
)
from .compliance import (
    class_limits,
    verify_filter_class,
    verify_weighting_class,
    weighting_class_limits,
)
from .core import FilterBankWarning, OctaveFilterBank
from .frequencies import (
    getansifrequencies,
    nominal_frequencies,
    normalized_frequencies,
    normalizedfreq,
)
from .intensity import (
    FieldIndicators,
    IntensityResult,
    dynamic_capability_index,
    field_indicators,
    sound_intensity,
)
from .levels import laeq, lc_peak, leq, lex_8h, ln_levels, sel, sound_exposure
from .loudness_zwicker import ZwickerLoudness, loudness_zwicker, loudness_zwicker_from_spectrum
from .loudness_ecma import EcmaLoudness, loudness_ecma
from .loudness_moore_glasberg import (
    MooreGlasbergLoudness,
    loudness_moore_glasberg,
    loudness_moore_glasberg_from_spectrum,
    loudness_moore_glasberg_from_third_octave,
)
from .loudness_moore_glasberg_time import (
    MooreGlasbergTimeVaryingLoudness,
    loudness_moore_glasberg_time,
)
from .loudness_contours import equal_loudness_contour, hearing_threshold, loudness_level
from .sharpness import sharpness_din, sharpness_din_from_specific
from .hearing import (
    AgeThresholdResult,
    age_threshold,
    reference_threshold,
)
from .noise_induced_hearing_loss import (
    HtlanResult,
    NiptsResult,
    htlan,
    nipts,
)
from .enclosed_space_absorption import (
    ReverberationResult,
    air_absorption_area,
    enclosed_space_reverberation,
    equivalent_absorption_area,
    hard_object_absorption,
    object_fraction,
    reverberation_time,
)
from .reverberation_prediction import (
    ReverberationModelResult,
    arau_puchades_reverberation_time,
    eyring_reverberation_time,
    fitzroy_reverberation_time,
    mean_absorption,
    millington_sette_reverberation_time,
    reverberation_time_models,
    sabine_reverberation_time,
)
from .dynamic_stiffness import (
    DynamicStiffnessResult,
    DynamicStiffnessWarning,
    apparent_dynamic_stiffness,
    enclosed_gas_stiffness,
    floating_floor_resonance,
    installed_dynamic_stiffness,
    natural_frequency,
)
from .mechanical_mobility import (
    FRF_UNITS,
    MobilityResult,
    convert_frf,
    resonance_frequency,
    sdof_accelerance,
    sdof_mobility,
    sdof_mobility_result,
    sdof_receptance,
)
from .transfer_stiffness import (
    REFERENCE_STIFFNESS,
    TransferStiffnessResult,
    base_transmissibility,
    indirect_transfer_stiffness_result,
    loss_factor,
    transfer_stiffness_direct,
    transfer_stiffness_indirect,
    transfer_stiffness_level,
)
from .vibration_sound_power import (
    NORMALIZED_IMPEDANCE,
    REFERENCE_VELOCITY,
    VibrationSoundPowerResult,
    extraneous_velocity_correction,
    mean_velocity_level,
    radiated_sound_power_level,
    radiation_factor,
    sound_power_from_vibration,
    velocity_level,
    velocity_level_from_acceleration,
)
from .structure_borne_power import (
    StructureBornePowerResult,
    plate_loss_factor,
    reception_plate_power,
    spatial_mean_velocity_level,
    structure_borne_power_level,
)
from .installed_structure_borne import (
    REFERENCE_AREA,
    InstalledSourceResult,
    coupling_term,
    coupling_term_force_source,
    coupling_term_velocity_source,
    installed_source_prediction,
    installed_structure_borne_power_level,
    structure_borne_pressure_level_path,
    total_structure_borne_pressure_level,
)
from .multiple_shock_vibration import (
    MultipleShockResult,
    acceleration_dose,
    compression_dose,
    daily_dose,
    daily_dose_multi,
    dose_from_peaks,
    injury_probability,
    injury_risk,
    multiple_shock_assessment,
    response_peaks,
    seat_to_spine_transfer,
    spinal_response,
    static_stress,
    ultimate_strength,
)
from .impulse_prominence import (
    ImpulseProminenceResult,
    impulse_adjustment,
    impulse_prominence,
    predicted_prominence,
    rating_level,
)
from .uncertainty import (
    MonteCarloResult,
    Quantity,
    UncertaintyResult,
    combine_uncertainty,
    monte_carlo,
    rectangular,
    triangular,
    u_shaped,
)
from .room_noise import (
    NCResult,
    RCResult,
    nc_curve,
    noise_criterion,
    rc_curve,
    room_criterion,
)
from .sii import (
    SIIResult,
    speech_intelligibility_index,
    standard_speech_spectrum,
)
from .sti import (
    STIResult,
    STIWarning,
    sti_from_impulse_response,
    stipa,
    stipa_signal,
)
from .tonality_ecma import EcmaTonality, tonality_ecma
from .roughness_ecma import EcmaRoughness, roughness_ecma
from .tonality import (
    TonalityWarning,
    ToneAssessment,
    prominence_ratio,
    tone_to_noise_ratio,
)
from .parametric_filters import (
    TimeWeighting,
    WeightingFilter,
    linkwitz_riley,
    time_weighting,
    weighting_filter,
)
from .insulation import (
    AirborneInsulationResult,
    FacadeInsulationResult,
    ImpactInsulationResult,
    ImpactRatingResult,
    WeightedRatingResult,
    airborne_insulation,
    energy_average_level,
    facade_insulation,
    impact_insulation,
    weighted_impact_improvement,
    weighted_impact_rating,
    weighted_rating,
)
from .lab_insulation import (
    LabAirborneInsulationResult,
    LabImpactInsulationResult,
    LabInsulationWarning,
    background_correction,
    lab_airborne_insulation,
    lab_impact_insulation,
)
from .intensity_insulation import (
    IntensityElementNormalizedResult,
    IntensityReductionResult,
    adaptation_term_kc,
    combine_subareas,
    intensity_element_normalized_difference,
    intensity_sound_reduction,
    surface_pressure_intensity_indicator,
)
from .floor_covering_improvement import (
    FloorCoveringImprovementResult,
    acceleration_level,
    background_corrected_level,
    impact_improvement,
    improvement_octave_bands,
)
from .survey_insulation import (
    SurveyAirborneResult,
    SurveyFacadeResult,
    SurveyImpactResult,
    SurveyServiceEquipmentResult,
    estimate_reverberation_index,
    reverberation_index,
    survey_airborne_insulation,
    survey_facade_insulation,
    survey_impact_insulation,
    survey_service_equipment_level,
)
from .open_plan import OpenPlanResult, open_plan_metrics
from .sound_power import (
    MeteorologicalCorrection,
    PrecisionCriteria,
    PrecisionFieldIndicators,
    PrecisionIntensityResult,
    PrecisionSoundPowerResult,
    SoundPowerResult,
    SoundPowerWarning,
    background_noise_correction,
    environmental_correction,
    measurement_positions,
    meteorological_corrections,
    precision_background_correction,
    precision_field_indicators,
    precision_positions,
    precision_qualification,
    precision_uncertainty,
    sound_power_anechoic,
    sound_power_intensity_precision,
    sound_power_pressure,
)
from .sound_power_reverberation import (
    ReverberationSoundPowerResult,
    sound_power_comparison,
    sound_power_reverberation,
)
from .sound_power_intensity import (
    SoundPowerIntensityResult,
    sound_power_intensity,
)
from .scattering_diffusion import (
    BASE_PLATE_BANDS,
    BASE_PLATE_MAX_SCATTERING,
    TWO_DIMENSIONAL_SOURCE_WEIGHTS,
    DiffusionResult,
    ScatteringDiffusionWarning,
    ScatteringResult,
    ScatteringUncertainty,
    absorption_coefficient_uncertainty,
    air_attenuation_coefficient,
    area_factors,
    base_plate_scattering,
    check_base_plate_scattering,
    directional_diffusion,
    directional_diffusion_coefficient,
    normalized_diffusion_coefficient,
    random_incidence_absorption,
    random_incidence_diffusion,
    reverberation_time_uncertainty,
    scattering_coefficient,
    scattering_coefficient_spectrum,
    scattering_coefficient_uncertainty,
    specular_absorption_coefficient,
    speed_of_sound,
)
from .road_absorption import (
    DEFAULT_MIC_HEIGHT,
    DEFAULT_SOURCE_HEIGHT,
    DEFAULT_SPEED_OF_SOUND,
    PART1_FREQUENCY_RANGE,
    SPOT_FREQUENCY_RANGE,
    SPOT_NARROW_BAND_RANGE,
    InsituAbsorptionResult,
    RoadAbsorptionWarning,
    absorption_reference_corrected,
    adrienne_window,
    check_spot_frequency_range,
    geometric_spreading_factor,
    geometric_spreading_factor_angle,
    insitu_absorption_coefficient,
    insitu_absorption_from_reflection,
    insitu_absorption_spectrum,
    insitu_reflection_factor,
    max_sampled_area_radius,
    msa_major_axis,
    one_third_octave_absorption,
    power_reflection_coefficient,
    reflected_path_delay,
    spot_internal_loss_correction,
    spot_microphone_spacing_bounds,
    spot_tube_upper_frequency,
)
from .human_vibration import (
    HAV_EAV_A8,
    HAV_ELV_A8,
    REFERENCE_ACCELERATION,
    REFERENCE_DURATION_S,
    WBV_EAV_A8,
    WBV_EAV_VDV,
    WBV_ELV_A8,
    WBV_ELV_VDV,
    WEIGHTING_NAMES,
    DailyVibrationExposure,
    ExposureAssessment,
    HumanVibrationWarning,
    WeightedSpectrum,
    WeightingResponse,
    apply_weighting,
    combine_partial_exposures,
    crest_factor,
    daily_exposure,
    daily_vibration_exposure,
    energy_equivalent_acceleration,
    exposure_assessment,
    frequency_weighting,
    hav_daily_exposure,
    hav_vwf_lifetime_years,
    motion_sickness_dose_value,
    mtvv,
    partial_exposure,
    running_rms,
    vibration_dose_value,
    vibration_total_value,
    weighted_acceleration,
    weighting_factors,
)
from .room_acoustics import (
    DecayCurve,
    RoomAcousticsResult,
    decay_curve,
    room_parameters,
)
from .sound_absorption import (
    AbsorptionWarning,
    absorption_area,
    absorption_coefficient,
    attenuation_from_alpha,
)
from .air_absorption import (
    AtmosphericAbsorptionWarning,
    air_attenuation,
    air_attenuation_m,
)
from .impedance_tube import (
    ImpedanceTubeResult,
    ImpedanceTubeWarning,
    TransferMatrix,
    absorption_from_reflection,
    air_density_astm,
    air_density_iso,
    air_layer_transfer_matrix,
    apply_mic_calibration,
    characteristic_impedance,
    face_quantities,
    mic_calibration_factor,
    normalized_surface_admittance,
    normalized_surface_impedance,
    plane_wave_frequency_range,
    reflection_factor,
    speed_of_sound_astm,
    speed_of_sound_iso,
    standing_wave_absorption,
    standing_wave_normalized_impedance,
    standing_wave_ratio_from_level,
    standing_wave_reflection,
    standing_wave_reflection_magnitude,
    surface_impedance,
    transfer_matrix_one_load,
    transfer_matrix_two_load,
    tube_attenuation_constant,
    tube_wavenumber,
    two_microphone_impedance,
    wave_decomposition,
)
from .airflow_resistance import (
    AirflowResistanceWarning,
    StaticAirflowResult,
    airflow_resistance,
    airflow_resistivity,
    alternating_airflow_resistance,
    effective_kappa,
    linear_airflow_velocity,
    piston_volume_flow_rate,
    specific_airflow_resistance,
    static_airflow_resistance,
    thermal_boundary_layer_thickness,
)
from .absorption_rating import (
    OCTAVE_BANDS,
    REFERENCE_CURVE,
    THIRD_OCTAVE_BANDS,
    AbsorptionRatingResult,
    absorption_class,
    practical_absorption_coefficient,
    weighted_absorption,
)
from .outdoor_propagation import (
    DEFAULT_FREQUENCIES,
    Barrier,
    OutdoorAttenuation,
    atmospheric_absorption,
    barrier_attenuation,
    directivity_omega,
    geometric_divergence,
    ground_attenuation,
    ground_attenuation_alternative,
    meteorological_correction,
    outdoor_propagation_attenuation,
    predicted_receiver_level,
)
from .room_ir import (
    ImpulseResponseResult,
    ImpulseResponseWarning,
    impulse_response,
    inverse_filter,
    mls_impulse_response,
    mls_signal,
    sweep_signal,
)
from ._plotting import plot_excitation
from ._warnings import PhonometryWarning, _warn_renamed
from .building_prediction import (
    AirbornePredictionResult,
    FlankingPath,
    ImpactPredictionResult,
    PathContribution,
    combine_linings,
    equivalent_impact_level,
    flanking_element,
    flanking_path,
    impact_flanking_correction,
    junction_min_vibration_reduction,
    junction_vibration_reduction,
    predicted_airborne_insulation,
    predicted_impact_insulation,
    standardized_impact_level,
)
from .facade_prediction import (
    FacadeElement,
    FacadePredictionResult,
    RadiatedPowerResult,
    facade_sound_reduction,
    outdoor_attenuation,
    outdoor_level,
    radiated_sound_power,
)
from .flanking_transmission import (
    FlankingImpactLevelResult,
    FlankingLevelDifferenceResult,
    VibrationReductionResult,
    band_mode_count,
    critical_frequency,
    direction_averaged_level_difference,
    equivalent_absorption_length,
    modal_density,
    modal_overlap_factor,
    normalized_flanking_impact_level,
    normalized_flanking_level_difference,
    strong_coupling_satisfied,
    total_loss_factor,
    velocity_level_difference,
    vibration_reduction_index,
    vibration_reduction_index_from_flanking,
)
from .building_uncertainty import (
    COVERAGE_FACTORS,
    BandUncertainty,
    UncertainValue,
    band_uncertainty,
    combine_uncertainties,
    coverage_factor,
    expanded_uncertainty,
    insulation_coverage_factor,
    insulation_expanded_uncertainty,
    maximum_repeatability_standard_deviation,
    prediction_input_uncertainty,
    reduce_by_independent_measurements,
    satisfies_lower_requirement,
    satisfies_upper_requirement,
    single_number_uncertainty,
    single_number_uncertainty_uncorrelated,
    uncertain_value,
)
from .absorption_uncertainty import (
    AbsorptionUncertaintyResult,
    sound_absorption_coefficient_uncertainty,
    absorption_coverage_factor,
    equivalent_area_uncertainty,
    practical_coefficient_uncertainty,
    single_number_rating_uncertainty,
    weighted_coefficient_uncertainty,
)
from .occupational_exposure import (
    COVERAGE_FACTOR,
    INSTRUMENT_U2,
    ExposureResult,
    OccupationalExposureWarning,
    Task,
    TaskContribution,
    full_day_exposure,
    job_based_exposure,
    minimum_cumulative_duration_hours,
    table_c4_contribution,
    task_based_exposure,
)
from ._version import __version__

# Public methods
__all__ = [
    "PhonometryWarning",
    "__version__",
    "octave_filter",
    "nominal_frequencies",
    "normalized_frequencies",
    # Deprecated aliases (remove in 4.0)
    "octavefilter",
    "getansifrequencies",
    "normalizedfreq",
    "OctaveFilterBank",
    "FilterBankWarning",
    "WeightingFilter",
    "weighting_filter",
    "time_weighting",
    "TimeWeighting",
    "linkwitz_riley",
    "sensitivity",
    "calculate_sensitivity",  # deprecated alias (remove in 4.0)
    "leq",
    "laeq",
    "ln_levels",
    "lc_peak",
    "sel",
    "sound_exposure",
    "lex_8h",
    "loudness_zwicker",
    "loudness_zwicker_from_spectrum",
    "ZwickerLoudness",
    "loudness_ecma",
    "EcmaLoudness",
    "loudness_moore_glasberg",
    "loudness_moore_glasberg_from_spectrum",
    "loudness_moore_glasberg_from_third_octave",
    "MooreGlasbergLoudness",
    "loudness_moore_glasberg_time",
    "MooreGlasbergTimeVaryingLoudness",
    "tonality_ecma",
    "EcmaTonality",
    "roughness_ecma",
    "EcmaRoughness",
    "sharpness_din",
    "sharpness_din_from_specific",
    "equal_loudness_contour",
    "loudness_level",
    "hearing_threshold",
    "tone_to_noise_ratio",
    "prominence_ratio",
    "ToneAssessment",
    "TonalityWarning",
    "sound_intensity",
    "IntensityResult",
    "field_indicators",
    "FieldIndicators",
    "dynamic_capability_index",
    "sti_from_impulse_response",
    "stipa",
    "stipa_signal",
    "STIResult",
    "STIWarning",
    "speech_intelligibility_index",
    "standard_speech_spectrum",
    "SIIResult",
    "noise_criterion",
    "room_criterion",
    "nc_curve",
    "rc_curve",
    "NCResult",
    "RCResult",
    "age_threshold",
    "reference_threshold",
    "AgeThresholdResult",
    "nipts",
    "htlan",
    "NiptsResult",
    "HtlanResult",
    "predicted_prominence",
    "impulse_adjustment",
    "impulse_prominence",
    "rating_level",
    "ImpulseProminenceResult",
    "seat_to_spine_transfer",
    "spinal_response",
    "response_peaks",
    "dose_from_peaks",
    "acceleration_dose",
    "daily_dose",
    "daily_dose_multi",
    "compression_dose",
    "static_stress",
    "ultimate_strength",
    "injury_risk",
    "injury_probability",
    "multiple_shock_assessment",
    "MultipleShockResult",
    "object_fraction",
    "hard_object_absorption",
    "air_absorption_area",
    "equivalent_absorption_area",
    "reverberation_time",
    "enclosed_space_reverberation",
    "ReverberationResult",
    "mean_absorption",
    "sabine_reverberation_time",
    "eyring_reverberation_time",
    "millington_sette_reverberation_time",
    "fitzroy_reverberation_time",
    "arau_puchades_reverberation_time",
    "reverberation_time_models",
    "ReverberationModelResult",
    "apparent_dynamic_stiffness",
    "enclosed_gas_stiffness",
    "installed_dynamic_stiffness",
    "natural_frequency",
    "floating_floor_resonance",
    "DynamicStiffnessResult",
    "DynamicStiffnessWarning",
    "convert_frf",
    "sdof_receptance",
    "sdof_mobility",
    "sdof_accelerance",
    "sdof_mobility_result",
    "resonance_frequency",
    "MobilityResult",
    "FRF_UNITS",
    "transfer_stiffness_level",
    "loss_factor",
    "transfer_stiffness_direct",
    "transfer_stiffness_indirect",
    "base_transmissibility",
    "indirect_transfer_stiffness_result",
    "TransferStiffnessResult",
    "REFERENCE_STIFFNESS",
    "velocity_level",
    "velocity_level_from_acceleration",
    "mean_velocity_level",
    "radiation_factor",
    "radiated_sound_power_level",
    "extraneous_velocity_correction",
    "sound_power_from_vibration",
    "VibrationSoundPowerResult",
    "REFERENCE_VELOCITY",
    "NORMALIZED_IMPEDANCE",
    "spatial_mean_velocity_level",
    "plate_loss_factor",
    "structure_borne_power_level",
    "reception_plate_power",
    "StructureBornePowerResult",
    "coupling_term",
    "coupling_term_force_source",
    "coupling_term_velocity_source",
    "installed_structure_borne_power_level",
    "structure_borne_pressure_level_path",
    "total_structure_borne_pressure_level",
    "installed_source_prediction",
    "InstalledSourceResult",
    "REFERENCE_AREA",
    "combine_uncertainty",
    "monte_carlo",
    "rectangular",
    "triangular",
    "u_shaped",
    "Quantity",
    "UncertaintyResult",
    "MonteCarloResult",
    "sweep_signal",
    "inverse_filter",
    "impulse_response",
    "mls_signal",
    "mls_impulse_response",
    "ImpulseResponseResult",
    "ImpulseResponseWarning",
    "plot_excitation",
    "room_parameters",
    "decay_curve",
    "DecayCurve",
    "RoomAcousticsResult",
    "absorption_area",
    "absorption_coefficient",
    "attenuation_from_alpha",
    "AbsorptionWarning",
    "air_attenuation",
    "air_attenuation_m",
    "AtmosphericAbsorptionWarning",
    # Impedance tube (ISO 10534-1/-2, ASTM E2611)
    "ImpedanceTubeResult",
    "ImpedanceTubeWarning",
    "TransferMatrix",
    "absorption_from_reflection",
    "air_density_astm",
    "air_density_iso",
    "air_layer_transfer_matrix",
    "apply_mic_calibration",
    "characteristic_impedance",
    "face_quantities",
    "mic_calibration_factor",
    "normalized_surface_admittance",
    "normalized_surface_impedance",
    "plane_wave_frequency_range",
    "reflection_factor",
    "speed_of_sound_astm",
    "speed_of_sound_iso",
    "standing_wave_absorption",
    "standing_wave_normalized_impedance",
    "standing_wave_ratio_from_level",
    "standing_wave_reflection",
    "standing_wave_reflection_magnitude",
    "surface_impedance",
    "transfer_matrix_one_load",
    "transfer_matrix_two_load",
    "tube_attenuation_constant",
    "tube_wavenumber",
    "two_microphone_impedance",
    "wave_decomposition",
    # Airflow resistance (ISO 9053-1/-2)
    "AirflowResistanceWarning",
    "StaticAirflowResult",
    "airflow_resistance",
    "airflow_resistivity",
    "alternating_airflow_resistance",
    "effective_kappa",
    "linear_airflow_velocity",
    "piston_volume_flow_rate",
    "specific_airflow_resistance",
    "static_airflow_resistance",
    "thermal_boundary_layer_thickness",
    # Absorption rating (ISO 11654)
    "AbsorptionRatingResult",
    "OCTAVE_BANDS",
    "REFERENCE_CURVE",
    "THIRD_OCTAVE_BANDS",
    "absorption_class",
    "practical_absorption_coefficient",
    "weighted_absorption",
    "DEFAULT_FREQUENCIES",
    "Barrier",
    "OutdoorAttenuation",
    "atmospheric_absorption",
    "barrier_attenuation",
    "directivity_omega",
    "geometric_divergence",
    "ground_attenuation",
    "ground_attenuation_alternative",
    "meteorological_correction",
    "outdoor_propagation_attenuation",
    "predicted_receiver_level",
    "open_plan_metrics",
    "OpenPlanResult",
    "sound_power_pressure",
    "measurement_positions",
    "background_noise_correction",
    "environmental_correction",
    "SoundPowerResult",
    "SoundPowerWarning",
    "sound_power_reverberation",
    "sound_power_comparison",
    "ReverberationSoundPowerResult",
    "sound_power_intensity",
    "SoundPowerIntensityResult",
    # ISO 3745 / ISO 9614-3 precision sound power
    "sound_power_anechoic",
    "PrecisionSoundPowerResult",
    "precision_positions",
    "precision_background_correction",
    "meteorological_corrections",
    "MeteorologicalCorrection",
    "precision_uncertainty",
    "sound_power_intensity_precision",
    "PrecisionIntensityResult",
    "precision_field_indicators",
    "PrecisionFieldIndicators",
    "precision_qualification",
    "PrecisionCriteria",
    # ISO 17497-1/-2 scattering & diffusion
    "speed_of_sound",
    "air_attenuation_coefficient",
    "random_incidence_absorption",
    "specular_absorption_coefficient",
    "scattering_coefficient",
    "scattering_coefficient_spectrum",
    "ScatteringResult",
    "base_plate_scattering",
    "BASE_PLATE_BANDS",
    "BASE_PLATE_MAX_SCATTERING",
    "check_base_plate_scattering",
    "reverberation_time_uncertainty",
    "absorption_coefficient_uncertainty",
    "scattering_coefficient_uncertainty",
    "ScatteringUncertainty",
    "directional_diffusion",
    "directional_diffusion_coefficient",
    "DiffusionResult",
    "normalized_diffusion_coefficient",
    "area_factors",
    "random_incidence_diffusion",
    "TWO_DIMENSIONAL_SOURCE_WEIGHTS",
    "ScatteringDiffusionWarning",
    # ISO 13472-1/-2 in-situ road-surface absorption
    "adrienne_window",
    "geometric_spreading_factor",
    "geometric_spreading_factor_angle",
    "reflected_path_delay",
    "insitu_reflection_factor",
    "insitu_absorption_from_reflection",
    "power_reflection_coefficient",
    "insitu_absorption_coefficient",
    "insitu_absorption_spectrum",
    "InsituAbsorptionResult",
    "absorption_reference_corrected",
    "one_third_octave_absorption",
    "max_sampled_area_radius",
    "msa_major_axis",
    "spot_tube_upper_frequency",
    "spot_microphone_spacing_bounds",
    "check_spot_frequency_range",
    "spot_internal_loss_correction",
    "DEFAULT_SOURCE_HEIGHT",
    "DEFAULT_MIC_HEIGHT",
    "DEFAULT_SPEED_OF_SOUND",
    "PART1_FREQUENCY_RANGE",
    "SPOT_FREQUENCY_RANGE",
    "SPOT_NARROW_BAND_RANGE",
    "RoadAbsorptionWarning",
    # ISO 8041-1 / ISO 2631 / ISO 5349 / Directive 2002/44/EC human vibration
    "frequency_weighting",
    "weighting_factors",
    "apply_weighting",
    "WeightingResponse",
    "WEIGHTING_NAMES",
    "weighted_acceleration",
    "WeightedSpectrum",
    "running_rms",
    "mtvv",
    "vibration_dose_value",
    "motion_sickness_dose_value",
    "crest_factor",
    "vibration_total_value",
    "daily_exposure",
    "partial_exposure",
    "combine_partial_exposures",
    "hav_daily_exposure",
    "energy_equivalent_acceleration",
    "hav_vwf_lifetime_years",
    "exposure_assessment",
    "ExposureAssessment",
    "daily_vibration_exposure",
    "DailyVibrationExposure",
    "HumanVibrationWarning",
    "REFERENCE_ACCELERATION",
    "REFERENCE_DURATION_S",
    "HAV_EAV_A8",
    "HAV_ELV_A8",
    "WBV_EAV_A8",
    "WBV_ELV_A8",
    "WBV_EAV_VDV",
    "WBV_ELV_VDV",
    "airborne_insulation",
    "AirborneInsulationResult",
    "impact_insulation",
    "ImpactInsulationResult",
    "weighted_rating",
    "WeightedRatingResult",
    "weighted_impact_rating",
    "weighted_impact_improvement",
    # ISO 16251-1 floor-covering impact sound improvement
    "FloorCoveringImprovementResult",
    "acceleration_level",
    "background_corrected_level",
    "impact_improvement",
    "improvement_octave_bands",
    "ImpactRatingResult",
    "facade_insulation",
    "FacadeInsulationResult",
    "energy_average_level",
    "lab_airborne_insulation",
    "LabAirborneInsulationResult",
    "lab_impact_insulation",
    "LabImpactInsulationResult",
    "background_correction",
    "LabInsulationWarning",
    "intensity_sound_reduction",
    "IntensityReductionResult",
    "intensity_element_normalized_difference",
    "IntensityElementNormalizedResult",
    "adaptation_term_kc",
    "surface_pressure_intensity_indicator",
    "combine_subareas",
    "reverberation_index",
    "estimate_reverberation_index",
    "survey_airborne_insulation",
    "SurveyAirborneResult",
    "survey_impact_insulation",
    "SurveyImpactResult",
    "survey_facade_insulation",
    "SurveyFacadeResult",
    "survey_service_equipment_level",
    "SurveyServiceEquipmentResult",
    "lden",
    "ldn",
    "composite_rating_level",
    "critical_bandwidth",
    "tonal_audibility",
    "tonal_adjustment",
    "assess_tonal_audibility",
    "tonal_adjustment_from_mean_audibility",
    "tonal_seeking_survey",
    "residual_sound_correction",
    "gaussian_residual_level",
    "combined_standard_uncertainty",
    "environmental_expanded_uncertainty",
    "residual_correction_uncertainty",
    "uncertainty_from_repeated_measurements",
    "TonalAssessmentResult",
    "ResidualCorrectionResult",
    "RepeatedMeasurementResult",
    "CalibrationWarning",
    "verify_filter_class",
    "class_limits",
    "verify_weighting_class",
    "weighting_class_limits",
    "predicted_airborne_insulation",
    "AirbornePredictionResult",
    "predicted_impact_insulation",
    "ImpactPredictionResult",
    "junction_vibration_reduction",
    "junction_min_vibration_reduction",
    "flanking_path",
    "flanking_element",
    "FlankingPath",
    "PathContribution",
    "combine_linings",
    "equivalent_impact_level",
    "impact_flanking_correction",
    "standardized_impact_level",
    "facade_sound_reduction",
    "radiated_sound_power",
    "outdoor_attenuation",
    "outdoor_level",
    "FacadeElement",
    "FacadePredictionResult",
    "RadiatedPowerResult",
    "velocity_level_difference",
    "direction_averaged_level_difference",
    "total_loss_factor",
    "equivalent_absorption_length",
    "vibration_reduction_index",
    "vibration_reduction_index_from_flanking",
    "normalized_flanking_level_difference",
    "normalized_flanking_impact_level",
    "critical_frequency",
    "strong_coupling_satisfied",
    "modal_density",
    "modal_overlap_factor",
    "band_mode_count",
    "VibrationReductionResult",
    "FlankingLevelDifferenceResult",
    "FlankingImpactLevelResult",
    "band_uncertainty",
    "single_number_uncertainty",
    "single_number_uncertainty_uncorrelated",
    "maximum_repeatability_standard_deviation",
    "coverage_factor",
    "expanded_uncertainty",
    "insulation_coverage_factor",
    "insulation_expanded_uncertainty",
    "uncertain_value",
    "combine_uncertainties",
    "prediction_input_uncertainty",
    "reduce_by_independent_measurements",
    "satisfies_lower_requirement",
    "satisfies_upper_requirement",
    "BandUncertainty",
    "UncertainValue",
    "COVERAGE_FACTORS",
    # ISO 12999-2 absorption uncertainty
    "AbsorptionUncertaintyResult",
    "absorption_coverage_factor",
    "sound_absorption_coefficient_uncertainty",
    "equivalent_area_uncertainty",
    "practical_coefficient_uncertainty",
    "weighted_coefficient_uncertainty",
    "single_number_rating_uncertainty",
    "task_based_exposure",
    "job_based_exposure",
    "full_day_exposure",
    "minimum_cumulative_duration_hours",
    "table_c4_contribution",
    "Task",
    "TaskContribution",
    "ExposureResult",
    "OccupationalExposureWarning",
    "COVERAGE_FACTOR",
    "INSTRUMENT_U2",
]


@lru_cache(maxsize=32)
def _cached_filter_bank(
    fs: int,
    fraction: float,
    order: int,
    limits: Tuple[float, ...] | None,
    filter_type: str,
    ripple: float,
    attenuation: float,
    calibration_factor: float,
    dbfs: bool,
) -> OctaveFilterBank:
    """Design (or reuse) a stateless filter bank for octave_filter()."""
    return OctaveFilterBank(
        fs=fs,
        fraction=fraction,
        order=order,
        limits=list(limits) if limits is not None else None,
        filter_type=filter_type,
        ripple=ripple,
        attenuation=attenuation,
        calibration_factor=calibration_factor,
        dbfs=dbfs,
    )


@overload
def octave_filter(
    x: List[float] | np.ndarray,  # NOSONAR - public API
    fs: int,
    fraction: float = 1,
    order: int = 6,
    limits: List[float] | None = None,
    show: bool = False,
    sigbands: Literal[False] = False,
    plot_file: str | None = None,
    detrend: bool = True,
    filter_type: str = "butter",
    ripple: float = 0.1,
    attenuation: float = 72.0,
    calibration_factor: float = 1.0,
    dbfs: bool = False,
    mode: str = "rms",
    nominal: Literal[False] = False,
) -> Tuple[np.ndarray, List[float]]: ...


@overload
def octave_filter(
    x: List[float] | np.ndarray,  # NOSONAR - public API
    fs: int,
    fraction: float = 1,
    order: int = 6,
    limits: List[float] | None = None,
    show: bool = False,
    sigbands: Literal[True] = True,
    plot_file: str | None = None,
    detrend: bool = True,
    filter_type: str = "butter",
    ripple: float = 0.1,
    attenuation: float = 72.0,
    calibration_factor: float = 1.0,
    dbfs: bool = False,
    mode: str = "rms",
    nominal: Literal[False] = False,
) -> Tuple[np.ndarray, List[float], List[np.ndarray]]: ...


@overload
def octave_filter(
    x: List[float] | np.ndarray,  # NOSONAR - public API
    fs: int,
    fraction: float = 1,
    order: int = 6,
    limits: List[float] | None = None,
    show: bool = False,
    sigbands: Literal[False] = False,
    plot_file: str | None = None,
    detrend: bool = True,
    filter_type: str = "butter",
    ripple: float = 0.1,
    attenuation: float = 72.0,
    calibration_factor: float = 1.0,
    dbfs: bool = False,
    mode: str = "rms",
    nominal: Literal[True] = ...,
) -> Tuple[np.ndarray, List[str]]: ...


@overload
def octave_filter(
    x: List[float] | np.ndarray,  # NOSONAR - public API
    fs: int,
    fraction: float = 1,
    order: int = 6,
    limits: List[float] | None = None,
    show: bool = False,
    sigbands: Literal[True] = True,
    plot_file: str | None = None,
    detrend: bool = True,
    filter_type: str = "butter",
    ripple: float = 0.1,
    attenuation: float = 72.0,
    calibration_factor: float = 1.0,
    dbfs: bool = False,
    mode: str = "rms",
    nominal: Literal[True] = ...,
) -> Tuple[np.ndarray, List[str], List[np.ndarray]]: ...


def octave_filter(
    x: List[float] | np.ndarray,  # NOSONAR - public API
    fs: int,
    fraction: float = 1,
    order: int = 6,
    limits: List[float] | None = None,
    show: bool = False,
    sigbands: bool = False,
    plot_file: str | None = None,
    detrend: bool = True,
    filter_type: str = "butter",
    ripple: float = 0.1,
    attenuation: float = 72.0,
    calibration_factor: float = 1.0,
    dbfs: bool = False,
    mode: str = "rms",
    nominal: bool = False,
) -> Tuple[np.ndarray, List[float]] | Tuple[np.ndarray, List[str]] | Tuple[np.ndarray, List[float], List[np.ndarray]] | Tuple[np.ndarray, List[str], List[np.ndarray]]:
    """
    Filter a signal with octave or fractional octave filter bank.

    This method uses a filter bank with Second-Order Sections (SOS) coefficients.
    To obtain the correct coefficients, automatic subsampling is applied to the
    signal in each filtered band.

    Multichannel support: If x is 2D (channels, samples), each channel is filtered.

    :param x: Input signal (1D array or 2D array [channels, samples]).
    :type x: Union[List[float], np.ndarray]
    :param fs: Sample rate in Hz.
    :type fs: int
    :param fraction: Bandwidth 'b'. Examples: 1/3-octave b=3, 1-octave b=1, 2/3-octave b=1.5. Default: 1.
    :type fraction: float
    :param order: Order of the filter. Default: 6.
    :type order: int
    :param limits: Minimum and maximum limit frequencies [f_min, f_max]. Default [12, 20000].
    :type limits: Optional[List[float]]
    :param show: If True, plot and show the filter response.
    :type show: bool
    :param sigbands: If True, also return the signal in the time domain divided into bands.
    :type sigbands: bool
    :param plot_file: Path to save the filter response plot.
    :type plot_file: Optional[str]
    :param detrend: If True, remove DC offset before filtering. Default: True.
    :type detrend: bool
    :param filter_type: Type of filter ('butter', 'cheby1', 'cheby2', 'ellip', 'bessel').
        Default: 'butter' (the only type that meets IEC 61260-1 class 1 with the
        default parameters).
    :param ripple: Passband ripple in dB (for cheby1, ellip). Default: 0.1.
    :param attenuation: Stopband attenuation in dB (for cheby2, ellip). Default: 72.0.
        For ``cheby2`` scipy pins the deep-stopband floor at exactly this value,
        so it must be >= 70 dB to clear the IEC 61260-1 class 1 limit (matches
        :class:`OctaveFilterBank`).
    :param calibration_factor: Calibration factor for SPL calculation. Default: 1.0.
    :param dbfs: If True, return results in dBFS. Default: False.
    :param mode: 'rms' or 'peak'. Default: 'rms'.
    :param nominal: If True, return IEC 61260-1 nominal frequency labels (List[str]) instead of exact floats.
    :return: A tuple containing (SPL_array, Frequencies_list) or (SPL_array, Frequencies_list, signals).
        When *nominal=True*, the frequency list contains ``List[str]`` labels instead of floats.
    :rtype: Union[Tuple[np.ndarray, List[float]], Tuple[np.ndarray, List[str]],
        Tuple[np.ndarray, List[float], List[np.ndarray]],
        Tuple[np.ndarray, List[str], List[np.ndarray]]]
    """
    
    if show or plot_file:
        # Plotting has side effects: bypass the cache.
        filter_bank = OctaveFilterBank(
            fs=fs,
            fraction=fraction,
            order=order,
            limits=limits,
            filter_type=filter_type,
            ripple=ripple,
            attenuation=attenuation,
            show=show,
            plot_file=plot_file,
            calibration_factor=calibration_factor,
            dbfs=dbfs,
        )
    else:
        # The bank is immutable in non-stateful mode: reuse the design.
        # Pass limits through as-is (tuple for hashability); the bank
        # constructor is the single place that validates them and owns
        # the default when None.
        limits_key = tuple(map(float, limits)) if limits is not None else None
        filter_bank = _cached_filter_bank(
            fs, fraction, order, limits_key, filter_type,
            ripple, attenuation, calibration_factor, dbfs,
        )

    return filter_bank.filter(x, sigbands=sigbands, mode=mode, detrend=detrend, nominal=nominal)  # type: ignore[call-overload,no-any-return]


def octavefilter(*args: Any, **kwargs: Any) -> Any:
    """Deprecated alias of :func:`octave_filter`."""
    _warn_renamed("octavefilter()", "octave_filter()")
    return octave_filter(*args, **kwargs)


#: Deprecated root-level name -> canonical name (phonometry 3.1 renames).
_RENAMED_ATTRIBUTES: dict[str, str] = {
    "OCTAVE_BANDS_HZ": "OCTAVE_BANDS",
    "THIRD_OCTAVE_BANDS_HZ": "THIRD_OCTAVE_BANDS",
    "BASE_PLATE_BANDS_HZ": "BASE_PLATE_BANDS",
    "ExposureWarning": "OccupationalExposureWarning",
}


def __getattr__(name: str) -> Any:
    """PEP 562 shim warning for names renamed in phonometry 3.1.

    Constants cannot warn through a wrapper, so the deprecated names live
    here (and in their home modules) as module ``__getattr__`` aliases.
    Remove in 4.0.
    """
    try:
        canonical = _RENAMED_ATTRIBUTES[name]
    except KeyError:
        raise AttributeError(
            f"module 'phonometry' has no attribute {name!r}"
        ) from None
    _warn_renamed(name, canonical)
    return globals()[canonical]
