0001 function map_sm = s2let_smoothing(map, L, fwhm, varargin) 0002 0003 p = inputParser; 0004 p.addParamValue('Reality', false, @islogical); 0005 p.addParamValue('Method', 'MW', @ischar); 0006 p.parse(varargin{:}); 0007 args = p.Results; 0008 0009 map_lm = ssht_forward(map, L, 'Reality', true, 'Method', args.Method); 0010 map_sm_lm = zeros(size(map_lm)); 0011 sigma = fwhm / 2.355; 0012 0013 ind = 1; 0014 for el = 0:L-1 0015 for m = -el:el 0016 map_sm_lm(ind) = map_lm(ind) * exp(-el*(el+1) * sigma^2 / 2); 0017 ind = ind+1; 0018 end 0019 end 0020 0021 map_sm = ssht_inverse(map_sm_lm, L, 'Reality', args.Reality, 'Method', args.Method);