spectrally positive K = one big atom + 3 residual atoms below m; theta = m
C(q5) rho^2 >= 0.00: max E|Z| - max(g(Lam), C) = +0.000000 at a = 1.3859 (w 1.0000), residual 0.0000@0.7205, 0.0000@0.7208, 0.0000@0.0144, m = 0.7215, Lam = 0.5206, E|Z| = 0.857409
Traceback (most recent call last):
  File "numerics/sp_idstar_scan.py", line 112, in <module>
    main()
  File "numerics/sp_idstar_scan.py", line 99, in main
    res = minimize(objective, z0, args=(J, C, rho0), method="Nelder-Mead",
          ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "site-packages/scipy/optimize/_minimize.py", line 772, in minimize
    res = _minimize_neldermead(fun, x0, args, callback, bounds=bounds,
          ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "site-packages/scipy/optimize/_optimize.py", line 877, in _minimize_neldermead
    fxr = func(xr)
          ^^^^^^^^
  File "site-packages/scipy/optimize/_optimize.py", line 561, in function_wrapper
    fx = function(np.copy(x), *(wrapper_args + args))
         ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "numerics/sp_idstar_scan.py", line 76, in objective
    a, wbig, xs, ws = unpack(z, J)
                      ^^^^^^^^^^^^
  File "numerics/sp_idstar_scan.py", line 63, in unpack
    e = np.exp(np.clip(z[1:J + 2], -30, 30))
               ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "site-packages/numpy/_core/fromnumeric.py", line 2334, in clip
    return _wrapfunc(a, 'clip', a_min, a_max, out=out, **kwargs)
           ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "site-packages/numpy/_core/fromnumeric.py", line 54, in _wrapfunc
    return bound(*args, **kwds)
           ^^^^^^^^^^^^^^^^^^^^
  File "site-packages/numpy/_core/_methods.py", line 96, in _clip
    def _clip(a, min=None, max=None, out=None, **kwargs):
    
KeyboardInterrupt
