s2let_check_cur_tiling - Checks the admissibity condition is satisfied. ----------------------------------------------------------- S2LET package to perform wavelets transform on the Sphere. Copyright (C) 2012-2016 Boris Leistedt, Jennifer Chan & Jason McEwen See LICENSE.txt for license details -----------------------------------------------------------
0001 function error_on_cur_tiling = s2let_check_cur_tiling(cur_lm, scal_l, L, Spin, J, J_min) 0002 % s2let_check_cur_tiling - Checks the admissibity condition is satisfied. 0003 % ----------------------------------------------------------- 0004 % S2LET package to perform wavelets transform on the Sphere. 0005 % Copyright (C) 2012-2016 Boris Leistedt, Jennifer Chan & Jason McEwen 0006 % See LICENSE.txt for license details 0007 % ----------------------------------------------------------- 0008 0009 % Scaling function 0010 identity = zeros(1,L); 0011 for el=abs(Spin):L-1 0012 identity(1,el+1) = identity(1,el+1)+4.*pi/(2*el+1)*scal_l(el^2+el+1,1)*conj(scal_l(el^2+el+1,1)); 0013 end 0014 0015 0016 % Curvelet functions 0017 for j= J_min: J 0018 ind = Spin*Spin + 1; 0019 for el=abs(Spin):L-1 0020 for m= -el:el 0021 identity(1,el+1) = identity(1,el+1)+(8.*pi^2/(2*el+1))* cur_lm{j-J_min+1}(ind)*conj(cur_lm{j-J_min+1}(ind)); 0022 ind = ind + 1; 0023 end 0024 end 0025 end 0026 0027 error_on_cur_tiling = 0; 0028 for el=abs(Spin):L-1 0029 error_on_cur_tiling = error_on_cur_tiling+(identity(1,el+1) - 1.0); 0030 end 0031 0032 0033 end 0034