Home > src > main > matlab > s2let_check_cur_tiling.m

s2let_check_cur_tiling

PURPOSE ^

s2let_check_cur_tiling - Checks the admissibity condition is satisfied.

SYNOPSIS ^

function error_on_cur_tiling = s2let_check_cur_tiling(cur_lm, scal_l, L, Spin, J, J_min)

DESCRIPTION ^

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

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

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

Generated on Fri 11-Nov-2016 11:50:36 by m2html © 2005