.\"
.de Id
..
.de Sp
.if n .sp
.if t .sp 0.4
..
.TH prism 7rheolef "rheolef-6.7" "rheolef-6.7" "rheolef-6.7"
.\" label: /*Class:prism
.SH NAME
\fBprism\fP - Prism reference element
.\" skip: @cindex prism
.\" skip: @cindex reference element
.SH DESCRIPTION

The prism reference element is
.\" begin_example
.Sp
.nf
        K = { 0 < x < 1 and 0 < y < 1-x and -1 < z < 1 }
.Sp
.fi
.\" end_example
.SH NUMBERING

The orientation is such that triedra (01, 02, 03) is direct
and all faces, see from exterior, are in the direct sens.  References:
P. L. Georges, "Generation automatique de maillages",
page 24-, coll RMA, 16, Masson, 1994.
.\" begin_example
.Sp
.nf
            z
.PP
            3 ---- 5
         '  | .    |
      4 .  '|      |
      |     |      |
      |     |      |
      |     0 . . .2     y
      |  '    .
      1 .  '
.PP
    x
.Sp
.fi
.\" end_example
.\" END
.SH IMPLEMENTATION
.\" begin_example
.Sp
.nf
const size_t dimension = 3;
const Float  measure = 1;
const size_t n_vertex = 6;
const point vertex [n_vertex] = {
        point( 0, 0,-1 ),
        point( 1, 0,-1 ),
        point( 0, 1,-1 ),
        point( 0, 0, 1 ),
        point( 1, 0, 1 ),
        point( 0, 1, 1 ) };
const size_t  n_face = 5;
const size_t face [n_face][4] = {
        { 0, 2, 1, size_t(-1) },
        { 3, 4, 5, size_t(-1) },
        { 0, 1, 4, 3 },
        { 1, 2, 5, 4 },
        { 0, 3, 5, 2 } };
const size_t  n_edge = 9;
const size_t edge [n_edge][2] = {
        { 0, 1 },
        { 1, 2 },
        { 2, 0 },
        { 0, 3 },
        { 1, 4 },
        { 2, 5 },
        { 3, 4 },
        { 4, 5 },
        { 5, 3 } };
.Sp
.fi
.\" end_example