[交流] 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 ]

回复此楼

» 收录本帖的淘帖专辑推荐

数学物理计算

» 猜你喜欢

• 上海工程技术大学张培磊教授团队招收博士生 已经有3人回复
• 上海工程技术大学【激光智能制造】课题组招收硕士 已经有5人回复
• 求助院士们，这个如何合成呀 已经有4人回复
• 临港实验室与上科大联培博士招生1名 已经有9人回复
• 想换工作。大多数高校都是 评职称时 认可5年内在原单位取得的成果吗？ 已经有7人回复
• 需要合成515-64-0，50g，能接单的留言 已经有4人回复
• 自荐读博 已经有4人回复
• 写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？ 已经有6人回复
• 带资进组求博导收留 已经有10人回复
• 最近几年招的学生写论文不引自己组发的文章 已经有11人回复

高级回复