diff --git a/handoff/NOTES_LEVEL.md b/handoff/NOTES_LEVEL.md index 4d7eeda..0113ff1 100644 --- a/handoff/NOTES_LEVEL.md +++ b/handoff/NOTES_LEVEL.md @@ -2746,3 +2746,61 @@ Our application: mask × audio via complex multiply (WITH 1/N from IFFT) 1. Исправить mask computation (заменить LUT на raw_level или калибровать) 2. Реализовать FIR convolution (OLA) вместо per-bin multiply 3. Оба исправления требуются для bit-exact + +## ============ UPDATE 2026-08-24j (24j): ПРИМЕНЁННЫЙ ОТКЛИК = ЖИВОЙ FIR ПОТЕНЦИАЛЬНО; ×1.805 = ФАНТОМ STALE-КЕРНЕЛА ============ + +### Метод (новый канал, скрипты rendersnap2.py + генератор зондов) +1. МЯГКИЙ снаппер `scripts/rendersnap2.py`: НЕ убивает reaper (частичный-wav + hazard 24c обойдён — рендер доестественного выхода), дополнительно дампит + банк скаляров ctx+0x540860..0x5408c0 каждый фазой. +2. Мультитон-проект: клоны dual_b1q_0.5.rpp с входом «2 главных тона 500+2000 + @0.4248 (= амплитуды dual.wav) + 20 зондов −30 дБ на сетке 100..7200 Гц». + Детектор в том же адаптивном состоянии, H измерима Goertzel'ем на всех + зондах. Генерация inline (/tmp/opencode/multi_in.wav + multi.rpp). +3. Свип НЕ годится для Y/X-отношения (энергия на гребне) — только мультитон. + +### Результат multi_b1q_0.5 (живой FIR ph069: m@43=0.320, m@171=0.267) +| f, Гц | cut дБ | gain | m(FIR) | gain/m | +|-------|--------|------|--------|--------| +| 500 | 10.32 | 0.3048 | 0.3234 | 0.943 | +| 2000 | 11.84 | 0.2558 | 0.2698 | 0.948 | +| 700..7200 (18 зондов) | 1.9–4.6 | — | — | **0.98–1.00** | +(450/550 аномальны — утечка боковых лепестков окна Гёрцеля от главного +тона 500; не в счёт.) Контрольные числа воспроизведены: 10.32/11.84 ✓✓. + +### ВЫВОДЫ (три, все жёсткие) +1. **Применённый отклик = mag(живого FIR@540668) побиновно** (±2%; главные + тона −0.5 дБ — недозахват глубины между фазами). НИКАКОЙ дополнительной + свёртки/OLA/×1.805 поверх per-bin multiply НЕ СУЩЕСТВУЕТ. +2. **«Свёртка ×1.805» (23e) и «OLA-нормализация» (24i) ОТЗЫВАЮТСЯ**: источник — + сравнение cut живого рендера с init-time кернелом. Старые phase*.npz из + rendersnap-сессий держат m@43=0.588 СЕКУНДАМИ — это ядро после param-rebuild + при ЗАГРУЗКЕ (52e9b0 строит FIR из init-state bands[]); живой рендер + перестраивает FIR КАЖДЫЙ кадр и мгновенно потребляет → STOP-семпл ловит + identity/reset (все 110 фаз sweep-рендера ≈ 1.0). 0.305/0.32≈0.95, + 0.524(init)/0.32(live)=1.64~«1.805» — вот вся «константа». +3. **R-слоты 0x5407f8/0x540788 — СТАТИКА**: r43=3.9807/r43=4.0285 ИДЕНТИЧНЫ во + всех трёх проектах (dual/sweep/multi, разный контент). Это twin-шаблон + полосы (геометрия), НЕ динамическая маска. Интерпретация 22y «R=1/mask + живая кривая редукции» — снята (кривая была от инициализации). +4. **Скаляр [ctx+0x540888]=1.000000** весь рендер (и s874/s87c/s88c=1.0; + s870=440.955=sens-scale). Wet-скаляр BLOCKMAP 23b НЕ вносит масштаба при + дефолтных параметрах. ×N-нормировка FFT (24g) тоже не подтверждается: + применённый gain совпадает с mag(FIR) без поправок. + +### Что это меняет в архитектуре транскрипции +- RT_FIRCONV как «особое применение» НЕ нужен: достаточно per-bin multiply + спектра кадра на маску-кривую (как default-путь spectral.cpp). +- ВСЕ усилия — в ФОРМУ кривой: bands[i] → (log→ops A-D→exp0.984→hannfall→ops) + → FIR. Наш structural chain должен выдать кривую, чья mag совпадёт с живым + FIR фазы. Эталонные кривые теперь ЕСТЬ: rendersnap2/ph*.npz (fir_re/fir_im). + +### NEXT (24k+) +1. Дискретизация формулы: снять ph-серии с ОДИНАКОВЫм статичным контентом при + вариации одного параметра (q-серия уже есть рефами; прогнать rendersnap2 по + dual_b1q_0.1..10.0) → восстановить зависимость m(fir) от параметров. +2. Сравнить живой FIR с нашим structural mask: где именно расходится форма + (хвост? ширина? пол?) — починка detector-цепи по живым эталонам. +3. Проверить hannfall-мультипликацию: fit exp(0.984·log(bands))·w против + захваченного fir_im≠0 (фаза!) — complex ops A-D дают мнимость. +4. Починить wav-path баг rendersnap2 (проверяет rpp.wav вместо RENDER_FILE). diff --git a/scripts/rendersnap2.py b/scripts/rendersnap2.py new file mode 100644 index 0000000..f35413f --- /dev/null +++ b/scripts/rendersnap2.py @@ -0,0 +1,182 @@ +#!/usr/bin/env python3 +"""rendersnap2.py — non-destructive variant of rendersnap.py. + +Differences (24j): +- does NOT kill the host: lets -renderproject finish naturally so the output + wav is complete (partial-wav hazard, NOTES 24c); +- additionally dumps the ctx SCALAR bank 0x540860..0x5408c0 every phase + (goal: wet/FIR scalar [ctx+0x540888], BLOCKMAP 23b line 52b832); +- keeps sampling until host exit or phase cap, dedupe by FIR+scalar md5. +""" +import hashlib +import os +import signal +import struct +import subprocess +import sys +import time + +import numpy as np + +SLOTS = [0x540668] +NARR = 8194 +SCAL_OFF = 0x540860 +SCAL_N = 24 # floats -> 0x540860..0x5408c0 +OUT = '/tmp/opencode/rendersnap2' + + +def find_host(): + import glob + for p in glob.glob('/proc/[0-9]*'): + pid = int(os.path.basename(p)) + try: + cmd = open(f'/proc/{pid}/cmdline', 'rb').read().replace(b'\0', b' ').decode('utf8', 'replace') + maps = open(f'/proc/{pid}/maps').read() + except Exception: + continue + if 'soothe2' in maps and 'reaper' not in cmd: + return pid + return None + + +def main(): + rpp = sys.argv[1] if len(sys.argv) > 1 else '/home/m/soothe-bt/dual_b1q_0.5.rpp' + cap = int(sys.argv[2]) if len(sys.argv) > 2 else 200 + os.makedirs(OUT, exist_ok=True) + subprocess.run("pkill -9 -x reaper; pkill -9 -f '[y]abridge'; " + "rm -rf /run/user/1000/yabridge-soothe2_x64-*; sleep 1", shell=True) + wav = rpp.replace('.rpp', '.wav') + if os.path.exists(wav): + os.remove(wav) + proc = subprocess.Popen(['/usr/bin/reaper', '-nosplash', '-ignoreerrors', + '-renderproject', rpp], + stdout=open('/dev/null', 'w'), stderr=subprocess.STDOUT) + t0 = time.time() + host = None + while time.time() - t0 < 30 and not host: + host = find_host() + time.sleep(0.001) + if not host: + print('NO HOST') + return 1 + print('host %d at %.3fs' % (host, time.time() - t0), flush=True) + + fd = os.open(f'/proc/{host}/mem', os.O_RDONLY) + + def rd(a, n): + try: + return os.pread(fd, n, a) + except OSError: + return None + + vt = struct.pack('= 0: + cand = a + j + sb = rd(cand + 0x540870, 4) + if sb and struct.unpack(' 100: + return cand + j = d.find(vt, j + 1) + j = d.find(m48) + while j >= 0: + cand = a + j - 0x24 + sb = rd(cand + 0x540870, 4) + if sb and struct.unpack(' 100: + return cand + j = d.find(m48, j + 1) + a += CH + return None + + ctx = None + while ctx is None and time.time() - t0 < 25: + try: + os.kill(host, signal.SIGSTOP) + except ProcessLookupError: + break + try: + ctx = scan_ctx() + finally: + if ctx is None: + try: + os.kill(host, signal.SIGCONT) + except ProcessLookupError: + break + if ctx is None: + time.sleep(0.004) + if not ctx: + print('NO CTX') + return 1 + print('ctx %#x' % ctx, flush=True) + + rng = np.random.default_rng(11) + prev = None + saved = 0 + while True: + try: + os.kill(host, signal.SIGSTOP) + except ProcessLookupError: + print('host exited at %.2fs' % (time.time() - t0), flush=True) + break + try: + pb = rd(ctx + 0x540668, 8) + if not pb: + continue + p = struct.unpack('= cap: + break + finally: + try: + os.kill(host, signal.SIGCONT) + except ProcessLookupError: + pass + time.sleep(float(rng.uniform(0.0005, 0.006))) + try: + os.kill(host, 0) + except ProcessLookupError: + print('host exited at %.2fs' % (time.time() - t0), flush=True) + break + # wait for natural finish so wav completes + for _ in range(300): + if proc.poll() is not None: + break + time.sleep(0.1) + print('saved=%d reaper_rc=%s wav=%s' % (saved, proc.poll(), + os.path.getsize(wav) if os.path.exists(wav) else 'NONE')) + return 0 + + +if __name__ == '__main__': + sys.exit(main())