tritet 3.1.0

Triangle and tetrahedron mesh generators
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163

use crate::StrError;

pub(crate) const TRITET_SUCCESS: i32 = 0;

pub(crate) const TRITET_ERROR_NULL_DATA: i32 = 10;
pub(crate) const TRITET_ERROR_STRING_CONCAT: i32 = 20;

pub(crate) const TRITET_ERROR_NULL_POINT_LIST: i32 = 100;
pub(crate) const TRITET_ERROR_NULL_SEGMENT_LIST: i32 = 200;
pub(crate) const TRITET_ERROR_NULL_FACET_LIST: i32 = 300;
pub(crate) const TRITET_ERROR_NULL_FACET_POLYGON_LIST: i32 = 400;
pub(crate) const TRITET_ERROR_NULL_REGION_LIST: i32 = 500;
pub(crate) const TRITET_ERROR_NULL_HOLE_LIST: i32 = 600;

pub(crate) const TRITET_ERROR_INVALID_POINT_INDEX: i32 = 1000;
pub(crate) const TRITET_ERROR_INVALID_SEGMENT_INDEX: i32 = 2000;
pub(crate) const TRITET_ERROR_INVALID_SEGMENT_POINT_ID: i32 = 3000;
pub(crate) const TRITET_ERROR_INVALID_FACET_INDEX: i32 = 4000;
pub(crate) const TRITET_ERROR_INVALID_FACET_NUM_POLYGON: i32 = 5000;
pub(crate) const TRITET_ERROR_INVALID_FACET_POINT_INDEX: i32 = 6000;
pub(crate) const TRITET_ERROR_INVALID_FACET_POINT_ID: i32 = 7000;
pub(crate) const TRITET_ERROR_INVALID_REGION_INDEX: i32 = 8000;
pub(crate) const TRITET_ERROR_INVALID_HOLE_INDEX: i32 = 9000;

pub(crate) const TRITET_ERROR_TETGEN_FAIL: i32 = 10000;

pub(crate) fn handle_status(status: i32) -> Result<(), StrError> {
    if status == TRITET_SUCCESS {
        return Ok(());
    }

    if status == TRITET_ERROR_NULL_DATA {
        return Err("INTERNAL ERROR: found NULL data");
    }
    if status == TRITET_ERROR_STRING_CONCAT {
        return Err("INTERNAL ERROR: cannot write string with commands for TetGen");
    }

    if status == TRITET_ERROR_NULL_POINT_LIST {
        return Err("INTERNAL ERROR: found NULL point list");
    }
    if status == TRITET_ERROR_NULL_SEGMENT_LIST {
        return Err("INTERNAL ERROR: found NULL segment list");
    }
    if status == TRITET_ERROR_NULL_FACET_LIST {
        return Err("INTERNAL ERROR: found NULL facet list");
    }
    if status == TRITET_ERROR_NULL_FACET_POLYGON_LIST {
        return Err("INTERNAL ERROR: found NULL facet polygon list");
    }
    if status == TRITET_ERROR_NULL_REGION_LIST {
        return Err("INTERNAL ERROR: found NULL region list");
    }
    if status == TRITET_ERROR_NULL_HOLE_LIST {
        return Err("INTERNAL ERROR: found NULL hole list");
    }

    if status == TRITET_ERROR_INVALID_POINT_INDEX {
        return Err("index of point is out of bounds");
    }
    if status == TRITET_ERROR_INVALID_SEGMENT_INDEX {
        return Err("index of segment is out of bounds");
    }
    if status == TRITET_ERROR_INVALID_SEGMENT_POINT_ID {
        return Err("id of segment point is out of bounds");
    }
    if status == TRITET_ERROR_INVALID_FACET_INDEX {
        return Err("index of facet is out of bounds");
    }
    if status == TRITET_ERROR_INVALID_FACET_NUM_POLYGON {
        return Err("INTERNAL ERROR: found invalid facet number of polygon");
    }
    if status == TRITET_ERROR_INVALID_FACET_POINT_INDEX {
        return Err("index of facet point is out of bounds");
    }
    if status == TRITET_ERROR_INVALID_FACET_POINT_ID {
        return Err("id of facet point is out of bounds");
    }
    if status == TRITET_ERROR_INVALID_REGION_INDEX {
        return Err("index of region is out of bounds");
    }
    if status == TRITET_ERROR_INVALID_HOLE_INDEX {
        return Err("index of hole is out of bounds");
    }

    if status == TRITET_ERROR_TETGEN_FAIL {
        return Err("TetGen failed: points are probably coplanar");
    }

    return Err("INTERNAL ERROR: some error occurred");
}

/// Maps indices used in this library (tritet) to indices used in Triangle
///
/// ```text
/// This library (tritet)      Triangle
///         NODES               CORNERS
///           2                    2
///          / \                  / \
///         /   \                /   \
///        5     4              4     3
///       /       \            /       \
///      /         \          /         \
///     0-----3-----1        0-----5-----1
/// ```
pub(crate) const TRITET_TO_TRIANGLE: [usize; 6] = [0, 1, 2, 5, 3, 4];

/// Maps indices used in this library (tritet) to indices used in Tetgen
///
/// ```text
///       This library (tritet)                          Tetgen
///               NODES                                  CORNERS
///             |                                         |
///             3                                         3
///            /|`.                                      /|`.
///            ||  `,                                    ||  `,
///           / |    ',                                 / |    ',
///           | |      \                                | |      \
///          /  |       `.                             /  |       `.
///          |  |         `,                           |  |         `,
///         /   7            9                        /   5            4
///         |   |             \                       |   |             \
///        /    |              `.                    /    |              `.
///        |    |                ',                  |    |                ',
///       8     |                  \                8     |                  \
///       |     0 ,,_               `.              |     0 ,,_               `.
///      |     /     ``'-., 6         `.           |     /     ``'-., 9         `.
///      |    /               `''-.,,_  ',         |    /               `''-.,,_  ',
///     |    /                        ``'2 ,,     |    /                        ``'2 ,,
///     |   '                       ,.-``         |   '                       ,.-``
///    |   4                   _,-'`             |   6                   _,-'`
///    ' /                 ,.'`                  ' /                 ,.'`
///   | /             _ 5 `                     | /             _ 7 `
///   '/          ,-'`                          '/          ,-'`
///  |/      ,.-``                             |/      ,.-``
///  /  _,-``                                  /  _,-``
/// 1 '`                                      1 '`
/// ```
pub(crate) const TRITET_TO_TETGEN: [usize; 10] = [0, 1, 2, 3, 6, 7, 9, 5, 8, 4];

/// Defines a set of "light" colors
pub(crate) const LIGHT_COLORS: [&'static str; 17] = [
    "#cbe4f9", "#cdf5f6", "#eff9da", "#f9ebdf", "#f9d8d6", "#d6cdea", "#acddde", "#caf1de", "#e1f8dc", "#fef8dd",
    "#ffe7c7", "#f7d8ba", "#d0fffe", "#fffddb", "#e4ffde", "#ffd3fd", "#ffe7d3",
];

/// Defines a set of "dark" colors
pub(crate) const DARK_COLORS: [&'static str; 12] = [
    "#2e3d7c", "#282528", "#ba292e", "#e15d3a", "#ffa73c", "#780000", "#540d4d", "#214b14", "#a36032", "#0f4539",
    "#2f3b22", "#152d32",
];

/// Sets the VTK code corresponding to triangles
pub(crate) const VTK_TRIANGLE: i32 = 5;

/// Sets the VTK code corresponding to quadratic triangles
pub(crate) const VTK_QUADRATIC_TRIANGLE: i32 = 22;

/// Sets the VTK code corresponding to tetrahedra
pub(crate) const VTK_TETRA: i32 = 10;

/// Sets the VTK code corresponding to quadratic tetrahedra
pub(crate) const VTK_QUADRATIC_TETRA: i32 = 24;