Polytope for wedge about gridap.jl

[fverdugo]

Perhaps this is not very important right now, but it seems that Polytope is not giving the right information on nfaces for wedges.

wedge = Polytope(2,2,1)
@show wedge.nf_nfs
# wedge.nf_nfs = Array{Int64,1}[[1], [2], [3], [4], [5], [6], [1, 4, 7], [2, 5, 8], [3, 6, 9], [1, 2, 10], [2, 4, 11], [1, 3, 12], [2, 3, 13], [3, 4, 14], [3, 5, 15], [1, 2, 4, 5, 7, 8, 10, 11, 16], [1, 3, 4, 6, 7, 9, 12, 14, 17], [2, 3, 5, 6, 8, 9, 13, 15, 18], [1, 2, 3, 10, 12, 13, 19], [2, 3, 4, 11, 13, 14, 20], [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21]]

I would say that info about nface number 11 is wrong. Is this a bug or I am missing something @santiagobadia ?

It can also be clearly seen by writing the wedge into vtu:

writevtk(wedge,"wedge")

For quads, tris, hexs, tets, and pyramids works OK. E.g. for pyramids:

pyramid = Polytope(1,1,2)
writevtk(pyramid,"pyramid")

[santiagobadia]

I will take a look at this.

[santiagobadia]

Fixed

[fverdugo]

Perfect! Nice pictures! 😎

Closing issue