time   calls  line
                  1 function C = convn(A,B,shape)
                  2 %CONVN  N-dimensional convolution.
                  3 %   C = CONVN(A, B) performs the N-dimensional convolution of
                  4 %   matrices A and B. If nak = size(A,k) and nbk = size(B,k), then
                  5 %   size(C,k) = max([nak+nbk-1,nak,nbk]);
                  6 %
                  7 %   C = CONVN(A, B, 'shape') controls the size of the answer C:
                  8 %     'full'   - (default) returns the full N-D convolution
                  9 %     'same'   - returns the central part of the convolution that
                 10 %                is the same size as A.
                 11 %     'valid'  - returns only the part of the result that can be
                 12 %                computed without assuming zero-padded arrays.
                 13 %                size(C,k) = max([nak-max(0,nbk-1)],0).
                 14 %
                 15 %   Class support for inputs A,B:
                 16 %      float: double, single
                 17 %
                 18 %   See also CONV, CONV2.
                 19 
                 20 %   Copyright 1984-2009 The MathWorks, Inc. 
                 21 %   $Revision: 1.10.4.6 $  $Date: 2009/10/24 19:17:38 $
                 22 
             1   23 error(nargchk(2,3,nargin,'struct')); 
                 24 
             1   25 if (nargin < 3) 
                 26   shape = 'full';
                 27 end
                 28 
             1   29 if ~isfloat(A) 
                 30     A = double(A);
                 31 end
             1   32 if ~isfloat(B) 
                 33     B = double(B);
                 34 end
                 35 
  0.02       1   36 C = convnc(A,B,shape);