Select Git revision
PixelBuffer.h
discreteEdge.cpp 18.34 KiB
// Gmsh - Copyright (C) 1997-2017 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@onelab.info>.
#include <vector>
#include <list>
#include "GmshConfig.h"
#include "GmshMessage.h"
#include "discreteEdge.h"
#include "MLine.h"
#include "Numeric.h"
#include "MPoint.h"
#include "MTriangle.h"
#include "MQuadrangle.h"
#include "MPrism.h"
#include "MTetrahedron.h"
#include "MHexahedron.h"
#include "MPyramid.h"
#include "GModelIO_GEO.h"
#include "Geo.h"
#include "OS.h"
#if defined(HAVE_MESH)
#include "meshGEdge.h"
#include "Context.h"
#endif
discreteEdge::discreteEdge(GModel *model, int num, GVertex *_v0, GVertex *_v1)
: GEdge(model, num, _v0, _v1)
{
createdTopo = false;
Curve *c = CreateCurve(num, MSH_SEGM_DISCRETE, 0, 0, 0, -1, -1, 0., 1.);
Tree_Add(model->getGEOInternals()->Curves, &c);
CreateReversedCurve(c);
}
// topology is already set
void discreteEdge::setTopo(std::vector<MLine*> mlines)
{
createdTopo = true;
lines = mlines;
_orientation = std::vector<int>(lines.size(),1);
}
void discreteEdge::createTopo()
{
if(!createdTopo){
orderMLines();
setBoundVertices();
createdTopo = true;
}
}
// FULL OF BUGS !!!!!!
void discreteEdge::orderMLines()
{
//printf("ordering line %d (%d,%d)\n", tag(),getBeginVertex()->mesh_vertices[0]->getNum(),getEndVertex()->mesh_vertices[0]->getNum());
//if(lines.size() <= 1) return;
std::vector<MLine*> _m;
std::list<MLine*> segments;
// store all lines in a list : segments
for (unsigned int i = 0; i < lines.size(); i++) segments.push_back(lines[i]);
// find a lonly MLine
for (std::list<MLine*>::iterator it = segments.begin();
it != segments.end(); ++it){
MVertex *vL = (*it)->getVertex(0);