Skip to content
Snippets Groups Projects
Select Git revision
  • fcecb101c342200368734ba25499b15bffcf4b34
  • 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

F_Python.cpp

Blame
  • test.py 1.61 KiB
    import gmsh
    import sys
    
    gmsh.initialize(sys.argv)
    print(gmsh.option.getNumber("Mesh.Algorithm"))
    
    gmsh.open("square.msh")
    
    gmsh.model.add("square")
    gmsh.model.geo.addPoint(0, 0, 0, 0.1, 1)
    gmsh.model.geo.addPoint(1, 0, 0, 0.1, 2)
    gmsh.model.geo.addPoint(1, 1, 0, 0.1, 3)
    gmsh.model.geo.addPoint(0, 1, 0, 0.1, 4)
    gmsh.model.geo.addLine(1, 2, 1)
    gmsh.model.geo.addLine(2, 3, 2)
    gmsh.model.geo.addLine(3, 4, 3)
    line4 = gmsh.model.geo.addLine(4, 1)
    print("line4 received tag ", line4)
    gmsh.model.geo.addCurveLoop([1, 2, 3, line4], 1)
    gmsh.model.geo.addPlaneSurface([1], 6)
    gmsh.model.geo.synchronize()
    
    ptag = gmsh.model.addPhysicalGroup(1, [1, 2, 3, 4])
    ent = gmsh.model.getEntitiesForPhysicalGroup(1, ptag)
    print("new physical group ", ptag, ":", ent, type(ent))
    
    gmsh.model.addPhysicalGroup(2, [6])
    
    print(gmsh.option.getString("General.BuildOptions"))
    print(gmsh.option.getNumber("Mesh.Algorithm"))
    gmsh.option.setNumber("Mesh.Algorithm", 3.0)
    print(gmsh.option.getNumber("Mesh.Algorithm"))
    gmsh.model.mesh.generate(2)
    
    gmsh.write("square.msh")
    
    print("Entities")
    entities = gmsh.model.getEntities()
    for e in entities:
        print("entity ", e)
        types, tags, nodes = gmsh.model.mesh.getElements(e[0], e[1])
        for i in range(len(types)):
            print("type ", types[i])
            print("tags : ", list(tags[i]))
            print("nodes : ", list(nodes[i]))
        if e[0] == [2] and e[1] == 6:
            gmsh.model.mesh.addElements(e[0], e[1], types, [tags[0][:10]],
                                        [nodes[0][:30]])
    
    gmsh.write("mesh_truncated.msh")
    print("Nodes")
    tags, coord, _ = gmsh.model.mesh.getNodes(2, 6)
    print(tags)
    print(coord)
    gmsh.finalize()