/usr/include/ug/fvgeom.h is in libug-dev 3.12.1-2.
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 | // -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
/****************************************************************************/
/* */
/* File: fvgeom.h */
/* */
/* Purpose: geometry evaluation for Finite-Volume discretization */
/* dimension independent, general elmement */
/* */
/* Author: Peter Bastian */
/* Institut fuer Computeranwendungen III */
/* Universitaet Stuttgart */
/* Pfaffenwaldring 27 */
/* 70569 Stuttgart */
/* email: peter@ica3.uni-stuttgart.de */
/* fon: 0049-(0)711-685-7003 */
/* fax: 0049-(0)711-685-7000 */
/* */
/* History: 06.05.96 begin, ug version 3.2 */
/* 06.08.96 modification and extension (Henrik Rentz-Reichert) */
/* */
/* Remarks: */
/* */
/****************************************************************************/
/* RCS_ID
$Header$
*/
/****************************************************************************/
/* */
/* include files */
/* system include files */
/* application include files */
/* */
/****************************************************************************/
#ifndef __FVGEOM__
#define __FVGEOM__
#include "gm.h"
#include "namespace.h"
START_UGDIM_NAMESPACE
/****************************************************************************/
/* */
/* defines in the following order */
/* */
/* compile time constants defining static data size (i.e. arrays) */
/* other constants */
/* macros */
/* */
/****************************************************************************/
#ifndef MAXNC
#define MAXNC MAX_CORNERS_OF_ELEM /* just to make it more readable*/
#endif
#ifndef MAXE
#define MAXE MAX_EDGES_OF_ELEM /* just to make it more readable*/
#endif
#ifndef MAXS
#define MAXS MAX_SIDES_OF_ELEM /* just to make it more readable*/
#endif
#define MAXF MAX_EDGES_OF_ELEM /* max # scvf */
#define MAXBF (MAX_SIDES_OF_ELEM*MAX_CORNERS_OF_SIDE)
/* max # boundary faces */
/* defines to specify data to be filled by EvaluateShapesAndDerivatives */
#define FILL_CORNER_DATA (1<<0)
#define FILL_SHAPES (1<<1)
#define FILL_DERIV (1<<2)
#define FILL_GRAD (1<<3)
#define FILL_J (1<<4)
#define FILL_COE (1<<5)
#define FILL_SHAPES_AND_GRAD (FILL_SHAPES | FILL_GRAD)
#define FILL_SHAPES_GRAD_AND_DERIV (FILL_SHAPES | FILL_GRAD | FILL_DERIV)
#define FILL_ALL (FILL_CORNER_DATA | FILL_SHAPES | FILL_GRAD | FILL_DERIV | FILL_J | FILL_COE)
/* macros for access to exported data structures */
#define FVG_ELEM(p) ((p)->e)
#define FVG_TAG(p) ((p)->tag)
#define FVG_NSCV(p) ((p)->n_scv)
#define FVG_NSCVF(p) ((p)->n_scvf)
#define FVG_NSCVBF(p) ((p)->n_bf)
#define FVG_GCO(p,i) ((p)->co_global[i])
#define FVG_GCOPTR(p) ((p)->co_global)
#define FVG_LCO(p,i) ((p)->co_local[i])
#define FVG_GEM(p,i) ((p)->em_global[i])
#define FVG_LEM(p,i) ((p)->em_local[i])
#define FVG_GSM(p,i) ((p)->sm_global[i])
#define FVG_LSM(p,i) ((p)->sm_local[i])
#define FVG_GCM(p) ((p)->s_global)
#define FVG_LCM(p) ((p)->s_local)
#define FVG_SCV(p,i) ((p)->scv+(i))
#define FVG_SCVF(p,i) ((p)->scvf+(i))
#define FVG_SCVBF(p,i) ((p)->bf+(i))
#define FVG_SDV(p) (&((p)->coe_sdv))
#define FVG_COSDV(p,i) ((p)->co_sdv+(i))
#define FVG_IPSDV(p,i) (&(((p)->scvf[i]).sdv))
#define FVG_BIPSDV(p,i) (&(((p)->bf[i]).sdv))
#define SCV_CO(p) ((p)->co)
#define SCV_GCO(p) ((p)->center)
#define SCV_VOL(p) ((p)->volume)
#define SCV_NDPROP(p) ((p)->node_property)
#define SCVF_FROM(p) ((p)->i)
#define SCVF_TO(p) ((p)->j)
#define SCVF_LIP(p) ((p)->ip_local)
#define SCVF_GIP(p) ((p)->ip_global)
#define SCVF_NORMAL(p) ((p)->normal)
#define SCVF_SDV(p) ((p)->sdv)
#define SCVBF_FROM(p) ((p)->co)
#define SCVBF_SIDE(p) ((p)->side)
#define SCVBF_LIP(p) ((p)->ip_local)
#define SCVBF_PARAMPTR(p) ((p)->param)
#define SCVBF_PARAM(p,i) ((p)->param[i])
#define SCVBF_NORMAL(p) ((p)->normal)
#define SCVBF_AREA(p) ((p)->area)
#define SCVBF_SDV(p) ((p)->sdv)
#define SDV_SHAPEPTR(p) ((p)->shape)
#define SDV_SHAPE(p,i) ((p)->shape[i])
#define SDV_GRADPTRPTR(p) ((p)->grad)
#define SDV_GRADPTR(p,i) ((p)->grad[i])
#define SDV_GRAD(p,i,j) ((p)->grad[i][j])
#define SDV_JINV(p) ((p)->Jinv)
#define SDV_J(p) ((p)->J)
#define SDV_DETJ(p) ((p)->detJ)
/****************************************************************************/
/* */
/* definition of exported data structures */
/* */
/****************************************************************************/
/* shape functions and derivatives etc. */
typedef struct {
DOUBLE shape[MAXNC]; /* values of shape functions at ip */
DOUBLE_VECTOR grad[MAXNC]; /* derivatives of shape functions at ip */
DOUBLE J[DIM][DIM]; /* jacobian at ip */
DOUBLE Jinv[DIM][DIM]; /* inverse of jacobian at ip */
DOUBLE detJ; /* det of jacobian at ip */
} SD_VALUES;
/* geometry related data */
typedef struct {
INT co; /* # of corner */
DOUBLE_VECTOR center; /* node position */
DOUBLE volume; /* volume (area) of scv */
INT node_property; /* subdomain info */
} SubControlVolume; /* FV intersected with element */
typedef struct {
INT i,j; /* scvf separates corner i and j of elem*/
DOUBLE_VECTOR ip_local; /* integration point in local coords */
DOUBLE_VECTOR ip_global; /* integration point in global coords */
DOUBLE_VECTOR normal; /* normal on face at ip pointing to CV j*/
SD_VALUES sdv; /* shape fcts, deriv. etc. at scv-faces */
} SubControlVolumeFace;
typedef struct {
INT co; /* corresponding corner */
INT side; /* boundary side of element */
DOUBLE_VECTOR ip_local; /* integration point in local coords */
DOUBLE param[DIM-1]; /* local side coordinates */
DOUBLE_VECTOR normal; /* normal on face at ip pointing to CV j*/
DOUBLE area; /* area of boundary face */
SD_VALUES sdv; /* shape fcts, deriv. etc. at b-faces */
} BoundaryFace;
typedef struct {
const ELEMENT *e; /* data for this element */
INT tag; /* element type */
INT n_scv; /* # sub control volumes (==corners) */
INT n_scvf; /* # sub control volume faces (==ip's) */
INT n_bf; /* # boundary faces */
DOUBLE_VECTOR co_global[MAXNC]; /* points in global space, corners */
DOUBLE_VECTOR co_local[MAXNC]; /* points in local space, corners */
DOUBLE_VECTOR em_global[MAXE]; /* points in global space, edge midpoint*/
DOUBLE_VECTOR em_local[MAXE]; /* points in local space, edge midpoints*/
DOUBLE_VECTOR sm_global[MAXS]; /* points in global space, side midpoint*/
DOUBLE_VECTOR sm_local[MAXS]; /* points in local space, side midpoints*/
DOUBLE_VECTOR s_global; /* points in global space, center */
DOUBLE_VECTOR s_local; /* points in local space, center */
SD_VALUES coe_sdv; /* shape fcts, deriv. etc. at coe */
SD_VALUES co_sdv[MAXNC]; /* shape fcts, deriv. etc. at corners */
SubControlVolume scv[MAXNC]; /* sub control volumes */
SubControlVolumeFace scvf[MAXF]; /* sub control volume faces */
BoundaryFace bf[MAXBF]; /* boundary faces */
} FVElementGeometry; /* geometry data for a general element */
/****************************************************************************/
/* */
/* definition of exported functions */
/* */
/****************************************************************************/
/* finite volume geometry */
INT EvaluateFVGeometry (const ELEMENT *e, FVElementGeometry *geo);
INT EvaluateAFVGeometry (const ELEMENT *e, const DOUBLE *conv, FVElementGeometry *geo);
/* upwinding procedures */
INT GetFullUpwindShapes (const FVElementGeometry *geo, const DOUBLE_VECTOR IPVel[MAXF], DOUBLE Shape[MAXF][MAXNC]);
INT GetSkewedUpwindShapes (const FVElementGeometry *geo, const DOUBLE_VECTOR IPVel[MAXF], DOUBLE Shape[MAXF][MAXNC]);
INT GetLPSUpwindShapes (const FVElementGeometry *geo, const DOUBLE_VECTOR IPVel[MAXF], DOUBLE Shape[MAXF][MAXNC]);
INT GetMWSUpwindShapes (const FVElementGeometry *geo, const DOUBLE_VECTOR IPVel[MAXF], DOUBLE Shape[MAXF][MAXNC]);
INT GetMJRawRegularUpwindShapes (const FVElementGeometry *geo, const DOUBLE_VECTOR IPVel[MAXF],
DOUBLE NodalShape[MAXF][MAXNC], DOUBLE IPShape[MAXF][MAXF]);
INT GetMJRawPositiveUpwindShapes (const FVElementGeometry *geo, const DOUBLE_VECTOR IPVel[MAXF],
DOUBLE NodalShape[MAXF][MAXNC], DOUBLE IPShape[MAXF][MAXF]);
/* */
INT SideIsCut (INT tag, const DOUBLE_VECTOR *x, const DOUBLE_VECTOR ip, const DOUBLE_VECTOR vel, INT side, DOUBLE_VECTOR y);
/* aligned finite volumes */
INT AFVGeometry (const ELEMENT *theElement, FVElementGeometry *geo, DOUBLE_VECTOR Convection);
/* shape functions and their derivatives */
INT EvaluateShapesAndDerivatives (FVElementGeometry *geo, INT flags);
/* intersect polygon with line */
INT Intersect2d (INT nco, const DOUBLE_VECTOR *x, const DOUBLE_VECTOR vel, const DOUBLE_VECTOR pt,
INT *Side, DOUBLE lambda[DIM_OF_BND]);
/* init */
INT InitFiniteVolumeGeom (void);
END_UGDIM_NAMESPACE
#endif
|