diff --git a/CMakeLists.txt b/CMakeLists.txt
index 8f2d2a0e0742f290e5229fcf376d49be249c6366..384ff887ee269a99394871d70fb3932efb5414c7 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -261,6 +261,31 @@ macro(set_config_option VARNAME STRING)
message(STATUS "Found " ${STRING})
endmacro(set_config_option)
+
+if (ENABLE_PARSER)
+ find_package(BISON)
+ find_package(FLEX)
+
+ if(BISON_FOUND AND FLEX_FOUND)
+ bison_target(ParserBison ${CMAKE_CURRENT_SOURCE_DIR}/Parser/Gmsh.y
+ ${CMAKE_CURRENT_BINARY_DIR}/Parser/Gmsh.tab.cpp
+ COMPILE_FLAGS -pgmsh_yy)
+ flex_target(ParserFlex ${CMAKE_CURRENT_SOURCE_DIR}/Parser/Gmsh.l
+ ${CMAKE_CURRENT_BINARY_DIR}/Parser/Gmsh.yy.cpp
+ COMPILE_FLAGS -Pgmsh_yy)
+ include_directories(${CMAKE_CURRENT_BINARY_DIR})
+ add_flex_bison_dependency(ParserFlex ParserBison)
+ set(GMSH_SRC ${GMSH_SRC} ${BISON_ParserBison_OUTPUTS} ${FLEX_ParserFlex_OUTPUTS})
+ message(STATUS "Gmsh sources are ${GMSH_SRC}")
+ else()
+ message(STATUS "Could not find bison or flex")
+ endif(BISON_FOUND AND FLEX_FOUND)
+
+ add_subdirectory(Parser)
+ set_config_option(HAVE_PARSER "Parser")
+endif(ENABLE_PARSER)
+
+
# check if the machine is 64 bits (this is more reliable than using
# CMAKE_SIZEOF_VOID_P, which does not seem to work e.g. on some Suse
# machines)
@@ -445,11 +470,6 @@ if(ENABLE_POST)
endif(ENABLE_PLUGINS)
endif(ENABLE_POST)
-if(ENABLE_PARSER)
- add_subdirectory(Parser)
- set_config_option(HAVE_PARSER "Parser")
-endif(ENABLE_PARSER)
-
if(ENABLE_FLTK)
# first, try to use fltk-config for fltk >= 1.3 (FindFLTK is buggy on Unix,
# where e.g. xft and xinerama options are not dealt with)
@@ -1479,16 +1499,6 @@ elseif(MSVC)
"/STACK:16777216 /SAFESEH:NO")
endif(WIN32 AND NOT MSVC OR CYGWIN)
-# parser target
-find_program(BISON bison)
-find_program(FLEX flex)
-if(BISON AND FLEX)
- add_custom_target(parser
- COMMAND ${BISON} -p gmsh_yy --output Gmsh.tab.cpp -d Gmsh.y
- COMMAND ${FLEX} -P gmsh_yy -o Gmsh.yy.cpp Gmsh.l
- WORKING_DIRECTORY ${CMAKE_CURRENT_SOURCE_DIR}/Parser)
-endif(BISON AND FLEX)
-
if(UNIX)
# cannot use cmake's file search functions here (they would only find files
# existing at configuration time)
diff --git a/Geo/GEdge.cpp b/Geo/GEdge.cpp
index 80604e799a1bcdeee813606e18d82778d5d1a1d0..271f460147d97c612b50620b9de01a935684b7c9 100644
--- a/Geo/GEdge.cpp
+++ b/Geo/GEdge.cpp
@@ -45,6 +45,79 @@ void GEdge::deleteMesh()
model()->destroyMeshCaches();
}
+void GEdge::setMeshMaster(GEdge* ge,int ori) {
+
+#warning missing:computation of affine transformation during setMeshMaster
+
+ GEntity::setMeshMaster(ge);
+ masterOrientation = ori > 0 ? 1:-1;
+
+ if (ori < 0) {
+ vertexCounterparts[getBeginVertex()] = ge->getEndVertex();
+ vertexCounterparts[getEndVertex()] = ge->getBeginVertex();
+ getBeginVertex()->setMeshMaster(ge->getEndVertex());
+ getEndVertex() ->setMeshMaster(ge->getBeginVertex());
+ }
+ else {
+ vertexCounterparts[getBeginVertex()] = ge->getBeginVertex();
+ vertexCounterparts[getEndVertex()] = ge->getEndVertex();
+ getBeginVertex()->setMeshMaster(ge->getBeginVertex());
+ getEndVertex() ->setMeshMaster(ge->getEndVertex());
+ }
+}
+
+void GEdge::setMeshMaster(GEdge* ge,const std::vector<double>& tfo) {
+
+ SPoint3 oriXYZ0 = ge->getBeginVertex()->xyz();
+ SPoint3 oriXYZ1 = ge->getEndVertex()->xyz();
+
+ SPoint3 tfoXYZ0(0,0,0);
+ SPoint3 tfoXYZ1(0,0,0);
+
+ int idx = 0;
+ for (int i=0;i<3;i++,idx++) {
+ for (int j=0;j<3;j++,idx++) {
+ tfoXYZ0[i] += tfo[idx] * oriXYZ0[j];
+ tfoXYZ1[i] += tfo[idx] * oriXYZ1[j];
+ }
+ tfoXYZ0[i] += tfo[idx];
+ tfoXYZ1[i] += tfo[idx];
+ }
+
+ SPoint3 locXYZ0 = getBeginVertex()->xyz();
+ SPoint3 locXYZ1 = getEndVertex()->xyz();
+
+ SVector3 d00 = locXYZ0 - tfoXYZ0;
+ SVector3 d10 = locXYZ1 - tfoXYZ0;
+ SVector3 d01 = locXYZ0 - tfoXYZ1;
+ SVector3 d11 = locXYZ1 - tfoXYZ1;
+
+ double tol = CTX::instance()->geom.tolerance;
+
+ if ((d00.norm() < tol) && (d11.norm() < tol)) {
+ GEntity::setMeshMaster(ge,tfo);
+ masterOrientation = 1;
+ vertexCounterparts[getBeginVertex()] = ge->getBeginVertex();
+ vertexCounterparts[getEndVertex()] = ge->getEndVertex();
+ getBeginVertex()->setMeshMaster(ge->getBeginVertex(),tfo);
+ getEndVertex() ->setMeshMaster(ge->getEndVertex() ,tfo);
+ return;
+ }
+ if ((d01.norm() < tol) && (d10.norm() < tol)) {
+ GEntity::setMeshMaster(ge,tfo);
+ masterOrientation = -1;
+ vertexCounterparts[getBeginVertex()] = ge->getEndVertex();
+ vertexCounterparts[getEndVertex()] = ge->getBeginVertex();
+ getBeginVertex()->setMeshMaster(ge->getEndVertex() ,tfo);
+ getEndVertex() ->setMeshMaster(ge->getBeginVertex(),tfo);
+ return;
+ }
+
+ Msg::Error("Transformation from edge %d (%d-%d) to %d (%d-%d) is incorrect",
+ ge->tag(),ge->getBeginVertex()->tag(),ge->getEndVertex()->tag(),
+ this->tag(),this->getBeginVertex()->tag(),this->getEndVertex()->tag());
+}
+
void GEdge::reverse()
{
GVertex* tmp = v0;
diff --git a/Geo/GEdge.h b/Geo/GEdge.h
index 5260301d8859d92abc17e8d7209b6deea64e7458..8aee20ffba78633595542b0106f64e78fb4d0088 100644
--- a/Geo/GEdge.h
+++ b/Geo/GEdge.h
@@ -26,6 +26,7 @@ class closestPointFinder;
// A model edge.
class GEdge : public GEntity {
+
private:
double _length;
bool _tooSmall;
@@ -38,9 +39,9 @@ class GEdge : public GEntity {
GEdgeCompound *compound; // this model edge belongs to a compound
std::list<GFace *> l_faces;
// for specific solid modelers that need to re-do the internal curve
- // if a topological change ending points is done (gluing)
+ // if a topological change ending points is done (glueing)
virtual void replaceEndingPointsInternals(GVertex *, GVertex *) {}
-
+
public:
GEdge(GModel *model, int tag, GVertex *_v0, GVertex *_v1);
virtual ~GEdge();
@@ -52,6 +53,15 @@ class GEdge : public GEntity {
virtual GVertex *getBeginVertex() const { return v0; }
virtual GVertex *getEndVertex() const { return v1; }
+ // same or opposite direction to the master
+ int masterOrientation;
+
+ // specify mesh master with transformation, deduce edgeCounterparts
+ virtual void setMeshMaster(GEdge* master,const std::vector<double>&);
+
+ // specify mesh master and edgeCounterparts, deduce transformation
+ virtual void setMeshMaster(GEdge* master,int sign);
+
void reverse();
// add/delete a face bounded by this edge
diff --git a/Geo/GEntity.cpp b/Geo/GEntity.cpp
index 6b898c0da651a92a98157363fddf5c053dd9cfc4..392a32f6eec5fdfb9fa448c30ae8303ee093be05 100644
--- a/Geo/GEntity.cpp
+++ b/Geo/GEntity.cpp
@@ -14,8 +14,8 @@
#include "GFace.h"
#include "GRegion.h"
-GEntity::GEntity(GModel *m, int t)
- : _model(m), _tag(t), _meshMaster(t), _visible(1), _selection(0),
+GEntity::GEntity(GModel *m,int t)
+ : _model(m), _tag(t),_meshMaster(this),_visible(1), _selection(0),
_allElementsVisible(1), _obb(0), va_lines(0), va_triangles(0)
{
_color = CTX::instance()->packColor(0, 0, 255, 0);
@@ -70,56 +70,40 @@ GFace *GEntity::cast2Face() { return dynamic_cast<GFace*>(this); }
GRegion *GEntity::cast2Region() { return dynamic_cast<GRegion*>(this); }
// sets the entity m from which the mesh will be copied
-void GEntity::setMeshMaster(int m_signed)
+void GEntity::setMeshMaster(GEntity* gMaster)
{
- if(m_signed == tag()){ _meshMaster = m_signed; return; }
-
- GEntity *gMaster = 0;
- int m = abs(m_signed);
- switch(dim()){
- case 0 : gMaster = model()->getVertexByTag(m); break;
- case 1 : gMaster = model()->getEdgeByTag(m); break;
- case 2 : gMaster = model()->getFaceByTag(m); break;
- case 3 : gMaster = model()->getRegionByTag(m); break;
- }
- if (!gMaster){
- Msg::Error("Model entity %d of dimension %d cannot be the mesh master of entity %d",
- m, dim(), tag());
+ if (gMaster->dim() != dim()){
+ Msg::Error("Model entity %d of dimension %d cannot"
+ "be the mesh master of entity %d of dimension %d",
+ gMaster->tag(),gMaster->dim(),tag(),dim());
return;
}
- int masterOfMaster = gMaster->meshMaster();
-
- if (masterOfMaster == gMaster->tag()){
- _meshMaster = m_signed;
- }
- else {
- setMeshMaster ( masterOfMaster * ((m_signed > 0) ? 1 : -1));
- }
+ _meshMaster = gMaster;
}
-// gets the entity from which the mesh will be copied
-int GEntity::meshMaster() const
+void GEntity::setMeshMaster(GEntity* gMaster,const std::vector<double>& tfo)
{
- if (_meshMaster == tag()) return tag();
-
- GEntity *gMaster = 0;
- switch(dim()){
- case 0 : gMaster = model()->getVertexByTag(abs(_meshMaster)); break;
- case 1 : gMaster = model()->getEdgeByTag(abs(_meshMaster)); break;
- case 2 : gMaster = model()->getFaceByTag(abs(_meshMaster)); break;
- case 3 : gMaster = model()->getRegionByTag(abs(_meshMaster)); break;
- }
- if (!gMaster){
- Msg::Error("Could not find mesh master entity %d",_meshMaster);
- return tag();
+ if (gMaster->dim() != dim()){
+ Msg::Error("Model entity %d of dimension %d cannot"
+ "be the mesh master of entity %d of dimension %d",
+ gMaster->tag(),gMaster->dim(),tag(),dim());
+ return;
}
- int masterOfMaster = gMaster->meshMaster();
- if (masterOfMaster == gMaster->tag()){
- return _meshMaster ;
- }
- else {
- return gMaster->meshMaster() * ((_meshMaster > 0) ? 1 : -1);
+ if (tfo.size() != 16) {
+ Msg::Error("Periodicity transformation from entity %d to %d (dim %d) has %d components"
+ ", while 16 are required",
+ gMaster->tag(),tag(),gMaster->dim(),tfo.size());
+ return;
}
+
+ affineTransform = tfo;
+ _meshMaster = gMaster;
}
+
+
+
+// gets the entity from which the mesh will be copied
+GEntity* GEntity::meshMaster() const {return _meshMaster;}
+
diff --git a/Geo/GEntity.h b/Geo/GEntity.h
index 9b935c0c6f5f11f702d920e3d1b6a64ebd2b3236..2146e85bc95fca0a8c28029e3405e47f82c92c89 100644
--- a/Geo/GEntity.h
+++ b/Geo/GEntity.h
@@ -34,7 +34,8 @@ class GEntity {
// gives the number of the master entity in periodic mesh, gives _tag
// if non-periodic
- int _meshMaster;
+ GEntity* _meshMaster;
+
// the visibility and the selection flag
char _visible, _selection;
@@ -47,11 +48,14 @@ class GEntity {
protected:
SOrientedBoundingBox *_obb;
-
+
public: // these will become protected at some point
// the mesh vertices uniquely owned by the entity
std::vector<MVertex*> mesh_vertices;
+ // corresponding principal vertices
+ std::map<GVertex*,GVertex*> vertexCounterparts;
+
// the physical entitites (if any) that contain this entity
std::vector<int> physicals;
@@ -167,7 +171,7 @@ class GEntity {
return name[type];
}
- GEntity(GModel *m, int t);
+ GEntity(GModel *m,int t);
virtual ~GEntity(){}
@@ -261,9 +265,11 @@ class GEntity {
return physicals;
}
- // returns the tag of the entity that its master entity (for mesh)
- int meshMaster() const;
- void setMeshMaster(int m);
+ // returns the master entity (for mesh)
+ GEntity* meshMaster() const;
+ virtual void setMeshMaster(GEntity*);
+
+ virtual void setMeshMaster(GEntity*,const std::vector<double>&);
// get the bounding box
virtual SBoundingBox3d bounds() const { return SBoundingBox3d(); }
@@ -315,7 +321,7 @@ class GEntity {
// get the number of mesh vertices in the entity
unsigned int getNumMeshVertices() { return (int)mesh_vertices.size(); }
-
+
// get the mesh vertex at the given index
MVertex *getMeshVertex(unsigned int index) { return mesh_vertices[index]; }
@@ -332,7 +338,7 @@ class GEntity {
GRegion *cast2Region();
// periodic data
- double periodicTransformation[4][4];
+ std::vector<double> affineTransform;
std::map<MVertex*,MVertex*> correspondingVertices;
};
diff --git a/Geo/GFace.cpp b/Geo/GFace.cpp
index 56755edc5f32ea05be4484c201bb0eb596ab8dd6..bf687799606ecb4b938b19e8197f1e0efe62ed9e 100644
--- a/Geo/GFace.cpp
+++ b/Geo/GFace.cpp
@@ -1479,4 +1479,404 @@ void GFace::relocateMeshVertices()
}
}
+void GFace::setMeshMaster(GFace* master,const std::vector<double>& tfo) {
+
+ std::list<GEdge*>::iterator eIter;
+ std::list<GVertex*>::iterator vIter;
+
+ // list all vertices and construct vertex to edge correspondence for local edge
+
+ std::set<GVertex*> l_vertices;
+ std::map<std::pair<GVertex*,GVertex*>,GEdge* > l_vtxToEdge;
+
+ for (eIter=l_edges.begin();eIter!=l_edges.end();++eIter) {
+ GVertex* v0 = (*eIter)->getBeginVertex();
+ GVertex* v1 = (*eIter)->getEndVertex();
+ l_vertices.insert(v0);
+ l_vertices.insert(v1);
+ l_vtxToEdge[std::make_pair(v0,v1)] = (*eIter);
+ }
+
+ for (eIter=embedded_edges.begin();eIter!=embedded_edges.end();++eIter) {
+ GVertex* v0 = (*eIter)->getBeginVertex();
+ GVertex* v1 = (*eIter)->getEndVertex();
+ l_vertices.insert(v0);
+ l_vertices.insert(v1);
+ l_vtxToEdge[std::make_pair(v0,v1)] = (*eIter);
+ }
+
+ l_vertices.insert(embedded_vertices.begin(),embedded_vertices.end());
+
+ // list all vertices and vertex to edge correspondence for remote edge
+
+ std::list<GEdge*> m_edges = master->edges();
+ std::set<GVertex*> m_vertices;
+ std::map<std::pair<GVertex*,GVertex*>,GEdge* > m_vtxToEdge;
+ for (eIter=m_edges.begin();eIter!=m_edges.end();++eIter) {
+ GVertex* v0 = (*eIter)->getBeginVertex();
+ GVertex* v1 = (*eIter)->getEndVertex();
+ m_vertices.insert(v0);
+ m_vertices.insert(v1);
+ m_vtxToEdge[std::make_pair(v0,v1)] = (*eIter);
+ }
+
+ std::list<GEdge*> m_embedded_edges = master->embeddedEdges();
+
+ for (eIter=m_embedded_edges.begin();eIter!=m_embedded_edges.end();eIter++) {
+ GVertex* v0 = (*eIter)->getBeginVertex();
+ GVertex* v1 = (*eIter)->getEndVertex();
+ m_vertices.insert(v0);
+ m_vertices.insert(v1);
+ m_vtxToEdge[std::make_pair(v0,v1)] = (*eIter);
+ }
+
+ std::list<GVertex*> m_embedded_vertices = master->embeddedVertices();
+ m_vertices.insert(m_embedded_vertices.begin(),m_embedded_vertices.end());
+
+ // check topological correspondence
+
+ if (l_vertices.size() != m_vertices.size()) {
+ Msg::Error("Periodic connection specified between topologically "
+ "incompatible surfaces %d and %d (that have %d vs %d model vertices)",
+ master->tag(),tag(),l_vertices.size(),m_vertices.size());
+ return;
+ }
+
+ if (l_vtxToEdge.size() != m_vtxToEdge.size()) {
+ Msg::Error("Periodic connection specified between topologically "
+ "incompatible surfaces %d and %d (that have %d vs %d model edges)",
+ master->tag(),tag(),l_vtxToEdge.size(),m_vtxToEdge.size());
+ return;
+ }
+
+ // compute corresponding vertices
+
+ std::map<GVertex*,GVertex*> gVertexCounterparts;
+
+ std::set<GVertex*>::iterator mvIter;
+ for (mvIter=m_vertices.begin();mvIter!=m_vertices.end();++mvIter) {
+
+ GVertex* m_vertex = *mvIter;
+
+ double xyzOri[4] = {m_vertex->x(),
+ m_vertex->y(),
+ m_vertex->z(),1};
+ SPoint3 xyzTfo(0,0,0);
+
+ for (size_t i=0,ij=0;i<3;i++) {
+ for (size_t j=0;j<4;j++,ij++) {
+ xyzTfo[i] += tfo[ij] * xyzOri[j];
+ }
+ }
+
+ GVertex* l_vertex = NULL;
+
+ std::set<GVertex*>::iterator lvIter = l_vertices.begin();
+ for (;lvIter!=l_vertices.end();++lvIter) {
+
+ SPoint3 xyz((*lvIter)->x(),(*lvIter)->y(),(*lvIter)->z());
+ SVector3 dist = xyz - xyzTfo;
+
+ if (dist.norm() < CTX::instance()->geom.tolerance) {
+ l_vertex = *lvIter;
+ break;
+ }
+ }
+
+ if (l_vertex==NULL) {
+ Msg::Error("Was not able to find corresponding node for %d "
+ "for periodic connection of surface %d to %d "
+ "using the specified transformation",
+ m_vertex->tag(),master->tag(),tag());
+ return;
+ }
+ gVertexCounterparts[l_vertex] = m_vertex;
+ }
+
+ if (gVertexCounterparts.size() != m_vertices.size()) {
+ Msg::Error("Could not find all node correspondances "
+ "for the periodic connection from surface %d to %d",
+ master->tag(),tag());
+ return;
+ }
+
+ // construct edge correspondence and update the edge masters
+
+ std::map<GEdge*,std::pair<GEdge*,int> > gEdgeCounterparts;
+
+ std::map<std::pair<GVertex*,GVertex*>,GEdge*>::iterator lv2eIter;
+ for (lv2eIter=l_vtxToEdge.begin();lv2eIter!=l_vtxToEdge.end();lv2eIter++) {
+
+ std::pair<GVertex*,GVertex*> lPair = lv2eIter->first;
+ GEdge* localEdge = lv2eIter->second;
+
+ std::pair<GVertex*,GVertex*> mPair(gVertexCounterparts[lPair.first],
+ gVertexCounterparts[lPair.second]);
+ int sign = 1;
+ std::map<std::pair<GVertex*,GVertex*>,GEdge*>::iterator mv2eIter = m_vtxToEdge.find(mPair);
+ if (mv2eIter == m_vtxToEdge.end()) {
+ sign *= -1;
+ std::pair<GVertex*,GVertex*> backward(mPair.second,mPair.first);
+ mv2eIter = m_vtxToEdge.find(backward);
+ }
+
+ if (mv2eIter == m_vtxToEdge.end()) {
+ Msg::Error("Could not find periodic copy of edge %d-%d "
+ "(corresponding to vertices %d %d) in face %d",
+ lPair.first->tag(),lPair.second->tag(),
+ mPair.first->tag(),mPair.second->tag(),
+ master->tag());
+ return;
+ }
+ GEdge* masterEdge = mv2eIter->second;
+ localEdge->setMeshMaster(masterEdge,tfo);
+ gEdgeCounterparts[localEdge] = std::make_pair(masterEdge,sign);
+ }
+
+ // complete the information at the edge level
+
+ edgeCounterparts = gEdgeCounterparts;
+ vertexCounterparts = gVertexCounterparts;
+ GEntity::setMeshMaster(master,tfo);
+}
+
+inline double myAngle(const SVector3 &a, const SVector3 &b, const SVector3 &d)
+{
+ double cosTheta = dot(a, b);
+ double sinTheta = dot(crossprod(a, b), d);
+ return atan2(sinTheta, cosTheta);
+}
+struct myPlane {
+ SPoint3 p;
+ SVector3 n;
+ double a;
+ // nx x + ny y + nz z + a = 0
+ myPlane(SPoint3 _p, SVector3 _n) : p(_p),n(_n)
+ {
+ n.normalize();
+ a = -(n.x()*p.x()+n.y()*p.y()+n.z()*p.z());
+ }
+ double eval (double x, double y, double z)
+ {
+ return n.x() * x + n.y() * y + n.z() * z + a;
+ }
+};
+struct myLine {
+ SPoint3 p;
+ SVector3 t;
+ myLine() : p(0,0,0) , t (0,0,1) {}
+ myLine(myPlane &p1, myPlane &p2)
+ {
+ t = crossprod(p1.n, p2.n);
+ if (t.norm() == 0.0){
+ Msg::Error("parallel planes do not intersect");
+ }
+ else
+ t.normalize();
+ // find a point, assume z = 0
+ double a[2][2] = {{p1.n.x(), p1.n.y()}, {p2.n.x(), p2.n.y()}};
+ double b[2] = {-p1.a, -p2.a}, x[2];
+ if (!sys2x2(a, b, x)){
+ // assume x = 0
+ double az[2][2] = {{p1.n.y(), p1.n.z()}, {p2.n.y(), p2.n.z()}};
+ double bz[2] = {-p1.a, -p2.a};
+ if (!sys2x2(az, bz, x)){
+ // assume y = 0
+ double ay[2][2] = {{p1.n.x(), p1.n.z()}, {p2.n.x(), p2.n.z()}};
+ double by[2] = {-p1.a, -p2.a};
+ if (!sys2x2(ay,by,x)){
+ Msg::Error("parallel planes do not intersect");
+ }
+ else {
+ p = SPoint3(x[0], 0., x[1]);
+ }
+ }
+ else{
+ p = SPoint3(0., x[0], x[1]);
+ }
+ }
+ else{
+ p = SPoint3(x[0], x[1], 0.);
+ }
+ }
+ SPoint3 orthogonalProjection (SPoint3 &a)
+ {
+ // (x - a) * t = 0 -->
+ // x = p + u t --> (p + ut - a) * t = 0 --> u = (a-p) * t
+ const double u = dot(a - p, t);
+ return SPoint3(p.x() + t.x() * u,p.y() + t.y() * u,p.z() + t.z() * u);
+ }
+};
+
+
+
+void GFace::setMeshMaster(GFace* master,const std::map<int,int>& edgeCopies) {
+
+ std::map<GVertex*,GVertex*> vs2vt;
+
+ for (std::list<GEdge*>::iterator it=l_edges.begin();it!=l_edges.end();++it){
+
+ // slave edge
+ GEdge* le = *it;
+
+ int sign = 1;
+ std::map<int,int>::const_iterator adnksd = edgeCopies.find(le->tag());
+ int source_e;
+ if (adnksd != edgeCopies.end()) source_e = adnksd->second;
+ else{
+ sign = -1;
+ adnksd = edgeCopies.find(-(*it)->tag());
+ if(adnksd != edgeCopies.end()) source_e = adnksd->second;
+ else{
+ Msg::Error("Could not find edge counterpart %d in slave surface %d",
+ (*it)->tag(), master->tag());
+ return;
+ }
+ }
+
+ // master edge
+ GEdge *me = master->model()->getEdgeByTag(abs(source_e));
+
+ if (source_e * sign > 0){
+ vs2vt[me->getBeginVertex()] = le->getBeginVertex();
+ vs2vt[me->getEndVertex()] = le->getEndVertex();
+ }
+ else {
+ vs2vt[me->getBeginVertex()] = le->getEndVertex();
+ vs2vt[me->getEndVertex()] = le->getBeginVertex();
+ }
+ }
+
+ // --- find out the transformation
+
+ bool translation = true;
+ SVector3 DX;
+
+ int count = 0;
+ for (std::map<GVertex*, GVertex*>::iterator it = vs2vt.begin();
+ it != vs2vt.end() ; ++it){
+ GVertex *vs = it->first;
+ GVertex *vt = it->second;
+ if (count == 0)
+ DX = SVector3(vt->x()-vs->x(),vt->y()-vs->y(),vt->z()-vs->z());
+ else {
+ SVector3 DX2(vt->x()-vs->x(),vt->y()-vs->y(),vt->z()-vs->z());
+ SVector3 DDX(DX2 - DX);
+ if (DDX.norm() > DX.norm() * 1.e-5) translation = false;
+ }
+ count ++;
+ }
+
+ std::vector<double> tfo(16);
+
+ if (translation) {
+ Msg::Info("Periodic mesh translation found: dx = (%g,%g,%g)",
+ DX.x(), DX.y(), DX.z());
+
+ for (size_t i=0;i<16;i++) tfo[i] = 0;
+ for (size_t i=0;i<3;i++) tfo[i*4+i] = 1;
+ tfo[3] = DX.x();
+ tfo[7] = DX.y();
+ tfo[11] = DX.z();
+ }
+
+ else {
+
+ bool rotation = false;
+ myLine LINE;
+ double ANGLE=0;
+
+ count = 0;
+ rotation = true;
+ std::vector<SPoint3> mps, mpt;
+ for (std::map<GVertex*, GVertex*>::iterator it = vs2vt.begin();
+ it != vs2vt.end() ; ++it){
+ GVertex *vs = it->first;
+ GVertex *vt = it->second;
+ mps.push_back(SPoint3(vs->x(), vs->y(), vs->z()));
+ mpt.push_back(SPoint3(vt->x(), vt->y(), vt->z()));
+ }
+ mean_plane mean_source, mean_target;
+ computeMeanPlaneSimple(mps, mean_source);
+ computeMeanPlaneSimple(mpt, mean_target);
+
+ myPlane PLANE_SOURCE(SPoint3(mean_source.x,mean_source.y,mean_source.z),
+ SVector3(mean_source.a,mean_source.b,mean_source.c));
+ myPlane PLANE_TARGET(SPoint3(mean_target.x,mean_target.y,mean_target.z),
+ SVector3(mean_target.a,mean_target.b,mean_target.c));
+
+ LINE = myLine(PLANE_SOURCE, PLANE_TARGET);
+
+ // LINE is the axis of rotation
+ // let us compute the angle of rotation
+ count = 0;
+ for (std::map<GVertex*, GVertex*>::iterator it = vs2vt.begin();
+ it != vs2vt.end(); ++it){
+ GVertex *vs = it->first;
+ GVertex *vt = it->second;
+ // project both points on the axis: that should be the same point !
+ SPoint3 ps = SPoint3(vs->x(), vs->y(), vs->z());
+ SPoint3 pt = SPoint3(vt->x(), vt->y(), vt->z());
+ SPoint3 p_ps = LINE.orthogonalProjection(ps);
+ SPoint3 p_pt = LINE.orthogonalProjection(pt);
+ SVector3 dist1 = ps - pt;
+ SVector3 dist2 = p_ps - p_pt;
+ if (dist1.norm() > CTX::instance()->geom.tolerance){
+ if (dist2.norm() > 1.e-8 * dist1.norm()){
+ rotation = false;
+ }
+ SVector3 t1 = ps - p_ps;
+ SVector3 t2 = pt - p_pt;
+ if (t1.norm() > 1.e-8 * dist1.norm()){
+ if (count == 0)
+ ANGLE = myAngle(t1, t2, LINE.t);
+ else {
+ double ANGLE2 = myAngle(t1, t2, LINE.t);
+ if (fabs (ANGLE2 - ANGLE) > 1.e-8){
+ rotation = false;
+ }
+ }
+ count++;
+ }
+ }
+ }
+
+ if (rotation){
+ Msg::Info("Periodic mesh rotation found: axis (%g,%g,%g) point (%g %g %g) angle %g",
+ LINE.t.x(), LINE.t.y(), LINE.t.z(), LINE.p.x(), LINE.p.y(), LINE.p.z(),
+ ANGLE * 180 / M_PI);
+
+ double ux = LINE.t.x();
+ double uy = LINE.t.y();
+ double uz = LINE.t.z();
+
+ tfo[0*4+0] = cos (ANGLE) + ux*ux*(1.-cos(ANGLE));
+ tfo[0*4+1] = ux*uy*(1.-cos(ANGLE)) - uz * sin(ANGLE);
+ tfo[0*4+2] = ux*uz*(1.-cos(ANGLE)) + uy * sin(ANGLE);
+
+ tfo[1*4+0] = ux*uy*(1.-cos(ANGLE)) + uz * sin(ANGLE);
+ tfo[1*4+1] = cos (ANGLE) + uy*uy*(1.-cos(ANGLE));
+ tfo[1*4+2] = uy*uz*(1.-cos(ANGLE)) - ux * sin(ANGLE);
+
+ tfo[2*4+0] = ux*uz*(1.-cos(ANGLE)) - uy * sin(ANGLE);
+ tfo[2*4+1] = uy*uz*(1.-cos(ANGLE)) + ux * sin(ANGLE);
+ tfo[2*4+2] = cos (ANGLE) + uz*uz*(1.-cos(ANGLE));
+
+ tfo[3] = tfo[7] = tfo[11] = 0;
+ for (int i=0;i<4;i++) tfo[12+i] = 0;
+
+ }
+ else {
+ Msg::Error("Only rotations or translations can currently be computed "
+ "automatically for periodic faces: face %d not meshed",
+ tag());
+ return;
+ }
+ }
+
+ // --- now check and encode the transformation
+ // --- including for edges and vertices
+
+ setMeshMaster(master,tfo);
+}
diff --git a/Geo/GFace.h b/Geo/GFace.h
index a4432a81ecc48f8425b9be3156b3b8c52958e5f3..c8e9428b90f5e2314c09f826cf546bb312e6882d 100644
--- a/Geo/GFace.h
+++ b/Geo/GFace.h
@@ -54,13 +54,16 @@ class GFace : public GEntity
public: // this will become protected or private
std::list<GEdgeLoop> edgeLoops;
-
+
// periodic counterparts of edges
- std::map<int,int> edgeCounterparts;
-
- // encoding of an explicit affine transformation for period meshing
- std::vector<double> affineTransform;
-
+ std::map<GEdge*,std::pair<GEdge*,int> > edgeCounterparts;
+
+ // specify mesh master with transformation, deduce edgeCounterparts
+ virtual void setMeshMaster(GFace* master,const std::vector<double>&);
+
+ // specify mesh master and edgeCounterparts, deduce transformation
+ virtual void setMeshMaster(GFace* master,const std::map<int,int>&);
+
// an array with additional vertices that are supposed to exist in
// the final mesh of the model face. This can be used for boundary
// layer meshes or when using Lloyd-like smoothing algorithms those
diff --git a/Geo/GModelFactory.cpp b/Geo/GModelFactory.cpp
index 54bfcc0d6cc193b6f23213f2718cd59cf5bfa04d..109af03ba801c2f415d2b30ea2411cae7c2e988c 100644
--- a/Geo/GModelFactory.cpp
+++ b/Geo/GModelFactory.cpp
@@ -1220,31 +1220,35 @@ void OCCFactory::setPeriodicAllFaces(GModel *gm, std::vector<double> FaceTransla
/* setPeriodicPairOfFaces: set periodic given a Slave/Master pair of
numFace-Edgelist */
-void OCCFactory::setPeriodicPairOfFaces(GModel *gm, int numFaceMaster, std::vector<int> EdgeListMaster,
- int numFaceSlave, std::vector<int> EdgeListSlave)
+void OCCFactory::setPeriodicPairOfFaces(GModel *gm, int numFaceMaster,
+ std::vector<int> EdgeListMaster,
+ int numFaceSlave,
+ std::vector<int> EdgeListSlave)
{
int NEdges=EdgeListMaster.size();
if (EdgeListMaster.size() != EdgeListSlave.size()){
Msg::Error("Slave/Master faces don't have the same number of edges!");
}
- else{
+ else {
+
+ std::map<int,int> edgeCounterparts;
+ std::vector<int>::iterator siter = EdgeListSlave.begin();
+ std::vector<int>::iterator miter = EdgeListMaster.begin();
+
+ for (;siter!=EdgeListSlave.end();++siter,++miter) {
+ edgeCounterparts[*siter] = *miter;
+ }
+
Surface *s_slave = FindSurface(abs(numFaceSlave));
if(s_slave){
- GModel::current()->getGEOInternals()->periodicFaces[numFaceSlave] = numFaceMaster;
- for (int i = 0; i < NEdges; i++){
- GModel::current()->getGEOInternals()->periodicEdges[EdgeListSlave[i]] = EdgeListMaster[i];
- s_slave->edgeCounterparts[EdgeListSlave[i]] = EdgeListMaster[i];
- }
+ s_slave->master = numFaceMaster;
+ s_slave->edgeCounterparts = edgeCounterparts;
}
else{
GFace *gf = GModel::current()->getFaceByTag(abs(numFaceSlave));
- if(gf){
- gf->setMeshMaster(numFaceMaster);
- for (int i = 0; i < NEdges; i++){
- gf->edgeCounterparts[EdgeListSlave[i]] = EdgeListMaster[i];
- GEdge *ges = GModel::current()->getEdgeByTag(abs(EdgeListSlave[i]));
- ges->setMeshMaster(EdgeListMaster[i]);
- }
+ if (gf) {
+ GFace *master = GModel::current()->getFaceByTag(abs(numFaceMaster));
+ gf->setMeshMaster(master,edgeCounterparts);
}
else Msg::Error("Slave surface %d not found", numFaceSlave);
}
diff --git a/Geo/GModelIO_GEO.cpp b/Geo/GModelIO_GEO.cpp
index c86089b72707297f450e339b8dc958479750fae1..bdbab02a040f65c252065ca2eddb454b449ffdf5 100644
--- a/Geo/GModelIO_GEO.cpp
+++ b/Geo/GModelIO_GEO.cpp
@@ -331,32 +331,41 @@ int GModel::importGEOInternals()
ge->physicals.push_back(pnum);
}
}
-
- // create periodic mesh relationships
- for(std::map<int,int>::iterator it = _geo_internals->periodicEdges.begin();
- it != _geo_internals->periodicEdges.end(); ++it){
- GEdge *ge = getEdgeByTag(abs(it->first));
- if (ge){
- int MASTER = it->second * (it->first > 0 ? 1 : -1);
- ge->setMeshMaster(MASTER);
- }
+
+ std::map<int,GEO_Internals::MasterEdge>::iterator peIter = _geo_internals->periodicEdges.begin();
+ for (;peIter!=_geo_internals->periodicEdges.end();++peIter) {
+
+ int iTarget = peIter->first;
+ GEO_Internals::MasterEdge& me = peIter->second;
+ int iSource = me.tag;
+
+ GEdge* target = getEdgeByTag(iTarget);
+ GEdge* source = getEdgeByTag(iSource);
+
+ if (!target) Msg::Error("Unknown target line for periodic connection from %d to %d",iTarget,iSource);
+ if (!source) Msg::Error("Unknown source line for periodic connection from %d to %d",iTarget,iSource);
+
+ if (me.affineTransform.size()==16) target->setMeshMaster(source,me.affineTransform);
+ else target->setMeshMaster(source,me.tag > 0 ? 1 : -1);
+
}
- for(std::map<int,int>::iterator it = _geo_internals->periodicFaces.begin();
- it != _geo_internals->periodicFaces.end(); ++it){
- GFace *gf = getFaceByTag(abs(it->first));
- if (gf)gf->setMeshMaster(it->second * (it->first > 0 ? 1 : -1));
- }
- for(eiter it = firstEdge() ; it != lastEdge() ; ++it){
- int meshMaster = (*it)->meshMaster();
- if (meshMaster != (*it)->tag()){
- GEdge *ge_master = getEdgeByTag(abs(meshMaster));
- if(ge_master)
- (*it)->getBeginVertex()->setMeshMaster((meshMaster > 0) ? ge_master->getBeginVertex()->tag() :
- ge_master->getEndVertex()->tag());
- if(ge_master)
- (*it)->getEndVertex()->setMeshMaster((meshMaster < 0) ? ge_master->getBeginVertex()->tag() :
- ge_master->getEndVertex()->tag());
- }
+
+
+ std::map<int,GEO_Internals::MasterFace>::iterator pfIter = _geo_internals->periodicFaces.begin();
+ for (;pfIter!=_geo_internals->periodicFaces.end();++pfIter) {
+
+ int iTarget = pfIter->first;
+ GEO_Internals::MasterFace& mf = pfIter->second;
+ int iSource = mf.tag;
+
+ GFace* target = getFaceByTag(iTarget);
+ GFace* source = getFaceByTag(iSource);
+
+ if (!target) Msg::Error("Unknown target surface for periodic connection from %d to %d",iTarget,iSource);
+ if (!source) Msg::Error("Unknown source surface for periodic connection from %d to %d",iTarget,iSource);
+
+ if (mf.affineTransform.size()==16) target->setMeshMaster(source,mf.affineTransform);
+ else target->setMeshMaster(source,mf.edgeCounterparts);
}
Msg::Debug("Gmsh model (GModel) imported:");
diff --git a/Geo/GModelIO_MSH.cpp b/Geo/GModelIO_MSH.cpp
index 832ce58c158ced2bb479df2cbdc0ce75021afd2e..ee696a558d9c1c4beefa2bb97a72218519b00e05 100644
--- a/Geo/GModelIO_MSH.cpp
+++ b/Geo/GModelIO_MSH.cpp
@@ -21,28 +21,28 @@ void writeMSHPeriodicNodes(FILE *fp, std::vector<GEntity*> &entities)
{
int count = 0;
for (unsigned int i = 0; i < entities.size(); i++)
- if (entities[i]->meshMaster() != entities[i]->tag()) count++;
+ if (entities[i]->meshMaster() != entities[i]) count++;
if (!count) return;
fprintf(fp, "$Periodic\n");
fprintf(fp, "%d\n", count);
for(unsigned int i = 0; i < entities.size(); i++){
- GEntity *g_slave = entities[i];
- int meshMaster = g_slave->meshMaster();
- if (g_slave->tag() != meshMaster){
- GEntity *g_master = 0;
- switch(g_slave->dim()){
- case 0 : g_master = g_slave->model()->getVertexByTag(abs(meshMaster));break;
- case 1 : g_master = g_slave->model()->getEdgeByTag(abs(meshMaster));break;
- case 2 : g_master = g_slave->model()->getFaceByTag(abs(meshMaster));break;
- case 3 : g_master = g_slave->model()->getRegionByTag(abs(meshMaster));break;
- }
+ GEntity *g_slave = entities[i];
+ GEntity *g_master = g_slave->meshMaster();
+ if (g_slave != g_master){
fprintf(fp,"%d %d %d\n", g_slave->dim(), g_slave->tag(), g_master->tag());
- fprintf(fp,"%d\n", (int)g_slave->correspondingVertices.size());
+
+ if (g_slave->affineTransform.size() == 16) {
+ fprintf(fp,"Affine");
+ for (int i=0;i<16;i++) fprintf(fp," %.16g",g_slave->affineTransform[i]);
+ fprintf(fp,"\n");
+ }
+
+ fprintf(fp,"%d\n", (int) g_slave->correspondingVertices.size());
for (std::map<MVertex*,MVertex*>::iterator it = g_slave->correspondingVertices.begin();
it != g_slave->correspondingVertices.end(); it++){
- MVertex *v1 = it->first;
- MVertex *v2 = it->second;
- fprintf(fp,"%d %d\n", v1->getIndex(), v2->getIndex());
+ MVertex *v1 = it->first;
+ MVertex *v2 = it->second;
+ fprintf(fp,"%d %d\n", v1->getIndex(), v2->getIndex());
}
}
}
@@ -59,19 +59,33 @@ void readMSHPeriodicNodes(FILE *fp, GModel *gm)
GEntity *s = 0, *m = 0;
switch(dim){
case 0 : s = gm->getVertexByTag(slave); m = gm->getVertexByTag(master); break;
- case 1 : s = gm->getEdgeByTag(slave); m = gm->getEdgeByTag(master); break;
- case 2 : s = gm->getFaceByTag(slave); m = gm->getFaceByTag(master); break;
+ case 1 : s = gm->getEdgeByTag(slave); m = gm->getEdgeByTag(master); break;
+ case 2 : s = gm->getFaceByTag(slave); m = gm->getFaceByTag(master); break;
}
if (s && m){
- s->setMeshMaster(m->tag());
+
+ char token[6];
+ fpos_t pos;
+ fgetpos(fp,&pos);
+ fscanf(fp,"%s",token);
+ if (strcmp(token,"Affine") == 0) {
+ std::vector<double> tfo;
+ for (int i=0;i<16;i++) fscanf(fp,"%lf",&tfo[i]);
+ s->setMeshMaster(m,tfo);
+ }
+ else {
+ fsetpos(fp,&pos);
+ s->setMeshMaster(m);
+ }
+
int numv;
if(fscanf(fp, "%d", &numv) != 1) numv = 0;
for(int j = 0; j < numv; j++){
- int v1, v2;
- if(fscanf(fp,"%d %d", &v1, &v2) != 2) continue;
- MVertex *mv1 = gm->getMeshVertexByTag(v1);
- MVertex *mv2 = gm->getMeshVertexByTag(v2);
- s->correspondingVertices[mv1] = mv2;
+ int v1, v2;
+ if(fscanf(fp,"%d %d", &v1, &v2) != 2) continue;
+ MVertex *mv1 = gm->getMeshVertexByTag(v1);
+ MVertex *mv2 = gm->getMeshVertexByTag(v2);
+ s->correspondingVertices[mv1] = mv2;
}
}
}
diff --git a/Geo/GModelIO_MSH2.cpp b/Geo/GModelIO_MSH2.cpp
index 80bda67c8ebe14d6240267389837485845907a56..560fbd83c43cffeb25e41e0c571d0db803da82b1 100644
--- a/Geo/GModelIO_MSH2.cpp
+++ b/Geo/GModelIO_MSH2.cpp
@@ -435,11 +435,11 @@ int GModel::_readMSH2(const std::string &name)
bool own = false;
// search parent element
- if (parent){
+ if (parent != 0){
#if (FAST_ELEMENTS == 1)
std::map<int, MElement* >::iterator ite = elems.find(parent);
if (ite == elems.end())
- Msg::Error("Parent element %d not found for element %d", parent, num);
+ Msg::Error("Parent element (ascii) %d not found for element %d of type %d", parent, num, type);
else{
p = ite->second;
parents[parent] = p;
@@ -456,7 +456,7 @@ int GModel::_readMSH2(const std::string &name)
p = getParent(parent, type, elements);
if (p)
parents[parent] = p;
- else Msg::Error("Parent element %d not found for element %d", parent, num);
+ else Msg::Error("Parent element %d not found for element %d dedju", parent, num);
}
else
p = itp->second;
@@ -564,7 +564,7 @@ int GModel::_readMSH2(const std::string &name)
#if (FAST_ELEMENTS == 1)
std::map<int, MElement* >::iterator ite = elems.find(parent);
if (ite == elems.end())
- Msg::Error("Parent element %d not found for element %d", parent, num);
+ Msg::Error("Parent (binary) element %d not found for element %d", parent, num);
else{
p = ite->second;
parents[parent] = p;
@@ -580,7 +580,7 @@ int GModel::_readMSH2(const std::string &name)
if(itp == parents.end()){
p = getParent(parent, type, elements);
if(p) parents[parent] = p;
- else Msg::Error("Parent element %d not found", parent);
+ else Msg::Error("Parent (binary) element %d not found", parent);
}
else p = itp->second;
std::set<MElement* >::iterator itpo = parentsOwned.find(p);
diff --git a/Geo/Geo.h b/Geo/Geo.h
index 01eb42b3a9bb73103506c5a7cba1c2680b6be438..7e22189f2ab1ef56b77da3b3af15a2adfbf5cb67 100644
--- a/Geo/Geo.h
+++ b/Geo/Geo.h
@@ -8,6 +8,7 @@
#include <math.h>
#include <vector>
+#include <list>
#include "GmshDefines.h"
#include "gmshSurface.h"
#include "ListUtils.h"
@@ -137,6 +138,8 @@ class Curve{
gmshSurface *geometry;
std::vector<int> compound;
int ReverseMesh;
+ int master;
+ std::list<double> affineTransformation;
void SetVisible(int value, bool recursive)
{
Visible = value;
@@ -191,6 +194,8 @@ class Surface{
// should be the only one in gmsh, so parameter "Type" should
// disappear from the class Surface.
gmshSurface *geometry;
+ //
+ int master;
// the mesh master surface
std::map<int,int> edgeCounterparts;
// prescribed affine transform for periodic meshing
@@ -283,6 +288,7 @@ class GEO_Internals{
void alloc_all();
void free_all();
public:
+
Tree_T *Points;
Tree_T *Curves;
Tree_T *Surfaces;
@@ -293,12 +299,28 @@ class GEO_Internals{
List_T *PhysicalGroups;
int MaxPointNum, MaxLineNum, MaxLineLoopNum, MaxSurfaceNum;
int MaxSurfaceLoopNum, MaxVolumeNum, MaxPhysicalNum;
+
+
+ struct MasterEdge {
+ int tag; // signed
+ std::vector<double> affineTransform;
+ };
+
+ std::map<int,MasterEdge> periodicEdges;
+
+ struct MasterFace {
+ int tag;
+ // map from slave to master edges
+ std::map<int,int> edgeCounterparts;
+ std::vector<double> affineTransform;
+ };
+
+ std::map<int,MasterFace> periodicFaces;
+
GEO_Internals(){ alloc_all(); }
~GEO_Internals(){ free_all(); }
void destroy(){ free_all(); alloc_all(); }
void reset_physicals();
- std::map<int,int> periodicFaces;
- std::map<int,int> periodicEdges;
};
class Shape{
diff --git a/Geo/gmshFace.cpp b/Geo/gmshFace.cpp
index 5762b8c0b0d837ebd0f6c824c452f37d98fd55ab..efdc745c864315b97034f8c177f188454d4964e9 100644
--- a/Geo/gmshFace.cpp
+++ b/Geo/gmshFace.cpp
@@ -25,8 +25,8 @@ gmshFace::gmshFace(GModel *m, Surface *face)
{
resetMeshAttributes();
- edgeCounterparts = s->edgeCounterparts;
- affineTransform = s->affineTransform;
+ // edgeCounterparts = s->edgeCounterparts;
+ // affineTransform = s->affineTransform;
std::vector<GEdge*> eds;
std::vector<int> nums;
diff --git a/Mesh/Generator.cpp b/Mesh/Generator.cpp
index 30ab454e94596facad10cb56cb0c6b58decd4b91..779b4885ef9f71449497886f2674a778d8b2d4dd 100644
--- a/Mesh/Generator.cpp
+++ b/Mesh/Generator.cpp
@@ -219,14 +219,13 @@ static void Mesh0D(GModel *m)
}
for(GModel::viter it = m->firstVertex(); it != m->lastVertex(); ++it){
GVertex *gv = *it;
- if (gv->meshMaster() != gv->tag()){
+ if (gv->meshMaster() != gv){
if (gv->correspondingVertices.empty()){
- GVertex *master = m->getVertexByTag(abs(gv->meshMaster()));
+ GVertex *master = dynamic_cast<GVertex*> (gv->meshMaster());
if(master)gv->correspondingVertices[gv->mesh_vertices[0]] = (master->mesh_vertices[0]);
}
}
}
-
}
static void Mesh1D(GModel *m)
diff --git a/Mesh/meshGEdge.cpp b/Mesh/meshGEdge.cpp
index 8ca50b2b48580a0db25be165edc862447843fb33..7e11458416a9d5019c8c2d1082e829b2e7ff3cbe 100644
--- a/Mesh/meshGEdge.cpp
+++ b/Mesh/meshGEdge.cpp
@@ -267,8 +267,8 @@ static double Integration(GEdge *ge, double t1, double t2,
return Points.back().p;
}
-static void copyMesh(GEdge *from, GEdge *to, int direction)
-{
+void copyMesh(GEdge *from, GEdge *to, int direction)
+{
Range<double> u_bounds = from->parBounds(0);
double u_min = u_bounds.low();
double u_max = u_bounds.high();
@@ -351,12 +351,12 @@ void meshGEdge::operator() (GEdge *ge)
if(MeshExtrudedCurve(ge)) return;
- if (ge->meshMaster() != ge->tag()){
- GEdge *gef = ge->model()->getEdgeByTag(abs(ge->meshMaster()));
+ if (ge->meshMaster() != ge){
+ GEdge *gef = dynamic_cast<GEdge*> (ge->meshMaster());
if (gef->meshStatistics.status == GEdge::PENDING) return;
Msg::Info("Meshing curve %d (%s) as a copy of %d", ge->tag(),
- ge->getTypeString().c_str(), ge->meshMaster());
- copyMesh(gef, ge, ge->meshMaster());
+ ge->getTypeString().c_str(), ge->meshMaster()->tag());
+ copyMesh(gef, ge, ge->masterOrientation);
ge->meshStatistics.status = GEdge::DONE;
return;
}
diff --git a/Mesh/meshGFace.cpp b/Mesh/meshGFace.cpp
index 2e700563277b98622850e55692f857bd07bf9ca4..a40c39c3b32580ad0a28b7bf8fbfff8b90e84205 100644
--- a/Mesh/meshGFace.cpp
+++ b/Mesh/meshGFace.cpp
@@ -49,6 +49,8 @@
// define this to use the old initial delaunay
#define OLD_CODE_DELAUNAY 1
+void copyMesh(GEdge*,GEdge*,int);
+
static void computeElementShapes(GFace *gf, double &worst, double &avg,
double &best, int &nT, int &greaterThan)
{
@@ -69,12 +71,6 @@ static void computeElementShapes(GFace *gf, double &worst, double &avg,
}
-inline double myAngle(const SVector3 &a, const SVector3 &b, const SVector3 &d)
-{
- double cosTheta = dot(a, b);
- double sinTheta = dot(crossprod(a, b), d);
- return atan2(sinTheta, cosTheta);
-}
class quadMeshRemoveHalfOfOneDMesh
{
@@ -221,257 +217,119 @@ public:
};
-struct myPlane {
- SPoint3 p;
- SVector3 n;
- double a;
- // nx x + ny y + nz z + a = 0
- myPlane(SPoint3 _p, SVector3 _n) : p(_p),n(_n)
- {
- n.normalize();
- a = -(n.x()*p.x()+n.y()*p.y()+n.z()*p.z());
- }
- double eval (double x, double y, double z)
- {
- return n.x() * x + n.y() * y + n.z() * z + a;
- }
-};
-
-struct myLine {
- SPoint3 p;
- SVector3 t;
- myLine() : p(0,0,0) , t (0,0,1) {}
- myLine(myPlane &p1, myPlane &p2)
- {
- t = crossprod(p1.n, p2.n);
- if (t.norm() == 0.0){
- Msg::Error("parallel planes do not intersect");
- }
- else
- t.normalize();
- // find a point, assume z = 0
- double a[2][2] = {{p1.n.x(), p1.n.y()}, {p2.n.x(), p2.n.y()}};
- double b[2] = {-p1.a, -p2.a}, x[2];
- if (!sys2x2(a, b, x)){
- // assume x = 0
- double az[2][2] = {{p1.n.y(), p1.n.z()}, {p2.n.y(), p2.n.z()}};
- double bz[2] = {-p1.a, -p2.a};
- if (!sys2x2(az, bz, x)){
- // assume y = 0
- double ay[2][2] = {{p1.n.x(), p1.n.z()}, {p2.n.x(), p2.n.z()}};
- double by[2] = {-p1.a, -p2.a};
- if (!sys2x2(ay,by,x)){
- Msg::Error("parallel planes do not intersect");
- }
- else {
- p = SPoint3(x[0], 0., x[1]);
- }
- }
- else{
- p = SPoint3(0., x[0], x[1]);
- }
- }
- else{
- p = SPoint3(x[0], x[1], 0.);
- }
- }
- SPoint3 orthogonalProjection (SPoint3 &a)
- {
- // (x - a) * t = 0 -->
- // x = p + u t --> (p + ut - a) * t = 0 --> u = (a-p) * t
- const double u = dot(a - p, t);
- return SPoint3(p.x() + t.x() * u,p.y() + t.y() * u,p.z() + t.z() * u);
- }
-};
static void copyMesh(GFace *source, GFace *target)
{
- std::map<MVertex*, MVertex*> vs2vt;
- std::list<GEdge*> edges = target->edges();
- {
- for (std::list<GEdge*>::iterator it = edges.begin(); it != edges.end(); ++it){
- int sign = 1;
- std::map<int, int>::iterator adnksd = target->edgeCounterparts.find((*it)->tag());
- int source_e;
- if(adnksd != target->edgeCounterparts.end())
- source_e = adnksd->second;
- else{
- sign = -1;
- adnksd = target->edgeCounterparts.find(-(*it)->tag());
- if(adnksd != target->edgeCounterparts.end())
- source_e = adnksd->second;
- else{
- Msg::Error("Could not find edge counterpart %d in slave surface %d",
- (*it)->tag(), target->tag());
- return;
- }
- }
- GEdge *se = source->model()->getEdgeByTag(abs(source_e));
- GEdge *te = *it;
- if (source_e * sign > 0){
- vs2vt[se->getBeginVertex()->mesh_vertices[0]] =
- te->getBeginVertex()->mesh_vertices[0];
- vs2vt[se->getEndVertex()->mesh_vertices[0]] =
- te->getEndVertex()->mesh_vertices[0];
- }
- else {
- vs2vt[se->getBeginVertex()->mesh_vertices[0]] =
- te->getEndVertex()->mesh_vertices[0];
- vs2vt[se->getEndVertex()->mesh_vertices[0]] =
- te->getBeginVertex()->mesh_vertices[0];
- }
- // iterate on source vertices
- for (unsigned i = 0; i < te->mesh_vertices.size(); i++){
- MVertex *vt = te->mesh_vertices[i];
- MVertex *vs = se->mesh_vertices[source_e * sign > 0 ? i :
- te->mesh_vertices.size() - i - 1];
- vs2vt[vs] = vt;
- }
- }
- }
-
- bool translation = true;
- SVector3 DX;
+
+#warning periodic face mesh does not take embedded vertices/edges into account
- bool rotation = false;
- double rot[3][3] ;
- myLine LINE;
- double ANGLE=0;
-
- bool affine = false;
- double mat[4][4];
-
- if(target->affineTransform.size() == 16){
- Msg::Info("Affine transformation specified");
- affine = true;
- for(int i = 0; i < 4; i++)
- for(int j = 0; j < 4; j++)
- mat[i][j] = target->affineTransform[4 * i + j];
- }
- else{
- int count = 0;
- for (std::map<MVertex*, MVertex*>::iterator it = vs2vt.begin();
- it != vs2vt.end() ; ++it){
- MVertex *vs = it->first;
- MVertex *vt = it->second;
- if (count == 0)
- DX = SVector3(vt->x() - vs->x(), vt->y() - vs->y(), vt->z() - vs->z());
- else {
- SVector3 DX2 = DX - SVector3(vt->x() - vs->x(), vt->y() - vs->y(),
- vt->z() - vs->z());
- if (DX2.norm() > DX.norm() * 1.e-5) translation = false;
- }
- count ++;
- }
-
- if (!translation){
- count = 0;
- rotation = true;
- std::vector<SPoint3> mps, mpt;
- for (std::map<MVertex*, MVertex*>::iterator it = vs2vt.begin();
- it != vs2vt.end() ; ++it){
- MVertex *vs = it->first;
- MVertex *vt = it->second;
- mps.push_back(SPoint3(vs->x(), vs->y(), vs->z()));
- mpt.push_back(SPoint3(vt->x(), vt->y(), vt->z()));
- }
- mean_plane mean_source, mean_target;
- computeMeanPlaneSimple(mps, mean_source);
- computeMeanPlaneSimple(mpt, mean_target);
- myPlane PLANE_SOURCE(SPoint3(mean_source.x,mean_source.y,mean_source.z),
- SVector3(mean_source.a,mean_source.b,mean_source.c));
- myPlane PLANE_TARGET(SPoint3(mean_target.x,mean_target.y,mean_target.z),
- SVector3(mean_target.a,mean_target.b,mean_target.c));
- LINE = myLine(PLANE_SOURCE, PLANE_TARGET);
-
- // LINE is the axis of rotation
- // let us compute the angle of rotation
- count = 0;
- for (std::map<MVertex*, MVertex*>::iterator it = vs2vt.begin();
- it != vs2vt.end(); ++it){
- MVertex *vs = it->first;
- MVertex *vt = it->second;
- // project both points on the axis: that should be the same point !
- SPoint3 ps = SPoint3(vs->x(), vs->y(), vs->z());
- SPoint3 pt = SPoint3(vt->x(), vt->y(), vt->z());
- SPoint3 p_ps = LINE.orthogonalProjection(ps);
- SPoint3 p_pt = LINE.orthogonalProjection(pt);
- SVector3 dist1 = ps - pt;
- SVector3 dist2 = p_ps - p_pt;
- if (dist1.norm() > CTX::instance()->geom.tolerance){
- if (dist2.norm() > 1.e-8 * dist1.norm()){
- rotation = false;
- }
- SVector3 t1 = ps - p_ps;
- SVector3 t2 = pt - p_pt;
- if (t1.norm() > 1.e-8 * dist1.norm()){
- if (count == 0)
- ANGLE = myAngle(t1, t2, LINE.t);
- else {
- double ANGLE2 = myAngle(t1, t2, LINE.t);
- if (fabs (ANGLE2 - ANGLE) > 1.e-8){
- rotation = false;
- }
- }
- count++;
- }
- }
- }
+ std::map<MVertex*, MVertex*> vs2vt;
+
+ // add principal vertex pairs
+
+ std::list<GVertex*> s_vtcs = source->vertices();
+ std::list<GVertex*> t_vtcs = target->vertices();
+
+
+ if (s_vtcs.size() != t_vtcs.size()) {
+ Msg::Info("Periodicity imposed on topologically incompatible surfaces"
+ "(%d vs %d bounding vertices)",s_vtcs.size(),t_vtcs.size());
+ }
+
+
+ std::set<GVertex*> checkVtcs;
+ checkVtcs.insert(s_vtcs.begin(),s_vtcs.end());
+
+ for (std::list<GVertex*>::iterator tvIter=t_vtcs.begin();tvIter!=t_vtcs.end();++tvIter) {
+
+ GVertex* gvt = *tvIter;
+ std::map<GVertex*,GVertex*>::iterator gvsIter = target->vertexCounterparts.find(gvt);
+
+
+ if (gvsIter == target->vertexCounterparts.end()) {
+ Msg::Info("Error during periodic meshing of surface %d with surface %d:"
+ "vertex %d has no periodic counterpart",
+ target->tag(),source->tag(),gvt->tag());
+ }
+
+ GVertex* gvs = gvsIter->second;
+ if (checkVtcs.find(gvs) == checkVtcs.end()) {
+ if (gvs) Msg::Info("Error during periodic meshing of surface %d with surface %d:"
+ "vertex %d has periodic counterpart %d outside of source surface",
+ target->tag(),source->tag(),gvt->tag(),gvs->tag());
+
+ else Msg::Info("Error during periodic meshing of surface %d with surface %d:"
+ "vertex %d has no periodic counterpart",
+ target->tag(),source->tag(),gvt->tag());
+ }
+
+ vs2vt[gvs->mesh_vertices[0]] = gvt->mesh_vertices[0];
+ }
+
+
+ // add corresponding edge nodes assuming edges were correctly meshed already
+
+ std::list<GEdge*> s_edges = source->edges();
+ std::list<GEdge*> t_edges = target->edges();
+
+ std::set<GEdge*> checkEdges;
+ checkEdges.insert(s_edges.begin(),s_edges.end());
+
+ for (std::list<GEdge*>::iterator te_iter = t_edges.begin();te_iter!=t_edges.end();++te_iter) {
+
+ GEdge* get = *te_iter;
+ std::map<GEdge*,std::pair<GEdge*,int> >::iterator gesIter = target->edgeCounterparts.find(get);
+
+ if (gesIter == target->edgeCounterparts.end()) {
+ Msg::Info("Error during periodic meshing of surface %d with surface %d:"
+ "edge %d has no periodic counterpart",
+ target->tag(),source->tag(),get->tag());
+ }
+
+ GEdge* ges = gesIter->second.first;
+
+ if (checkEdges.find(ges) == checkEdges.end()) {
+ Msg::Info("Error during periodic meshing of surface %d with surface %d:"
+ "edge %d has periodic counterpart %d outside of get surface",
+ target->tag(),source->tag(),get->tag(),ges->tag());
+ }
+
+
+ if (get->mesh_vertices.size() != ges->mesh_vertices.size()) {
+ Msg::Info("Error during periodic meshing of surface %d with surface %d:"
+ "edge %d has %d vertices, whereas correspondant %d has %d",
+ target->tag(),source->tag(),
+ get->tag(),get->mesh_vertices.size(),
+ ges->tag(),ges->mesh_vertices.size());
+ }
+
+ int orientation = gesIter->second.second;
+ int is = orientation == 1 ? 0 : get->mesh_vertices.size()-1;
+
+ for (unsigned it=0;it<get->mesh_vertices.size();it++,is+=orientation) {
+ vs2vt[ges->mesh_vertices[is]] = get->mesh_vertices[it];
+ }
+ }
+
+ // now transform
+
+ std::vector<double>& tfo = target->affineTransform;
- if (rotation){
- Msg::Info("Periodic mesh rotation found: axis (%g,%g,%g) point (%g %g %g) angle %g",
- LINE.t.x(), LINE.t.y(), LINE.t.z(), LINE.p.x(), LINE.p.y(), LINE.p.z(),
- ANGLE * 180 / M_PI);
- double ux = LINE.t.x();
- double uy = LINE.t.y();
- double uz = LINE.t.z();
- rot[0][0] = cos (ANGLE) + ux*ux*(1.-cos(ANGLE));
- rot[0][1] = ux*uy*(1.-cos(ANGLE)) - uz * sin(ANGLE);
- rot[0][2] = ux*uz*(1.-cos(ANGLE)) + uy * sin(ANGLE);
- rot[1][0] = ux*uy*(1.-cos(ANGLE)) + uz * sin(ANGLE);
- rot[1][1] = cos (ANGLE) + uy*uy*(1.-cos(ANGLE));
- rot[1][2] = uy*uz*(1.-cos(ANGLE)) - ux * sin(ANGLE);
- rot[2][0] = ux*uz*(1.-cos(ANGLE)) - uy * sin(ANGLE);
- rot[2][1] = uy*uz*(1.-cos(ANGLE)) + ux * sin(ANGLE);
- rot[2][2] = cos (ANGLE) + uz*uz*(1.-cos(ANGLE));
- }
- else {
- Msg::Error("Only rotations or translations can be currently taken into account "
- "automatically for periodic faces: face %d not meshed", target->tag());
- return;
- }
- }
- else{
- Msg::Info("Periodic mesh translation found: dx = (%g,%g,%g)",
- DX.x(), DX.y(), DX.z());
- }
- }
- // now transform
for(unsigned int i = 0; i < source->mesh_vertices.size(); i++){
MVertex *vs = source->mesh_vertices[i];
SPoint2 XXX;
- if (affine) {
- double ps[4] = {vs->x(), vs->y(), vs->z(), 1.};
- double res[4] = {0., 0., 0., 0.};
- for(int i = 0; i < 4; i++)
- for(int j = 0; j < 4; j++)
- res[i] += mat[i][j] * ps[j];
- SPoint3 tp (res[0], res[1], res[2]);
- XXX = target->parFromPoint(tp);
- }
- else if (translation) {
- SPoint3 tp (vs->x() + DX.x(), vs->y() + DX.y(), vs->z() + DX.z());
- XXX = target->parFromPoint(tp);
- }
- else if (rotation){
- SPoint3 ps = SPoint3(vs->x(),vs->y(),vs->z());
- SPoint3 p_ps = LINE.orthogonalProjection(ps);
- SPoint3 P = ps - p_ps, res;
- matvec(rot, P, res);
- res += p_ps;
- XXX = target->parFromPoint(res);
- }
+
+ double ps[4] = {vs->x(), vs->y(), vs->z(), 1.};
+ double res[4] = {0., 0., 0., 0.};
+ int idx = 0;
+ for(int i = 0; i < 4; i++)
+ for(int j = 0; j < 4; j++)
+ res[i] += tfo[idx++] * ps[j];
+
+ SPoint3 tp (res[0], res[1], res[2]);
+ XXX = target->parFromPoint(tp);
+
GPoint gp = target->point(XXX);
MVertex *vt = new MFaceVertex(gp.x(), gp.y(), gp.z(), target, gp.u(), gp.v());
target->mesh_vertices.push_back(vt);
@@ -2534,21 +2392,21 @@ void meshGFace::operator() (GFace *gf, bool print)
if(MeshTransfiniteSurface(gf)) return;
if(MeshExtrudedSurface(gf)) return;
- if(gf->meshMaster() != gf->tag()){
- GFace *gff = gf->model()->getFaceByTag(abs(gf->meshMaster()));
+ if(gf->meshMaster() != gf){
+ GFace *gff = dynamic_cast<GFace*> (gf->meshMaster());
if(gff){
if (gff->meshStatistics.status != GFace::DONE){
gf->meshStatistics.status = GFace::PENDING;
return;
}
Msg::Info("Meshing face %d (%s) as a copy of %d", gf->tag(),
- gf->getTypeString().c_str(), gf->meshMaster());
+ gf->getTypeString().c_str(), gf->meshMaster()->tag());
copyMesh(gff, gf);
gf->meshStatistics.status = GFace::DONE;
return;
}
else
- Msg::Warning("Unknown mesh master face %d", abs(gf->meshMaster()));
+ Msg::Warning("Unknown mesh master face %d", gf->meshMaster()->tag());
}
const char *algo = "Unknown";
diff --git a/Parser/CMakeLists.txt b/Parser/CMakeLists.txt
index 3795b90be4b3668e0893c8d9f82baad0c083b20f..c909e3b90d35605ae43d30f6ebd6114ae41d8622 100644
--- a/Parser/CMakeLists.txt
+++ b/Parser/CMakeLists.txt
@@ -4,8 +4,6 @@
# bugs and problems to the public mailing list <gmsh@geuz.org>.
set(SRC
- Gmsh.tab.cpp
- Gmsh.yy.cpp
FunctionManager.cpp
)
diff --git a/Parser/Gmsh.l b/Parser/Gmsh.l
index 304126ce8625bc48e6e2c9e4b7f544e28616e152..2ab84df4a9faa9d235c13a5c44ad048bda395ab4 100644
--- a/Parser/Gmsh.l
+++ b/Parser/Gmsh.l
@@ -196,6 +196,7 @@ Knots return tNurbsKnots;
Length return tLength;
Line return tLine;
+Lines return tLines;
List return tList;
Log return tLog;
Log10 return tLog10;
@@ -257,6 +258,7 @@ Sphere return tSphere;
Spline return tSpline;
Split return tSplit;
Surface return tSurface;
+Surfaces return tSurfaces;
Str return tStr;
Sprintf return tSprintf;
StrCat return tStrCat;
@@ -275,6 +277,7 @@ T3 return tText3D;
TextAttributes return tTextAttributes;
TIME return tTime;
Transfinite return tTransfinite;
+Transform return tTransform;
Transf[qQ]uad[tT]ri return tTransfQuadTri;
Translate return tTranslate;
Tanh return tTanh;
diff --git a/Parser/Gmsh.tab.hpp b/Parser/Gmsh.tab.hpp
deleted file mode 100644
index b80b1d82ce499dbc0b13e0755ddd874c41f51fd5..0000000000000000000000000000000000000000
--- a/Parser/Gmsh.tab.hpp
+++ /dev/null
@@ -1,267 +0,0 @@
-/* A Bison parser, made by GNU Bison 2.7.12-4996. */
-
-/* Bison interface for Yacc-like parsers in C
-
- Copyright (C) 1984, 1989-1990, 2000-2013 Free Software Foundation, Inc.
-
- This program is free software: you can redistribute it and/or modify
- it under the terms of the GNU General Public License as published by
- the Free Software Foundation, either version 3 of the License, or
- (at your option) any later version.
-
- This program is distributed in the hope that it will be useful,
- but WITHOUT ANY WARRANTY; without even the implied warranty of
- MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
- GNU General Public License for more details.
-
- You should have received a copy of the GNU General Public License
- along with this program. If not, see <http://www.gnu.org/licenses/>. */
-
-/* As a special exception, you may create a larger work that contains
- part or all of the Bison parser skeleton and distribute that work
- under terms of your choice, so long as that work isn't itself a
- parser generator using the skeleton or a modified version thereof
- as a parser skeleton. Alternatively, if you modify or redistribute
- the parser skeleton itself, you may (at your option) remove this
- special exception, which will cause the skeleton and the resulting
- Bison output files to be licensed under the GNU General Public
- License without this special exception.
-
- This special exception was added by the Free Software Foundation in
- version 2.2 of Bison. */
-
-#ifndef YY_GMSH_YY_GMSH_TAB_HPP_INCLUDED
-# define YY_GMSH_YY_GMSH_TAB_HPP_INCLUDED
-/* Enabling traces. */
-#ifndef YYDEBUG
-# define YYDEBUG 0
-#endif
-#if YYDEBUG
-extern int gmsh_yydebug;
-#endif
-
-/* Tokens. */
-#ifndef YYTOKENTYPE
-# define YYTOKENTYPE
- /* Put the tokens into the symbol table, so that GDB and other debuggers
- know about them. */
- enum yytokentype {
- tDOUBLE = 258,
- tSTRING = 259,
- tBIGSTR = 260,
- tEND = 261,
- tAFFECT = 262,
- tDOTS = 263,
- tPi = 264,
- tMPI_Rank = 265,
- tMPI_Size = 266,
- tEuclidian = 267,
- tCoordinates = 268,
- tExp = 269,
- tLog = 270,
- tLog10 = 271,
- tSqrt = 272,
- tSin = 273,
- tAsin = 274,
- tCos = 275,
- tAcos = 276,
- tTan = 277,
- tRand = 278,
- tAtan = 279,
- tAtan2 = 280,
- tSinh = 281,
- tCosh = 282,
- tTanh = 283,
- tFabs = 284,
- tFloor = 285,
- tCeil = 286,
- tRound = 287,
- tFmod = 288,
- tModulo = 289,
- tHypot = 290,
- tList = 291,
- tPrintf = 292,
- tError = 293,
- tStr = 294,
- tSprintf = 295,
- tStrCat = 296,
- tStrPrefix = 297,
- tStrRelative = 298,
- tStrReplace = 299,
- tFind = 300,
- tStrFind = 301,
- tStrCmp = 302,
- tStrChoice = 303,
- tTextAttributes = 304,
- tBoundingBox = 305,
- tDraw = 306,
- tSetChanged = 307,
- tToday = 308,
- tFixRelativePath = 309,
- tCurrentDirectory = 310,
- tSyncModel = 311,
- tNewModel = 312,
- tOnelabAction = 313,
- tOnelabRun = 314,
- tCpu = 315,
- tMemory = 316,
- tTotalMemory = 317,
- tCreateTopology = 318,
- tCreateTopologyNoHoles = 319,
- tDistanceFunction = 320,
- tDefineConstant = 321,
- tUndefineConstant = 322,
- tDefineNumber = 323,
- tDefineString = 324,
- tSetNumber = 325,
- tSetString = 326,
- tPoint = 327,
- tCircle = 328,
- tEllipse = 329,
- tLine = 330,
- tSphere = 331,
- tPolarSphere = 332,
- tSurface = 333,
- tSpline = 334,
- tVolume = 335,
- tCharacteristic = 336,
- tLength = 337,
- tParametric = 338,
- tElliptic = 339,
- tRefineMesh = 340,
- tAdaptMesh = 341,
- tRelocateMesh = 342,
- tPlane = 343,
- tRuled = 344,
- tTransfinite = 345,
- tComplex = 346,
- tPhysical = 347,
- tCompound = 348,
- tPeriodic = 349,
- tUsing = 350,
- tPlugin = 351,
- tDegenerated = 352,
- tRecursive = 353,
- tRotate = 354,
- tTranslate = 355,
- tSymmetry = 356,
- tDilate = 357,
- tExtrude = 358,
- tLevelset = 359,
- tRecombine = 360,
- tSmoother = 361,
- tSplit = 362,
- tDelete = 363,
- tCoherence = 364,
- tIntersect = 365,
- tMeshAlgorithm = 366,
- tReverse = 367,
- tLayers = 368,
- tScaleLast = 369,
- tHole = 370,
- tAlias = 371,
- tAliasWithOptions = 372,
- tCopyOptions = 373,
- tQuadTriAddVerts = 374,
- tQuadTriNoNewVerts = 375,
- tQuadTriSngl = 376,
- tQuadTriDbl = 377,
- tRecombLaterals = 378,
- tTransfQuadTri = 379,
- tText2D = 380,
- tText3D = 381,
- tInterpolationScheme = 382,
- tTime = 383,
- tCombine = 384,
- tBSpline = 385,
- tBezier = 386,
- tNurbs = 387,
- tNurbsOrder = 388,
- tNurbsKnots = 389,
- tColor = 390,
- tColorTable = 391,
- tFor = 392,
- tIn = 393,
- tEndFor = 394,
- tIf = 395,
- tEndIf = 396,
- tExit = 397,
- tAbort = 398,
- tField = 399,
- tReturn = 400,
- tCall = 401,
- tMacro = 402,
- tShow = 403,
- tHide = 404,
- tGetValue = 405,
- tGetEnv = 406,
- tGetString = 407,
- tGetNumber = 408,
- tHomology = 409,
- tCohomology = 410,
- tBetti = 411,
- tSetOrder = 412,
- tExists = 413,
- tFileExists = 414,
- tGMSH_MAJOR_VERSION = 415,
- tGMSH_MINOR_VERSION = 416,
- tGMSH_PATCH_VERSION = 417,
- tGmshExecutableName = 418,
- tSetPartition = 419,
- tAFFECTDIVIDE = 420,
- tAFFECTTIMES = 421,
- tAFFECTMINUS = 422,
- tAFFECTPLUS = 423,
- tOR = 424,
- tAND = 425,
- tNOTEQUAL = 426,
- tEQUAL = 427,
- tGREATEROREQUAL = 428,
- tLESSOREQUAL = 429,
- UNARYPREC = 430,
- tMINUSMINUS = 431,
- tPLUSPLUS = 432
- };
-#endif
-
-
-#if ! defined YYSTYPE && ! defined YYSTYPE_IS_DECLARED
-typedef union YYSTYPE
-{
-/* Line 2053 of yacc.c */
-#line 99 "Gmsh.y"
-
- char *c;
- int i;
- unsigned int u;
- double d;
- double v[5];
- Shape s;
- List_T *l;
-
-
-/* Line 2053 of yacc.c */
-#line 245 "Gmsh.tab.hpp"
-} YYSTYPE;
-# define YYSTYPE_IS_TRIVIAL 1
-# define yystype YYSTYPE /* obsolescent; will be withdrawn */
-# define YYSTYPE_IS_DECLARED 1
-#endif
-
-extern YYSTYPE gmsh_yylval;
-
-#ifdef YYPARSE_PARAM
-#if defined __STDC__ || defined __cplusplus
-int gmsh_yyparse (void *YYPARSE_PARAM);
-#else
-int gmsh_yyparse ();
-#endif
-#else /* ! YYPARSE_PARAM */
-#if defined __STDC__ || defined __cplusplus
-int gmsh_yyparse (void);
-#else
-int gmsh_yyparse ();
-#endif
-#endif /* ! YYPARSE_PARAM */
-
-#endif /* !YY_GMSH_YY_GMSH_TAB_HPP_INCLUDED */
diff --git a/Parser/Gmsh.y b/Parser/Gmsh.y
index e1affe4680d7ff3442c6d7eacce0099f06e08614..8d1e662afe14a87fdb9ad5641be1de1ca70d1b76 100644
--- a/Parser/Gmsh.y
+++ b/Parser/Gmsh.y
@@ -5,6 +5,7 @@
// bugs and problems to the public mailing list <gmsh@geuz.org>.
#include <sstream>
+#include <map>
#include <string.h>
#include <stdarg.h>
#include <time.h>
@@ -89,6 +90,11 @@ int PrintListOfDouble(char *format, List_T *list, char *buffer);
void PrintParserSymbols(std::vector<std::string> &vec);
fullMatrix<double> ListOfListOfDouble2Matrix(List_T *list);
+void addPeriodicEdge(int,int,const std::vector<double>&);
+void addPeriodicFace(int,int,const std::map<int,int>&);
+void addPeriodicFace(int,int,const std::vector<double>&);
+void computeAffineTransformation(SPoint3&,SPoint3&,double,SPoint3&,std::vector<double>&);
+
struct doubleXstring{
double d;
char *s;
@@ -124,9 +130,10 @@ struct doubleXstring{
%token tDistanceFunction tDefineConstant tUndefineConstant
%token tDefineNumber tDefineString tSetNumber tSetString
%token tPoint tCircle tEllipse tLine tSphere tPolarSphere tSurface tSpline tVolume
+%token tLines tSurfaces
%token tCharacteristic tLength tParametric tElliptic tRefineMesh tAdaptMesh
%token tRelocateMesh
-%token tPlane tRuled tTransfinite tComplex tPhysical tCompound tPeriodic
+%token tPlane tRuled tTransfinite tComplex tPhysical tCompound tPeriodic tTransform
%token tUsing tPlugin tDegenerated tRecursive
%token tRotate tTranslate tSymmetry tDilate tExtrude tLevelset
%token tRecombine tSmoother tSplit tDelete tCoherence
@@ -4019,99 +4026,201 @@ Constraints :
| tSmoother tSurface ListOfDouble tAFFECT FExpr tEND
{
for(int i = 0; i < List_Nbr($3); i++){
- double d;
- List_Read($3, i, &d);
- int j = (int)d;
- Surface *s = FindSurface(j);
- if(s){
- s->TransfiniteSmoothing = (int)$5;
- }
- else{
- GFace *gf = GModel::current()->getFaceByTag(j);
- if(gf)
- gf->meshAttributes.transfiniteSmoothing = (int)$5;
- else
- yymsg(1, "Unknown surface %d", (int)$5);
- }
+ double d;
+ List_Read($3, i, &d);
+ int j = (int)d;
+ Surface *s = FindSurface(j);
+ if(s){
+ s->TransfiniteSmoothing = (int)$5;
+ }
+ else{
+ GFace *gf = GModel::current()->getFaceByTag(j);
+ if(gf)
+ gf->meshAttributes.transfiniteSmoothing = (int)$5;
+ else
+ yymsg(1, "Unknown surface %d", (int)$5);
+ }
}
List_Delete($3);
}
- | tPeriodic tLine ListOfDouble tAFFECT ListOfDouble tEND
- {
- if(List_Nbr($5) != List_Nbr($3)){
- yymsg(0, "Number of master (%d) different from number of slave (%d) lines",
- List_Nbr($5), List_Nbr($3));
- }
- else{
- for(int i = 0; i < List_Nbr($3); i++){
- double d_master, d_slave;
- List_Read($5, i, &d_master);
- List_Read($3, i, &d_slave);
- int j_master = (int)d_master;
- int j_slave = (int)d_slave;
- Curve *c_slave = FindCurve(abs(j_slave));
- if(c_slave){
- GModel::current()->getGEOInternals()->periodicEdges[j_slave] = j_master;
- }
- else{
- GEdge *ge = GModel::current()->getEdgeByTag(abs(j_slave));
- if(ge) ge->setMeshMaster(j_master);
- else yymsg(0, "Unknown line %d", j_slave);
- }
- }
- }
- List_Delete($3);
- List_Delete($5);
- }
- | tPeriodic tSurface FExpr '{' RecursiveListOfDouble '}' tAFFECT FExpr
- '{' RecursiveListOfDouble '}' PeriodicTransform tEND
- {
- if (List_Nbr($5) != List_Nbr($10)){
- yymsg(0, "Number of master surface edges (%d) different from number of "
- "slave (%d) edges", List_Nbr($10), List_Nbr($5));
- }
- else{
- int j_master = (int)$8;
- int j_slave = (int)$3;
- Surface *s_slave = FindSurface(abs(j_slave));
- std::vector<double> transfo;
- for(int i = 0; i < List_Nbr($12); i++){
- double d;
- List_Read($12, i, &d);
- transfo.push_back(d);
+ | tPeriodic tLine '{' RecursiveListOfDouble '}' tAFFECT '{' RecursiveListOfDouble '}' tTransform PeriodicTransform tEND
+ {
+ if (List_Nbr($4) != List_Nbr($8)){
+ yymsg(0, "Number of master lines (%d) different from number of "
+ "slaves (%d) ", List_Nbr($8), List_Nbr($4));
+ }
+ else{
+ if (List_Nbr($11) < 12){
+ yymsg(0, "Affine transformation requires at least 12 entries");
+ }
+ else {
+
+ std::vector<double> transfo(16,0);
+ for(int i = 0; i < List_Nbr($11); i++) List_Read($11, i, &transfo[i]);
+
+ for(int i = 0; i < List_Nbr($4); i++){
+ double d_master, d_slave;
+ List_Read($8, i, &d_master);
+ List_Read($4, i, &d_slave);
+ int j_master = (int)d_master;
+ int j_slave = (int)d_slave;
+
+ addPeriodicEdge(j_slave,j_master,transfo);
+ }
+ }
+ }
+ List_Delete($8);
+ List_Delete($4);
+ }
+ | tPeriodic tSurface '{' RecursiveListOfDouble '}' tAFFECT '{' RecursiveListOfDouble '}' tTransform PeriodicTransform tEND
+ {
+ if (List_Nbr($4) != List_Nbr($8)){
+ yymsg(0, "Number of master faces (%d) different from number of "
+ "slaves (%d) ", List_Nbr($8), List_Nbr($4));
+ }
+ else{
+ if (List_Nbr($11) < 12){
+ yymsg(0, "Affine transformation requires at least 12 entries");
+ }
+ else {
+ std::vector<double> transfo(16,0);
+ for(int i = 0; i < List_Nbr($11); i++) List_Read($11, i, &transfo[i]);
+
+ for(int i = 0; i < List_Nbr($4); i++){
+ double d_master, d_slave;
+ List_Read($8, i, &d_master);
+ List_Read($4, i, &d_slave);
+ addPeriodicFace(d_slave,d_master,transfo);
+ }
+ }
+ }
+ List_Delete($4);
+ List_Delete($8);
+ }
+ | tPeriodic tLine '{' RecursiveListOfDouble '}' tAFFECT '{' RecursiveListOfDouble '}' tRotate '{' VExpr ',' VExpr ',' FExpr '}' tEND
+ {
+ if (List_Nbr($4) != List_Nbr($8)){
+ yymsg(0, "Number of master edges (%d) different from number of "
+ "slaves (%d) ", List_Nbr($8), List_Nbr($4));
+ }
+ else{
+ SPoint3 origin($12[0],$12[1],$12[2]);
+ SPoint3 axis($14[0],$14[1],$14[2]);
+ double angle($16);
+ SPoint3 translation(0,0,0);
+
+ std::vector<double> transfo;
+ computeAffineTransformation(origin,axis,angle,translation,transfo);
+
+ for(int i = 0; i < List_Nbr($4); i++){
+ double d_master, d_slave;
+ List_Read($8, i, &d_master);
+ List_Read($4, i, &d_slave);
+ addPeriodicEdge(d_slave,d_master,transfo);
+ }
+ }
+ List_Delete($4);
+ List_Delete($8);
+ }
+ | tPeriodic tSurface '{' RecursiveListOfDouble '}' tAFFECT '{' RecursiveListOfDouble '}' tRotate '{' VExpr ',' VExpr ',' FExpr '}' tEND
+ {
+ if (List_Nbr($4) != List_Nbr($8)){
+ yymsg(0, "Number of master faces (%d) different from number of "
+ "slaves (%d) ", List_Nbr($8), List_Nbr($4));
+ }
+ else{
+
+ SPoint3 origin($12[0],$12[1],$12[2]);
+ SPoint3 axis($14[0],$14[1],$14[2]);
+ double angle($16);
+ SPoint3 translation(0,0,0);
+
+ std::vector<double> transfo;
+ computeAffineTransformation(origin,axis,angle,translation,transfo);
+
+ for(int i = 0; i < List_Nbr($4); i++){
+ double d_master, d_slave;
+ List_Read($8, i, &d_master);
+ List_Read($4, i, &d_slave);
+ addPeriodicFace(d_slave,d_master,transfo);
+ }
+ }
+ List_Delete($4);
+ List_Delete($8);
+ }
+ | tPeriodic tLine '{' RecursiveListOfDouble '}' tAFFECT '{' RecursiveListOfDouble '}' tTranslate VExpr tEND
+ {
+ if (List_Nbr($4) != List_Nbr($8)){
+ yymsg(0, "Number of master edges (%d) different from number of "
+ "slaves (%d) ", List_Nbr($8), List_Nbr($4));
+ }
+ else{
+ SPoint3 origin(0,0,0);
+ SPoint3 axis(0,0,0);
+ double angle(0);
+ SPoint3 translation($11[0],$11[1],$11[2]);
+
+ std::vector<double> transfo;
+ computeAffineTransformation(origin,axis,angle,translation,transfo);
+
+ for(int i = 0; i < List_Nbr($4); i++){
+ double d_master, d_slave;
+ List_Read($8, i, &d_master);
+ List_Read($4, i, &d_slave);
+ addPeriodicEdge(d_slave,d_master,transfo);
+ }
+ }
+ List_Delete($4);
+ List_Delete($8);
+ }
+ | tPeriodic tSurface '{' RecursiveListOfDouble '}' tAFFECT '{' RecursiveListOfDouble '}' tTranslate VExpr tEND
+ {
+ if (List_Nbr($4) != List_Nbr($8)){
+ yymsg(0, "Number of master faces (%d) different from number of "
+ "slaves (%d) ", List_Nbr($8), List_Nbr($4));
+ }
+ else{
+ SPoint3 origin(0,0,0);
+ SPoint3 axis(0,0,0);
+ double angle(0);
+ SPoint3 translation($11[0],$11[1],$11[2]);
+
+ std::vector<double> transfo;
+ computeAffineTransformation(origin,axis,angle,translation,transfo);
+
+ for(int i = 0; i < List_Nbr($4); i++){
+ double d_master, d_slave;
+ List_Read($8, i, &d_master);
+ List_Read($4, i, &d_slave);
+ addPeriodicFace(d_slave,d_master,transfo);
+ }
+ }
+ List_Delete($4);
+ List_Delete($8);
+ }
+ | tPeriodic tSurface FExpr '{' RecursiveListOfDouble '}' tAFFECT FExpr '{' RecursiveListOfDouble '}' tEND
+ {
+ if (List_Nbr($5) != List_Nbr($10)){
+ yymsg(0, "Number of master surface edges (%d) different from number of "
+ "slave (%d) edges", List_Nbr($10), List_Nbr($5));
+ }
+ else{
+ int j_master = (int)$8;
+ int j_slave = (int)$3;
+
+ std::map<int,int> edgeCounterParts;
+ for (int i = 0; i < List_Nbr($5); i++){
+ double ds,dm;
+ List_Read($5,i,&ds);
+ List_Read($10,i,&dm);
+ edgeCounterParts[(int) ds] = (int) dm;
+ std::cout << "edge " << ds << " to " << dm << std::endl;
+ }
+ addPeriodicFace(j_slave,j_master,edgeCounterParts);
+ }
+ List_Delete($5);
+ List_Delete($10);
}
- if(s_slave){
- GModel::current()->getGEOInternals()->periodicFaces[j_slave] = j_master;
- for (int i = 0; i < List_Nbr($5); i++){
- double dm, ds;
- List_Read($5, i, &ds);
- List_Read($10, i, &dm);
- GModel::current()->getGEOInternals()->periodicEdges[(int)ds] = (int)dm;
- s_slave->edgeCounterparts[(int)ds] = (int)dm;
- s_slave->affineTransform = transfo;
- }
- }
- else{
- GFace *gf = GModel::current()->getFaceByTag(abs(j_slave));
- if(gf){
- gf->setMeshMaster(j_master);
- for (int i = 0; i < List_Nbr($5); i++){
- double dm, ds;
- List_Read($5, i, &ds);
- List_Read($10, i, &dm);
- gf->edgeCounterparts[(int)ds] = (int)dm;
- gf->affineTransform = transfo;
- GEdge *ges = GModel::current()->getEdgeByTag(abs((int)ds));
- ges->setMeshMaster((int)dm);
- }
- }
- else yymsg(0, "Unknown surface %d", j_slave);
- }
- }
- List_Delete($5);
- List_Delete($10);
- List_Delete($12);
- }
| tPoint '{' RecursiveListOfDouble '}' tIn tSurface '{' FExpr '}' tEND
{
Surface *s = FindSurface((int)$8);
@@ -5882,3 +5991,124 @@ void yymsg(int level, const char *fmt, ...)
Msg::Warning("'%s', line %d : %s", gmsh_yyname.c_str(), gmsh_yylineno - 1, tmp);
}
}
+
+
+void addPeriodicFace(int iTarget,int iSource,
+ const std::vector<double>& affineTransform) {
+
+ Surface *target = FindSurface(abs(iTarget));
+
+ if (target) {
+ GEO_Internals::MasterFace& mf =
+ GModel::current()->getGEOInternals()->periodicFaces[iTarget];
+ mf.tag = iSource;
+ mf.edgeCounterparts.clear();
+ mf.affineTransform = affineTransform;
+ }
+ else{
+ GFace *target = GModel::current()->getFaceByTag(abs(iTarget));
+ GFace *source = GModel::current()->getFaceByTag(abs(iSource));
+ if (!target) Msg::Error("Could not find edge %d for periodic copy from %d",
+ iTarget,iSource);
+ target->setMeshMaster(source,affineTransform);
+ }
+}
+
+
+void addPeriodicFace(int iTarget,int iSource,
+ const std::map<int,int>& edgeCounterparts) {
+
+ Surface *target = FindSurface(abs(iTarget));
+
+ std::cout << "Encoding periodic connection between " << iTarget << " and " << iSource << std::endl;
+ std::map<int,int>::const_iterator sIter = edgeCounterparts.begin();
+ for (;sIter!=edgeCounterparts.end();++sIter) {
+ std::cout << sIter->first << " - " << sIter->second << std::endl;
+ }
+
+
+ if (target) {
+ GEO_Internals::MasterFace& mf =
+ GModel::current()->getGEOInternals()->periodicFaces[iTarget];
+ mf.tag = iSource;
+ mf.edgeCounterparts = edgeCounterparts;
+ mf.affineTransform.clear();
+ }
+ else{
+ GFace *target = GModel::current()->getFaceByTag(abs(iTarget));
+ GFace *source = GModel::current()->getFaceByTag(abs(iSource));
+ if (!target || !source) Msg::Error("Could not find surface %d or %d for periodic copy",
+ iTarget,iSource);
+ target->setMeshMaster(source,edgeCounterparts);
+ }
+}
+
+void addPeriodicEdge(int iTarget,int iSource,
+ const std::vector<double>& affineTransform) {
+
+ Curve *target = FindCurve(abs(iTarget));
+ if (target) {
+ GEO_Internals::MasterEdge& me = GModel::current()->getGEOInternals()->periodicEdges[iTarget];
+ me.tag = iSource;
+ me.affineTransform = affineTransform;
+
+ GEO_Internals::MasterEdge& test = GModel::current()->getGEOInternals()->periodicEdges[iTarget];
+ }
+ else{
+ GEdge *target = GModel::current()->getEdgeByTag(abs(iTarget));
+ GEdge *source = GModel::current()->getEdgeByTag(abs(iSource));
+ if (!target || !source) Msg::Error("Could not find surface %d or %d for periodic copy",
+ iTarget,iSource);
+ if (affineTransform.size() == 16) {
+ target->setMeshMaster(source,affineTransform);
+ }
+ else {
+ target->setMeshMaster(source,iSource*iTarget < 0 ? -1:1);
+ }
+ }
+}
+
+void computeAffineTransformation(SPoint3& origin,
+ SPoint3& axis,
+ double angle,
+ SPoint3& translation,
+ std::vector<double>& tfo) {
+
+ tfo.resize(16,0.0);
+
+ angle *= acos(-1.)/180.;
+
+ double ca = cos(angle);
+ double sa = sin(angle);
+
+ double ux = axis.x();
+ double uy = axis.y();
+ double uz = axis.z();
+
+ tfo.resize(16);
+
+ tfo[0*4+0] = ca + ux*ux*(1.-ca);
+ tfo[0*4+1] = ux*uy*(1.-ca) - uz * sa;
+ tfo[0*4+2] = ux*uz*(1.-ca) + uy * sa;
+
+ tfo[1*4+0] = ux*uy*(1.-ca) + uz * sa;
+ tfo[1*4+1] = ca + uy*uy*(1.-ca);
+ tfo[1*4+2] = uy*uz*(1.-ca) - ux * sa;
+
+ tfo[2*4+0] = ux*uz*(1.-ca) - uy * sa;
+ tfo[2*4+1] = uy*uz*(1.-ca) + ux * sa;
+ tfo[2*4+2] = ca + uz*uz*(1.-ca);
+
+ int idx = 0;
+ for (size_t i=0;i<3;i++) {
+ int tIdx = i*4+3;
+ tfo[tIdx] = translation[i];
+ for (int j=0;j<3;j++,idx++) {
+ tfo[tIdx] -= tfo[idx] * origin[j];
+ }
+ idx++;
+ }
+
+ for (int i=0;i<4;i++) tfo[12+i] = 0;
+ tfo[15] = 1;
+}
diff --git a/Parser/Gmsh.yy.cpp b/Parser/Gmsh.yy.cpp
deleted file mode 100644
index 6090e94da1895a9f916016f6d0f7c9fda078112e..0000000000000000000000000000000000000000
--- a/Parser/Gmsh.yy.cpp
+++ /dev/null
@@ -1,3509 +0,0 @@
-#line 2 "Gmsh.yy.cpp"
-
-#line 4 "Gmsh.yy.cpp"
-
-#define YY_INT_ALIGNED short int
-
-/* A lexical scanner generated by flex */
-
-#define yy_create_buffer gmsh_yy_create_buffer
-#define yy_delete_buffer gmsh_yy_delete_buffer
-#define yy_flex_debug gmsh_yy_flex_debug
-#define yy_init_buffer gmsh_yy_init_buffer
-#define yy_flush_buffer gmsh_yy_flush_buffer
-#define yy_load_buffer_state gmsh_yy_load_buffer_state
-#define yy_switch_to_buffer gmsh_yy_switch_to_buffer
-#define yyin gmsh_yyin
-#define yyleng gmsh_yyleng
-#define yylex gmsh_yylex
-#define yylineno gmsh_yylineno
-#define yyout gmsh_yyout
-#define yyrestart gmsh_yyrestart
-#define yytext gmsh_yytext
-#define yywrap gmsh_yywrap
-#define yyalloc gmsh_yyalloc
-#define yyrealloc gmsh_yyrealloc
-#define yyfree gmsh_yyfree
-
-#define FLEX_SCANNER
-#define YY_FLEX_MAJOR_VERSION 2
-#define YY_FLEX_MINOR_VERSION 5
-#define YY_FLEX_SUBMINOR_VERSION 37
-#if YY_FLEX_SUBMINOR_VERSION > 0
-#define FLEX_BETA
-#endif
-
-/* First, we deal with platform-specific or compiler-specific issues. */
-
-/* begin standard C headers. */
-#include <stdio.h>
-#include <string.h>
-#include <errno.h>
-#include <stdlib.h>
-
-/* end standard C headers. */
-
-/* flex integer type definitions */
-
-#ifndef FLEXINT_H
-#define FLEXINT_H
-
-/* C99 systems have <inttypes.h>. Non-C99 systems may or may not. */
-
-#if defined (__STDC_VERSION__) && __STDC_VERSION__ >= 199901L
-
-/* C99 says to define __STDC_LIMIT_MACROS before including stdint.h,
- * if you want the limit (max/min) macros for int types.
- */
-#ifndef __STDC_LIMIT_MACROS
-#define __STDC_LIMIT_MACROS 1
-#endif
-
-#include <inttypes.h>
-typedef int8_t flex_int8_t;
-typedef uint8_t flex_uint8_t;
-typedef int16_t flex_int16_t;
-typedef uint16_t flex_uint16_t;
-typedef int32_t flex_int32_t;
-typedef uint32_t flex_uint32_t;
-#else
-typedef signed char flex_int8_t;
-typedef short int flex_int16_t;
-typedef int flex_int32_t;
-typedef unsigned char flex_uint8_t;
-typedef unsigned short int flex_uint16_t;
-typedef unsigned int flex_uint32_t;
-
-/* Limits of integral types. */
-#ifndef INT8_MIN
-#define INT8_MIN (-128)
-#endif
-#ifndef INT16_MIN
-#define INT16_MIN (-32767-1)
-#endif
-#ifndef INT32_MIN
-#define INT32_MIN (-2147483647-1)
-#endif
-#ifndef INT8_MAX
-#define INT8_MAX (127)
-#endif
-#ifndef INT16_MAX
-#define INT16_MAX (32767)
-#endif
-#ifndef INT32_MAX
-#define INT32_MAX (2147483647)
-#endif
-#ifndef UINT8_MAX
-#define UINT8_MAX (255U)
-#endif
-#ifndef UINT16_MAX
-#define UINT16_MAX (65535U)
-#endif
-#ifndef UINT32_MAX
-#define UINT32_MAX (4294967295U)
-#endif
-
-#endif /* ! C99 */
-
-#endif /* ! FLEXINT_H */
-
-#ifdef __cplusplus
-
-/* The "const" storage-class-modifier is valid. */
-#define YY_USE_CONST
-
-#else /* ! __cplusplus */
-
-/* C99 requires __STDC__ to be defined as 1. */
-#if defined (__STDC__)
-
-#define YY_USE_CONST
-
-#endif /* defined (__STDC__) */
-#endif /* ! __cplusplus */
-
-#ifdef YY_USE_CONST
-#define yyconst const
-#else
-#define yyconst
-#endif
-
-/* Returned upon end-of-file. */
-#define YY_NULL 0
-
-/* Promotes a possibly negative, possibly signed char to an unsigned
- * integer for use as an array index. If the signed char is negative,
- * we want to instead treat it as an 8-bit unsigned char, hence the
- * double cast.
- */
-#define YY_SC_TO_UI(c) ((unsigned int) (unsigned char) c)
-
-/* Enter a start condition. This macro really ought to take a parameter,
- * but we do it the disgusting crufty way forced on us by the ()-less
- * definition of BEGIN.
- */
-#define BEGIN (yy_start) = 1 + 2 *
-
-/* Translate the current start state into a value that can be later handed
- * to BEGIN to return to the state. The YYSTATE alias is for lex
- * compatibility.
- */
-#define YY_START (((yy_start) - 1) / 2)
-#define YYSTATE YY_START
-
-/* Action number for EOF rule of a given start state. */
-#define YY_STATE_EOF(state) (YY_END_OF_BUFFER + state + 1)
-
-/* Special action meaning "start processing a new file". */
-#define YY_NEW_FILE gmsh_yyrestart(gmsh_yyin )
-
-#define YY_END_OF_BUFFER_CHAR 0
-
-/* Size of default input buffer. */
-#ifndef YY_BUF_SIZE
-#define YY_BUF_SIZE 16384
-#endif
-
-/* The state buf must be large enough to hold one state per character in the main buffer.
- */
-#define YY_STATE_BUF_SIZE ((YY_BUF_SIZE + 2) * sizeof(yy_state_type))
-
-#ifndef YY_TYPEDEF_YY_BUFFER_STATE
-#define YY_TYPEDEF_YY_BUFFER_STATE
-typedef struct yy_buffer_state *YY_BUFFER_STATE;
-#endif
-
-#ifndef YY_TYPEDEF_YY_SIZE_T
-#define YY_TYPEDEF_YY_SIZE_T
-typedef size_t yy_size_t;
-#endif
-
-extern yy_size_t gmsh_yyleng;
-
-extern FILE *gmsh_yyin, *gmsh_yyout;
-
-#define EOB_ACT_CONTINUE_SCAN 0
-#define EOB_ACT_END_OF_FILE 1
-#define EOB_ACT_LAST_MATCH 2
-
- #define YY_LESS_LINENO(n)
-
-/* Return all but the first "n" matched characters back to the input stream. */
-#define yyless(n) \
- do \
- { \
- /* Undo effects of setting up gmsh_yytext. */ \
- int yyless_macro_arg = (n); \
- YY_LESS_LINENO(yyless_macro_arg);\
- *yy_cp = (yy_hold_char); \
- YY_RESTORE_YY_MORE_OFFSET \
- (yy_c_buf_p) = yy_cp = yy_bp + yyless_macro_arg - YY_MORE_ADJ; \
- YY_DO_BEFORE_ACTION; /* set up gmsh_yytext again */ \
- } \
- while ( 0 )
-
-#define unput(c) yyunput( c, (yytext_ptr) )
-
-#ifndef YY_STRUCT_YY_BUFFER_STATE
-#define YY_STRUCT_YY_BUFFER_STATE
-struct yy_buffer_state
- {
- FILE *yy_input_file;
-
- char *yy_ch_buf; /* input buffer */
- char *yy_buf_pos; /* current position in input buffer */
-
- /* Size of input buffer in bytes, not including room for EOB
- * characters.
- */
- yy_size_t yy_buf_size;
-
- /* Number of characters read into yy_ch_buf, not including EOB
- * characters.
- */
- yy_size_t yy_n_chars;
-
- /* Whether we "own" the buffer - i.e., we know we created it,
- * and can realloc() it to grow it, and should free() it to
- * delete it.
- */
- int yy_is_our_buffer;
-
- /* Whether this is an "interactive" input source; if so, and
- * if we're using stdio for input, then we want to use getc()
- * instead of fread(), to make sure we stop fetching input after
- * each newline.
- */
- int yy_is_interactive;
-
- /* Whether we're considered to be at the beginning of a line.
- * If so, '^' rules will be active on the next match, otherwise
- * not.
- */
- int yy_at_bol;
-
- int yy_bs_lineno; /**< The line count. */
- int yy_bs_column; /**< The column count. */
-
- /* Whether to try to fill the input buffer when we reach the
- * end of it.
- */
- int yy_fill_buffer;
-
- int yy_buffer_status;
-
-#define YY_BUFFER_NEW 0
-#define YY_BUFFER_NORMAL 1
- /* When an EOF's been seen but there's still some text to process
- * then we mark the buffer as YY_EOF_PENDING, to indicate that we
- * shouldn't try reading from the input source any more. We might
- * still have a bunch of tokens to match, though, because of
- * possible backing-up.
- *
- * When we actually see the EOF, we change the status to "new"
- * (via gmsh_yyrestart()), so that the user can continue scanning by
- * just pointing gmsh_yyin at a new input file.
- */
-#define YY_BUFFER_EOF_PENDING 2
-
- };
-#endif /* !YY_STRUCT_YY_BUFFER_STATE */
-
-/* Stack of input buffers. */
-static size_t yy_buffer_stack_top = 0; /**< index of top of stack. */
-static size_t yy_buffer_stack_max = 0; /**< capacity of stack. */
-static YY_BUFFER_STATE * yy_buffer_stack = 0; /**< Stack as an array. */
-
-/* We provide macros for accessing buffer states in case in the
- * future we want to put the buffer states in a more general
- * "scanner state".
- *
- * Returns the top of the stack, or NULL.
- */
-#define YY_CURRENT_BUFFER ( (yy_buffer_stack) \
- ? (yy_buffer_stack)[(yy_buffer_stack_top)] \
- : NULL)
-
-/* Same as previous macro, but useful when we know that the buffer stack is not
- * NULL or when we need an lvalue. For internal use only.
- */
-#define YY_CURRENT_BUFFER_LVALUE (yy_buffer_stack)[(yy_buffer_stack_top)]
-
-/* yy_hold_char holds the character lost when gmsh_yytext is formed. */
-static char yy_hold_char;
-static yy_size_t yy_n_chars; /* number of characters read into yy_ch_buf */
-yy_size_t gmsh_yyleng;
-
-/* Points to current character in buffer. */
-static char *yy_c_buf_p = (char *) 0;
-static int yy_init = 0; /* whether we need to initialize */
-static int yy_start = 0; /* start state number */
-
-/* Flag which is used to allow gmsh_yywrap()'s to do buffer switches
- * instead of setting up a fresh gmsh_yyin. A bit of a hack ...
- */
-static int yy_did_buffer_switch_on_eof;
-
-void gmsh_yyrestart (FILE *input_file );
-void gmsh_yy_switch_to_buffer (YY_BUFFER_STATE new_buffer );
-YY_BUFFER_STATE gmsh_yy_create_buffer (FILE *file,int size );
-void gmsh_yy_delete_buffer (YY_BUFFER_STATE b );
-void gmsh_yy_flush_buffer (YY_BUFFER_STATE b );
-void gmsh_yypush_buffer_state (YY_BUFFER_STATE new_buffer );
-void gmsh_yypop_buffer_state (void );
-
-static void gmsh_yyensure_buffer_stack (void );
-static void gmsh_yy_load_buffer_state (void );
-static void gmsh_yy_init_buffer (YY_BUFFER_STATE b,FILE *file );
-
-#define YY_FLUSH_BUFFER gmsh_yy_flush_buffer(YY_CURRENT_BUFFER )
-
-YY_BUFFER_STATE gmsh_yy_scan_buffer (char *base,yy_size_t size );
-YY_BUFFER_STATE gmsh_yy_scan_string (yyconst char *yy_str );
-YY_BUFFER_STATE gmsh_yy_scan_bytes (yyconst char *bytes,yy_size_t len );
-
-void *gmsh_yyalloc (yy_size_t );
-void *gmsh_yyrealloc (void *,yy_size_t );
-void gmsh_yyfree (void * );
-
-#define yy_new_buffer gmsh_yy_create_buffer
-
-#define yy_set_interactive(is_interactive) \
- { \
- if ( ! YY_CURRENT_BUFFER ){ \
- gmsh_yyensure_buffer_stack (); \
- YY_CURRENT_BUFFER_LVALUE = \
- gmsh_yy_create_buffer(gmsh_yyin,YY_BUF_SIZE ); \
- } \
- YY_CURRENT_BUFFER_LVALUE->yy_is_interactive = is_interactive; \
- }
-
-#define yy_set_bol(at_bol) \
- { \
- if ( ! YY_CURRENT_BUFFER ){\
- gmsh_yyensure_buffer_stack (); \
- YY_CURRENT_BUFFER_LVALUE = \
- gmsh_yy_create_buffer(gmsh_yyin,YY_BUF_SIZE ); \
- } \
- YY_CURRENT_BUFFER_LVALUE->yy_at_bol = at_bol; \
- }
-
-#define YY_AT_BOL() (YY_CURRENT_BUFFER_LVALUE->yy_at_bol)
-
-/* Begin user sect3 */
-
-typedef unsigned char YY_CHAR;
-
-FILE *gmsh_yyin = (FILE *) 0, *gmsh_yyout = (FILE *) 0;
-
-typedef int yy_state_type;
-
-extern int gmsh_yylineno;
-
-int gmsh_yylineno = 1;
-
-extern char *gmsh_yytext;
-#define yytext_ptr gmsh_yytext
-
-static yy_state_type yy_get_previous_state (void );
-static yy_state_type yy_try_NUL_trans (yy_state_type current_state );
-static int yy_get_next_buffer (void );
-static void yy_fatal_error (yyconst char msg[] );
-
-/* Done after the current pattern has been matched and before the
- * corresponding action - sets up gmsh_yytext.
- */
-#define YY_DO_BEFORE_ACTION \
- (yytext_ptr) = yy_bp; \
- gmsh_yyleng = (size_t) (yy_cp - yy_bp); \
- (yy_hold_char) = *yy_cp; \
- *yy_cp = '\0'; \
- (yy_c_buf_p) = yy_cp;
-
-#define YY_NUM_RULES 200
-#define YY_END_OF_BUFFER 201
-/* This struct is not used in this scanner,
- but its presence is necessary. */
-struct yy_trans_info
- {
- flex_int32_t yy_verify;
- flex_int32_t yy_nxt;
- };
-static yyconst flex_int16_t yy_accept[1002] =
- { 0,
- 0, 0, 201, 199, 1, 1, 199, 5, 199, 6,
- 199, 199, 199, 199, 199, 194, 21, 2, 199, 16,
- 199, 198, 198, 198, 198, 198, 198, 198, 198, 198,
- 198, 198, 198, 198, 198, 198, 198, 198, 198, 198,
- 198, 198, 198, 198, 199, 28, 24, 19, 25, 17,
- 26, 18, 0, 196, 3, 4, 20, 195, 194, 0,
- 29, 27, 30, 198, 198, 198, 198, 198, 198, 198,
- 198, 198, 198, 198, 198, 198, 198, 198, 198, 198,
- 198, 198, 198, 198, 198, 198, 198, 198, 198, 198,
- 198, 198, 198, 198, 198, 198, 198, 198, 198, 198,
-
- 198, 198, 198, 109, 108, 198, 198, 198, 198, 198,
- 198, 198, 198, 198, 198, 198, 198, 198, 198, 198,
- 198, 133, 198, 198, 198, 198, 198, 198, 198, 198,
- 198, 198, 198, 198, 198, 198, 198, 198, 198, 198,
- 180, 181, 198, 198, 198, 198, 198, 198, 198, 198,
- 198, 23, 22, 0, 195, 0, 0, 197, 198, 198,
- 198, 198, 198, 198, 198, 198, 198, 198, 198, 198,
- 198, 198, 198, 198, 198, 198, 198, 198, 198, 54,
- 65, 198, 198, 198, 198, 198, 198, 198, 198, 198,
- 198, 198, 198, 76, 198, 198, 198, 198, 198, 198,
-
- 198, 198, 94, 198, 198, 198, 198, 198, 198, 198,
- 198, 198, 198, 198, 198, 198, 198, 198, 198, 116,
- 198, 198, 198, 198, 198, 198, 198, 198, 198, 198,
- 198, 198, 198, 198, 198, 198, 198, 198, 198, 198,
- 198, 198, 198, 198, 198, 198, 198, 198, 198, 198,
- 162, 198, 198, 198, 198, 198, 168, 198, 198, 198,
- 198, 188, 198, 198, 198, 198, 198, 198, 198, 198,
- 0, 196, 0, 0, 195, 198, 32, 198, 198, 198,
- 198, 198, 37, 39, 198, 198, 198, 198, 62, 198,
- 47, 198, 198, 198, 198, 198, 198, 198, 198, 198,
-
- 53, 198, 198, 198, 198, 198, 198, 74, 198, 198,
- 198, 198, 198, 198, 85, 198, 87, 198, 198, 89,
- 198, 198, 93, 198, 198, 198, 198, 198, 198, 198,
- 104, 105, 198, 198, 198, 198, 198, 198, 198, 198,
- 114, 115, 198, 198, 198, 198, 198, 198, 198, 198,
- 198, 198, 198, 198, 198, 198, 198, 198, 198, 198,
- 198, 150, 198, 198, 198, 198, 198, 198, 198, 198,
- 198, 198, 198, 198, 198, 198, 198, 177, 163, 198,
- 198, 198, 198, 161, 198, 198, 198, 198, 198, 198,
- 198, 183, 187, 198, 198, 198, 198, 198, 198, 198,
-
- 10, 15, 9, 8, 198, 12, 14, 0, 195, 31,
- 198, 35, 198, 198, 198, 41, 198, 43, 198, 198,
- 198, 198, 198, 198, 198, 57, 198, 198, 198, 198,
- 198, 198, 198, 198, 198, 198, 198, 198, 198, 82,
- 84, 198, 198, 198, 88, 198, 198, 92, 198, 198,
- 198, 198, 198, 198, 198, 198, 107, 198, 198, 112,
- 198, 198, 198, 117, 198, 198, 120, 198, 198, 198,
- 198, 127, 198, 130, 198, 198, 198, 134, 198, 135,
- 198, 198, 198, 198, 198, 198, 198, 198, 198, 198,
- 148, 149, 198, 198, 198, 198, 198, 198, 198, 198,
-
- 198, 166, 198, 198, 198, 198, 198, 198, 198, 198,
- 198, 198, 198, 189, 198, 198, 198, 192, 198, 11,
- 198, 13, 198, 198, 34, 38, 40, 198, 44, 198,
- 198, 198, 48, 198, 198, 198, 198, 198, 198, 198,
- 198, 198, 198, 198, 198, 71, 73, 198, 198, 81,
- 198, 86, 198, 198, 198, 198, 198, 198, 96, 198,
- 198, 198, 198, 198, 198, 198, 118, 113, 198, 198,
- 198, 125, 198, 122, 198, 198, 198, 198, 198, 139,
- 198, 138, 198, 198, 198, 198, 198, 153, 198, 147,
- 198, 198, 198, 198, 198, 198, 198, 164, 165, 198,
-
- 170, 198, 175, 198, 198, 198, 198, 198, 198, 198,
- 198, 198, 198, 198, 198, 193, 7, 198, 198, 42,
- 45, 198, 198, 198, 198, 198, 198, 52, 56, 198,
- 198, 198, 198, 198, 198, 198, 198, 198, 78, 198,
- 198, 198, 79, 198, 198, 198, 198, 198, 198, 198,
- 198, 198, 198, 198, 198, 198, 198, 198, 198, 198,
- 198, 198, 198, 198, 198, 198, 198, 198, 198, 198,
- 198, 198, 198, 154, 198, 198, 198, 198, 198, 198,
- 198, 169, 198, 174, 198, 198, 198, 167, 198, 198,
- 198, 198, 198, 198, 198, 198, 198, 198, 198, 198,
-
- 198, 198, 198, 198, 59, 198, 198, 198, 198, 198,
- 198, 198, 198, 77, 80, 198, 198, 198, 95, 198,
- 198, 198, 198, 198, 99, 198, 106, 198, 198, 119,
- 123, 124, 198, 126, 198, 198, 198, 131, 132, 198,
- 198, 198, 198, 198, 198, 198, 198, 198, 198, 198,
- 198, 198, 157, 198, 198, 156, 198, 198, 198, 198,
- 178, 198, 198, 198, 198, 198, 198, 198, 33, 198,
- 198, 198, 198, 50, 198, 198, 198, 198, 198, 198,
- 198, 198, 198, 198, 83, 198, 198, 198, 198, 198,
- 98, 97, 198, 198, 110, 198, 198, 129, 198, 198,
-
- 198, 198, 198, 198, 198, 144, 146, 198, 198, 198,
- 198, 158, 198, 159, 176, 172, 198, 198, 179, 198,
- 198, 198, 198, 186, 198, 198, 198, 61, 198, 51,
- 58, 198, 198, 198, 67, 198, 198, 198, 198, 90,
- 198, 198, 198, 198, 198, 198, 198, 198, 136, 198,
- 198, 142, 198, 198, 198, 151, 198, 198, 75, 198,
- 198, 171, 198, 198, 198, 198, 198, 198, 46, 198,
- 60, 49, 198, 198, 198, 198, 198, 72, 198, 198,
- 198, 198, 198, 198, 198, 198, 137, 198, 198, 143,
- 198, 198, 198, 198, 173, 198, 190, 198, 184, 198,
-
- 198, 198, 198, 198, 198, 69, 70, 198, 198, 198,
- 198, 198, 198, 198, 128, 198, 198, 198, 152, 198,
- 160, 198, 198, 198, 198, 198, 198, 198, 198, 198,
- 198, 198, 198, 198, 198, 121, 198, 198, 198, 198,
- 198, 185, 198, 198, 55, 63, 198, 68, 198, 198,
- 198, 198, 198, 198, 198, 198, 145, 155, 182, 198,
- 198, 198, 198, 91, 198, 198, 198, 198, 198, 140,
- 198, 198, 36, 198, 66, 198, 198, 198, 198, 198,
- 198, 191, 198, 198, 198, 198, 198, 198, 141, 198,
- 100, 101, 102, 103, 198, 198, 198, 198, 111, 64,
-
- 0
- } ;
-
-static yyconst flex_int32_t yy_ec[256] =
- { 0,
- 1, 1, 1, 1, 1, 1, 1, 1, 2, 3,
- 1, 2, 2, 1, 1, 1, 1, 1, 1, 1,
- 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
- 1, 2, 4, 5, 1, 1, 1, 6, 7, 1,
- 1, 8, 9, 1, 10, 11, 12, 13, 14, 15,
- 16, 17, 17, 17, 17, 17, 17, 18, 19, 20,
- 21, 22, 1, 1, 23, 24, 25, 26, 27, 28,
- 29, 30, 31, 32, 33, 34, 35, 36, 37, 38,
- 39, 40, 41, 42, 43, 44, 45, 46, 46, 46,
- 1, 1, 1, 1, 47, 1, 48, 49, 50, 51,
-
- 52, 53, 54, 55, 56, 46, 57, 58, 59, 60,
- 61, 62, 63, 64, 65, 66, 67, 68, 69, 70,
- 71, 72, 1, 73, 1, 1, 1, 1, 1, 1,
- 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
- 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
- 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
- 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
- 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
- 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
- 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
-
- 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
- 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
- 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
- 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
- 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
- 1, 1, 1, 1, 1
- } ;
-
-static yyconst flex_int32_t yy_meta[74] =
- { 0,
- 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
- 1, 1, 2, 2, 2, 2, 2, 1, 1, 1,
- 1, 1, 2, 2, 2, 2, 2, 2, 2, 2,
- 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
- 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
- 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
- 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
- 2, 2, 1
- } ;
-
-static yyconst flex_int16_t yy_base[1003] =
- { 0,
- 0, 0, 1134, 1135, 1135, 1135, 1112, 1135, 1126, 1135,
- 1110, 65, 66, 64, 76, 78, 1135, 1135, 1109, 1108,
- 1107, 49, 49, 56, 69, 64, 79, 50, 80, 90,
- 0, 1067, 96, 101, 80, 94, 108, 1059, 107, 121,
- 163, 105, 1064, 1072, 1050, 1135, 1135, 1135, 1135, 1135,
- 1135, 1135, 1111, 182, 1135, 1135, 1135, 187, 205, 227,
- 1135, 1135, 1135, 0, 1060, 1059, 1071, 1062, 1067, 1060,
- 1067, 1052, 95, 1046, 123, 1056, 1063, 1046, 190, 1042,
- 1056, 1043, 132, 1048, 1057, 1046, 1052, 1038, 1051, 151,
- 1051, 198, 1038, 1037, 1033, 1036, 1054, 1028, 1028, 1041,
-
- 24, 1029, 1048, 0, 1023, 1027, 1016, 163, 145, 1032,
- 1054, 1034, 147, 1032, 1013, 1017, 1028, 1028, 1014, 1013,
- 1005, 0, 180, 135, 1019, 1026, 1013, 209, 159, 1014,
- 1023, 1004, 1008, 1008, 1006, 205, 1002, 1001, 1000, 170,
- 0, 0, 1028, 1002, 991, 187, 1012, 1008, 1002, 999,
- 987, 1135, 1135, 269, 274, 283, 257, 289, 991, 989,
- 991, 1004, 239, 991, 990, 991, 982, 991, 986, 987,
- 985, 985, 978, 991, 255, 979, 54, 975, 967, 982,
- 0, 988, 971, 978, 981, 980, 983, 961, 973, 114,
- 967, 969, 200, 0, 962, 960, 966, 971, 971, 981,
-
- 959, 968, 0, 968, 987, 281, 961, 963, 962, 952,
- 951, 984, 958, 943, 956, 953, 954, 953, 938, 989,
- 955, 937, 939, 944, 931, 962, 947, 937, 942, 945,
- 936, 926, 930, 935, 928, 939, 926, 934, 933, 248,
- 927, 921, 914, 928, 931, 918, 925, 918, 293, 906,
- 919, 912, 920, 915, 914, 903, 295, 915, 908, 916,
- 938, 909, 897, 914, 913, 900, 907, 898, 890, 274,
- 297, 330, 339, 344, 349, 890, 0, 889, 889, 892,
- 896, 903, 0, 935, 893, 892, 895, 895, 0, 878,
- 0, 896, 885, 878, 882, 876, 883, 279, 887, 900,
-
- 0, 870, 883, 874, 873, 866, 865, 0, 868, 868,
- 875, 863, 870, 859, 0, 857, 0, 872, 895, 0,
- 869, 856, 0, 853, 871, 857, 849, 849, 866, 886,
- 0, 0, 854, 845, 870, 845, 843, 843, 840, 847,
- 0, 0, 891, 254, 842, 838, 878, 842, 838, 833,
- 849, 832, 836, 833, 837, 840, 835, 824, 825, 822,
- 308, 0, 828, 822, 825, 834, 819, 818, 815, 829,
- 828, 826, 822, 809, 811, 826, 807, 0, 0, 806,
- 807, 307, 810, 0, 320, 813, 804, 815, 818, 813,
- 829, 0, 0, 840, 791, 803, 795, 806, 804, 798,
-
- 0, 0, 798, 0, 803, 796, 0, 367, 372, 0,
- 818, 807, 786, 790, 789, 0, 788, 0, 783, 790,
- 787, 794, 791, 790, 780, 798, 779, 786, 770, 780,
- 773, 782, 773, 780, 779, 778, 777, 304, 764, 0,
- 0, 776, 761, 774, 0, 754, 765, 0, 766, 355,
- 753, 761, 755, 760, 747, 755, 0, 777, 749, 0,
- 748, 757, 746, 0, 762, 753, 0, 737, 749, 745,
- 754, 0, 755, 0, 751, 751, 751, 0, 740, 0,
- 758, 745, 733, 747, 730, 742, 745, 732, 726, 738,
- 0, 0, 755, 740, 728, 735, 721, 720, 728, 730,
-
- 729, 0, 714, 713, 717, 715, 716, 723, 314, 724,
- 707, 711, 705, 0, 735, 338, 713, 0, 716, 0,
- 713, 0, 714, 709, 0, 0, 749, 711, 0, 702,
- 703, 694, 0, 699, 700, 709, 704, 685, 694, 693,
- 686, 709, 684, 367, 685, 0, 0, 342, 692, 0,
- 691, 0, 694, 689, 696, 682, 223, 719, 0, 692,
- 684, 672, 686, 683, 699, 683, 0, 0, 682, 673,
- 660, 0, 677, 0, 678, 355, 663, 672, 679, 0,
- 664, 0, 669, 343, 668, 688, 656, 0, 669, 0,
- 672, 659, 669, 665, 650, 659, 662, 0, 0, 660,
-
- 0, 656, 0, 660, 657, 661, 650, 655, 642, 654,
- 638, 651, 358, 654, 641, 0, 0, 635, 633, 0,
- 0, 644, 657, 644, 645, 633, 644, 0, 0, 641,
- 643, 634, 628, 662, 626, 619, 619, 636, 0, 618,
- 632, 633, 0, 615, 613, 618, 645, 640, 633, 622,
- 613, 620, 621, 599, 635, 618, 601, 609, 613, 603,
- 605, 612, 594, 596, 609, 600, 602, 379, 608, 595,
- 586, 601, 600, 0, 586, 596, 597, 584, 591, 586,
- 581, 0, 594, 0, 587, 576, 593, 0, 569, 587,
- 574, 578, 569, 36, 40, 67, 74, 120, 158, 147,
-
- 169, 183, 207, 218, 0, 253, 260, 266, 285, 291,
- 312, 336, 338, 0, 0, 346, 341, 353, 0, 373,
- 374, 387, 349, 362, 0, 350, 0, 395, 353, 0,
- 0, 0, 358, 0, 357, 364, 369, 0, 0, 374,
- 377, 380, 370, 372, 367, 382, 383, 371, 402, 372,
- 388, 377, 0, 376, 389, 0, 393, 376, 391, 398,
- 0, 391, 394, 390, 404, 397, 402, 430, 0, 419,
- 396, 399, 403, 0, 389, 409, 410, 403, 403, 401,
- 401, 418, 412, 417, 0, 405, 403, 432, 433, 444,
- 0, 0, 409, 434, 0, 421, 422, 0, 429, 416,
-
- 430, 424, 447, 430, 433, 0, 0, 431, 435, 454,
- 438, 0, 434, 0, 0, 0, 423, 440, 0, 444,
- 430, 444, 430, 0, 436, 436, 429, 0, 435, 0,
- 0, 436, 437, 445, 452, 439, 454, 448, 458, 0,
- 458, 464, 465, 466, 466, 484, 450, 456, 0, 466,
- 475, 0, 468, 463, 458, 0, 458, 476, 0, 464,
- 474, 0, 460, 457, 487, 478, 471, 466, 0, 467,
- 0, 0, 473, 485, 488, 473, 484, 0, 501, 496,
- 497, 498, 495, 508, 491, 487, 0, 496, 480, 0,
- 502, 496, 481, 494, 0, 489, 0, 492, 0, 492,
-
- 502, 503, 506, 495, 502, 0, 0, 515, 537, 538,
- 539, 509, 532, 510, 0, 506, 527, 514, 0, 521,
- 0, 522, 519, 510, 516, 528, 508, 519, 515, 516,
- 543, 544, 545, 534, 540, 0, 522, 537, 525, 527,
- 527, 0, 545, 534, 0, 560, 533, 0, 543, 558,
- 559, 560, 566, 562, 539, 541, 0, 0, 0, 546,
- 542, 547, 538, 0, 579, 580, 581, 565, 589, 0,
- 549, 550, 0, 587, 0, 581, 582, 583, 562, 592,
- 558, 0, 563, 589, 590, 591, 576, 602, 0, 572,
- 0, 0, 0, 0, 596, 580, 606, 569, 0, 0,
-
- 1135, 633
- } ;
-
-static yyconst flex_int16_t yy_def[1003] =
- { 0,
- 1001, 1, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001,
- 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001,
- 1001, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1001, 1001, 1001, 1001, 1001, 1001,
- 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001,
- 1001, 1001, 1001, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
-
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
-
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1001, 1001, 1001, 1001, 1001, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
-
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
-
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1001, 1001, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
-
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
-
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
-
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
-
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
-
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
- 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002, 1002,
-
- 0, 1001
- } ;
-
-static yyconst flex_int16_t yy_nxt[1209] =
- { 0,
- 4, 5, 6, 7, 8, 9, 10, 11, 12, 13,
- 14, 15, 16, 16, 16, 16, 16, 17, 18, 19,
- 20, 21, 22, 23, 24, 25, 26, 27, 28, 29,
- 30, 31, 32, 33, 34, 35, 36, 37, 38, 39,
- 40, 41, 42, 43, 31, 31, 31, 31, 31, 31,
- 31, 31, 31, 31, 31, 31, 31, 31, 31, 44,
- 31, 31, 31, 31, 31, 31, 31, 31, 31, 31,
- 31, 31, 45, 49, 53, 51, 54, 54, 54, 54,
- 54, 209, 210, 55, 97, 50, 52, 56, 58, 72,
- 59, 59, 59, 59, 59, 766, 57, 65, 66, 67,
-
- 73, 98, 297, 75, 60, 767, 68, 76, 99, 74,
- 77, 78, 69, 70, 71, 298, 79, 80, 768, 81,
- 83, 86, 82, 87, 84, 103, 91, 88, 769, 60,
- 89, 115, 85, 90, 92, 100, 93, 94, 111, 95,
- 101, 310, 104, 107, 311, 96, 116, 108, 112, 105,
- 102, 109, 113, 117, 127, 119, 110, 118, 128, 120,
- 167, 114, 121, 122, 148, 123, 168, 129, 124, 149,
- 131, 125, 132, 130, 770, 133, 134, 141, 142, 135,
- 170, 771, 136, 137, 184, 185, 138, 139, 171, 186,
- 235, 140, 236, 143, 54, 54, 54, 54, 54, 155,
-
- 155, 155, 155, 155, 218, 223, 193, 772, 154, 219,
- 144, 224, 194, 156, 145, 58, 195, 59, 59, 59,
- 59, 59, 216, 146, 245, 246, 147, 233, 259, 260,
- 217, 60, 773, 154, 774, 157, 157, 264, 156, 158,
- 158, 158, 158, 158, 175, 647, 234, 176, 177, 197,
- 178, 179, 265, 648, 180, 198, 60, 199, 240, 253,
- 775, 241, 254, 280, 314, 315, 242, 200, 255, 158,
- 158, 158, 158, 158, 243, 776, 244, 271, 271, 281,
- 282, 272, 272, 272, 272, 272, 155, 155, 155, 155,
- 155, 274, 274, 465, 466, 275, 275, 275, 275, 275,
-
- 273, 158, 158, 158, 158, 158, 294, 326, 363, 272,
- 272, 272, 272, 272, 364, 295, 327, 373, 777, 385,
- 778, 328, 386, 401, 329, 273, 402, 779, 374, 375,
- 376, 403, 387, 377, 388, 404, 428, 405, 406, 429,
- 780, 407, 272, 272, 272, 272, 272, 408, 408, 483,
- 781, 409, 409, 409, 409, 409, 275, 275, 275, 275,
- 275, 275, 275, 275, 275, 275, 501, 504, 548, 549,
- 782, 606, 502, 483, 505, 607, 669, 662, 506, 409,
- 409, 409, 409, 409, 409, 409, 409, 409, 409, 557,
- 613, 635, 558, 639, 663, 614, 693, 640, 670, 783,
-
- 669, 741, 636, 784, 742, 785, 786, 637, 787, 788,
- 789, 790, 791, 694, 743, 792, 793, 794, 795, 744,
- 693, 796, 797, 798, 799, 800, 741, 801, 802, 742,
- 803, 804, 805, 806, 807, 808, 809, 810, 743, 811,
- 812, 813, 814, 744, 815, 816, 817, 818, 819, 820,
- 821, 822, 823, 824, 825, 826, 827, 828, 829, 830,
- 831, 832, 833, 834, 835, 836, 837, 838, 839, 840,
- 841, 842, 843, 844, 845, 846, 847, 848, 849, 850,
- 851, 852, 853, 854, 855, 856, 857, 858, 859, 860,
- 861, 862, 863, 864, 865, 866, 867, 868, 869, 870,
-
- 871, 872, 873, 874, 875, 876, 853, 877, 878, 879,
- 880, 881, 882, 883, 884, 885, 886, 887, 888, 889,
- 890, 891, 892, 893, 894, 895, 896, 897, 898, 899,
- 900, 901, 902, 903, 904, 905, 906, 907, 908, 909,
- 910, 911, 888, 912, 913, 914, 915, 916, 917, 918,
- 919, 920, 898, 921, 922, 923, 924, 925, 926, 927,
- 928, 929, 930, 931, 932, 933, 934, 935, 936, 937,
- 938, 939, 940, 941, 942, 943, 944, 945, 946, 947,
- 948, 949, 950, 951, 952, 953, 954, 955, 956, 957,
- 958, 959, 960, 961, 938, 962, 963, 964, 965, 966,
-
- 967, 968, 969, 970, 971, 972, 973, 974, 975, 976,
- 977, 978, 979, 980, 981, 982, 983, 984, 985, 986,
- 987, 988, 989, 990, 991, 992, 993, 994, 995, 996,
- 997, 998, 999, 1000, 64, 765, 764, 763, 762, 761,
- 760, 759, 758, 757, 756, 755, 754, 753, 752, 751,
- 750, 749, 748, 747, 746, 745, 740, 739, 738, 737,
- 736, 735, 734, 733, 732, 731, 730, 729, 728, 727,
- 726, 725, 724, 723, 722, 721, 720, 719, 718, 717,
- 716, 715, 714, 713, 712, 711, 710, 709, 708, 707,
- 706, 705, 704, 703, 702, 701, 700, 699, 698, 697,
-
- 696, 695, 692, 691, 690, 689, 688, 687, 686, 685,
- 684, 683, 682, 681, 680, 679, 678, 677, 676, 675,
- 674, 673, 672, 671, 668, 667, 666, 665, 664, 661,
- 660, 659, 658, 657, 656, 655, 654, 653, 652, 651,
- 650, 649, 646, 645, 644, 643, 642, 641, 638, 634,
- 633, 632, 631, 630, 629, 628, 627, 626, 625, 624,
- 623, 622, 621, 620, 619, 618, 617, 616, 615, 612,
- 611, 610, 609, 608, 605, 604, 603, 602, 601, 600,
- 599, 598, 597, 596, 595, 594, 593, 592, 591, 590,
- 589, 588, 587, 586, 585, 584, 583, 582, 581, 580,
-
- 579, 578, 577, 576, 575, 574, 573, 572, 571, 570,
- 569, 568, 567, 566, 565, 564, 563, 562, 561, 560,
- 559, 556, 555, 554, 553, 552, 551, 550, 547, 546,
- 545, 544, 543, 542, 541, 540, 539, 538, 537, 536,
- 535, 534, 533, 532, 531, 530, 529, 528, 527, 526,
- 525, 524, 523, 522, 521, 520, 519, 518, 517, 516,
- 515, 514, 513, 512, 511, 510, 509, 508, 507, 503,
- 500, 499, 498, 497, 496, 495, 494, 493, 492, 491,
- 490, 489, 488, 487, 486, 485, 484, 482, 481, 480,
- 479, 478, 477, 476, 475, 474, 473, 472, 471, 470,
-
- 469, 468, 467, 464, 463, 462, 461, 460, 459, 458,
- 457, 456, 455, 454, 453, 452, 451, 450, 449, 448,
- 447, 446, 445, 444, 443, 442, 441, 440, 439, 438,
- 437, 436, 435, 434, 433, 432, 431, 430, 427, 426,
- 425, 424, 423, 422, 421, 420, 419, 418, 417, 416,
- 415, 414, 413, 412, 411, 410, 400, 399, 398, 397,
- 396, 395, 394, 393, 392, 391, 390, 389, 384, 383,
- 382, 381, 380, 379, 378, 372, 371, 370, 369, 368,
- 367, 366, 365, 362, 361, 360, 359, 358, 357, 356,
- 355, 354, 353, 352, 351, 350, 349, 348, 347, 346,
-
- 345, 344, 343, 342, 341, 340, 339, 338, 337, 336,
- 335, 334, 333, 332, 331, 330, 325, 324, 323, 322,
- 321, 320, 319, 318, 317, 316, 313, 312, 309, 308,
- 307, 306, 305, 304, 303, 302, 301, 300, 299, 296,
- 293, 292, 291, 290, 289, 288, 287, 286, 285, 284,
- 283, 279, 278, 277, 276, 270, 269, 268, 267, 266,
- 263, 262, 261, 258, 257, 256, 252, 251, 250, 249,
- 248, 247, 239, 238, 237, 232, 231, 230, 229, 228,
- 227, 226, 225, 222, 221, 220, 215, 214, 213, 212,
- 211, 208, 207, 206, 205, 204, 203, 202, 201, 196,
-
- 192, 191, 190, 189, 188, 187, 183, 182, 181, 174,
- 173, 172, 169, 166, 165, 164, 163, 162, 161, 160,
- 159, 153, 152, 151, 150, 126, 106, 63, 62, 61,
- 48, 47, 46, 1001, 3, 1001, 1001, 1001, 1001, 1001,
- 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001,
- 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001,
- 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001,
- 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001,
- 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001,
- 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001,
-
- 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001
- } ;
-
-static yyconst flex_int16_t yy_chk[1209] =
- { 0,
- 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
- 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
- 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
- 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
- 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
- 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
- 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
- 1, 1, 1, 12, 14, 13, 14, 14, 14, 14,
- 14, 101, 101, 15, 28, 12, 13, 15, 16, 23,
- 16, 16, 16, 16, 16, 694, 15, 22, 22, 22,
-
- 23, 28, 177, 24, 16, 695, 22, 24, 28, 23,
- 24, 24, 22, 22, 22, 177, 24, 24, 696, 24,
- 25, 26, 24, 26, 25, 30, 27, 26, 697, 16,
- 26, 35, 25, 26, 27, 29, 27, 27, 34, 27,
- 29, 190, 30, 33, 190, 27, 35, 33, 34, 30,
- 29, 33, 34, 36, 39, 37, 33, 36, 39, 37,
- 73, 34, 37, 37, 42, 37, 73, 39, 37, 42,
- 40, 37, 40, 39, 698, 40, 40, 41, 41, 40,
- 75, 699, 40, 40, 83, 83, 40, 40, 75, 83,
- 124, 40, 124, 41, 54, 54, 54, 54, 54, 58,
-
- 58, 58, 58, 58, 109, 113, 90, 700, 54, 109,
- 41, 113, 90, 58, 41, 59, 90, 59, 59, 59,
- 59, 59, 108, 41, 129, 129, 41, 123, 140, 140,
- 108, 59, 701, 54, 702, 60, 60, 146, 58, 60,
- 60, 60, 60, 60, 79, 557, 123, 79, 79, 92,
- 79, 79, 146, 557, 79, 92, 59, 92, 128, 136,
- 703, 128, 136, 163, 193, 193, 128, 92, 136, 157,
- 157, 157, 157, 157, 128, 704, 128, 154, 154, 163,
- 163, 154, 154, 154, 154, 154, 155, 155, 155, 155,
- 155, 156, 156, 344, 344, 156, 156, 156, 156, 156,
-
- 155, 158, 158, 158, 158, 158, 175, 206, 240, 271,
- 271, 271, 271, 271, 240, 175, 206, 249, 706, 257,
- 707, 206, 257, 270, 206, 155, 270, 708, 249, 249,
- 249, 270, 257, 249, 257, 270, 298, 270, 270, 298,
- 709, 270, 272, 272, 272, 272, 272, 273, 273, 361,
- 710, 273, 273, 273, 273, 273, 274, 274, 274, 274,
- 274, 275, 275, 275, 275, 275, 382, 385, 438, 438,
- 711, 509, 382, 361, 385, 509, 584, 576, 385, 408,
- 408, 408, 408, 408, 409, 409, 409, 409, 409, 450,
- 516, 544, 450, 548, 576, 516, 613, 548, 584, 712,
-
- 584, 668, 544, 713, 668, 716, 717, 544, 718, 720,
- 721, 722, 723, 613, 668, 724, 726, 728, 729, 668,
- 613, 733, 735, 736, 737, 740, 668, 741, 742, 668,
- 743, 744, 745, 746, 747, 748, 749, 750, 668, 751,
- 752, 754, 755, 668, 757, 758, 759, 760, 762, 763,
- 764, 765, 766, 767, 768, 770, 771, 772, 773, 775,
- 776, 777, 778, 779, 780, 781, 782, 783, 784, 786,
- 787, 788, 789, 790, 793, 794, 796, 797, 799, 800,
- 801, 802, 803, 804, 805, 808, 809, 810, 811, 813,
- 817, 818, 820, 821, 822, 823, 825, 826, 827, 829,
-
- 832, 833, 834, 835, 836, 837, 803, 838, 839, 841,
- 842, 843, 844, 845, 846, 847, 848, 850, 851, 853,
- 854, 855, 857, 858, 860, 861, 863, 864, 865, 866,
- 867, 868, 870, 873, 874, 875, 876, 877, 879, 880,
- 881, 882, 851, 883, 884, 885, 886, 888, 889, 891,
- 892, 893, 865, 894, 896, 898, 900, 901, 902, 903,
- 904, 905, 908, 909, 910, 911, 912, 913, 914, 916,
- 917, 918, 920, 922, 923, 924, 925, 926, 927, 928,
- 929, 930, 931, 932, 933, 934, 935, 937, 938, 939,
- 940, 941, 943, 944, 917, 946, 947, 949, 950, 951,
-
- 952, 953, 954, 955, 956, 960, 961, 962, 963, 965,
- 966, 967, 968, 969, 971, 972, 974, 976, 977, 978,
- 979, 980, 981, 983, 984, 985, 986, 987, 988, 990,
- 995, 996, 997, 998, 1002, 693, 692, 691, 690, 689,
- 687, 686, 685, 683, 681, 680, 679, 678, 677, 676,
- 675, 673, 672, 671, 670, 669, 667, 666, 665, 664,
- 663, 662, 661, 660, 659, 658, 657, 656, 655, 654,
- 653, 652, 651, 650, 649, 648, 647, 646, 645, 644,
- 642, 641, 640, 638, 637, 636, 635, 634, 633, 632,
- 631, 630, 627, 626, 625, 624, 623, 622, 619, 618,
-
- 615, 614, 612, 611, 610, 609, 608, 607, 606, 605,
- 604, 602, 600, 597, 596, 595, 594, 593, 592, 591,
- 589, 587, 586, 585, 583, 581, 579, 578, 577, 575,
- 573, 571, 570, 569, 566, 565, 564, 563, 562, 561,
- 560, 558, 556, 555, 554, 553, 551, 549, 545, 543,
- 542, 541, 540, 539, 538, 537, 536, 535, 534, 532,
- 531, 530, 528, 527, 524, 523, 521, 519, 517, 515,
- 513, 512, 511, 510, 508, 507, 506, 505, 504, 503,
- 501, 500, 499, 498, 497, 496, 495, 494, 493, 490,
- 489, 488, 487, 486, 485, 484, 483, 482, 481, 479,
-
- 477, 476, 475, 473, 471, 470, 469, 468, 466, 465,
- 463, 462, 461, 459, 458, 456, 455, 454, 453, 452,
- 451, 449, 447, 446, 444, 443, 442, 439, 437, 436,
- 435, 434, 433, 432, 431, 430, 429, 428, 427, 426,
- 425, 424, 423, 422, 421, 420, 419, 417, 415, 414,
- 413, 412, 411, 406, 405, 403, 400, 399, 398, 397,
- 396, 395, 394, 391, 390, 389, 388, 387, 386, 383,
- 381, 380, 377, 376, 375, 374, 373, 372, 371, 370,
- 369, 368, 367, 366, 365, 364, 363, 360, 359, 358,
- 357, 356, 355, 354, 353, 352, 351, 350, 349, 348,
-
- 347, 346, 345, 343, 340, 339, 338, 337, 336, 335,
- 334, 333, 330, 329, 328, 327, 326, 325, 324, 322,
- 321, 319, 318, 316, 314, 313, 312, 311, 310, 309,
- 307, 306, 305, 304, 303, 302, 300, 299, 297, 296,
- 295, 294, 293, 292, 290, 288, 287, 286, 285, 284,
- 282, 281, 280, 279, 278, 276, 269, 268, 267, 266,
- 265, 264, 263, 262, 261, 260, 259, 258, 256, 255,
- 254, 253, 252, 251, 250, 248, 247, 246, 245, 244,
- 243, 242, 241, 239, 238, 237, 236, 235, 234, 233,
- 232, 231, 230, 229, 228, 227, 226, 225, 224, 223,
-
- 222, 221, 220, 219, 218, 217, 216, 215, 214, 213,
- 212, 211, 210, 209, 208, 207, 205, 204, 202, 201,
- 200, 199, 198, 197, 196, 195, 192, 191, 189, 188,
- 187, 186, 185, 184, 183, 182, 180, 179, 178, 176,
- 174, 173, 172, 171, 170, 169, 168, 167, 166, 165,
- 164, 162, 161, 160, 159, 151, 150, 149, 148, 147,
- 145, 144, 143, 139, 138, 137, 135, 134, 133, 132,
- 131, 130, 127, 126, 125, 121, 120, 119, 118, 117,
- 116, 115, 114, 112, 111, 110, 107, 106, 105, 103,
- 102, 100, 99, 98, 97, 96, 95, 94, 93, 91,
-
- 89, 88, 87, 86, 85, 84, 82, 81, 80, 78,
- 77, 76, 74, 72, 71, 70, 69, 68, 67, 66,
- 65, 53, 45, 44, 43, 38, 32, 21, 20, 19,
- 11, 9, 7, 3, 1001, 1001, 1001, 1001, 1001, 1001,
- 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001,
- 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001,
- 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001,
- 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001,
- 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001,
- 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001,
-
- 1001, 1001, 1001, 1001, 1001, 1001, 1001, 1001
- } ;
-
-static yy_state_type yy_last_accepting_state;
-static char *yy_last_accepting_cpos;
-
-extern int gmsh_yy_flex_debug;
-int gmsh_yy_flex_debug = 0;
-
-/* The intent behind this definition is that it'll catch
- * any uses of REJECT which flex missed.
- */
-#define REJECT reject_used_but_not_detected
-#define yymore() yymore_used_but_not_detected
-#define YY_MORE_ADJ 0
-#define YY_RESTORE_YY_MORE_OFFSET
-char *gmsh_yytext;
-#line 1 "Gmsh.l"
-#line 2 "Gmsh.l"
-// Gmsh - Copyright (C) 1997-2015 C. Geuzaine, J.-F. Remacle
-//
-// See the LICENSE.txt file for license information. Please report all
-// bugs and problems to the public mailing list <gmsh@geuz.org>.
-
-#include <stdio.h>
-#include <stdlib.h>
-#include <string.h>
-#include <math.h>
-#include "GmshConfig.h"
-#include "GmshMessage.h"
-#include "Geo.h"
-#include "Gmsh.tab.hpp"
-#include "GmshIO.h"
-
-void parsestring(char endchar);
-char *strsave(char *ptr);
-void skipcomments(void);
-void skipline(void);
-
-#if defined(HAVE_COMPRESSED_IO) && defined(HAVE_ZLIB)
-#define YY_INPUT(buf,result,max_size) \
- { \
- int c = '*', n; \
- for ( n = 0; n < (int) max_size && \
- (c = gzgetc( gmsh_yyin )) != EOF && c != '\n'; ++n ) \
- buf[n] = (char) c; \
- if ( c == '\n' ){ \
- buf[n++] = (char) c; \
- gmsh_yylineno++; \
- } \
- if ( c == EOF ) { \
- int ernum; \
- const char *msg=gzerror(gmsh_yyin,&ernum); \
- if (ernum) \
- Msg::Fatal("Input in flex scanner failed"); \
- } \
- result = n; \
- }
-#else
-#define YY_INPUT(buf,result,max_size) \
- { \
- int c = '*', n; \
- for ( n = 0; n < (int) max_size && \
- (c = fgetc( gmsh_yyin )) != EOF && c != '\n'; ++n ) \
- buf[n] = (char) c; \
- if ( c == '\n' ){ \
- buf[n++] = (char) c; \
- gmsh_yylineno++; \
- } \
- if ( c == EOF && ferror( gmsh_yyin ) ) \
- Msg::Fatal("Input in flex scanner failed"); \
- result = n; \
- }
-#endif
-
-#if defined(WIN32)
-#define isatty(arg) -1
-#define YY_NO_UNISTD_H
-#endif
-
-#line 1133 "Gmsh.yy.cpp"
-
-#define INITIAL 0
-
-#ifndef YY_NO_UNISTD_H
-/* Special case for "unistd.h", since it is non-ANSI. We include it way
- * down here because we want the user's section 1 to have been scanned first.
- * The user has a chance to override it with an option.
- */
-#include <unistd.h>
-#endif
-
-#ifndef YY_EXTRA_TYPE
-#define YY_EXTRA_TYPE void *
-#endif
-
-static int yy_init_globals (void );
-
-/* Accessor methods to globals.
- These are made visible to non-reentrant scanners for convenience. */
-
-int gmsh_yylex_destroy (void );
-
-int gmsh_yyget_debug (void );
-
-void gmsh_yyset_debug (int debug_flag );
-
-YY_EXTRA_TYPE gmsh_yyget_extra (void );
-
-void gmsh_yyset_extra (YY_EXTRA_TYPE user_defined );
-
-FILE *gmsh_yyget_in (void );
-
-void gmsh_yyset_in (FILE * in_str );
-
-FILE *gmsh_yyget_out (void );
-
-void gmsh_yyset_out (FILE * out_str );
-
-yy_size_t gmsh_yyget_leng (void );
-
-char *gmsh_yyget_text (void );
-
-int gmsh_yyget_lineno (void );
-
-void gmsh_yyset_lineno (int line_number );
-
-/* Macros after this point can all be overridden by user definitions in
- * section 1.
- */
-
-#ifndef YY_SKIP_YYWRAP
-#ifdef __cplusplus
-extern "C" int gmsh_yywrap (void );
-#else
-extern int gmsh_yywrap (void );
-#endif
-#endif
-
- static void yyunput (int c,char *buf_ptr );
-
-#ifndef yytext_ptr
-static void yy_flex_strncpy (char *,yyconst char *,int );
-#endif
-
-#ifdef YY_NEED_STRLEN
-static int yy_flex_strlen (yyconst char * );
-#endif
-
-#ifndef YY_NO_INPUT
-
-#ifdef __cplusplus
-static int yyinput (void );
-#else
-static int input (void );
-#endif
-
-#endif
-
-/* Amount of stuff to slurp up with each read. */
-#ifndef YY_READ_BUF_SIZE
-#define YY_READ_BUF_SIZE 8192
-#endif
-
-/* Copy whatever the last rule matched to the standard output. */
-#ifndef ECHO
-/* This used to be an fputs(), but since the string might contain NUL's,
- * we now use fwrite().
- */
-#define ECHO do { if (fwrite( gmsh_yytext, gmsh_yyleng, 1, gmsh_yyout )) {} } while (0)
-#endif
-
-/* Gets input and stuffs it into "buf". number of characters read, or YY_NULL,
- * is returned in "result".
- */
-#ifndef YY_INPUT
-#define YY_INPUT(buf,result,max_size) \
- if ( YY_CURRENT_BUFFER_LVALUE->yy_is_interactive ) \
- { \
- int c = '*'; \
- size_t n; \
- for ( n = 0; n < max_size && \
- (c = getc( gmsh_yyin )) != EOF && c != '\n'; ++n ) \
- buf[n] = (char) c; \
- if ( c == '\n' ) \
- buf[n++] = (char) c; \
- if ( c == EOF && ferror( gmsh_yyin ) ) \
- YY_FATAL_ERROR( "input in flex scanner failed" ); \
- result = n; \
- } \
- else \
- { \
- errno=0; \
- while ( (result = fread(buf, 1, max_size, gmsh_yyin))==0 && ferror(gmsh_yyin)) \
- { \
- if( errno != EINTR) \
- { \
- YY_FATAL_ERROR( "input in flex scanner failed" ); \
- break; \
- } \
- errno=0; \
- clearerr(gmsh_yyin); \
- } \
- }\
-\
-
-#endif
-
-/* No semi-colon after return; correct usage is to write "yyterminate();" -
- * we don't want an extra ';' after the "return" because that will cause
- * some compilers to complain about unreachable statements.
- */
-#ifndef yyterminate
-#define yyterminate() return YY_NULL
-#endif
-
-/* Number of entries by which start-condition stack grows. */
-#ifndef YY_START_STACK_INCR
-#define YY_START_STACK_INCR 25
-#endif
-
-/* Report a fatal error. */
-#ifndef YY_FATAL_ERROR
-#define YY_FATAL_ERROR(msg) yy_fatal_error( msg )
-#endif
-
-/* end tables serialization structures and prototypes */
-
-/* Default declaration of generated scanner - a define so the user can
- * easily add parameters.
- */
-#ifndef YY_DECL
-#define YY_DECL_IS_OURS 1
-
-extern int gmsh_yylex (void);
-
-#define YY_DECL int gmsh_yylex (void)
-#endif /* !YY_DECL */
-
-/* Code executed at the beginning of each rule, after gmsh_yytext and gmsh_yyleng
- * have been set up.
- */
-#ifndef YY_USER_ACTION
-#define YY_USER_ACTION
-#endif
-
-/* Code executed at the end of each rule. */
-#ifndef YY_BREAK
-#define YY_BREAK break;
-#endif
-
-#define YY_RULE_SETUP \
- YY_USER_ACTION
-
-/** The main scanner function which does all the work.
- */
-YY_DECL
-{
- register yy_state_type yy_current_state;
- register char *yy_cp, *yy_bp;
- register int yy_act;
-
-#line 72 "Gmsh.l"
-
-
-#line 1318 "Gmsh.yy.cpp"
-
- if ( !(yy_init) )
- {
- (yy_init) = 1;
-
-#ifdef YY_USER_INIT
- YY_USER_INIT;
-#endif
-
- if ( ! (yy_start) )
- (yy_start) = 1; /* first start state */
-
- if ( ! gmsh_yyin )
- gmsh_yyin = stdin;
-
- if ( ! gmsh_yyout )
- gmsh_yyout = stdout;
-
- if ( ! YY_CURRENT_BUFFER ) {
- gmsh_yyensure_buffer_stack ();
- YY_CURRENT_BUFFER_LVALUE =
- gmsh_yy_create_buffer(gmsh_yyin,YY_BUF_SIZE );
- }
-
- gmsh_yy_load_buffer_state( );
- }
-
- while ( 1 ) /* loops until end-of-file is reached */
- {
- yy_cp = (yy_c_buf_p);
-
- /* Support of gmsh_yytext. */
- *yy_cp = (yy_hold_char);
-
- /* yy_bp points to the position in yy_ch_buf of the start of
- * the current run.
- */
- yy_bp = yy_cp;
-
- yy_current_state = (yy_start);
-yy_match:
- do
- {
- register YY_CHAR yy_c = yy_ec[YY_SC_TO_UI(*yy_cp)];
- if ( yy_accept[yy_current_state] )
- {
- (yy_last_accepting_state) = yy_current_state;
- (yy_last_accepting_cpos) = yy_cp;
- }
- while ( yy_chk[yy_base[yy_current_state] + yy_c] != yy_current_state )
- {
- yy_current_state = (int) yy_def[yy_current_state];
- if ( yy_current_state >= 1002 )
- yy_c = yy_meta[(unsigned int) yy_c];
- }
- yy_current_state = yy_nxt[yy_base[yy_current_state] + (unsigned int) yy_c];
- ++yy_cp;
- }
- while ( yy_base[yy_current_state] != 1135 );
-
-yy_find_action:
- yy_act = yy_accept[yy_current_state];
- if ( yy_act == 0 )
- { /* have to back up */
- yy_cp = (yy_last_accepting_cpos);
- yy_current_state = (yy_last_accepting_state);
- yy_act = yy_accept[yy_current_state];
- }
-
- YY_DO_BEFORE_ACTION;
-
-do_action: /* This label is used only to access EOF actions. */
-
- switch ( yy_act )
- { /* beginning of action switch */
- case 0: /* must back up */
- /* undo the effects of YY_DO_BEFORE_ACTION */
- *yy_cp = (yy_hold_char);
- yy_cp = (yy_last_accepting_cpos);
- yy_current_state = (yy_last_accepting_state);
- goto yy_find_action;
-
-case 1:
-/* rule 1 can match eol */
-YY_RULE_SETUP
-#line 74 "Gmsh.l"
-/* none */;
- YY_BREAK
-case 2:
-YY_RULE_SETUP
-#line 75 "Gmsh.l"
-return tEND;
- YY_BREAK
-case 3:
-YY_RULE_SETUP
-#line 76 "Gmsh.l"
-skipcomments();
- YY_BREAK
-case 4:
-YY_RULE_SETUP
-#line 77 "Gmsh.l"
-skipline();
- YY_BREAK
-case 5:
-YY_RULE_SETUP
-#line 78 "Gmsh.l"
-{ parsestring('\"'); return tBIGSTR; }
- YY_BREAK
-case 6:
-YY_RULE_SETUP
-#line 79 "Gmsh.l"
-{ parsestring('\''); return tBIGSTR; }
- YY_BREAK
-case 7:
-YY_RULE_SETUP
-#line 80 "Gmsh.l"
-{ gmsh_yylval.d = NEWREG(); return tDOUBLE; }
- YY_BREAK
-case 8:
-YY_RULE_SETUP
-#line 81 "Gmsh.l"
-{ gmsh_yylval.d = NEWPOINT(); return tDOUBLE; }
- YY_BREAK
-case 9:
-YY_RULE_SETUP
-#line 82 "Gmsh.l"
-{ gmsh_yylval.d = NEWLINE(); return tDOUBLE; }
- YY_BREAK
-case 10:
-YY_RULE_SETUP
-#line 83 "Gmsh.l"
-{ gmsh_yylval.d = NEWLINE(); return tDOUBLE; }
- YY_BREAK
-case 11:
-YY_RULE_SETUP
-#line 84 "Gmsh.l"
-{ gmsh_yylval.d = NEWLINELOOP(); return tDOUBLE; }
- YY_BREAK
-case 12:
-YY_RULE_SETUP
-#line 85 "Gmsh.l"
-{ gmsh_yylval.d = NEWSURFACE(); return tDOUBLE; }
- YY_BREAK
-case 13:
-YY_RULE_SETUP
-#line 86 "Gmsh.l"
-{ gmsh_yylval.d = NEWSURFACELOOP(); return tDOUBLE; }
- YY_BREAK
-case 14:
-YY_RULE_SETUP
-#line 87 "Gmsh.l"
-{ gmsh_yylval.d = NEWVOLUME(); return tDOUBLE; }
- YY_BREAK
-case 15:
-YY_RULE_SETUP
-#line 88 "Gmsh.l"
-{ gmsh_yylval.d = NEWFIELD(); return tDOUBLE; }
- YY_BREAK
-case 16:
-YY_RULE_SETUP
-#line 89 "Gmsh.l"
-return tAFFECT;
- YY_BREAK
-case 17:
-YY_RULE_SETUP
-#line 90 "Gmsh.l"
-return tAFFECTPLUS;
- YY_BREAK
-case 18:
-YY_RULE_SETUP
-#line 91 "Gmsh.l"
-return tAFFECTMINUS;
- YY_BREAK
-case 19:
-YY_RULE_SETUP
-#line 92 "Gmsh.l"
-return tAFFECTTIMES;
- YY_BREAK
-case 20:
-YY_RULE_SETUP
-#line 93 "Gmsh.l"
-return tAFFECTDIVIDE;
- YY_BREAK
-case 21:
-YY_RULE_SETUP
-#line 94 "Gmsh.l"
-return tDOTS;
- YY_BREAK
-case 22:
-YY_RULE_SETUP
-#line 95 "Gmsh.l"
-return tDOTS;
- YY_BREAK
-case 23:
-YY_RULE_SETUP
-#line 96 "Gmsh.l"
-return tOR;
- YY_BREAK
-case 24:
-YY_RULE_SETUP
-#line 97 "Gmsh.l"
-return tAND;
- YY_BREAK
-case 25:
-YY_RULE_SETUP
-#line 98 "Gmsh.l"
-return tPLUSPLUS;
- YY_BREAK
-case 26:
-YY_RULE_SETUP
-#line 99 "Gmsh.l"
-return tMINUSMINUS;
- YY_BREAK
-case 27:
-YY_RULE_SETUP
-#line 100 "Gmsh.l"
-return tEQUAL;
- YY_BREAK
-case 28:
-YY_RULE_SETUP
-#line 101 "Gmsh.l"
-return tNOTEQUAL;
- YY_BREAK
-case 29:
-YY_RULE_SETUP
-#line 102 "Gmsh.l"
-return tLESSOREQUAL;
- YY_BREAK
-case 30:
-YY_RULE_SETUP
-#line 103 "Gmsh.l"
-return tGREATEROREQUAL;
- YY_BREAK
-case 31:
-YY_RULE_SETUP
-#line 105 "Gmsh.l"
-return tAbort;
- YY_BREAK
-case 32:
-YY_RULE_SETUP
-#line 106 "Gmsh.l"
-return tAcos;
- YY_BREAK
-case 33:
-YY_RULE_SETUP
-#line 107 "Gmsh.l"
-return tAdaptMesh;
- YY_BREAK
-case 34:
-YY_RULE_SETUP
-#line 108 "Gmsh.l"
-return tAcos;
- YY_BREAK
-case 35:
-YY_RULE_SETUP
-#line 109 "Gmsh.l"
-return tAlias;
- YY_BREAK
-case 36:
-YY_RULE_SETUP
-#line 110 "Gmsh.l"
-return tAliasWithOptions;
- YY_BREAK
-case 37:
-YY_RULE_SETUP
-#line 111 "Gmsh.l"
-return tAsin;
- YY_BREAK
-case 38:
-YY_RULE_SETUP
-#line 112 "Gmsh.l"
-return tAsin;
- YY_BREAK
-case 39:
-YY_RULE_SETUP
-#line 113 "Gmsh.l"
-return tAtan;
- YY_BREAK
-case 40:
-YY_RULE_SETUP
-#line 114 "Gmsh.l"
-return tAtan;
- YY_BREAK
-case 41:
-YY_RULE_SETUP
-#line 115 "Gmsh.l"
-return tAtan2;
- YY_BREAK
-case 42:
-YY_RULE_SETUP
-#line 116 "Gmsh.l"
-return tAtan2;
- YY_BREAK
-case 43:
-YY_RULE_SETUP
-#line 118 "Gmsh.l"
-return tBetti;
- YY_BREAK
-case 44:
-YY_RULE_SETUP
-#line 119 "Gmsh.l"
-return tBezier;
- YY_BREAK
-case 45:
-YY_RULE_SETUP
-#line 120 "Gmsh.l"
-return tBSpline;
- YY_BREAK
-case 46:
-YY_RULE_SETUP
-#line 121 "Gmsh.l"
-return tBoundingBox;
- YY_BREAK
-case 47:
-YY_RULE_SETUP
-#line 123 "Gmsh.l"
-return tCeil;
- YY_BREAK
-case 48:
-YY_RULE_SETUP
-#line 124 "Gmsh.l"
-return tCircle;
- YY_BREAK
-case 49:
-YY_RULE_SETUP
-#line 125 "Gmsh.l"
-return tCopyOptions;
- YY_BREAK
-case 50:
-YY_RULE_SETUP
-#line 126 "Gmsh.l"
-return tCoherence;
- YY_BREAK
-case 51:
-YY_RULE_SETUP
-#line 127 "Gmsh.l"
-return tCohomology;
- YY_BREAK
-case 52:
-YY_RULE_SETUP
-#line 128 "Gmsh.l"
-return tCombine;
- YY_BREAK
-case 53:
-YY_RULE_SETUP
-#line 129 "Gmsh.l"
-return tCosh;
- YY_BREAK
-case 54:
-YY_RULE_SETUP
-#line 130 "Gmsh.l"
-return tCos;
- YY_BREAK
-case 55:
-YY_RULE_SETUP
-#line 131 "Gmsh.l"
-return tCharacteristic;
- YY_BREAK
-case 56:
-YY_RULE_SETUP
-#line 132 "Gmsh.l"
-return tComplex;
- YY_BREAK
-case 57:
-YY_RULE_SETUP
-#line 133 "Gmsh.l"
-return tColor;
- YY_BREAK
-case 58:
-YY_RULE_SETUP
-#line 134 "Gmsh.l"
-return tColorTable;
- YY_BREAK
-case 59:
-YY_RULE_SETUP
-#line 135 "Gmsh.l"
-return tCompound;
- YY_BREAK
-case 60:
-YY_RULE_SETUP
-#line 136 "Gmsh.l"
-return tCoordinates;
- YY_BREAK
-case 61:
-YY_RULE_SETUP
-#line 137 "Gmsh.l"
-return tSpline;
- YY_BREAK
-case 62:
-YY_RULE_SETUP
-#line 138 "Gmsh.l"
-return tCall;
- YY_BREAK
-case 63:
-YY_RULE_SETUP
-#line 139 "Gmsh.l"
-return tCreateTopology;
- YY_BREAK
-case 64:
-YY_RULE_SETUP
-#line 140 "Gmsh.l"
-return tCreateTopologyNoHoles;
- YY_BREAK
-case 65:
-YY_RULE_SETUP
-#line 141 "Gmsh.l"
-return tCpu;
- YY_BREAK
-case 66:
-YY_RULE_SETUP
-#line 142 "Gmsh.l"
-return tCurrentDirectory;
- YY_BREAK
-case 67:
-YY_RULE_SETUP
-#line 143 "Gmsh.l"
-return tCurrentDirectory;
- YY_BREAK
-case 68:
-YY_RULE_SETUP
-#line 145 "Gmsh.l"
-return tDefineConstant;
- YY_BREAK
-case 69:
-YY_RULE_SETUP
-#line 146 "Gmsh.l"
-return tDefineNumber;
- YY_BREAK
-case 70:
-YY_RULE_SETUP
-#line 147 "Gmsh.l"
-return tDefineString;
- YY_BREAK
-case 71:
-YY_RULE_SETUP
-#line 148 "Gmsh.l"
-return tDelete;
- YY_BREAK
-case 72:
-YY_RULE_SETUP
-#line 149 "Gmsh.l"
-return tDegenerated;
- YY_BREAK
-case 73:
-YY_RULE_SETUP
-#line 150 "Gmsh.l"
-return tDilate;
- YY_BREAK
-case 74:
-YY_RULE_SETUP
-#line 151 "Gmsh.l"
-return tDraw;
- YY_BREAK
-case 75:
-YY_RULE_SETUP
-#line 152 "Gmsh.l"
-return tSetChanged;
- YY_BREAK
-case 76:
-YY_RULE_SETUP
-#line 154 "Gmsh.l"
-return tExp;
- YY_BREAK
-case 77:
-YY_RULE_SETUP
-#line 155 "Gmsh.l"
-return tEllipse;
- YY_BREAK
-case 78:
-YY_RULE_SETUP
-#line 156 "Gmsh.l"
-return tEllipse;
- YY_BREAK
-case 79:
-YY_RULE_SETUP
-#line 157 "Gmsh.l"
-return tExtrude;
- YY_BREAK
-case 80:
-YY_RULE_SETUP
-#line 158 "Gmsh.l"
-return tElliptic;
- YY_BREAK
-case 81:
-YY_RULE_SETUP
-#line 159 "Gmsh.l"
-return tEndFor;
- YY_BREAK
-case 82:
-YY_RULE_SETUP
-#line 160 "Gmsh.l"
-return tEndIf;
- YY_BREAK
-case 83:
-YY_RULE_SETUP
-#line 161 "Gmsh.l"
-return tEuclidian;
- YY_BREAK
-case 84:
-YY_RULE_SETUP
-#line 162 "Gmsh.l"
-return tError;
- YY_BREAK
-case 85:
-YY_RULE_SETUP
-#line 163 "Gmsh.l"
-return tExit;
- YY_BREAK
-case 86:
-YY_RULE_SETUP
-#line 164 "Gmsh.l"
-return tExists;
- YY_BREAK
-case 87:
-YY_RULE_SETUP
-#line 166 "Gmsh.l"
-return tFabs;
- YY_BREAK
-case 88:
-YY_RULE_SETUP
-#line 167 "Gmsh.l"
-return tField;
- YY_BREAK
-case 89:
-YY_RULE_SETUP
-#line 168 "Gmsh.l"
-return tFind;
- YY_BREAK
-case 90:
-YY_RULE_SETUP
-#line 169 "Gmsh.l"
-return tFileExists;
- YY_BREAK
-case 91:
-YY_RULE_SETUP
-#line 170 "Gmsh.l"
-return tFixRelativePath;
- YY_BREAK
-case 92:
-YY_RULE_SETUP
-#line 171 "Gmsh.l"
-return tFloor;
- YY_BREAK
-case 93:
-YY_RULE_SETUP
-#line 172 "Gmsh.l"
-return tFmod;
- YY_BREAK
-case 94:
-YY_RULE_SETUP
-#line 173 "Gmsh.l"
-return tFor;
- YY_BREAK
-case 95:
-YY_RULE_SETUP
-#line 174 "Gmsh.l"
-return tMacro;
- YY_BREAK
-case 96:
-YY_RULE_SETUP
-#line 176 "Gmsh.l"
-return tGetEnv;
- YY_BREAK
-case 97:
-YY_RULE_SETUP
-#line 177 "Gmsh.l"
-return tGetString;
- YY_BREAK
-case 98:
-YY_RULE_SETUP
-#line 178 "Gmsh.l"
-return tGetNumber;
- YY_BREAK
-case 99:
-YY_RULE_SETUP
-#line 179 "Gmsh.l"
-return tGetValue;
- YY_BREAK
-case 100:
-YY_RULE_SETUP
-#line 180 "Gmsh.l"
-return tGMSH_MAJOR_VERSION;
- YY_BREAK
-case 101:
-YY_RULE_SETUP
-#line 181 "Gmsh.l"
-return tGMSH_MINOR_VERSION;
- YY_BREAK
-case 102:
-YY_RULE_SETUP
-#line 182 "Gmsh.l"
-return tGMSH_PATCH_VERSION;
- YY_BREAK
-case 103:
-YY_RULE_SETUP
-#line 183 "Gmsh.l"
-return tGmshExecutableName;
- YY_BREAK
-case 104:
-YY_RULE_SETUP
-#line 185 "Gmsh.l"
-return tHide;
- YY_BREAK
-case 105:
-YY_RULE_SETUP
-#line 186 "Gmsh.l"
-return tHole;
- YY_BREAK
-case 106:
-YY_RULE_SETUP
-#line 187 "Gmsh.l"
-return tHomology;
- YY_BREAK
-case 107:
-YY_RULE_SETUP
-#line 188 "Gmsh.l"
-return tHypot;
- YY_BREAK
-case 108:
-YY_RULE_SETUP
-#line 190 "Gmsh.l"
-return tIn;
- YY_BREAK
-case 109:
-YY_RULE_SETUP
-#line 191 "Gmsh.l"
-return tIf;
- YY_BREAK
-case 110:
-YY_RULE_SETUP
-#line 192 "Gmsh.l"
-return tIntersect;
- YY_BREAK
-case 111:
-YY_RULE_SETUP
-#line 193 "Gmsh.l"
-return tInterpolationScheme;
- YY_BREAK
-case 112:
-YY_RULE_SETUP
-#line 195 "Gmsh.l"
-return tNurbsKnots;
- YY_BREAK
-case 113:
-YY_RULE_SETUP
-#line 197 "Gmsh.l"
-return tLength;
- YY_BREAK
-case 114:
-YY_RULE_SETUP
-#line 198 "Gmsh.l"
-return tLine;
- YY_BREAK
-case 115:
-YY_RULE_SETUP
-#line 199 "Gmsh.l"
-return tList;
- YY_BREAK
-case 116:
-YY_RULE_SETUP
-#line 200 "Gmsh.l"
-return tLog;
- YY_BREAK
-case 117:
-YY_RULE_SETUP
-#line 201 "Gmsh.l"
-return tLog10;
- YY_BREAK
-case 118:
-YY_RULE_SETUP
-#line 202 "Gmsh.l"
-return tLayers;
- YY_BREAK
-case 119:
-YY_RULE_SETUP
-#line 203 "Gmsh.l"
-return tLevelset;
- YY_BREAK
-case 120:
-YY_RULE_SETUP
-#line 205 "Gmsh.l"
-return tMacro;
- YY_BREAK
-case 121:
-YY_RULE_SETUP
-#line 206 "Gmsh.l"
-return tMeshAlgorithm;
- YY_BREAK
-case 122:
-YY_RULE_SETUP
-#line 207 "Gmsh.l"
-return tModulo;
- YY_BREAK
-case 123:
-YY_RULE_SETUP
-#line 208 "Gmsh.l"
-return tMPI_Rank;
- YY_BREAK
-case 124:
-YY_RULE_SETUP
-#line 209 "Gmsh.l"
-return tMPI_Size;
- YY_BREAK
-case 125:
-YY_RULE_SETUP
-#line 210 "Gmsh.l"
-return tMemory;
- YY_BREAK
-case 126:
-YY_RULE_SETUP
-#line 212 "Gmsh.l"
-return tNewModel;
- YY_BREAK
-case 127:
-YY_RULE_SETUP
-#line 213 "Gmsh.l"
-return tNurbs;
- YY_BREAK
-case 128:
-YY_RULE_SETUP
-#line 215 "Gmsh.l"
-return tOnelabAction;
- YY_BREAK
-case 129:
-YY_RULE_SETUP
-#line 216 "Gmsh.l"
-return tOnelabRun;
- YY_BREAK
-case 130:
-YY_RULE_SETUP
-#line 217 "Gmsh.l"
-return tNurbsOrder;
- YY_BREAK
-case 131:
-YY_RULE_SETUP
-#line 219 "Gmsh.l"
-return tPeriodic;
- YY_BREAK
-case 132:
-YY_RULE_SETUP
-#line 220 "Gmsh.l"
-return tPhysical;
- YY_BREAK
-case 133:
-YY_RULE_SETUP
-#line 221 "Gmsh.l"
-return tPi;
- YY_BREAK
-case 134:
-YY_RULE_SETUP
-#line 222 "Gmsh.l"
-return tPlane;
- YY_BREAK
-case 135:
-YY_RULE_SETUP
-#line 223 "Gmsh.l"
-return tPoint;
- YY_BREAK
-case 136:
-YY_RULE_SETUP
-#line 224 "Gmsh.l"
-return tParametric;
- YY_BREAK
-case 137:
-YY_RULE_SETUP
-#line 225 "Gmsh.l"
-return tPolarSphere;
- YY_BREAK
-case 138:
-YY_RULE_SETUP
-#line 226 "Gmsh.l"
-return tPrintf;
- YY_BREAK
-case 139:
-YY_RULE_SETUP
-#line 227 "Gmsh.l"
-return tPlugin;
- YY_BREAK
-case 140:
-YY_RULE_SETUP
-#line 229 "Gmsh.l"
-return tQuadTriAddVerts;
- YY_BREAK
-case 141:
-YY_RULE_SETUP
-#line 230 "Gmsh.l"
-return tQuadTriNoNewVerts;
- YY_BREAK
-case 142:
-YY_RULE_SETUP
-#line 232 "Gmsh.l"
-return tQuadTriDbl;
- YY_BREAK
-case 143:
-YY_RULE_SETUP
-#line 233 "Gmsh.l"
-return tQuadTriSngl;
- YY_BREAK
-case 144:
-YY_RULE_SETUP
-#line 235 "Gmsh.l"
-return tRecombine;
- YY_BREAK
-case 145:
-YY_RULE_SETUP
-#line 236 "Gmsh.l"
-return tRecombLaterals;
- YY_BREAK
-case 146:
-YY_RULE_SETUP
-#line 237 "Gmsh.l"
-return tRecursive;
- YY_BREAK
-case 147:
-YY_RULE_SETUP
-#line 238 "Gmsh.l"
-return tRotate;
- YY_BREAK
-case 148:
-YY_RULE_SETUP
-#line 239 "Gmsh.l"
-return tRound;
- YY_BREAK
-case 149:
-YY_RULE_SETUP
-#line 240 "Gmsh.l"
-return tRuled;
- YY_BREAK
-case 150:
-YY_RULE_SETUP
-#line 241 "Gmsh.l"
-return tRand;
- YY_BREAK
-case 151:
-YY_RULE_SETUP
-#line 242 "Gmsh.l"
-return tRefineMesh;
- YY_BREAK
-case 152:
-YY_RULE_SETUP
-#line 243 "Gmsh.l"
-return tRelocateMesh;
- YY_BREAK
-case 153:
-YY_RULE_SETUP
-#line 244 "Gmsh.l"
-return tReturn;
- YY_BREAK
-case 154:
-YY_RULE_SETUP
-#line 245 "Gmsh.l"
-return tReverse;
- YY_BREAK
-case 155:
-YY_RULE_SETUP
-#line 247 "Gmsh.l"
-return tScaleLast;
- YY_BREAK
-case 156:
-YY_RULE_SETUP
-#line 248 "Gmsh.l"
-return tSmoother;
- YY_BREAK
-case 157:
-YY_RULE_SETUP
-#line 249 "Gmsh.l"
-return tSetOrder;
- YY_BREAK
-case 158:
-YY_RULE_SETUP
-#line 250 "Gmsh.l"
-return tSetNumber;
- YY_BREAK
-case 159:
-YY_RULE_SETUP
-#line 251 "Gmsh.l"
-return tSetString;
- YY_BREAK
-case 160:
-YY_RULE_SETUP
-#line 252 "Gmsh.l"
-return tSetPartition;
- YY_BREAK
-case 161:
-YY_RULE_SETUP
-#line 253 "Gmsh.l"
-return tSqrt;
- YY_BREAK
-case 162:
-YY_RULE_SETUP
-#line 254 "Gmsh.l"
-return tSin;
- YY_BREAK
-case 163:
-YY_RULE_SETUP
-#line 255 "Gmsh.l"
-return tSinh;
- YY_BREAK
-case 164:
-YY_RULE_SETUP
-#line 256 "Gmsh.l"
-return tSphere;
- YY_BREAK
-case 165:
-YY_RULE_SETUP
-#line 257 "Gmsh.l"
-return tSpline;
- YY_BREAK
-case 166:
-YY_RULE_SETUP
-#line 258 "Gmsh.l"
-return tSplit;
- YY_BREAK
-case 167:
-YY_RULE_SETUP
-#line 259 "Gmsh.l"
-return tSurface;
- YY_BREAK
-case 168:
-YY_RULE_SETUP
-#line 260 "Gmsh.l"
-return tStr;
- YY_BREAK
-case 169:
-YY_RULE_SETUP
-#line 261 "Gmsh.l"
-return tSprintf;
- YY_BREAK
-case 170:
-YY_RULE_SETUP
-#line 262 "Gmsh.l"
-return tStrCat;
- YY_BREAK
-case 171:
-YY_RULE_SETUP
-#line 263 "Gmsh.l"
-return tStrReplace;
- YY_BREAK
-case 172:
-YY_RULE_SETUP
-#line 264 "Gmsh.l"
-return tStrPrefix;
- YY_BREAK
-case 173:
-YY_RULE_SETUP
-#line 265 "Gmsh.l"
-return tStrRelative;
- YY_BREAK
-case 174:
-YY_RULE_SETUP
-#line 266 "Gmsh.l"
-return tStrFind;
- YY_BREAK
-case 175:
-YY_RULE_SETUP
-#line 267 "Gmsh.l"
-return tStrCmp;
- YY_BREAK
-case 176:
-YY_RULE_SETUP
-#line 268 "Gmsh.l"
-return tStrChoice;
- YY_BREAK
-case 177:
-YY_RULE_SETUP
-#line 269 "Gmsh.l"
-return tShow;
- YY_BREAK
-case 178:
-YY_RULE_SETUP
-#line 270 "Gmsh.l"
-return tSymmetry;
- YY_BREAK
-case 179:
-YY_RULE_SETUP
-#line 271 "Gmsh.l"
-return tSyncModel;
- YY_BREAK
-case 180:
-YY_RULE_SETUP
-#line 273 "Gmsh.l"
-return tText2D;
- YY_BREAK
-case 181:
-YY_RULE_SETUP
-#line 274 "Gmsh.l"
-return tText3D;
- YY_BREAK
-case 182:
-YY_RULE_SETUP
-#line 275 "Gmsh.l"
-return tTextAttributes;
- YY_BREAK
-case 183:
-YY_RULE_SETUP
-#line 276 "Gmsh.l"
-return tTime;
- YY_BREAK
-case 184:
-YY_RULE_SETUP
-#line 277 "Gmsh.l"
-return tTransfinite;
- YY_BREAK
-case 185:
-YY_RULE_SETUP
-#line 278 "Gmsh.l"
-return tTransfQuadTri;
- YY_BREAK
-case 186:
-YY_RULE_SETUP
-#line 279 "Gmsh.l"
-return tTranslate;
- YY_BREAK
-case 187:
-YY_RULE_SETUP
-#line 280 "Gmsh.l"
-return tTanh;
- YY_BREAK
-case 188:
-YY_RULE_SETUP
-#line 281 "Gmsh.l"
-return tTan;
- YY_BREAK
-case 189:
-YY_RULE_SETUP
-#line 282 "Gmsh.l"
-return tToday;
- YY_BREAK
-case 190:
-YY_RULE_SETUP
-#line 283 "Gmsh.l"
-return tTotalMemory;
- YY_BREAK
-case 191:
-YY_RULE_SETUP
-#line 285 "Gmsh.l"
-return tUndefineConstant;
- YY_BREAK
-case 192:
-YY_RULE_SETUP
-#line 286 "Gmsh.l"
-return tUsing;
- YY_BREAK
-case 193:
-YY_RULE_SETUP
-#line 288 "Gmsh.l"
-return tVolume;
- YY_BREAK
-case 194:
-#line 291 "Gmsh.l"
-case 195:
-#line 292 "Gmsh.l"
-case 196:
-#line 293 "Gmsh.l"
-case 197:
-YY_RULE_SETUP
-#line 293 "Gmsh.l"
-{ gmsh_yylval.d = atof((char *)gmsh_yytext); return tDOUBLE; }
- YY_BREAK
-case 198:
-YY_RULE_SETUP
-#line 295 "Gmsh.l"
-{ gmsh_yylval.c = strsave((char*)gmsh_yytext); return tSTRING; }
- YY_BREAK
-case 199:
-YY_RULE_SETUP
-#line 297 "Gmsh.l"
-return gmsh_yytext[0];
- YY_BREAK
-case 200:
-YY_RULE_SETUP
-#line 299 "Gmsh.l"
-ECHO;
- YY_BREAK
-#line 2393 "Gmsh.yy.cpp"
-case YY_STATE_EOF(INITIAL):
- yyterminate();
-
- case YY_END_OF_BUFFER:
- {
- /* Amount of text matched not including the EOB char. */
- int yy_amount_of_matched_text = (int) (yy_cp - (yytext_ptr)) - 1;
-
- /* Undo the effects of YY_DO_BEFORE_ACTION. */
- *yy_cp = (yy_hold_char);
- YY_RESTORE_YY_MORE_OFFSET
-
- if ( YY_CURRENT_BUFFER_LVALUE->yy_buffer_status == YY_BUFFER_NEW )
- {
- /* We're scanning a new file or input source. It's
- * possible that this happened because the user
- * just pointed gmsh_yyin at a new source and called
- * gmsh_yylex(). If so, then we have to assure
- * consistency between YY_CURRENT_BUFFER and our
- * globals. Here is the right place to do so, because
- * this is the first action (other than possibly a
- * back-up) that will match for the new input source.
- */
- (yy_n_chars) = YY_CURRENT_BUFFER_LVALUE->yy_n_chars;
- YY_CURRENT_BUFFER_LVALUE->yy_input_file = gmsh_yyin;
- YY_CURRENT_BUFFER_LVALUE->yy_buffer_status = YY_BUFFER_NORMAL;
- }
-
- /* Note that here we test for yy_c_buf_p "<=" to the position
- * of the first EOB in the buffer, since yy_c_buf_p will
- * already have been incremented past the NUL character
- * (since all states make transitions on EOB to the
- * end-of-buffer state). Contrast this with the test
- * in input().
- */
- if ( (yy_c_buf_p) <= &YY_CURRENT_BUFFER_LVALUE->yy_ch_buf[(yy_n_chars)] )
- { /* This was really a NUL. */
- yy_state_type yy_next_state;
-
- (yy_c_buf_p) = (yytext_ptr) + yy_amount_of_matched_text;
-
- yy_current_state = yy_get_previous_state( );
-
- /* Okay, we're now positioned to make the NUL
- * transition. We couldn't have
- * yy_get_previous_state() go ahead and do it
- * for us because it doesn't know how to deal
- * with the possibility of jamming (and we don't
- * want to build jamming into it because then it
- * will run more slowly).
- */
-
- yy_next_state = yy_try_NUL_trans( yy_current_state );
-
- yy_bp = (yytext_ptr) + YY_MORE_ADJ;
-
- if ( yy_next_state )
- {
- /* Consume the NUL. */
- yy_cp = ++(yy_c_buf_p);
- yy_current_state = yy_next_state;
- goto yy_match;
- }
-
- else
- {
- yy_cp = (yy_c_buf_p);
- goto yy_find_action;
- }
- }
-
- else switch ( yy_get_next_buffer( ) )
- {
- case EOB_ACT_END_OF_FILE:
- {
- (yy_did_buffer_switch_on_eof) = 0;
-
- if ( gmsh_yywrap( ) )
- {
- /* Note: because we've taken care in
- * yy_get_next_buffer() to have set up
- * gmsh_yytext, we can now set up
- * yy_c_buf_p so that if some total
- * hoser (like flex itself) wants to
- * call the scanner after we return the
- * YY_NULL, it'll still work - another
- * YY_NULL will get returned.
- */
- (yy_c_buf_p) = (yytext_ptr) + YY_MORE_ADJ;
-
- yy_act = YY_STATE_EOF(YY_START);
- goto do_action;
- }
-
- else
- {
- if ( ! (yy_did_buffer_switch_on_eof) )
- YY_NEW_FILE;
- }
- break;
- }
-
- case EOB_ACT_CONTINUE_SCAN:
- (yy_c_buf_p) =
- (yytext_ptr) + yy_amount_of_matched_text;
-
- yy_current_state = yy_get_previous_state( );
-
- yy_cp = (yy_c_buf_p);
- yy_bp = (yytext_ptr) + YY_MORE_ADJ;
- goto yy_match;
-
- case EOB_ACT_LAST_MATCH:
- (yy_c_buf_p) =
- &YY_CURRENT_BUFFER_LVALUE->yy_ch_buf[(yy_n_chars)];
-
- yy_current_state = yy_get_previous_state( );
-
- yy_cp = (yy_c_buf_p);
- yy_bp = (yytext_ptr) + YY_MORE_ADJ;
- goto yy_find_action;
- }
- break;
- }
-
- default:
- YY_FATAL_ERROR(
- "fatal flex scanner internal error--no action found" );
- } /* end of action switch */
- } /* end of scanning one token */
-} /* end of gmsh_yylex */
-
-/* yy_get_next_buffer - try to read in a new buffer
- *
- * Returns a code representing an action:
- * EOB_ACT_LAST_MATCH -
- * EOB_ACT_CONTINUE_SCAN - continue scanning from current position
- * EOB_ACT_END_OF_FILE - end of file
- */
-static int yy_get_next_buffer (void)
-{
- register char *dest = YY_CURRENT_BUFFER_LVALUE->yy_ch_buf;
- register char *source = (yytext_ptr);
- register int number_to_move, i;
- int ret_val;
-
- if ( (yy_c_buf_p) > &YY_CURRENT_BUFFER_LVALUE->yy_ch_buf[(yy_n_chars) + 1] )
- YY_FATAL_ERROR(
- "fatal flex scanner internal error--end of buffer missed" );
-
- if ( YY_CURRENT_BUFFER_LVALUE->yy_fill_buffer == 0 )
- { /* Don't try to fill the buffer, so this is an EOF. */
- if ( (yy_c_buf_p) - (yytext_ptr) - YY_MORE_ADJ == 1 )
- {
- /* We matched a single character, the EOB, so
- * treat this as a final EOF.
- */
- return EOB_ACT_END_OF_FILE;
- }
-
- else
- {
- /* We matched some text prior to the EOB, first
- * process it.
- */
- return EOB_ACT_LAST_MATCH;
- }
- }
-
- /* Try to read more data. */
-
- /* First move last chars to start of buffer. */
- number_to_move = (int) ((yy_c_buf_p) - (yytext_ptr)) - 1;
-
- for ( i = 0; i < number_to_move; ++i )
- *(dest++) = *(source++);
-
- if ( YY_CURRENT_BUFFER_LVALUE->yy_buffer_status == YY_BUFFER_EOF_PENDING )
- /* don't do the read, it's not guaranteed to return an EOF,
- * just force an EOF
- */
- YY_CURRENT_BUFFER_LVALUE->yy_n_chars = (yy_n_chars) = 0;
-
- else
- {
- yy_size_t num_to_read =
- YY_CURRENT_BUFFER_LVALUE->yy_buf_size - number_to_move - 1;
-
- while ( num_to_read <= 0 )
- { /* Not enough room in the buffer - grow it. */
-
- /* just a shorter name for the current buffer */
- YY_BUFFER_STATE b = YY_CURRENT_BUFFER_LVALUE;
-
- int yy_c_buf_p_offset =
- (int) ((yy_c_buf_p) - b->yy_ch_buf);
-
- if ( b->yy_is_our_buffer )
- {
- yy_size_t new_size = b->yy_buf_size * 2;
-
- if ( new_size <= 0 )
- b->yy_buf_size += b->yy_buf_size / 8;
- else
- b->yy_buf_size *= 2;
-
- b->yy_ch_buf = (char *)
- /* Include room in for 2 EOB chars. */
- gmsh_yyrealloc((void *) b->yy_ch_buf,b->yy_buf_size + 2 );
- }
- else
- /* Can't grow it, we don't own it. */
- b->yy_ch_buf = 0;
-
- if ( ! b->yy_ch_buf )
- YY_FATAL_ERROR(
- "fatal error - scanner input buffer overflow" );
-
- (yy_c_buf_p) = &b->yy_ch_buf[yy_c_buf_p_offset];
-
- num_to_read = YY_CURRENT_BUFFER_LVALUE->yy_buf_size -
- number_to_move - 1;
-
- }
-
- if ( num_to_read > YY_READ_BUF_SIZE )
- num_to_read = YY_READ_BUF_SIZE;
-
- /* Read in more data. */
- YY_INPUT( (&YY_CURRENT_BUFFER_LVALUE->yy_ch_buf[number_to_move]),
- (yy_n_chars), num_to_read );
-
- YY_CURRENT_BUFFER_LVALUE->yy_n_chars = (yy_n_chars);
- }
-
- if ( (yy_n_chars) == 0 )
- {
- if ( number_to_move == YY_MORE_ADJ )
- {
- ret_val = EOB_ACT_END_OF_FILE;
- gmsh_yyrestart(gmsh_yyin );
- }
-
- else
- {
- ret_val = EOB_ACT_LAST_MATCH;
- YY_CURRENT_BUFFER_LVALUE->yy_buffer_status =
- YY_BUFFER_EOF_PENDING;
- }
- }
-
- else
- ret_val = EOB_ACT_CONTINUE_SCAN;
-
- if ((yy_size_t) ((yy_n_chars) + number_to_move) > YY_CURRENT_BUFFER_LVALUE->yy_buf_size) {
- /* Extend the array by 50%, plus the number we really need. */
- yy_size_t new_size = (yy_n_chars) + number_to_move + ((yy_n_chars) >> 1);
- YY_CURRENT_BUFFER_LVALUE->yy_ch_buf = (char *) gmsh_yyrealloc((void *) YY_CURRENT_BUFFER_LVALUE->yy_ch_buf,new_size );
- if ( ! YY_CURRENT_BUFFER_LVALUE->yy_ch_buf )
- YY_FATAL_ERROR( "out of dynamic memory in yy_get_next_buffer()" );
- }
-
- (yy_n_chars) += number_to_move;
- YY_CURRENT_BUFFER_LVALUE->yy_ch_buf[(yy_n_chars)] = YY_END_OF_BUFFER_CHAR;
- YY_CURRENT_BUFFER_LVALUE->yy_ch_buf[(yy_n_chars) + 1] = YY_END_OF_BUFFER_CHAR;
-
- (yytext_ptr) = &YY_CURRENT_BUFFER_LVALUE->yy_ch_buf[0];
-
- return ret_val;
-}
-
-/* yy_get_previous_state - get the state just before the EOB char was reached */
-
- static yy_state_type yy_get_previous_state (void)
-{
- register yy_state_type yy_current_state;
- register char *yy_cp;
-
- yy_current_state = (yy_start);
-
- for ( yy_cp = (yytext_ptr) + YY_MORE_ADJ; yy_cp < (yy_c_buf_p); ++yy_cp )
- {
- register YY_CHAR yy_c = (*yy_cp ? yy_ec[YY_SC_TO_UI(*yy_cp)] : 1);
- if ( yy_accept[yy_current_state] )
- {
- (yy_last_accepting_state) = yy_current_state;
- (yy_last_accepting_cpos) = yy_cp;
- }
- while ( yy_chk[yy_base[yy_current_state] + yy_c] != yy_current_state )
- {
- yy_current_state = (int) yy_def[yy_current_state];
- if ( yy_current_state >= 1002 )
- yy_c = yy_meta[(unsigned int) yy_c];
- }
- yy_current_state = yy_nxt[yy_base[yy_current_state] + (unsigned int) yy_c];
- }
-
- return yy_current_state;
-}
-
-/* yy_try_NUL_trans - try to make a transition on the NUL character
- *
- * synopsis
- * next_state = yy_try_NUL_trans( current_state );
- */
- static yy_state_type yy_try_NUL_trans (yy_state_type yy_current_state )
-{
- register int yy_is_jam;
- register char *yy_cp = (yy_c_buf_p);
-
- register YY_CHAR yy_c = 1;
- if ( yy_accept[yy_current_state] )
- {
- (yy_last_accepting_state) = yy_current_state;
- (yy_last_accepting_cpos) = yy_cp;
- }
- while ( yy_chk[yy_base[yy_current_state] + yy_c] != yy_current_state )
- {
- yy_current_state = (int) yy_def[yy_current_state];
- if ( yy_current_state >= 1002 )
- yy_c = yy_meta[(unsigned int) yy_c];
- }
- yy_current_state = yy_nxt[yy_base[yy_current_state] + (unsigned int) yy_c];
- yy_is_jam = (yy_current_state == 1001);
-
- return yy_is_jam ? 0 : yy_current_state;
-}
-
- static void yyunput (int c, register char * yy_bp )
-{
- register char *yy_cp;
-
- yy_cp = (yy_c_buf_p);
-
- /* undo effects of setting up gmsh_yytext */
- *yy_cp = (yy_hold_char);
-
- if ( yy_cp < YY_CURRENT_BUFFER_LVALUE->yy_ch_buf + 2 )
- { /* need to shift things up to make room */
- /* +2 for EOB chars. */
- register yy_size_t number_to_move = (yy_n_chars) + 2;
- register char *dest = &YY_CURRENT_BUFFER_LVALUE->yy_ch_buf[
- YY_CURRENT_BUFFER_LVALUE->yy_buf_size + 2];
- register char *source =
- &YY_CURRENT_BUFFER_LVALUE->yy_ch_buf[number_to_move];
-
- while ( source > YY_CURRENT_BUFFER_LVALUE->yy_ch_buf )
- *--dest = *--source;
-
- yy_cp += (int) (dest - source);
- yy_bp += (int) (dest - source);
- YY_CURRENT_BUFFER_LVALUE->yy_n_chars =
- (yy_n_chars) = YY_CURRENT_BUFFER_LVALUE->yy_buf_size;
-
- if ( yy_cp < YY_CURRENT_BUFFER_LVALUE->yy_ch_buf + 2 )
- YY_FATAL_ERROR( "flex scanner push-back overflow" );
- }
-
- *--yy_cp = (char) c;
-
- (yytext_ptr) = yy_bp;
- (yy_hold_char) = *yy_cp;
- (yy_c_buf_p) = yy_cp;
-}
-
-#ifndef YY_NO_INPUT
-#ifdef __cplusplus
- static int yyinput (void)
-#else
- static int input (void)
-#endif
-
-{
- int c;
-
- *(yy_c_buf_p) = (yy_hold_char);
-
- if ( *(yy_c_buf_p) == YY_END_OF_BUFFER_CHAR )
- {
- /* yy_c_buf_p now points to the character we want to return.
- * If this occurs *before* the EOB characters, then it's a
- * valid NUL; if not, then we've hit the end of the buffer.
- */
- if ( (yy_c_buf_p) < &YY_CURRENT_BUFFER_LVALUE->yy_ch_buf[(yy_n_chars)] )
- /* This was really a NUL. */
- *(yy_c_buf_p) = '\0';
-
- else
- { /* need more input */
- yy_size_t offset = (yy_c_buf_p) - (yytext_ptr);
- ++(yy_c_buf_p);
-
- switch ( yy_get_next_buffer( ) )
- {
- case EOB_ACT_LAST_MATCH:
- /* This happens because yy_g_n_b()
- * sees that we've accumulated a
- * token and flags that we need to
- * try matching the token before
- * proceeding. But for input(),
- * there's no matching to consider.
- * So convert the EOB_ACT_LAST_MATCH
- * to EOB_ACT_END_OF_FILE.
- */
-
- /* Reset buffer status. */
- gmsh_yyrestart(gmsh_yyin );
-
- /*FALLTHROUGH*/
-
- case EOB_ACT_END_OF_FILE:
- {
- if ( gmsh_yywrap( ) )
- return EOF;
-
- if ( ! (yy_did_buffer_switch_on_eof) )
- YY_NEW_FILE;
-#ifdef __cplusplus
- return yyinput();
-#else
- return input();
-#endif
- }
-
- case EOB_ACT_CONTINUE_SCAN:
- (yy_c_buf_p) = (yytext_ptr) + offset;
- break;
- }
- }
- }
-
- c = *(unsigned char *) (yy_c_buf_p); /* cast for 8-bit char's */
- *(yy_c_buf_p) = '\0'; /* preserve gmsh_yytext */
- (yy_hold_char) = *++(yy_c_buf_p);
-
- return c;
-}
-#endif /* ifndef YY_NO_INPUT */
-
-/** Immediately switch to a different input stream.
- * @param input_file A readable stream.
- *
- * @note This function does not reset the start condition to @c INITIAL .
- */
- void gmsh_yyrestart (FILE * input_file )
-{
-
- if ( ! YY_CURRENT_BUFFER ){
- gmsh_yyensure_buffer_stack ();
- YY_CURRENT_BUFFER_LVALUE =
- gmsh_yy_create_buffer(gmsh_yyin,YY_BUF_SIZE );
- }
-
- gmsh_yy_init_buffer(YY_CURRENT_BUFFER,input_file );
- gmsh_yy_load_buffer_state( );
-}
-
-/** Switch to a different input buffer.
- * @param new_buffer The new input buffer.
- *
- */
- void gmsh_yy_switch_to_buffer (YY_BUFFER_STATE new_buffer )
-{
-
- /* TODO. We should be able to replace this entire function body
- * with
- * gmsh_yypop_buffer_state();
- * gmsh_yypush_buffer_state(new_buffer);
- */
- gmsh_yyensure_buffer_stack ();
- if ( YY_CURRENT_BUFFER == new_buffer )
- return;
-
- if ( YY_CURRENT_BUFFER )
- {
- /* Flush out information for old buffer. */
- *(yy_c_buf_p) = (yy_hold_char);
- YY_CURRENT_BUFFER_LVALUE->yy_buf_pos = (yy_c_buf_p);
- YY_CURRENT_BUFFER_LVALUE->yy_n_chars = (yy_n_chars);
- }
-
- YY_CURRENT_BUFFER_LVALUE = new_buffer;
- gmsh_yy_load_buffer_state( );
-
- /* We don't actually know whether we did this switch during
- * EOF (gmsh_yywrap()) processing, but the only time this flag
- * is looked at is after gmsh_yywrap() is called, so it's safe
- * to go ahead and always set it.
- */
- (yy_did_buffer_switch_on_eof) = 1;
-}
-
-static void gmsh_yy_load_buffer_state (void)
-{
- (yy_n_chars) = YY_CURRENT_BUFFER_LVALUE->yy_n_chars;
- (yytext_ptr) = (yy_c_buf_p) = YY_CURRENT_BUFFER_LVALUE->yy_buf_pos;
- gmsh_yyin = YY_CURRENT_BUFFER_LVALUE->yy_input_file;
- (yy_hold_char) = *(yy_c_buf_p);
-}
-
-/** Allocate and initialize an input buffer state.
- * @param file A readable stream.
- * @param size The character buffer size in bytes. When in doubt, use @c YY_BUF_SIZE.
- *
- * @return the allocated buffer state.
- */
- YY_BUFFER_STATE gmsh_yy_create_buffer (FILE * file, int size )
-{
- YY_BUFFER_STATE b;
-
- b = (YY_BUFFER_STATE) gmsh_yyalloc(sizeof( struct yy_buffer_state ) );
- if ( ! b )
- YY_FATAL_ERROR( "out of dynamic memory in gmsh_yy_create_buffer()" );
-
- b->yy_buf_size = size;
-
- /* yy_ch_buf has to be 2 characters longer than the size given because
- * we need to put in 2 end-of-buffer characters.
- */
- b->yy_ch_buf = (char *) gmsh_yyalloc(b->yy_buf_size + 2 );
- if ( ! b->yy_ch_buf )
- YY_FATAL_ERROR( "out of dynamic memory in gmsh_yy_create_buffer()" );
-
- b->yy_is_our_buffer = 1;
-
- gmsh_yy_init_buffer(b,file );
-
- return b;
-}
-
-/** Destroy the buffer.
- * @param b a buffer created with gmsh_yy_create_buffer()
- *
- */
- void gmsh_yy_delete_buffer (YY_BUFFER_STATE b )
-{
-
- if ( ! b )
- return;
-
- if ( b == YY_CURRENT_BUFFER ) /* Not sure if we should pop here. */
- YY_CURRENT_BUFFER_LVALUE = (YY_BUFFER_STATE) 0;
-
- if ( b->yy_is_our_buffer )
- gmsh_yyfree((void *) b->yy_ch_buf );
-
- gmsh_yyfree((void *) b );
-}
-
-/* Initializes or reinitializes a buffer.
- * This function is sometimes called more than once on the same buffer,
- * such as during a gmsh_yyrestart() or at EOF.
- */
- static void gmsh_yy_init_buffer (YY_BUFFER_STATE b, FILE * file )
-
-{
- int oerrno = errno;
-
- gmsh_yy_flush_buffer(b );
-
- b->yy_input_file = file;
- b->yy_fill_buffer = 1;
-
- /* If b is the current buffer, then gmsh_yy_init_buffer was _probably_
- * called from gmsh_yyrestart() or through yy_get_next_buffer.
- * In that case, we don't want to reset the lineno or column.
- */
- if (b != YY_CURRENT_BUFFER){
- b->yy_bs_lineno = 1;
- b->yy_bs_column = 0;
- }
-
- b->yy_is_interactive = file ? (isatty( fileno(file) ) > 0) : 0;
-
- errno = oerrno;
-}
-
-/** Discard all buffered characters. On the next scan, YY_INPUT will be called.
- * @param b the buffer state to be flushed, usually @c YY_CURRENT_BUFFER.
- *
- */
- void gmsh_yy_flush_buffer (YY_BUFFER_STATE b )
-{
- if ( ! b )
- return;
-
- b->yy_n_chars = 0;
-
- /* We always need two end-of-buffer characters. The first causes
- * a transition to the end-of-buffer state. The second causes
- * a jam in that state.
- */
- b->yy_ch_buf[0] = YY_END_OF_BUFFER_CHAR;
- b->yy_ch_buf[1] = YY_END_OF_BUFFER_CHAR;
-
- b->yy_buf_pos = &b->yy_ch_buf[0];
-
- b->yy_at_bol = 1;
- b->yy_buffer_status = YY_BUFFER_NEW;
-
- if ( b == YY_CURRENT_BUFFER )
- gmsh_yy_load_buffer_state( );
-}
-
-/** Pushes the new state onto the stack. The new state becomes
- * the current state. This function will allocate the stack
- * if necessary.
- * @param new_buffer The new state.
- *
- */
-void gmsh_yypush_buffer_state (YY_BUFFER_STATE new_buffer )
-{
- if (new_buffer == NULL)
- return;
-
- gmsh_yyensure_buffer_stack();
-
- /* This block is copied from gmsh_yy_switch_to_buffer. */
- if ( YY_CURRENT_BUFFER )
- {
- /* Flush out information for old buffer. */
- *(yy_c_buf_p) = (yy_hold_char);
- YY_CURRENT_BUFFER_LVALUE->yy_buf_pos = (yy_c_buf_p);
- YY_CURRENT_BUFFER_LVALUE->yy_n_chars = (yy_n_chars);
- }
-
- /* Only push if top exists. Otherwise, replace top. */
- if (YY_CURRENT_BUFFER)
- (yy_buffer_stack_top)++;
- YY_CURRENT_BUFFER_LVALUE = new_buffer;
-
- /* copied from gmsh_yy_switch_to_buffer. */
- gmsh_yy_load_buffer_state( );
- (yy_did_buffer_switch_on_eof) = 1;
-}
-
-/** Removes and deletes the top of the stack, if present.
- * The next element becomes the new top.
- *
- */
-void gmsh_yypop_buffer_state (void)
-{
- if (!YY_CURRENT_BUFFER)
- return;
-
- gmsh_yy_delete_buffer(YY_CURRENT_BUFFER );
- YY_CURRENT_BUFFER_LVALUE = NULL;
- if ((yy_buffer_stack_top) > 0)
- --(yy_buffer_stack_top);
-
- if (YY_CURRENT_BUFFER) {
- gmsh_yy_load_buffer_state( );
- (yy_did_buffer_switch_on_eof) = 1;
- }
-}
-
-/* Allocates the stack if it does not exist.
- * Guarantees space for at least one push.
- */
-static void gmsh_yyensure_buffer_stack (void)
-{
- yy_size_t num_to_alloc;
-
- if (!(yy_buffer_stack)) {
-
- /* First allocation is just for 2 elements, since we don't know if this
- * scanner will even need a stack. We use 2 instead of 1 to avoid an
- * immediate realloc on the next call.
- */
- num_to_alloc = 1;
- (yy_buffer_stack) = (struct yy_buffer_state**)gmsh_yyalloc
- (num_to_alloc * sizeof(struct yy_buffer_state*)
- );
- if ( ! (yy_buffer_stack) )
- YY_FATAL_ERROR( "out of dynamic memory in gmsh_yyensure_buffer_stack()" );
-
- memset((yy_buffer_stack), 0, num_to_alloc * sizeof(struct yy_buffer_state*));
-
- (yy_buffer_stack_max) = num_to_alloc;
- (yy_buffer_stack_top) = 0;
- return;
- }
-
- if ((yy_buffer_stack_top) >= ((yy_buffer_stack_max)) - 1){
-
- /* Increase the buffer to prepare for a possible push. */
- int grow_size = 8 /* arbitrary grow size */;
-
- num_to_alloc = (yy_buffer_stack_max) + grow_size;
- (yy_buffer_stack) = (struct yy_buffer_state**)gmsh_yyrealloc
- ((yy_buffer_stack),
- num_to_alloc * sizeof(struct yy_buffer_state*)
- );
- if ( ! (yy_buffer_stack) )
- YY_FATAL_ERROR( "out of dynamic memory in gmsh_yyensure_buffer_stack()" );
-
- /* zero only the new slots.*/
- memset((yy_buffer_stack) + (yy_buffer_stack_max), 0, grow_size * sizeof(struct yy_buffer_state*));
- (yy_buffer_stack_max) = num_to_alloc;
- }
-}
-
-/** Setup the input buffer state to scan directly from a user-specified character buffer.
- * @param base the character buffer
- * @param size the size in bytes of the character buffer
- *
- * @return the newly allocated buffer state object.
- */
-YY_BUFFER_STATE gmsh_yy_scan_buffer (char * base, yy_size_t size )
-{
- YY_BUFFER_STATE b;
-
- if ( size < 2 ||
- base[size-2] != YY_END_OF_BUFFER_CHAR ||
- base[size-1] != YY_END_OF_BUFFER_CHAR )
- /* They forgot to leave room for the EOB's. */
- return 0;
-
- b = (YY_BUFFER_STATE) gmsh_yyalloc(sizeof( struct yy_buffer_state ) );
- if ( ! b )
- YY_FATAL_ERROR( "out of dynamic memory in gmsh_yy_scan_buffer()" );
-
- b->yy_buf_size = size - 2; /* "- 2" to take care of EOB's */
- b->yy_buf_pos = b->yy_ch_buf = base;
- b->yy_is_our_buffer = 0;
- b->yy_input_file = 0;
- b->yy_n_chars = b->yy_buf_size;
- b->yy_is_interactive = 0;
- b->yy_at_bol = 1;
- b->yy_fill_buffer = 0;
- b->yy_buffer_status = YY_BUFFER_NEW;
-
- gmsh_yy_switch_to_buffer(b );
-
- return b;
-}
-
-/** Setup the input buffer state to scan a string. The next call to gmsh_yylex() will
- * scan from a @e copy of @a str.
- * @param yystr a NUL-terminated string to scan
- *
- * @return the newly allocated buffer state object.
- * @note If you want to scan bytes that may contain NUL values, then use
- * gmsh_yy_scan_bytes() instead.
- */
-YY_BUFFER_STATE gmsh_yy_scan_string (yyconst char * yystr )
-{
-
- return gmsh_yy_scan_bytes(yystr,strlen(yystr) );
-}
-
-/** Setup the input buffer state to scan the given bytes. The next call to gmsh_yylex() will
- * scan from a @e copy of @a bytes.
- * @param yybytes the byte buffer to scan
- * @param _yybytes_len the number of bytes in the buffer pointed to by @a bytes.
- *
- * @return the newly allocated buffer state object.
- */
-YY_BUFFER_STATE gmsh_yy_scan_bytes (yyconst char * yybytes, yy_size_t _yybytes_len )
-{
- YY_BUFFER_STATE b;
- char *buf;
- yy_size_t n;
- int i;
-
- /* Get memory for full buffer, including space for trailing EOB's. */
- n = _yybytes_len + 2;
- buf = (char *) gmsh_yyalloc(n );
- if ( ! buf )
- YY_FATAL_ERROR( "out of dynamic memory in gmsh_yy_scan_bytes()" );
-
- for ( i = 0; i < _yybytes_len; ++i )
- buf[i] = yybytes[i];
-
- buf[_yybytes_len] = buf[_yybytes_len+1] = YY_END_OF_BUFFER_CHAR;
-
- b = gmsh_yy_scan_buffer(buf,n );
- if ( ! b )
- YY_FATAL_ERROR( "bad buffer in gmsh_yy_scan_bytes()" );
-
- /* It's okay to grow etc. this buffer, and we should throw it
- * away when we're done.
- */
- b->yy_is_our_buffer = 1;
-
- return b;
-}
-
-#ifndef YY_EXIT_FAILURE
-#define YY_EXIT_FAILURE 2
-#endif
-
-static void yy_fatal_error (yyconst char* msg )
-{
- (void) fprintf( stderr, "%s\n", msg );
- exit( YY_EXIT_FAILURE );
-}
-
-/* Redefine yyless() so it works in section 3 code. */
-
-#undef yyless
-#define yyless(n) \
- do \
- { \
- /* Undo effects of setting up gmsh_yytext. */ \
- int yyless_macro_arg = (n); \
- YY_LESS_LINENO(yyless_macro_arg);\
- gmsh_yytext[gmsh_yyleng] = (yy_hold_char); \
- (yy_c_buf_p) = gmsh_yytext + yyless_macro_arg; \
- (yy_hold_char) = *(yy_c_buf_p); \
- *(yy_c_buf_p) = '\0'; \
- gmsh_yyleng = yyless_macro_arg; \
- } \
- while ( 0 )
-
-/* Accessor methods (get/set functions) to struct members. */
-
-/** Get the current line number.
- *
- */
-int gmsh_yyget_lineno (void)
-{
-
- return gmsh_yylineno;
-}
-
-/** Get the input stream.
- *
- */
-FILE *gmsh_yyget_in (void)
-{
- return gmsh_yyin;
-}
-
-/** Get the output stream.
- *
- */
-FILE *gmsh_yyget_out (void)
-{
- return gmsh_yyout;
-}
-
-/** Get the length of the current token.
- *
- */
-yy_size_t gmsh_yyget_leng (void)
-{
- return gmsh_yyleng;
-}
-
-/** Get the current token.
- *
- */
-
-char *gmsh_yyget_text (void)
-{
- return gmsh_yytext;
-}
-
-/** Set the current line number.
- * @param line_number
- *
- */
-void gmsh_yyset_lineno (int line_number )
-{
-
- gmsh_yylineno = line_number;
-}
-
-/** Set the input stream. This does not discard the current
- * input buffer.
- * @param in_str A readable stream.
- *
- * @see gmsh_yy_switch_to_buffer
- */
-void gmsh_yyset_in (FILE * in_str )
-{
- gmsh_yyin = in_str ;
-}
-
-void gmsh_yyset_out (FILE * out_str )
-{
- gmsh_yyout = out_str ;
-}
-
-int gmsh_yyget_debug (void)
-{
- return gmsh_yy_flex_debug;
-}
-
-void gmsh_yyset_debug (int bdebug )
-{
- gmsh_yy_flex_debug = bdebug ;
-}
-
-static int yy_init_globals (void)
-{
- /* Initialization is the same as for the non-reentrant scanner.
- * This function is called from gmsh_yylex_destroy(), so don't allocate here.
- */
-
- (yy_buffer_stack) = 0;
- (yy_buffer_stack_top) = 0;
- (yy_buffer_stack_max) = 0;
- (yy_c_buf_p) = (char *) 0;
- (yy_init) = 0;
- (yy_start) = 0;
-
-/* Defined in main.c */
-#ifdef YY_STDINIT
- gmsh_yyin = stdin;
- gmsh_yyout = stdout;
-#else
- gmsh_yyin = (FILE *) 0;
- gmsh_yyout = (FILE *) 0;
-#endif
-
- /* For future reference: Set errno on error, since we are called by
- * gmsh_yylex_init()
- */
- return 0;
-}
-
-/* gmsh_yylex_destroy is for both reentrant and non-reentrant scanners. */
-int gmsh_yylex_destroy (void)
-{
-
- /* Pop the buffer stack, destroying each element. */
- while(YY_CURRENT_BUFFER){
- gmsh_yy_delete_buffer(YY_CURRENT_BUFFER );
- YY_CURRENT_BUFFER_LVALUE = NULL;
- gmsh_yypop_buffer_state();
- }
-
- /* Destroy the stack itself. */
- gmsh_yyfree((yy_buffer_stack) );
- (yy_buffer_stack) = NULL;
-
- /* Reset the globals. This is important in a non-reentrant scanner so the next time
- * gmsh_yylex() is called, initialization will occur. */
- yy_init_globals( );
-
- return 0;
-}
-
-/*
- * Internal utility routines.
- */
-
-#ifndef yytext_ptr
-static void yy_flex_strncpy (char* s1, yyconst char * s2, int n )
-{
- register int i;
- for ( i = 0; i < n; ++i )
- s1[i] = s2[i];
-}
-#endif
-
-#ifdef YY_NEED_STRLEN
-static int yy_flex_strlen (yyconst char * s )
-{
- register int n;
- for ( n = 0; s[n]; ++n )
- ;
-
- return n;
-}
-#endif
-
-void *gmsh_yyalloc (yy_size_t size )
-{
- return (void *) malloc( size );
-}
-
-void *gmsh_yyrealloc (void * ptr, yy_size_t size )
-{
- /* The cast to (char *) in the following accommodates both
- * implementations that use char* generic pointers, and those
- * that use void* generic pointers. It works with the latter
- * because both ANSI C and C++ allow castless assignment from
- * any pointer type to void*, and deal with argument conversions
- * as though doing an assignment.
- */
- return (void *) realloc( (char *) ptr, size );
-}
-
-void gmsh_yyfree (void * ptr )
-{
- free( (char *) ptr ); /* see gmsh_yyrealloc() for (char *) cast */
-}
-
-#define YYTABLES_NAME "yytables"
-
-#line 299 "Gmsh.l"
-
-
-
-#undef gmsh_yywrap
-
-int gmsh_yywrap() {return 1;}
-
-void skipcomments(void)
-{
- int c;
-
- while (1) {
- while ((c = yyinput()) != '*'){
- if(gmsheof(gmsh_yyin)){
- Msg::Error("End of file in commented region");
- return;
- }
- }
- if ((c = yyinput()) == '/')
- return;
- unput(c);
- }
-}
-
-void parsestring(char endchar)
-{
- int c;
- char tmp[1024];
-
- // Note that we keep special characters (end-of-line \n, tabs \t,
- // etc.) "as is" in the output string: see yyinput() above
- int i = 0;
- while ((c = yyinput()) != endchar) {
- if(gmsheof(gmsh_yyin)){
- Msg::Error("End of file in string");
- break;
- }
- else if(i >= (int)sizeof(tmp)-1){
- Msg::Error("String too long");
- break;
- }
- else{
- tmp[i++] = (char)c;
- }
- }
- tmp[i] = '\0';
- gmsh_yylval.c = strsave(tmp);
-}
-
-char *strsave(char *ptr)
-{
- return((char*)strcpy((char*)malloc(strlen(ptr)+1),ptr));
-}
-
-void skipline()
-{
- int c;
- while ((c = yyinput()) != '\n'){
- if(gmsheof(gmsh_yyin)) return;
- }
-}
-
-void skip_until(const char *skip, const char *until)
-{
- int i, nb_skip;
- int l, l_skip, l_until;
- char chars[256];
-
- nb_skip = 0;
-
- if(skip)
- l_skip = strlen(skip);
- else
- l_skip = 0;
-
- l_until = strlen(until);
-
- while(1){
- while (1){
- chars[0] = yyinput();
- if(gmsheof(gmsh_yyin)){
- Msg::Error("Unexpected end of file");
- return;
- }
- if(chars[0] == until[0]) break;
- if(skip && chars[0] == skip[0]) break;
- }
-
- l = std::max(l_skip,l_until);
- if(l >= (int)sizeof(chars)){
- Msg::Error("Search pattern too long in skip_until");
- return;
- }
- for(i=1; i<l; i++){
- chars[i] = yyinput();
- if(gmsheof(gmsh_yyin)){
- l = i;
- break;
- }
- }
-
- if(!strncmp(chars,until,l_until)){
- if(!nb_skip){
- return;
- }
- else{
- nb_skip--;
- }
- }
- else if(skip && !strncmp(chars,skip,l_skip)){
- nb_skip++;
- }
- else{
- for(i=1;i<l-1;i++){
- unput(chars[l-i]);
- }
- }
-
- }
-}
-
-void gmsh_yyflush() { YY_FLUSH_BUFFER; }
-