bndfct.m

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

function r=bndfct(phi) % This defines domains by specifying the
                       % boundary radius r in polar coordinates.
                       % No other routine is needed for domain definition.
                       % Make sure that the domain is in [-1,1]
                       % and that it is star-shaped wrt. the origin.
global RBFdomain
if RBFdomain==1      % circle
    r=1.0;
    return;
end
if RBFdomain==2   % the square [-1,1]^2
    r=1.0./max(abs(cos(phi)),abs(sin(phi)));
    return; 
end
if RBFdomain==3   % cardioid
    r=0.4*(1.0+cos(phi));
    return; 
end
if RBFdomain==4   % five- pointed star
    r=0.3*(1.5+cos(5.0*phi));
    return; 
end
if RBFdomain==5   % five- pointed star
    r=0.3*(1.5+cos(2.0*phi));
    return; 
end
if RBFdomain==6   % five- pointed star
    r=0.3*(1.2+cos(6.0*phi));
    return; 
end

error('Undefined domain')