Home > src > main > matlab > flaglet_axisym_tiling.m

flaglet_axisym_tiling

PURPOSE ^

flaglet_axisym_tiling - Compute tiling in l-n harmonic space.

SYNOPSIS ^

function [kappa kappa0] = flaglet_axisym_tiling(B_l, B_p, L, P, J_min_l, J_min_p)

DESCRIPTION ^

 flaglet_axisym_tiling - Compute tiling in l-n harmonic space.
 -- Axisymmetric wavelets on the solid sphere.

 Default usage :

   [kappa kappa0] = flaglet_axisym_tiling(B_l, B_n, L, N, J_min_l, J_min_n)

 kappa is an array containing wavelet tiling functions,
 kappa0 contains the scaling function,
 B_l is the wavelet parameter for angular space,
 B_n is the wavelet parameter for radial space,
 L is the angular band-limit,
 N is the radial band-limit,
 J_min_l the first angular wavelet scale to use,
 J_min_n the first radial wavelet scale to use.

 B3LET package to perform Wavelets transform on the solid Sphere.
 Copyright (C) 2012  Boris Leistedt & Jason McEwen
 See LICENSE.txt for license details

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function [kappa kappa0] = flaglet_axisym_tiling(B_l, B_p, L, P, J_min_l, J_min_p)
0002 
0003 % flaglet_axisym_tiling - Compute tiling in l-n harmonic space.
0004 % -- Axisymmetric wavelets on the solid sphere.
0005 %
0006 % Default usage :
0007 %
0008 %   [kappa kappa0] = flaglet_axisym_tiling(B_l, B_n, L, N, J_min_l, J_min_n)
0009 %
0010 % kappa is an array containing wavelet tiling functions,
0011 % kappa0 contains the scaling function,
0012 % B_l is the wavelet parameter for angular space,
0013 % B_n is the wavelet parameter for radial space,
0014 % L is the angular band-limit,
0015 % N is the radial band-limit,
0016 % J_min_l the first angular wavelet scale to use,
0017 % J_min_n the first radial wavelet scale to use.
0018 %
0019 % B3LET package to perform Wavelets transform on the solid Sphere.
0020 % Copyright (C) 2012  Boris Leistedt & Jason McEwen
0021 % See LICENSE.txt for license details
0022 
0023 p = inputParser;
0024 p.addRequired('B_l', @isnumeric);       
0025 p.addRequired('B_p', @isnumeric);          
0026 p.addRequired('L', @isnumeric);   
0027 p.addRequired('P', @isnumeric);    
0028 p.addRequired('J_min_l', @isnumeric);   
0029 p.addRequired('J_min_p', @isnumeric); 
0030 p.parse(B_l, B_p, L, P, J_min_l, J_min_p);
0031 args = p.Results;
0032 
0033 [kappa_vec kappa0] = flaglet_axisym_tiling_mex(B_l, B_p, L, P, J_min_l, J_min_p);
0034 
0035 J_l = ceil(log(L) ./ log(B_l));
0036 J_p = ceil(log(P) ./ log(B_p));
0037 
0038 kappa = cell(J_l + 1, J_p + 1);
0039 for jl = 0:J_l
0040     for jp = 0:J_p
0041         temp = zeros(P,L);
0042         for l = 0:L-1
0043             for p = 0:P-1
0044                 temp(p+1,l+1) = kappa_vec(1,jp*(J_l+1)*L*P  + jl*L*P + p*L + l + 1);
0045             end
0046         end
0047         kappa{jl+1,jp+1} = temp;
0048     end   
0049 end
0050 
0051 end

Generated on Mon 24-Sep-2012 12:26:33 by m2html © 2005