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' '''