Skip to content
Snippets Groups Projects
Select Git revision
  • 4cc4b58d90556c85f48e602c7fe37506099be4eb
  • master default protected
  • dof-renumbering
  • test-dof-hash
  • gdemesy-master-patch-30528
  • eval-space-time
  • oscillating_multiharm
  • MH_movement
  • axisqu
  • write_vtu_and_ensight_formats
  • movingband
  • CP_1972_add_vtu_file_writing
  • mortar
  • fast_freq_sweep_Resolution
  • applyresolvent_again
  • marteaua-master-patch-54323
  • patch-1
  • binde-master-patch-08072
  • binde-master-patch-52461
  • BCGSL
  • resolvent
  • getdp_3_5_0
  • getdp_3_4_0
  • getdp_3_3_0
  • getdp_3_2_0
  • getdp_3_1_0
  • getdp_3_0_4
  • getdp_3_0_3
  • getdp_3_0_2
  • getdp_3_0_1
  • getdp_3_0_0
  • onelab_mobile_2.1.0
  • getdp_2_11_3 protected
  • getdp_2_11_2 protected
  • getdp_2_11_1 protected
  • getdp_2_11_0 protected
  • getdp_2_10_0 protected
  • getdp_2_9_2 protected
  • getdp_2_9_1 protected
  • getdp_2_9_0 protected
  • getdp_2_8_0 protected
41 results

CMakeLists.txt

Blame
  • gmshVertex.cpp 3.77 KiB
    // Gmsh - Copyright (C) 1997-2014 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 "GFace.h"
    #include "GEdge.h"
    #include "gmshVertex.h"
    #include "Geo.h"
    #include "GeoInterpolation.h"
    #include "GmshMessage.h"
    #include "MVertex.h"
    #include "MPoint.h"
    
    gmshVertex::gmshVertex(GModel *m, Vertex *_v)
      : GVertex(m, _v->Num, _v->lc), v(_v)
    {
    }
    
    void gmshVertex::resetMeshAttributes()
    {
      meshSize = v->lc;
    }
    
    void gmshVertex::setPosition(GPoint &p)
    {
      v->Pos.X = p.x();
      v->Pos.Y = p.y();
      v->Pos.Z = p.z();
      if(mesh_vertices.size()){
        mesh_vertices[0]->x() = p.x();
        mesh_vertices[0]->y() = p.y();
        mesh_vertices[0]->z() = p.z();
      }
    }
    
    GEntity::GeomType gmshVertex::geomType() const
    {
      if(v->Typ == MSH_POINT_BND_LAYER)
        return BoundaryLayerPoint;
      else
        return Point;
    }
    
    SPoint2 gmshVertex::reparamOnFace(const GFace *face, int dir) const
    {
    
      Surface *s = (Surface*)face->getNativePtr();
    
      if(s->geometry){
        // It is not always right if it is periodic.
        if(l_edges.size() == 1 &&
           (*l_edges.begin())->getBeginVertex() ==
           (*l_edges.begin())->getEndVertex()){
          Range<double> bb = (*l_edges.begin())->parBounds(0);
          return (*l_edges.begin())->reparamOnFace(face, bb.low(), dir);
        }
        return v->pntOnGeometry;
      }
    
      if(s->Typ ==  MSH_SURF_REGL){
        Curve *C[4];
        for(int i = 0; i < 4; i++)
          List_Read(s->Generatrices, i, &C[i]);
    
        double U, V;
        if ((C[0]->beg == v && C[3]->beg == v) ||
            (C[0]->end == v && C[3]->beg == v) ||
            (C[0]->beg == v && C[3]->end == v) ||
            (C[0]->end == v && C[3]->end == v)){
          U = V = 0;
        }
        else if ((C[0]->beg == v && C[1]->beg == v) ||
                 (C[0]->end == v && C[1]->beg == v) ||
                 (C[0]->beg == v && C[1]->end == v) ||
                 (C[0]->end == v && C[1]->end == v)){
          U = 1;
          V = 0;
        }
        else if ((C[2]->beg == v && C[1]->beg == v) ||
                 (C[2]->end == v && C[1]->beg == v) ||
                 (C[2]->beg == v && C[1]->end == v) ||
                 (C[2]->end == v && C[1]->end == v)){
          U = 1;
          V = 1;
        }
        else if ((C[2]->beg == v && C[3]->beg == v) ||
                 (C[2]->end == v && C[3]->beg == v) ||
                 (C[2]->beg == v && C[3]->end == v) ||
                 (C[2]->end == v && C[3]->end == v)){
          U = 0;
          V = 1;
        }
        else{
          Msg::Info("Reparameterizing point %d on face %d", v->Num, s->Num);
          return GVertex::reparamOnFace(face, dir);
        }
        return SPoint2(U, V);
      }
      else if(s->Typ ==  MSH_SURF_TRIC){
        Curve *C[3];
        for(int i = 0; i < 3; i++)
          List_Read(s->Generatrices, i, &C[i]);
    
        double U, V;
        if ((C[0]->beg == v && C[2]->beg == v) ||
            (C[0]->end == v && C[2]->beg == v) ||
            (C[0]->beg == v && C[2]->end == v) ||
            (C[0]->end == v && C[2]->end == v)){
          U = V = 0;
        }
        else if ((C[0]->beg == v && C[1]->beg == v) ||
                 (C[0]->end == v && C[1]->beg == v) ||
                 (C[0]->beg == v && C[1]->end == v) ||
                 (C[0]->end == v && C[1]->end == v)){
          U = 1;
          V = 0;
        }
        else if ((C[2]->beg == v && C[1]->beg == v) ||
                 (C[2]->end == v && C[1]->beg == v) ||
                 (C[2]->beg == v && C[1]->end == v) ||
                 (C[2]->end == v && C[1]->end == v)){
          U = 1;
          V = 1;
        }
        else{
          Msg::Info("Reparameterizing point %d on face %d", v->Num, s->Num);
          return GVertex::reparamOnFace(face, dir);
        }
        return SPoint2(U, V);
      }
      else{
        return GVertex::reparamOnFace(face, dir);
      }
    }
    
    void gmshVertex::writeGEO(FILE *fp, const std::string &meshSizeParameter)
    {
      fprintf(fp, "Point(%d) = {%.16g, %.16g, %.16g, %.16g};\n",
              v->Num, v->Pos.X, v->Pos.Y, v->Pos.Z, v->lc);
    }