This file is indexed.

/usr/include/gerris/fluid.h is in libgfs-mpi-dev 20131206+dfsg-7.

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
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
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
/* Gerris - The GNU Flow Solver
 * Copyright (C) 2001 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.  
 */

#ifndef __FLUID_H__
#define __FLUID_H__

#include <gts.h>

#ifdef __cplusplus
extern "C" {
#endif /* __cplusplus */

#include "ftt.h"

typedef struct _GfsVariable               GfsVariable;
typedef struct _GfsDerivedVariable        GfsDerivedVariable;
typedef struct _GfsDomain                 GfsDomain;

typedef struct _GfsStateVector     GfsStateVector;
typedef struct _GfsSolidVector     GfsSolidVector;
typedef struct _GfsFaceStateVector GfsFaceStateVector;

struct _GfsFaceStateVector {
  gdouble un;
  gdouble v;
};

struct _GfsStateVector {
  /* temporary face variables */
  GfsFaceStateVector f[FTT_NEIGHBORS];

  /* solid boundaries */
  GfsSolidVector * solid;

  gdouble place_holder;
};

struct _GfsSolidVector {
  gdouble s[FTT_NEIGHBORS];
  gdouble a, fv;
  FttCell * merged;
  FttVector cm, ca, v;
};

typedef enum {
  GFS_FLAG_USED =              1 <<  FTT_FLAG_USER,
  GFS_FLAG_BOUNDARY          = 1 << (FTT_FLAG_USER + 1),
  GFS_FLAG_DIRICHLET         = 1 << (FTT_FLAG_USER + 2),
  GFS_FLAG_GRADIENT_BOUNDARY = 1 << (FTT_FLAG_USER + 3),
  GFS_FLAG_PERMANENT         = 1 << (FTT_FLAG_USER + 4),
  GFS_FLAG_THIN              = 1 << (FTT_FLAG_USER + 5),
  GFS_FLAG_USER =                    FTT_FLAG_USER + 6 /* user flags start here */
} GfsFlags;

#define GFS_STATE(cell)               ((GfsStateVector *) (cell)->data)
#define GFS_VALUEI(cell, index)     ((&GFS_STATE (cell)->place_holder)[index])

#define GFS_FACE_NORMAL_VELOCITY(fa)\
  (GFS_STATE ((fa)->cell)->f[(fa)->d].un)
#define GFS_FACE_NORMAL_VELOCITY_LEFT(fa)\
  (GFS_STATE ((fa)->cell)->f[(fa)->d].un)
#define GFS_FACE_NORMAL_VELOCITY_RIGHT(fa)\
  (GFS_STATE ((fa)->neighbor)->f[FTT_OPPOSITE_DIRECTION ((fa)->d)].un)

#define GFS_FACE_FRACTION(fa) (GFS_IS_MIXED ((fa)->cell) ?\
                               GFS_STATE ((fa)->cell)->solid->s[(fa)->d] : 1.)
#define GFS_FACE_FRACTION_LEFT(fa) GFS_FACE_FRACTION(fa)
#define GFS_FACE_FRACTION_RIGHT(fa) (GFS_IS_MIXED ((fa)->neighbor) ?\
                 GFS_STATE ((fa)->neighbor)->solid->s[FTT_OPPOSITE_DIRECTION ((fa)->d)] : 1.)

#define GFS_IS_FLUID(cell)      ((cell) != NULL &&\
                                 GFS_STATE (cell)->solid == NULL)
#define GFS_IS_MIXED(cell)      ((cell) != NULL &&\
                                 GFS_STATE (cell)->solid != NULL)
#define GFS_CELL_IS_BOUNDARY(cell) (((cell)->flags & GFS_FLAG_BOUNDARY) != 0)
#define GFS_CELL_IS_PERMANENT(cell) (((cell)->flags & GFS_FLAG_PERMANENT) != 0)
#define GFS_CELL_IS_GRADIENT_BOUNDARY(cell) (((cell)->flags & GFS_FLAG_GRADIENT_BOUNDARY) != 0)

FttCellFace           gfs_cell_face                 (FttCell * cell,
						     FttDirection d);
void                  gfs_cell_cleanup              (FttCell * cell,
						     GfsDomain * domain);
void                  gfs_cell_reset                (FttCell * cell, 
						     GfsVariable * v);
void                  gfs_face_reset                (FttCellFace * face,
						     GfsVariable * v);
void                  gfs_get_from_below_intensive  (FttCell * cell, 
						     const GfsVariable * v);
void                  gfs_cell_coarse_fine          (FttCell * cell,
						     GfsVariable * v);
gdouble               gfs_neighbor_value            (const FttCellFace * face,
						     guint v,
						     gdouble * x);
gdouble               gfs_face_interpolated_value   (const FttCellFace * face,
						     guint v);
gdouble               gfs_face_interpolated_value_generic (const FttCellFace * face, 
							   const GfsVariable * v);
gdouble               gfs_face_weighted_interpolated_value (const FttCellFace * face,
							    guint v);
typedef gdouble    (* GfsCenterGradient)            (FttCell * cell,
						     FttComponent c,
						     guint v);
gdouble               gfs_center_gradient           (FttCell * cell,
						     FttComponent c,
						     guint v);
void                  gfs_cm_gradient               (FttCell * cell,
						     GfsVariable * v,
						     FttVector * g);
void                  gfs_center_gradient_stencil   (FttCell * cell,
						     FttComponent c,
						     guint v);
gdouble               gfs_center_van_leer_gradient  (FttCell * cell,
						     FttComponent c,
						     guint v);
gdouble               gfs_center_minmod_gradient    (FttCell * cell,
						     FttComponent c,
						     guint v);
gdouble               gfs_center_superbee_gradient  (FttCell * cell,
						     FttComponent c,
						     guint v);
gdouble               gfs_center_sweby_gradient     (FttCell * cell,
						     FttComponent c,
						     guint v);
gdouble               gfs_center_regular_gradient   (FttCell * cell,
						     FttComponent c,
						     GfsVariable * v);
gdouble               gfs_center_regular_2nd_derivative (FttCell * cell, 
							 FttComponent c, 
							 GfsVariable * v);

typedef struct _GfsGradient GfsGradient;

struct _GfsGradient {
  gdouble a, b;
};

void                  gfs_face_gradient              (const FttCellFace * face,
						      GfsGradient * g,
						      guint v,
						      gint max_level);
void                  gfs_face_weighted_gradient     (const FttCellFace * face,
						      GfsGradient * g,
						      guint v,
						      gint max_level);
void                  gfs_face_weighted_gradient_2D  (const FttCellFace * face,
						      GfsGradient * g,
						      guint v,
						      gint max_level);
void                  gfs_face_cm_gradient           (const FttCellFace * face,
						      GfsGradient * g,
						      guint v,
						      gint max_level);
void                  gfs_face_cm_weighted_gradient  (const FttCellFace * face,
						      GfsGradient * g,
						      guint v,
						      gint max_level);
void                  gfs_cell_dirichlet_gradient    (FttCell * cell,
						      guint v,
						      gint max_level,
						      gdouble v0,
						      FttVector * grad);
void                  gfs_mixed_cell_gradient        (FttCell * cell,
						      GfsVariable * v,
						      FttVector * g);
gdouble               gfs_cell_dirichlet_gradient_flux (FttCell * cell,
							guint v,
							gint max_level,
							gdouble v0);
gdouble               gfs_cell_dirichlet_value         (FttCell * cell,
							GfsVariable * v,
							gint max_level);
gdouble               gfs_mixed_cell_interpolate       (FttCell * cell,
							FttVector p,
							GfsVariable * v);

void                  gfs_normal_divergence          (FttCell * cell,
						      GfsVariable * v);
void                  gfs_normal_divergence_2D       (FttCell * cell,
						      GfsVariable * v);
gdouble               gfs_divergence                 (FttCell * cell,
						      GfsVariable ** v);
gdouble               gfs_vorticity                  (FttCell * cell,
						      GfsVariable ** v);
gdouble               gfs_vector_norm                (FttCell * cell,
						      GfsVariable ** v);
gdouble               gfs_vector_norm2               (FttCell * cell,
						      GfsVariable ** v);
gdouble               gfs_vector_lambda2             (FttCell * cell,
						      GfsVariable ** v);
void                  gfs_pressure_force             (FttCell * cell,
						      GfsVariable * p,
						      FttVector * f);
GtsRange              gfs_stats_variable             (FttCell * root, 
						      GfsVariable * v, 
						      FttTraverseFlags flags,
						      gint max_depth);

typedef struct _GfsNorm GfsNorm;

struct _GfsNorm {
  gdouble bias, first, second, infty, w;
};

void                  gfs_norm_init                 (GfsNorm * n);
void                  gfs_norm_reset                (GfsNorm * n);
void                  gfs_norm_add                  (GfsNorm * n, 
						     gdouble val,
						     gdouble weight);
void                  gfs_norm_update               (GfsNorm * n);

GfsNorm               gfs_norm_variable             (FttCell * root, 
						     GfsVariable * v, 
						     FttTraverseFlags flags,
						     gint max_depth);
  
void                  gfs_cell_traverse_mixed       (FttCell * root,
						     FttTraverseType order,
						     FttTraverseFlags flags,
						     FttCellTraverseFunc func,
						     gpointer data);
void                  gfs_cell_corner_values        (FttCell * cell, 
						     GfsVariable * v, 
						     gint max_level,
						     gdouble f[4*(FTT_DIMENSION - 1)]);
gdouble               gfs_interpolate_from_corners  (FttCell * cell,
						     FttVector p,
						     gdouble * f);
gdouble               gfs_interpolate               (FttCell * cell,
						     FttVector p,
						     GfsVariable * v);
void                  gfs_interpolate_stencil       (FttCell * cell,
						     GfsVariable * v);
void                  ftt_cell_refine_corners       (FttCell * cell,
						     FttCellInitFunc init,
						     gpointer data);
gdouble               gfs_center_curvature          (FttCell * cell,
						     FttComponent c,
						     guint v);
gdouble               gfs_streamline_curvature      (FttCell * cell,
						     GfsVariable ** v);
void                  gfs_shear_strain_rate_tensor  (FttCell * cell, 
						     GfsVariable ** u,
						     gdouble t[FTT_DIMENSION][FTT_DIMENSION]);
gdouble               gfs_2nd_principal_invariant   (FttCell * cell, 
						     GfsVariable ** u);

typedef struct {
#if FTT_2D
  FttCell * c[7];
  gdouble w[7];
#else  /* 3D */
  FttCell * c[29];
  gdouble w[29];
#endif /* 3D */
  guint n;  
} GfsInterpolator;

void                  gfs_cell_corner_interpolator  (FttCell * cell,
						     FttDirection d[FTT_DIMENSION],
						     gint max_level,
						     gboolean centered,
						     GfsInterpolator * inter);
gdouble               gfs_cell_corner_value         (FttCell * cell,
						     FttDirection * d,
						     GfsVariable * v,
						     gint max_level);

/* GfsStencil: Header */

typedef struct _GfsLinearProblem GfsLinearProblem;
typedef struct _GfsStencil GfsStencil;

struct _GfsStencil {
  GArray * coeff; /* doubles */
  GArray * id;    /* ints */
};

GfsStencil *  gfs_stencil_new                       (FttCell * cell, 
						     GfsLinearProblem * lp,
						     gdouble coeff);
void          gfs_stencil_add_element               (GfsStencil * stencil,
						     FttCell * cell,
						     GfsLinearProblem * lp,
						     gdouble coeff);
void          gfs_stencil_destroy                   (GfsStencil * stencil);
void          gfs_face_weighted_gradient_stencil    (const FttCellFace * face,
						     GfsGradient * g,
						     gint max_level,
						     GfsLinearProblem * lp,
						     GfsStencil * stencil);
void          gfs_face_cm_weighted_gradient_stencil (const FttCellFace * face,
						     GfsGradient * g,
						     gint max_level,
						     GfsLinearProblem * lp,
						     GfsStencil * stencil);
gdouble       gfs_cell_dirichlet_gradient_flux_stencil (FttCell * cell,
							gint max_level,
							gdouble v0,
							GfsLinearProblem * lp,
							GfsStencil * stencil);

#ifdef __cplusplus
}
#endif /* __cplusplus */

#endif /* __FLUID_H__ */