def test_axisym_adjoint(
px2wav, wav2px, rng: np.random.Generator, L=10, B=2, J_min=2
):
nwvlts = pys2let_j_max(B, L, J_min) - J_min + 1
x = random_lms(rng, L)
y_scal, y_wav = random_wavlet_lms(rng, L, nwvlts)
y = wav2px(y_wav, y_scal, B, L, J_min)
x_wav, x_scal = px2wav(x, B, L, J_min)
# y'Ax
yAx = y_wav.conj().T @ x_wav + y_scal.conj() @ x_scal
# (A'y)'x
Ayx = approx(y.conj().T @ x)
assert yAx == Ayx
del finalmaphalf
return 0
if __name__ == "__main__":
##Input
nprocess = 12
nmaps = 9 #No. maps (WMAP = 5) (Planck = 9)
ellmax = 2048 #S2LET parameters - actually band-limits to 1 less
wavparam = 2
ndir = 2 #No. directions for each wavelet scale
spin = 0 #0 for temp, 1 for spin signals
upsample = 0 #0 for multiresolution, 1 for all scales at full resolution
jmin = 6
jmax = ps.pys2let_j_max(wavparam, ellmax, jmin)
fitsdir = '/home/keir/s2let_ilc_data/' #'/Users/keir/Documents/s2let_ilc/simu_data/'
fitsroot = 'planck_deconv_' #'simu_dirty_beam_wmap_9yr_' #'wmap_deconv_nosource_smoothw_extrapolated_9yr_'
fitscode = ['30', '44', '70', '100', '143', '217', '353', '545',
'857'] #['k','ka','q','v','w']
scal_fits = [None] * nmaps
wav_fits = [None] * nmaps
for i in xrange(len(scal_fits)):
scal_fits[i] = fitsdir + fitsroot + fitscode[i] + '_scal_' + str(
ellmax) + '_' + str(wavparam) + '_' + str(jmin) + '_' + str(
ndir) + '.npy'
wav_fits[i] = fitsdir + fitsroot + fitscode[i] + '_wav_' + str(
ellmax) + '_' + str(wavparam) + '_' + str(jmin) + '_' + str(
ndir) + '.npy'
#Increase line thicknesses
plt.rc('axes', linewidth=2.0)
plt.rc('xtick.major', width=2.0)
plt.rc('xtick.minor', width=2.0)
plt.rc('ytick.major', width=2.0)
plt.rc('ytick.minor', width=2.0)
plt.rc('lines', linewidth=1.5)
##Input
L = 128 #255 #500 #3600
J_min = 3 #0
#Bs = np.array([14,2]) #60,2]) #,1.3,1.2])
B = 2
#l_transitions = np.array([15]) #61]) #,513,2017])
J = ps.pys2let_j_max(B,L,J_min)
print J
N = 1 #3
nside = 64 #256
spin = 0
upsample = 1
#Construct valid hybrid tiling
'''hybrid_scal_l, hybrid_wav_l, hybrid_scal_bandlimit, hybrid_wav_bandlimits, J, L_bounds = ps.construct_hybrid_tiling(L,J_min,Bs,l_transitions)
res = ps.verify_tiling(L, hybrid_scal_l, hybrid_wav_l.T.ravel(), hybrid_scal_bandlimit, hybrid_wav_bandlimits)
if res == False:
raise ValueError('\nAn invalid wavelet tiling has been chosen.\n')
else:
print '\nA valid wavelet tiling has been chosen.\n'
'''