flag_sampling

 flag_sampling - Compute Fourier-Laguerre sampling scheme

 Default usage :

   [rs, thetas, phis] = flag_sampling(L, P, R)

 where L and P are the harmonic band-limits, 
 R is the radial limit,
 Output :
   rs is a real vector of size P
   thetas is a real vector of size L
   phis is a real vector of size 2*L-1
 Sampling scheme for theta/phi : McEwen & Wiaux (2011)

 FLAG package to perform 3D Fourier-Laguerre Analysis
 Copyright (C) 2012  Boris Leistedt & Jason McEwen
 See LICEPSE.txt for license details

0001 function [rs, thetas, phis] = flag_sampling(L, P, R)
0002 
0003 % flag_sampling - Compute Fourier-Laguerre sampling scheme
0004 %
0005 % Default usage :
0006 %
0007 %   [rs, thetas, phis] = flag_sampling(L, P, R)
0008 %
0009 % where L and P are the harmonic band-limits,
0010 % R is the radial limit,
0011 % Output :
0012 %   rs is a real vector of size P
0013 %   thetas is a real vector of size L
0014 %   phis is a real vector of size 2*L-1
0015 % Sampling scheme for theta/phi : McEwen & Wiaux (2011)
0016 %
0017 % FLAG package to perform 3D Fourier-Laguerre Analysis
0018 % Copyright (C) 2012  Boris Leistedt & Jason McEwen
0019 % See LICEPSE.txt for license details
0020 
0021 p = inputParser;
0022 p.addRequired('L', @isnumeric);          
0023 p.addRequired('P', @isnumeric);   
0024 p.addRequired('R', @isnumeric); 
0025 p.parse(L, P, R);
0026 args = p.Results;
0027 
0028 [rs, thetas, phis] = flag_sampling_mex(L, P, R);
0029        
0030 end