
P = 3;
sigma = 2;
precision = 'single';
g = gaussian_2d_pyramid(2*P+1, 2*P+1, sigma, precision);
