Skip to content
Snippets Groups Projects
Commit d9091914 authored by Christophe Geuzaine's avatar Christophe Geuzaine
Browse files

fix some of the demos + clarify some comments
parent c4b176e7
No related branches found
No related tags found
No related merge requests found
// $Id: meshGFaceExtruded.cpp,v 1.21 2007-07-31 22:09:11 geuzaine Exp $ // $Id: meshGFaceExtruded.cpp,v 1.22 2007-08-02 16:16:19 geuzaine Exp $
// //
// Copyright (C) 1997-2007 C. Geuzaine, J.-F. Remacle // Copyright (C) 1997-2007 C. Geuzaine, J.-F. Remacle
// //
...@@ -193,7 +193,7 @@ int MeshExtrudedSurface(GFace *gf, ...@@ -193,7 +193,7 @@ int MeshExtrudedSurface(GFace *gf,
Msg(STATUS2, "Meshing surface %d (extruded)", gf->tag()); Msg(STATUS2, "Meshing surface %d (extruded)", gf->tag());
// build a set with all the vertices on the boundary of gf // build a set with all the vertices on the boundary of the face gf
double old_tol = MVertexLessThanLexicographic::tolerance; double old_tol = MVertexLessThanLexicographic::tolerance;
MVertexLessThanLexicographic::tolerance = 1.e-12 * CTX.lc; MVertexLessThanLexicographic::tolerance = 1.e-12 * CTX.lc;
...@@ -209,8 +209,8 @@ int MeshExtrudedSurface(GFace *gf, ...@@ -209,8 +209,8 @@ int MeshExtrudedSurface(GFace *gf,
++it; ++it;
} }
// if the edges are constrained, the vertices already exist on the // if the edges of the mesh are constrained, the vertices already
// face--so we add them to the set // exist on the face--so we add them to the set
if(constrainedEdges) if(constrainedEdges)
pos.insert(gf->mesh_vertices.begin(), gf->mesh_vertices.end()); pos.insert(gf->mesh_vertices.begin(), gf->mesh_vertices.end());
......
// $Id: meshGRegion.cpp,v 1.31 2007-05-24 17:34:04 geuzaine Exp $ // $Id: meshGRegion.cpp,v 1.32 2007-08-02 16:16:19 geuzaine Exp $
// //
// Copyright (C) 1997-2007 C. Geuzaine, J.-F. Remacle // Copyright (C) 1997-2007 C. Geuzaine, J.-F. Remacle
// //
...@@ -192,7 +192,7 @@ void MeshDelaunayVolume(std::vector<GRegion*> &regions) ...@@ -192,7 +192,7 @@ void MeshDelaunayVolume(std::vector<GRegion*> &regions)
#else #else
for(unsigned int i = 0; i < regions.size(); i++) for(unsigned int i = 0; i < regions.size(); i++)
Msg(STATUS2, "Meshing volume %d (Delaunay)", regions[i]->tag()); Msg(STATUS2, "Meshing volume %d (Tetgen+Delaunay)", regions[i]->tag());
// put all the faces in the same model // put all the faces in the same model
GRegion *gr = regions[0]; GRegion *gr = regions[0];
......
...@@ -227,8 +227,8 @@ Volume(vAir) = {news-2,news-1,surftotprobe};//Rest Air ...@@ -227,8 +227,8 @@ Volume(vAir) = {news-2,news-1,surftotprobe};//Rest Air
Characteristic Length{10,11,28,32}= pflaw ;//Watch out! Changing characteristic lenths of the flaw at the end Characteristic Length{10,11,28,32}= pflaw ;//Watch out! Changing characteristic lenths of the flaw at the end
//Characteristic Length{66,93}= ph3 ; //Characteristic Length{66,93}= ph3 ;
//Transfinite Line{cd[0],cd[1],cd[2]} = 5 ; Transfinite Line{cd[0],cd[1],cd[2]} = 5 ;
//Transfinite Line{15,6} = nDefect+1 ; Transfinite Line{15,6} = nDefect+1 ;
//------------------------------------------------------------------- //-------------------------------------------------------------------
......
...@@ -38,28 +38,28 @@ Line Loop(14) = {5,6,7,8,9,10,11,12}; ...@@ -38,28 +38,28 @@ Line Loop(14) = {5,6,7,8,9,10,11,12};
Plane Surface(15) = {13,14}; Plane Surface(15) = {13,14};
Extrude Line {8, {0,-1,0}, {0,0,0}, Pi/2}{ Extrude Line {8, {0,-1,0}, {0,0,0}, Pi/2}{
Layers {{R*Pi/2/mafil},{1},{1.0}}; Layers {{R*Pi/2/mafil},{1.0}};
}; };
Extrude Line {9, {0,-1,0}, {0,0,0}, Pi/2}{ Extrude Line {9, {0,-1,0}, {0,0,0}, Pi/2}{
Layers {{R*Pi/2/mafil},{1},{1.0}}; Layers {{R*Pi/2/mafil},{1.0}};
}; };
Extrude Line {10, {0,-1,0}, {0,0,0}, Pi/2}{ Extrude Line {10, {0,-1,0}, {0,0,0}, Pi/2}{
Layers {{R*Pi/2/mafil},{1},{1.0}}; Layers {{R*Pi/2/mafil},{1.0}};
}; };
Extrude Line {7, {0,-1,0}, {0,0,0}, Pi/2}{ Extrude Line {7, {0,-1,0}, {0,0,0}, Pi/2}{
Layers {{R*Pi/2/mafil},{1},{1.0}}; Layers {{R*Pi/2/mafil},{1.0}};
}; };
Extrude Line {6, {0,-1,0}, {0,0,0}, Pi/2}{ Extrude Line {6, {0,-1,0}, {0,0,0}, Pi/2}{
Layers {{R*Pi/2/mafil},{1},{1.0}}; Layers {{R*Pi/2/mafil},{1.0}};
}; };
Extrude Line {5, {0,-1,0}, {0,0,0}, Pi/2}{ Extrude Line {5, {0,-1,0}, {0,0,0}, Pi/2}{
Layers {{R*Pi/2/mafil},{1},{1.0}}; Layers {{R*Pi/2/mafil},{1.0}};
}; };
Extrude Line {12, {0,-1,0}, {0,0,0}, Pi/2}{ Extrude Line {12, {0,-1,0}, {0,0,0}, Pi/2}{
Layers {{R*Pi/2/mafil},{1},{1.0}}; Layers {{R*Pi/2/mafil},{1.0}};
}; };
Extrude Line {11, {0,-1,0}, {0,0,0}, Pi/2}{ Extrude Line {11, {0,-1,0}, {0,0,0}, Pi/2}{
Layers {{R*Pi/2/mafil},{1},{1.0}}; Layers {{R*Pi/2/mafil},{1.0}};
}; };
Point(51) = {0,-wsubstrat/2,0,maboite}; Point(51) = {0,-wsubstrat/2,0,maboite};
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment