1#ifndef DUNE_FEM_PYRAMIDPOINTS_HH
2#define DUNE_FEM_PYRAMIDPOINTS_HH
4#include <dune/common/fvector.hh>
5#include <dune/common/visibility.hh>
24 const FieldVector<double, 3>&
point(
int m,
int i)
const
56 G[m][0][0] =0.58541020;
57 G[m][0][1] =0.72819660;
58 G[m][0][2] =0.13819660;
60 G[m][1][0] =0.13819660;
61 G[m][1][1] =0.72819660;
62 G[m][1][2] =0.13819660;
64 G[m][2][0] =0.13819660;
65 G[m][2][1] =0.27630920;
66 G[m][2][2] =0.58541020;
68 G[m][3][0] =0.13819660;
69 G[m][3][1] =0.27630920;
70 G[m][3][2] =0.13819660;
72 G[m][4][0] =0.72819660;
73 G[m][4][1] =0.13819660;
74 G[m][4][2] =0.13819660;
76 G[m][5][0] =0.72819660;
77 G[m][5][1] =0.58541020;
78 G[m][5][2] =0.13819660;
80 G[m][6][0] =0.27630920;
81 G[m][6][1] =0.13819660;
82 G[m][6][2] =0.58541020;
84 G[m][7][0] =0.27630920;
85 G[m][7][1] =0.13819660;
86 G[m][7][2] =0.13819660;
88 W[m][0] = 0.041666667;
89 W[m][1] = 0.041666667;
90 W[m][2] = 0.041666667;
91 W[m][3] = 0.041666667;
92 W[m][4] = 0.041666667;
93 W[m][5] = 0.041666667;
94 W[m][6] = 0.041666667;
95 W[m][7] = 0.041666667;
Definition: bindguard.hh:11
Definition: pyramidpoints.hh:17
const FieldVector< double, 3 > & point(int m, int i) const
Access to the ith point of quadrature rule m.
Definition: pyramidpoints.hh:24
@ MAXP
Definition: pyramidpoints.hh:20
@ highest_order
Definition: pyramidpoints.hh:21
PyramidPoints()
initialize quadrature points on the interval for all orders
Definition: pyramidpoints.hh:48
int numPoints(int m) const
Number of points in the quadrature rule m.
Definition: pyramidpoints.hh:42
int order(int m) const
Actual order of quadrature rule m.
Definition: pyramidpoints.hh:36
double weight(int m, int i) const
Access to the ith weight of quadrature rule m.
Definition: pyramidpoints.hh:30
@ numQuads
Definition: pyramidpoints.hh:19