Skip to content
Snippets Groups Projects
Select Git revision
  • 9c646474c3c78a76957518976dc1b393bee3228d
  • master default protected
  • revert-ef4a3a4f
  • patch_releases_4_14
  • overlaps_tags_and_distributed_export
  • overlaps_tags_and_distributed_export_rebased
  • relaying
  • alphashapes
  • steplayer
  • bl
  • pluginMeshQuality
  • fixBugsAmaury
  • hierarchical-basis
  • new_export_boris
  • oras_vs_osm
  • reassign_partitions
  • distributed_fwi
  • rename-classes
  • fix/fortran-api-example-t4
  • robust_partitions
  • reducing_files
  • gmsh_4_14_0
  • gmsh_4_13_1
  • gmsh_4_13_0
  • gmsh_4_12_2
  • gmsh_4_12_1
  • gmsh_4_12_0
  • gmsh_4_11_1
  • gmsh_4_11_0
  • gmsh_4_10_5
  • gmsh_4_10_4
  • gmsh_4_10_3
  • gmsh_4_10_2
  • gmsh_4_10_1
  • gmsh_4_10_0
  • gmsh_4_9_5
  • gmsh_4_9_4
  • gmsh_4_9_3
  • gmsh_4_9_2
  • gmsh_4_9_1
  • gmsh_4_9_0
41 results

CutPlane.cpp

Blame
  • meshGRegionDelaunayInsertion.cpp 50.83 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 <set>
    #include <map>
    #include <algorithm>
    #include "GmshMessage.h"
    #include "robustPredicates.h"
    #include "OS.h"
    #include "meshGRegion.h"
    #include "meshGRegionLocalMeshMod.h"
    #include "meshGRegionDelaunayInsertion.h"
    #include "GModel.h"
    #include "GRegion.h"
    #include "MTriangle.h"
    #include "Numeric.h"
    #include "Context.h"
    #include "delaunay3d.h"
    #include "MEdge.h"
    #include "MLine.h"
    
    int MTet4::radiusNorm = 2;
    static double LIMIT_ = 1;
    
    static void createAllEmbeddedEdges (GRegion *gr, std::set<MEdge, Less_Edge> &allEmbeddedEdges)
    {
      std::list<GEdge*> e = gr->embeddedEdges();
      //  printf("=================> %d embedded GEdges\n",e.size());
      for (std::list<GEdge*>::iterator it = e.begin() ; it != e.end(); ++it){
        for (unsigned int i = 0; i < (*it)->lines.size(); i++){
          allEmbeddedEdges.insert (MEdge((*it)->lines[i]->getVertex(0),(*it)->lines[i]->getVertex(1)));
        }
      }
    }
    
    static void createAllEmbeddedFaces (GRegion *gr, std::set<MFace, Less_Face> &allEmbeddedFaces)
    {
      std::list<GFace*> f = gr->embeddedFaces();
      for (std::list<GFace*>::iterator it = f.begin() ; it != f.end(); ++it){
        for (unsigned int i = 0; i < (*it)->triangles.size(); i++){
          allEmbeddedFaces.insert ((*it)->triangles[i]->getFace(0));
        }
      }
    }
    
    static bool isActive(MTet4 *t, double limit_, int &active)
    {
      if (t->isDeleted()) return false;
      for (active = 0; active < 4; active++){
        MTet4 *neigh = t->getNeigh(active);
        if (!neigh || (neigh->getRadius() < limit_ && neigh->getRadius() > 0)) {
          return true;
        }
      }
      return false;
    }
    
    
    int MTet4::inCircumSphere(const double *p) const
    {
      double pa[3] = {base->getVertex(0)->x(),
                      base->getVertex(0)->y(),
                      base->getVertex(0)->z()};
      double pb[3] = {base->getVertex(1)->x(),
                      base->getVertex(1)->y(),
                      base->getVertex(1)->z()};
      double pc[3] = {base->getVertex(2)->x(),
                      base->getVertex(2)->y(),