24j: applied response == live FIR per-bin (x1.805 revoked as stale-kernel phantom); R-slots static twin template; rendersnap2 soft snapper + probe-multitone method
This commit is contained in:
@@ -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).
|
||||
|
||||
@@ -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('<Q', 0x1824AC210)
|
||||
m48 = struct.pack('<I', 0x473b8000)
|
||||
|
||||
def scan_ctx():
|
||||
for line in open(f'/proc/{host}/maps'):
|
||||
parts = line.split()
|
||||
if 'rw' not in parts[1]:
|
||||
continue
|
||||
lo, hi = (int(x, 16) for x in parts[0].split('-'))
|
||||
CH = 16 * 1024 * 1024
|
||||
a = lo
|
||||
while a < hi:
|
||||
d = rd(a, min(CH + 4096, hi - a))
|
||||
if not d:
|
||||
break
|
||||
j = d.find(vt)
|
||||
while j >= 0:
|
||||
cand = a + j
|
||||
sb = rd(cand + 0x540870, 4)
|
||||
if sb and struct.unpack('<f', sb)[0] > 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('<f', sb)[0] > 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('<Q', pb)[0]
|
||||
fb = rd(p, NARR * 4)
|
||||
if not fb:
|
||||
continue
|
||||
scal = rd(ctx + SCAL_OFF, SCAL_N * 4)
|
||||
if scal is None:
|
||||
continue
|
||||
key = hashlib.md5(fb[:2049 * 8] + scal).digest()
|
||||
if key != prev:
|
||||
prev = key
|
||||
arr = np.frombuffer(fb[:2049 * 8], dtype='<f4')
|
||||
mag = np.hypot(arr[0::2], arr[1::2])
|
||||
sv = np.frombuffer(scal, dtype='<f4')
|
||||
np.savez_compressed(f'{OUT}/ph{saved:03d}.npz',
|
||||
fir_re=arr[0::2].copy(), fir_im=arr[1::2].copy(),
|
||||
scal=sv.copy())
|
||||
print('PH%03d t=%.2f fir43=%.4f fir171=%.4f | s888=%.6f s88c=%.6f s874=%.6f s87c=%.6f s870=%.3f' %
|
||||
(saved, time.time() - t0, mag[43], mag[171],
|
||||
sv[10], sv[11], sv[5], sv[7], sv[4]), flush=True)
|
||||
saved += 1
|
||||
if saved >= 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())
|
||||
Reference in New Issue
Block a user