Skip to content
Snippets Groups Projects
Select Git revision
  • 82ef2ab118fb009c015bd2404bcb003f25d6d6d1
  • master default protected
  • relaying
  • overlaps_tags_and_distributed_export
  • patches-4.14
  • steplayer
  • bl
  • pluginMeshQuality
  • fixBugsAmaury
  • hierarchical-basis
  • alphashapes
  • new_export_boris
  • oras_vs_osm
  • reassign_partitions
  • distributed_fwi
  • rename-classes
  • fix/fortran-api-example-t4
  • robust_partitions
  • reducing_files
  • fix_overlaps
  • 3115-issue-fix
  • 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

explore.jl

Blame
  • explore.jl 1.22 KiB
    import gmsh
    
    if length(ARGS) < 1
        println("Usage: julia explore.jl file.msh")
        exit(0)
    end
    
    gmsh.initialize()
    gmsh.option.setNumber("General.Terminal", 1)
    gmsh.open(ARGS[1])
    
    # get all elementary entities in the model
    entities = gmsh.model.getEntities()
    
    for e in entities
        # get the mesh nodes for each elementary entity
        nodeTags, nodeCoords, nodeParams = gmsh.model.mesh.getNodes(e[1], e[2])
        # get the mesh elements for each elementary entity
        elemTypes, elemTags, elemNodeTags = gmsh.model.mesh.getElements(e[1], e[2])
        # report some statistics
        numElem = sum([length(t) for t in elemTags])
        println(length(nodeTags), " mesh nodes and ", numElem, " mesh elements ",
                "on entity ", e, " of type ", gmsh.model.getType(e[1], e[2]))
        partitions = gmsh.model.getPartitions(e[1], e[2])
        if length(partitions) > 0
            println(" - Partition tag(s): ", partitions, " - parent entity",
                    gmsh.model.getParent(e[1], e[2]))
        end
        for t in elemTypes
            name, dim, order, numv, parv = gmsh.model.mesh.getElementProperties(t)
            println(" - Element type: ", name, ", order ", order)
            println("   with ", numv, " nodes in param coord: ", parv)
        end
    end
    
    gmsh.finalize()