clear all;
close all;
clear functions;
clc;
disp('---------------------------------------------');
disp(sprintf('Example to demonstrate the use of Gauss Transform'));
disp('---------------------------------------------');
d=2;
disp(sprintf('Dimensionality d=%d\n',d));
N=100;
disp(sprintf('Number of source points N=%d\n',N));
M=100;
disp(sprintf('Number of target points M=%d\n',N));
X=rand(d,N);
Y=rand(d,M);
q=rand(1,N);
h=sqrt(2)*0.2*sqrt(d);
disp(sprintf('Bandwidth h=%f\n',h));
disp('---------------------------------------------');
disp(sprintf('Running the Gauss Transform.\n'));
disp('---------------------------------------------');
to=clock;
[G_direct]=GaussTransform(d,N,M,X,h,q,Y);
GT_time=etime(clock,to);
disp(sprintf('Time taken=%f secs\n',GT_time));
---------------------------------------------
Example to demonstrate the use of Gauss Transform
---------------------------------------------
Dimensionality d=2
Number of source points N=100
Number of target points M=100
Bandwidth h=0.400000
---------------------------------------------
Running the Gauss Transform.
---------------------------------------------
Time taken=0.000000 secs