sigma[current["key"]] = current["sigma"]
current = move_one_dyadic_step(current, Q, alpha=alpha)
// get the last key
last_key = max(xi.keys())n = last_key[1] + 1
last_xi = xi[last_key]
// fill num_interm wavelets between last_xi and 0, both excluded
num_intermediate = Q - 1
for q in range(1, num_intermediate + 1):
After Change
xi_max = compute_xi_max(Q)
sigma_max = compute_sigma_psi(xi_max, Q, r=r_psi)
xi = []
sigma = []
j = []
if sigma_max <= sigma_low: