quadrature_gauss.h
Go to the documentation of this file.
1 // The libMesh Finite Element Library.
2 // Copyright (C) 2002-2014 Benjamin S. Kirk, John W. Peterson, Roy H. Stogner
3 
4 // This library is free software; you can redistribute it and/or
5 // modify it under the terms of the GNU Lesser General Public
6 // License as published by the Free Software Foundation; either
7 // version 2.1 of the License, or (at your option) any later version.
8 
9 // This library is distributed in the hope that it will be useful,
10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 // Lesser General Public License for more details.
13 
14 // You should have received a copy of the GNU Lesser General Public
15 // License along with this library; if not, write to the Free Software
16 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
17 
18 
19 
20 #ifndef LIBMESH_QUADRATURE_GAUSS_H
21 #define LIBMESH_QUADRATURE_GAUSS_H
22 
23 // Local includes
24 #include "libmesh/quadrature.h"
25 
26 // C++ includes
27 
28 namespace libMesh
29 {
30 
31 
32 
33 
40 // ------------------------------------------------------------
41 // QGauss class definition
42 
43 class QGauss : public QBase
44 {
45  public:
46 
50  QGauss (const unsigned int _dim,
51  const Order _order=INVALID_ORDER);
52 
56  ~QGauss();
57 
61  QuadratureType type() const { return QGAUSS; }
62 
63 
64  private:
65 
66  void init_1D (const ElemType _type=INVALID_ELEM,
67  unsigned int p_level=0);
68  void init_2D (const ElemType _type=INVALID_ELEM,
69  unsigned int p_level=0);
70  void init_3D (const ElemType _type=INVALID_ELEM,
71  unsigned int p_level=0);
72 
73 
80  void dunavant_rule(const Real rule_data[][4],
81  const unsigned int n_pts);
82 
83  void dunavant_rule2(const Real* wts,
84  const Real* a,
85  const Real* b,
86  const unsigned int* permutation_ids,
87  const unsigned int n_wts);
88 
95  void keast_rule(const Real rule_data[][4],
96  const unsigned int n_pts);
97 
98 };
99 
100 
101 
102 // ------------------------------------------------------------
103 // QGauss class members
104 inline
105 QGauss::QGauss(const unsigned int d,
106  const Order o) : QBase(d,o)
107 {
108  // explicitly call the init function in 1D since the
109  // other tensor-product rules require this one.
110  // note that EDGE will not be used internally, however
111  // if we called the function with INVALID_ELEM it would try to
112  // be smart and return, thinking it had already done the work.
113  if (_dim == 1)
114  init(EDGE2);
115 }
116 
117 
118 
119 
120 inline
122 {
123 }
124 
125 
126 } // namespace libMesh
127 
128 
129 
130 #endif // LIBMESH_QUADRATURE_GAUSS_H

Site Created By: libMesh Developers
Last modified: February 07 2014 16:57:06 UTC

Hosted By:
SourceForge.net Logo