/usr/include/ug/quadrature.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 | // -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
/****************************************************************************/
/* */
/* File: quadrature.h */
/* */
/* Purpose: quadrature formulas (header file) */
/* */
/* Author: Christian Wieners */
/* Institut fuer Computeranwendungen III */
/* Universitaet Stuttgart */
/* Pfaffenwaldring 27 */
/* 70569 Stuttgart */
/* email: ug@ica3.uni-stuttgart.de */
/* */
/* History: Sep 25 95 begin */
/* */
/* Remarks: */
/* */
/****************************************************************************/
/* RCS_ID
$Header$
*/
/****************************************************************************/
/* */
/* auto include mechanism and other include files */
/* */
/****************************************************************************/
#ifndef __QUADRATURE__
#define __QUADRATURE__
#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 */
/* */
/****************************************************************************/
#ifdef __TWODIM__
#define MAX_INT_POINTS 33
#else
#define MAX_INT_POINTS 27
#endif
#define Q_NIP(p) ((p)->nip)
#define Q_LOCAL(p,i) ((DOUBLE*)((const DOUBLE_VECTOR_3D *)((p)->local))[i])
#define Q_WEIGHT(p,i) ((DOUBLE)((const DOUBLE *)((p)->weight))[i])
#define G_LOCAL(p) ((p)->local)
#define G_GLOBAL(p) ((p)->global)
#define G_WEIGHT(p) ((p)->weight)
#define G_JINV(p) ((p)->Jinv)
/****************************************************************************/
/* */
/* data structures exported by the corresponding source file */
/* */
/****************************************************************************/
typedef const struct quadrature {
INT nip; /* number of integration points */
INT order; /* NEW ! */ /* order of quadrature rule */
const DOUBLE_VECTOR_3D *local; /* array[nip] for local coordinates */
const DOUBLE *weight; /* array[nip] for weights */
} QUADRATURE;
typedef struct {
DOUBLE_VECTOR local;
DOUBLE_VECTOR global;
DOUBLE weight;
DOUBLE_VECTOR Jinv[DIM];
} GAUSS_POINT;
/****************************************************************************/
/* */
/* definition of exported global variables */
/* */
/****************************************************************************/
/****************************************************************************/
/* */
/* function declarations */
/* */
/****************************************************************************/
QUADRATURE *GetQuadrature(INT dim, INT n, INT order);
QUADRATURE *GetQuadratureRule(INT dim, INT n, INT order);
QUADRATURE *GetSymmetricQuadratureRule(INT dim, INT n, INT order);
INT GaussPoints(INT dim, INT n, INT order, DOUBLE_VECTOR *x, GAUSS_POINT *gp);
END_UGDIM_NAMESPACE
#endif
|