search

OptimalDesignLab / SummationByParts.jl

Other
6 stars 4 forks source link

Reference Element Topology #21

Open JaredCrean2 opened 6 years ago

JaredCrean2 commented 6 years ago

Not having a good way to get the reference element is a hindrance (although not blocking) supporting quadratic coordinate fields for FFD and MeshMovement2. What happens inside PDESolver is this:

      face_verts = SummationByParts.SymCubatures.getfacevertexindices(sbp.cub)
      edge_verts = [1 2 1 1 2 3;  # TODO: SBP should provide this
                             2 3 3 4 4 4]
      topo2 = ElementTopology2()  
      topo = ElementTopology{3}(face_verts, edge_verts, topo2=topo2)

But there should be some kind of API to get the edge_verts (or perhaps to supply the ElementTopology object itself, which is defined in ODLCommonTools?).