rbflaplace.m

本代码为基于matlab的求解偏微分方程的利用径向基函数的无网格方法的源码

function y =  rbfLaplace (rr)   % Laplace of rbf, see rbflap2.mws
global RBFscale;
global RBFpar;
global RBFtype;   % 1=MQ, 2=G, 3=IMQ, 4=W2, 5=TPS, 6=W4
r=rr/RBFscale;
if RBFtype==1
    y=(-(2+r.*r).*(1+r.*r).^(-3/2))/RBFscale^2;             % Multiquadric
    return;
end
if RBFtype==2
    y=(4*(r.*r-1).*exp(-r.*r))/RBFscale^2;                  % Gaussian
    return;
end
if RBFtype==3
    y=((r.*r-2).*(1+r.*r).^(-5/2))/RBFscale^2;            % inverse Multiquadric
    return;
end
if RBFtype==4
  y=(20*(5*r-2).*max(0,(1-r)).^2)/RBFscale^2;      % Wendland 2
   return;
end
if RBFtype==5
    y=(4*(1+log(r+eps)))/RBFscale^2;                % TPS
    return;
end
if RBFtype==6
    y=(112*(20*r.^2-4*r-1).*max(0,(1-r)).^4)/RBFscale^2;     % Wendland C4
    return;
end
if RBFtype==7
    y=2*RBFpar.^2./((r.^2+RBFpar.^2).^2.*RBFscale^2);      % regularized log = AFS of Laplace in R^2
    return;
end

error('RBF type not implemented')