-- status: library-dev — CORRECTED root cause (deep research): EP fits consume GLOBAL np.random via LaplaceOptimiser.n_refine=3 → refine_state → mean_field.sample → NormalMessage.sample (normal.py:274,278). Fit's fixed point (good vs sigma-collapse) depends on ambient np.random state at fit start; test seeds only before data gen, not the fit. #1351 added test_ep_statistics_fixes.py (Monte-Carlo KL tests) which shifts ambient RNG into a failing region. PROVEN: reseed-1 before fit = PASS 5/5; sweep pre-fit seed 0..7 = PASS×7, seed7 FAIL. Factor ordering ALREADY deterministic (my earlier hypothesis wrong). Fix candidates A(seed refine sampler)/B(deterministic refinement)/C(test seed) — awaiting user steer; approved "factor-schedule" approach SUPERSEDED
0 commit comments