Function details for flagRFI

flagRFI_std>padSmooth (12 calls, 0.645 sec)
Generated 05-Aug-2011 13:00:40 using cpu time.
subfunction in file /home/LeechJ/cbass_analysis/reduc/flagRFI_std.m
Copy to new window for comparing multiple runs

Parents (calling functions)

Function Name	Function Type	Calls
flagRFI_std	function	12

Lines where the most time was spent

Line Number	Code	Calls	Total Time	% Time	Time Plot
562	dataVecT = interp1(xVec(flagVe...	12	0.404 s	62.7%
590	dataVecFilt = real(ifft(fft(da...	12	0.164 s	25.4%
583	dataVecPad(1:padParams(2)) = ....	12	0.044 s	6.8%
586	dataVecPad(end-padParams(3)+1:...	12	0.022 s	3.4%
566	dataVecF = dataVec .* (1 - fla...	12	0.011 s	1.7%
All other lines			0 s	0%
Totals			0.645 s	100%

Children (called functions)

Function Name	Function Type	Calls	Total Time	% Time	Time Plot
interp1	function	12	0.262 s	40.7%
smooth	function	12	0.120 s	18.6%
repmat	function	24	0.033 s	5.1%
median	function	24	0.022 s	3.4%
substruct	function	2	0 s	0%
mean	function	2	0 s	0%
Self time (built-ins, overhead, etc.)			0.208 s	32.2%
Totals			0.645 s	100%

Code Analyzer results

Line number

Message

Coverage results
[ Show coverage for parent directory ]

Total lines in function	37
Non-code lines (comments, blank lines)	20
Code lines (lines that can run)	17
Code lines that did run	14
Code lines that did not run	3
Coverage (did run/can run)	82.35 %

Function listing

   time   calls  line
                558 function smoothVec = padSmooth(xVec,dataVec,kernelVec,flagVec,padParams)
                559 
                560 
                561 % Interpolate over the flagged regions.
  0.40      12  562 dataVecT = interp1(xVec(flagVec == 0), ... 
                563     smooth(dataVec(flagVec == 0),ceil(length(dataVec) / 100)), ...
                564     (1:length(dataVec)), 'pchip');
                565 
  0.01      12  566 dataVecF = dataVec .* (1 - flagVec) + flagVec .* dataVecT; 
                567 
            12  568 if flagVec(1) == 1 
             2  569     noFlagPos = find(flagVec == 0,1,'first'); 
                570     
             2  571     dataVecF(1:noFlagPos-1) = mean(dataVecF(noFlagPos:noFlagPos+99)); 
             2  572 end     
            12  573 if flagVec(end) == 1 
                574     noFlagPos = find(flagVec == 0,1,'last');
                575     
                576     dataVecF(noFlagPos+1:end) = mean(dataVecF(noFlagPos-99:noFlagPos));
                577 end    
                578     
                579 
                580 
                581 % Pad the data up to a power of two.
            12  582 dataVecPad = zeros(padParams(1),1); 
  0.04      12  583 dataVecPad(1:padParams(2)) = ... 
                584     repmat(median(dataVecF(1:10)),padParams(2),1);
            12  585 dataVecPad(padParams(2)+1:padParams(2)+length(dataVec)) = dataVecF; 
  0.02      12  586 dataVecPad(end-padParams(3)+1:end) = ... 
                587     repmat(median(dataVecF(end-9:end)),padParams(3),1);
                588 
                589 % Perform a Fourier filtering.
  0.16      12  590 dataVecFilt = real(ifft(fft(dataVecPad).*fft(kernelVec))); 
                591     
                592 % And then clip off the padding.
            12  593 smoothVec = dataVecFilt(padParams(2)+1:end - padParams(3)); 
            12  594 end

Other subfunctions in this file are not included in this listing.