/usr/include/gerris/isocube.h is in libgfs-mpi-dev 20131206+dfsg-5.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
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 | /* Gerris - The GNU Flow Solver
* Copyright (C) 2001-2004 National Institute of Water and Atmospheric
* Research
*
* This program is free software; you can redistribute it and/or
* modify it under the terms of the GNU General Public License as
* published by the Free Software Foundation; either version 2 of the
* License, or (at your option) any later version.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this program; if not, write to the Free Software
* Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA
* 02111-1307, USA.
*/
/* isocube adapted from GTS (see gts/src/iso.c and gts/doc/isocube.fig) */
static guint edge1[12][2] = {
{0, 4}, {1, 5}, {3, 7}, {2, 6},
{0, 2}, {1, 3}, {5, 7}, {4, 6},
{0, 1}, {4, 5}, {6, 7}, {2, 3}
};
static FttVector vertex[8] = {
{0.,0.,0.},{0.,0.,1.},{0.,1.,0.},{0.,1.,1.},
{1.,0.,0.},{1.,0.,1.},{1.,1.,0.},{1.,1.,1.}
};
static guint face[6][4][2] = {
{{7,0},{10,0},{6,1},{9,1}}, /* right */
{{4,0},{11,0},{5,1},{8,1}}, /* left */
{{3,0},{10,0},{2,1},{11,1}},/* top */
{{0,0},{9,0},{1,1},{8,1}}, /* bottom */
{{1,0},{6,0},{2,1},{5,1}}, /* front */
{{0,0},{7,0},{3,1},{4,1}} /* back */
};
static guint face_v[6][4] = {
{4,6,7,5},/* right */
{0,2,3,1},/* left */
{2,6,7,3},/* top */
{0,4,5,1},/* bottom */
{1,5,7,3},/* front */
{0,4,6,2} /* back */
};
static FttVector edge[12][2] = {
{{0.,0.,0.},{1.,0.,0.}},{{0.,0.,1.},{1.,0.,1.}},{{0.,1.,1.},{1.,1.,1.}},{{0.,1.,0.},{1.,1.,0.}},
{{0.,0.,0.},{0.,1.,0.}},{{0.,0.,1.},{0.,1.,1.}},{{1.,0.,1.},{1.,1.,1.}},{{1.,0.,0.},{1.,1.,0.}},
{{0.,0.,0.},{0.,0.,1.}},{{1.,0.,0.},{1.,0.,1.}},{{1.,1.,0.},{1.,1.,1.}},{{0.,1.,0.},{0.,1.,1.}}
};
static guint connectv[12][2][4] = {
{{4, 5, 1, 0}, {0, 2, 6, 4}}, /* 0 */
{{5, 7, 3, 1}, {1, 0, 4, 5}}, /* 1 */
{{7, 6, 2, 3}, {3, 1, 5, 7}}, /* 2 */
{{6, 4, 0, 2}, {2, 3, 7, 6}}, /* 3 */
{{2, 6, 4, 0}, {0, 1, 3, 2}}, /* 4 */
{{3, 2, 0, 1}, {1, 5, 7, 3}}, /* 5 */
{{7, 3, 1, 5}, {5, 4, 6, 7}}, /* 6 */
{{6, 7, 5, 4}, {4, 0, 2, 6}}, /* 7 */
{{1, 3, 2, 0}, {0, 4, 5, 1}}, /* 8 */
{{5, 1, 0, 4}, {4, 6, 7, 5}}, /* 9 */
{{7, 5, 4, 6}, {6, 2, 3, 7}}, /* 10 */
{{3, 7, 6, 2}, {2, 0, 1, 3}} /* 11 */
};
static FttVector cvertex[8] = {
{0., 0., 0.}, {0., 0., 1.}, {0., 1., 0.}, {0., 1., 1.},
{1., 0., 0.}, {1., 0., 1.}, {1., 1., 0.}, {1., 1., 1.}
};
#if FTT_2D
# define NCORNERS 4
static FttDirection corner[4][2] = {
{FTT_LEFT, FTT_BOTTOM}, {FTT_RIGHT, FTT_BOTTOM}, {FTT_RIGHT, FTT_TOP}, {FTT_LEFT, FTT_TOP}
};
#else /* 3D */
# define NCORNERS 8
/* first index is the edge number, second index is the edge orientation
(0 or 1), third index are the edges which this edge may connect to
in order and the corresponding face direction */
static guint connect[12][2][4] = {
{{9, 1, 8, FTT_BOTTOM}, {4, 3, 7, FTT_BACK}}, /* 0 */
{{6, 2, 5, FTT_FRONT}, {8, 0, 9, FTT_BOTTOM}}, /* 1 */
{{10, 3, 11, FTT_TOP}, {5, 1, 6, FTT_FRONT}}, /* 2 */
{{7, 0, 4, FTT_BACK}, {11, 2, 10, FTT_TOP}}, /* 3 */
{{3, 7, 0, FTT_BACK}, {8, 5, 11, FTT_LEFT}}, /* 4 */
{{11, 4, 8, FTT_LEFT}, {1, 6, 2, FTT_FRONT}}, /* 5 */
{{2, 5, 1, FTT_FRONT}, {9, 7, 10, FTT_RIGHT}}, /* 6 */
{{10, 6, 9, FTT_RIGHT}, {0, 4, 3, FTT_BACK}}, /* 7 */
{{5, 11, 4, FTT_LEFT}, {0, 9, 1, FTT_BOTTOM}}, /* 8 */
{{1, 8, 0, FTT_BOTTOM}, {7, 10, 6, FTT_RIGHT}}, /* 9 */
{{6, 9, 7, FTT_RIGHT}, {3, 11, 2, FTT_TOP}}, /* 10 */
{{2, 10, 3, FTT_TOP}, {4, 8, 5, FTT_LEFT}} /* 11 */
};
static FttDirection corner[8][3] = {
{FTT_LEFT, FTT_BOTTOM, FTT_BACK},
{FTT_LEFT, FTT_BOTTOM, FTT_FRONT},
{FTT_LEFT, FTT_TOP, FTT_BACK},
{FTT_LEFT, FTT_TOP, FTT_FRONT},
{FTT_RIGHT, FTT_BOTTOM, FTT_BACK},
{FTT_RIGHT, FTT_BOTTOM, FTT_FRONT},
{FTT_RIGHT, FTT_TOP, FTT_BACK},
{FTT_RIGHT, FTT_TOP, FTT_FRONT}
};
#endif /* 3D */
#define SLIGHTLY_LARGER 1.001
|