regintpointsro.m

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

function points=regpolintpointsRO(h) % generate regular interior data points
% with approximate fill distance h
n=ceil(sqrt(2)/h);
dh=2/n;
[X,Y]=meshgrid(-1:dh:1,-1:dh:1);
[nx,ny]=size(X);
xr=reshape(X,[(nx)*(nx) 1]);
yr=reshape(Y,[(nx)*(nx) 1]);
zr=domainfct(xr,yr);
intpoints=[xr yr];     
I=find(zr<-0.001);
points=[xr(I) yr(I)];