Visualization of meshes with quads and general polytopes

[paurierap]

[fverdugo]

We started to investigate the extension to quads and hex with @Kevin-Mattheus-Moerman

preliminar code in branch: https://github.com/gridap/GridapMakie.jl/tree/quads_and_hexs

For the record, I paste here the test we used:

using Gridap
using Gridap.Geometry
using GridapMakie
using GLMakie
using Gridap.Visualization

domain = (0,1,0,1,0,1)
cell_nums = 10 .* (1,1,1)
model = CartesianDiscreteModel(domain,cell_nums)
Ω = Triangulation(model)
f = CellField(x->sin(π*(x[1]+x[2])),Ω)
fig = plot(Ω,f,strokewidth=2,strokecolor=:black,shading=true,transparency=true)
display(fig)

Pending issues:

• [ ] Check that we have not destroyed the code for tris/tets
• [ ] Problems with shadings
• [ ] Simplification of function to_face_grid_with_map
• [ ] GeometryBasics.normal_mesh seems to convert to tris/tets

[Kevin-Mattheus-Moerman]

@fverdugo @ericneiva @paurierap seeing this extremely late but I've done quadrilateral visualisations here: https://github.com/Kevin-Mattheus-Moerman/JuliaAdventures/blob/main/gridap_makie_visualisation_01.jl. It requires Makie's flexible poly function.