| 查看: 3343 | 回复: 9 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[交流]
2D凸包Delaunay三角剖分程序 已有3人参与
|
||||
|
之前用Python写的网格程序，速度实在太慢。没有办法，花了几天时间用C++重写了一次，虽然还有很多瑕疵，不过用用已经可以了，一段实现如下剖分的代码和效果： 资料地址：http://persson.berkeley.edu/distmesh/ 程序代码： #include "nr.h" //简易向量矩阵类 #include "distmesh_2d.h" //参考persson的论文写的基于力平衡的网格程序 #include "delaunay.h" //随机增量算法Delaunay三角剖分（效果一般） #include "dfunction.h" //一般区域的符号距离函数 struct Ftor{ vector<Point<2>> poly; Ftor(vector<Point<2>> &polyy):poly(polyy){} Doub fd(const Point<2> &p){ Doub d1=dcircle(p,0.0,0.0,0.5); //圆的符号距离函数 Doub d2=dpoly(p,poly); //多边形的符号距离函数 return ddiff(d2,d1); //区域做差（挖空圆） } Doub fh(const Point<2> &p){ return 0.05+0.3*dcircle(p,0.0,0.0,0.5); //密度函数，内圆处加密 } }; int main(){ ofstream fp("p.txt"  ; //输出点分布ofstream ft("t.txt" ; //输出三角形节点编号Doub pi=3.1415926; Box<2> box; //包围区域初始化剖分点 box.lo.x[0]=-1.0;box.lo.x[1]=-1.0; box.hi.x[0]=1.0;box.hi.x[1]=1.0; vector<Point<2>> pfix(6); //生成多边形+固定剖分点 for(int i=0;i<6;i++){ //正六变形生成 Doub theta=pi/3.0; pfix.x[0]=cos(theta*i); pfix.x[1]=sin(theta*i); } Ftor Fun(pfix);; //距离函数和密度函数初始化 Distmesh2d<Ftor> mesh(Fun,0.1,box,pfix); //调用网格程序 //以下输出，Matlab数组是1偏的，三角形编号+1处理 fp<<fixed<<setprecision(4); ft<<fixed<<setprecision(6); int n=mesh.p.size(); for(int i=0;i<n;i++){ fp<<setw(4)<<mesh.p.x[0]<<setw(10)<<mesh.p.x[1]<<endl; } int k=mesh.tri.size(); for(int i=0;i<k;i++){ ft<<setw(10)<<mesh.tri.p[0]+1 <<setw(10)<<mesh.tri.p[1]+1 <<setw(10)<<mesh.tri.p[2]+1<<endl; } fp.close(); ft.close(); cin.get(); return 0; }  密度控制.jpg  1.jpg  2.jpg [ Last edited by nagami on 2015-1-16 at 21:30 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 数学物理计算 |
» 猜你喜欢
FileCode能看出啥？
已经有23人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
关于Filecode分析方法
已经有9人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有10人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
pippi6
铁杆木虫 (著名写手)
工程和科学数值计算咨询
- 数学EPI: 6
- 应助: 413 (硕士)
- 贵宾: 0.002
- 金币: 7116.5
- 散金: 15
- 红花: 63
- 帖子: 1639
- 在线: 798.9小时
- 虫号: 2469437
- 注册: 2013-05-14
- 专业: 计算数学与科学工程计算
7楼2015-01-10 10:52:56
2楼2015-01-08 20:30:02
3楼2015-01-08 21:57:52
人走茶不凉
银虫 (职业作家)
- 应助: 9 (幼儿园)
- 金币: 13053.8
- 红花: 3
- 帖子: 3748
- 在线: 708.9小时
- 虫号: 667594
- 注册: 2008-12-04
- 性别: GG
- 专业: 计算机应用技术
4楼2015-01-09 00:45:44
