function A=smoothG(B,s)

% function performing two-dimensional Gaussian smoothing with 
% sigma=s and a 2D circular Gaussian kernel

n1=ceil(3*s);sigma1=s;n2=n1;sigma2=sigma1;
filter1=d2gauss0(n1,sigma1,n2,sigma2);
A=conv2(B,filter1,'same');
