Skip to content
Snippets Groups Projects
Select Git revision
  • 1f58a04887353c62aee188496874e5046e0218f2
  • master default
  • cgnsUnstructured
  • partitioning
  • poppler
  • HighOrderBLCurving
  • gmsh_3_0_4
  • gmsh_3_0_3
  • gmsh_3_0_2
  • gmsh_3_0_1
  • gmsh_3_0_0
  • gmsh_2_16_0
  • gmsh_2_15_0
  • gmsh_2_14_1
  • gmsh_2_14_0
  • gmsh_2_13_2
  • gmsh_2_13_1
  • gmsh_2_12_0
  • gmsh_2_11_0
  • gmsh_2_10_1
  • gmsh_2_10_0
  • gmsh_2_9_3
  • gmsh_2_9_2
  • gmsh_2_9_1
  • gmsh_2_9_0
  • gmsh_2_8_6
26 results

gl2gif.cpp

Blame
  • Forked from gmsh / gmsh
    Source project has a limited visibility.
    GFace.cpp 21.71 KiB
    // Gmsh - Copyright (C) 1997-2008 C. Geuzaine, J.-F. Remacle
    //
    // See the LICENSE.txt file for license information. Please report all
    // bugs and problems to <gmsh@geuz.org>.
    
    #include "GModel.h"
    #include "GFace.h"
    #include "GEdge.h"
    #include "MElement.h"
    
    #if defined(HAVE_GMSH_EMBEDDED)
    #  include "GmshEmbedded.h"
    #else
    #  include "Message.h"
    #  include "Numeric.h"
    #  include "GaussLegendre1D.h"
    #  include "VertexArray.h"
    #  include "Context.h"
    #  if defined(HAVE_GSL)
    #    include <gsl/gsl_vector.h>
    #    include <gsl/gsl_linalg.h>
    #  else
    #    define NRANSI
    #    include "nrutil.h"
    void dsvdcmp(double **a, int m, int n, double w[], double **v);
    #  endif
    #endif
    
    extern Context_T CTX;
    
    #define SQU(a)      ((a)*(a))
    
    GFace::GFace(GModel *model, int tag)
      : GEntity(model, tag), r1(0), r2(0), va_geom_triangles(0)
    {
      meshStatistics.status = GFace::PENDING;
      resetMeshAttributes();
    }
    
    GFace::~GFace()
    {
      std::list<GEdge*>::iterator it = l_edges.begin();
    
      while (it != l_edges.end()){
        (*it)->delFace(this);
        ++it;
      }
    
      for(unsigned int i = 0; i < mesh_vertices.size(); i++)
        delete mesh_vertices[i];
    
      for(unsigned int i = 0; i < triangles.size(); i++)
        delete triangles[i];
    
      for(unsigned int i = 0; i < quadrangles.size(); i++)
        delete quadrangles[i];
    
      if(va_geom_triangles)
        delete va_geom_triangles;
    }
    
    unsigned int GFace::getNumMeshElements()
    { 
      return triangles.size() + quadrangles.size(); 
    }
    
    MElement *GFace::getMeshElement(unsigned int index)
    { 
      if(index < triangles.size())
        return triangles[index];