File size: 2,158 Bytes
debe354 | 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 | """Tests for Signal-to-Noise Ratio analysis."""
import numpy as np
import math
import sys
sys.path.insert(0, str(__file__).rsplit("/tests/", 1)[0])
from src.sovereign_shift import Q, N_ACTIVE, THETA, THETA_NUM
from src.metasum import compute, hallucination_bound
def test_snr_sweep():
"""SNR remains positive across hallucination sweep."""
rng = np.random.default_rng(42)
print(" Hallucination Sweep:")
for n_halluc in [0, 10, 100, 500, 1000]:
weights = np.zeros(Q)
active = rng.choice(Q, size=N_ACTIVE, replace=False)
weights[active] = 1.0
if n_halluc > 0:
inactive = np.setdiff1d(np.arange(Q), active)
halluc = rng.choice(inactive, size=min(n_halluc, len(inactive)), replace=False)
weights[halluc] = 1.0
displacements = np.arange(Q, dtype=float)
ms_mag = abs(compute(weights, displacements))
bound = hallucination_bound(n_halluc) if n_halluc > 0 else 0
print(f" n_halluc={n_halluc:5d}: |MetaSum|={ms_mag:.2f}, Weyl={bound:.2f}")
print(" PASS: SNR sweep complete")
def test_89_is_weyl_bound():
"""The Weyl bound √(N·log Q) ≈ 89 = Sovereign Shift numerator."""
bound = hallucination_bound(N_ACTIVE)
# Should be within 5% of 89
ratio = bound / THETA_NUM
assert 0.95 < ratio < 1.10, f"Weyl bound = {bound:.2f}, 89 * ratio = {ratio:.3f}"
print(f" PASS: Weyl bound = {bound:.2f}, ratio to 89 = {ratio:.3f}")
def test_impossible_hallucination_dominance():
"""Verify >134,000 agents needed to dominate (impossible with Q=2462)."""
# |halluc| > N² / log(Q) needed to dominate
needed = N_ACTIVE**2 / math.log(Q)
max_possible = Q - N_ACTIVE # 1438
assert max_possible < needed, \
f"max_possible={max_possible} should be < needed={needed:.0f}"
print(f" PASS: need {needed:.0f} halluc agents, max possible = {max_possible} (IMPOSSIBLE)")
if __name__ == "__main__":
print("SNR Tests:")
test_snr_sweep()
test_89_is_weyl_bound()
test_impossible_hallucination_dominance()
print("ALL PASS")
|