From ac9d14f30e085fbd56c75db35fd52366966da395 Mon Sep 17 00:00:00 2001
From: Jean-Francois Remacle <jean-francois.remacle@uclouvain.be>
Date: Fri, 27 Feb 2009 19:24:22 +0000
Subject: [PATCH] *** empty log message ***
---
Geo/GEdgeCompound.cpp | 131 ++++++++++++++++++++++++++++++++++++++++++
Geo/GEdgeCompound.h | 42 ++++++++++++++
2 files changed, 173 insertions(+)
create mode 100644 Geo/GEdgeCompound.cpp
create mode 100644 Geo/GEdgeCompound.h
diff --git a/Geo/GEdgeCompound.cpp b/Geo/GEdgeCompound.cpp
new file mode 100644
index 0000000000..aba627092c
--- /dev/null
+++ b/Geo/GEdgeCompound.cpp
@@ -0,0 +1,131 @@
+// Gmsh - Copyright (C) 1997-2009 C. Geuzaine, J.-F. Remacle
+//
+// See the LICENSE.txt file for license information. Please report all
+// bugs and problems to <gmsh@geuz.org>.
+
+#include "GmshConfig.h"
+#include "GEdgeCompound.h"
+#include "Numeric.h"
+
+GEdgeCompound::GEdgeCompound(GModel *m, int tag, std::vector<GEdge*> &compound)
+ : GEdge(m, tag, 0 , 0), _compound(compound)
+{
+ orderEdges ();
+ int N = _compound.size();
+ v0 = _orientation[0] ? _compound[0]->getBeginVertex() : _compound[0]->getEndVertex();
+ v1 = _orientation[N-1] ? _compound[N-1]->getEndVertex() : _compound[N-1]->getBeginVertex();
+ v0->addEdge(this);
+ v1->addEdge(this);
+ parametrize();
+}
+
+void GEdgeCompound::orderEdges()
+{
+ std::list<GEdge*> edges ;
+ std::vector<GEdge*> _c ;
+ for (int i=0;i<_compound.size();i++)edges.push_back(_compound[i]);
+ _c.push_back(*(edges.begin()));
+ edges.erase(edges.begin());
+ _orientation.push_back(true);
+ GVertex *first = _c[0]->getBeginVertex();
+ GVertex *last = _c[0]->getEndVertex();
+
+ while (first != last){
+ if (edges.empty())break;
+ for (std::list<GEdge*>::iterator it = edges.begin() ; it != edges.end() ; ++it){
+ GEdge *e = *it;
+ // printf("last %d edge %d %d\n",last->tag(),e->getBeginVertex()->tag(),
+ // e->getEndVertex()->tag());
+ if (e->getBeginVertex() == last){
+ _c.push_back(e);
+ edges.erase(it);
+ _orientation.push_back(true);
+ last = e->getEndVertex();
+ break;
+ }
+ else if (e->getEndVertex() == last){
+ _c.push_back(e);
+ edges.erase(it);
+ _orientation.push_back(false);
+ last = e->getBeginVertex();
+ break;
+ }
+ }
+ }
+ _compound = _c;
+}
+
+int GEdgeCompound::minimumMeshSegments() const{
+ int N = 0;
+ for (int i=0;i<_compound.size();i++)N +=_compound[i]->minimumMeshSegments();
+ return N;
+}
+int GEdgeCompound::minimumDrawSegments() const{
+ int N = 0;
+ for (int i=0;i<_compound.size();i++)N +=_compound[i]->minimumDrawSegments();
+ return N;
+}
+
+
+GEdgeCompound::~GEdgeCompound()
+{
+}
+
+Range<double> GEdgeCompound::parBounds(int i) const {
+ return Range<double>(0, _pars[_compound.size()-1]);
+}
+
+/*
+ +--------+-----------+----------- ... ----+
+ 0 _par[1] _par[2] _par[N-1]
+
+*/
+
+void GEdgeCompound::getLocalParameter ( const double &t,
+ int &iEdge,
+ double & tLoc) const
+{
+ for (iEdge=0 ; iEdge<_compound.size() ;iEdge++){
+ // printf("iEdge %d par %g\n",iEdge,_pars[iEdge]);
+ if (t <= _pars[iEdge]){
+ tLoc = _orientation[iEdge] ? t - _pars[iEdge] : _pars[iEdge+1] - t;
+ return;
+ }
+ }
+}
+
+void GEdgeCompound::parametrize()
+{
+ for (int i=0;i<_compound.size();i++){
+ Range<double> b = _compound[i]->parBounds(0);
+ if (!i)_pars.push_back(b.high() - b.low());
+ else _pars.push_back(_pars[i-1]+(b.high() - b.low()));
+ }
+}
+
+
+double GEdgeCompound::curvature(double par) const
+{
+ double tLoc;
+ int iEdge;
+ getLocalParameter(par,iEdge,tLoc);
+ return _compound[iEdge]->curvature(tLoc);
+}
+
+GPoint GEdgeCompound::point(double par) const
+{
+ double tLoc;
+ int iEdge;
+ getLocalParameter(par,iEdge,tLoc);
+ // printf("compound %d par %g edge %d loc %g\n",tag(),par,iEdge,tLoc);
+ return _compound[iEdge]->point(tLoc);
+}
+
+SVector3 GEdgeCompound::firstDer(double par) const
+{
+ double tLoc;
+ int iEdge;
+ getLocalParameter(par,iEdge,tLoc);
+ return _compound[iEdge]->firstDer(tLoc);
+}
+
diff --git a/Geo/GEdgeCompound.h b/Geo/GEdgeCompound.h
new file mode 100644
index 0000000000..8171e3ac83
--- /dev/null
+++ b/Geo/GEdgeCompound.h
@@ -0,0 +1,42 @@
+// Gmsh - Copyright (C) 1997-2009 C. Geuzaine, J.-F. Remacle
+//
+// See the LICENSE.txt file for license information. Please report all
+// bugs and problems to <gmsh@geuz.org>.
+
+#ifndef _GEDGE_COMPOUND_H_
+#define _GEDGE_COMPOUND_H_
+
+#include "GFace.h"
+#include "GEdge.h"
+
+/*
+A GEdgeCompound is a model edge that is the compound of model edges.
+*/
+
+class GEdgeCompound : public GEdge {
+ protected:
+ std::vector<GEdge*> _compound;
+ std::vector<bool> _orientation;
+ std::vector<double> _pars;
+ void parametrize() ;
+ void orderEdges() ;
+ void getLocalParameter ( const double &t,
+ int &iEdge,
+ double & tLoc) const;
+
+public:
+ GEdgeCompound(GModel *m, int tag, std::vector<GEdge*> &compound);
+ virtual ~GEdgeCompound();
+ Range<double> parBounds(int i) const;
+ virtual GPoint point(double par) const;
+ virtual SVector3 firstDer(double par) const;
+ virtual GEntity::GeomType geomType() const { return CompoundCurve; }
+ ModelType getNativeType() const { return GmshModel; }
+ void * getNativePtr() const { return 0; }
+ virtual double curvature(double t) const;
+ virtual int minimumMeshSegments() const;
+ virtual int minimumDrawSegments() const;
+
+};
+
+#endif
--
GitLab