Grid data access

You may want to use the this package to obtain grid data for your application. The ExtendableGrids package provides a lightweight grid data structure which in its core is a Dict{Type,Any} and the possibility to dispatch the return type on the key.

function extract_2d()
    builder=SimplexGridBuilder(Generator=Triangulate)

    p1=point!(builder,0,0)
    p2=point!(builder,1,0)
    p3=point!(builder,1,1)
    p4=point!(builder,0,1)

    facetregion!(builder,1)
    facet!(builder,p1,p2)
    facetregion!(builder,2)
    facet!(builder,p2,p3)
    facetregion!(builder,3)
    facet!(builder,p3,p4)
    facetregion!(builder,4)
    facet!(builder,p4,p1)

    grid=simplexgrid(builder,maxvolume=0.25)

    @show grid[Coordinates]
    @show grid[CellNodes]
    @show grid[CellRegions]
    @show grid[BFaceNodes]
    @show grid[BFaceRegions]
    grid
end

The output of this call is this:

grid[Coordinates] = [0.0 1.0 1.0 0.0 0.5; 0.0 0.0 1.0 1.0 0.5]
grid[CellNodes] = Int32[2 4 5 1; 3 1 3 2; 5 5 4 5]
grid[CellRegions] = Int32[1, 1, 1, 1]
grid[BFaceNodes] = Int32[2 3 4 1; 1 2 3 4]
grid[BFaceRegions] = Int32[1, 2, 3, 4]

Thus, the grid is described by five arrays:

• 2 x npoints coordinates,
• 3 x ntriangles connectivity,
• ntriangles region markers,
• 2 x nfacets boundary faces,
• nfacets boundary face markers.

This page was generated using Literate.jl.