fe_type.C
Go to the documentation of this file.
1 
2 // The libMesh Finite Element Library.
3 // Copyright (C) 2002-2014 Benjamin S. Kirk, John W. Peterson, Roy H. Stogner
4 
5 // This library is free software; you can redistribute it and/or
6 // modify it under the terms of the GNU Lesser General Public
7 // License as published by the Free Software Foundation; either
8 // version 2.1 of the License, or (at your option) any later version.
9 
10 // This library is distributed in the hope that it will be useful,
11 // but WITHOUT ANY WARRANTY; without even the implied warranty of
12 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 // Lesser General Public License for more details.
14 
15 // You should have received a copy of the GNU Lesser General Public
16 // License along with this library; if not, write to the Free Software
17 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
18 
19 // Local includes
20 #include "libmesh/fe_type.h"
23 
24 namespace libMesh
25 {
26 
27 // ---------------------------------------
28 // FEType class members
29 
30 AutoPtr<QBase>
32  const int extraorder) const
33 {
34 
35  // Clough elements have at least piecewise cubic functions
36  if (family == CLOUGH)
37  {
38  // this seems ridiculous but for some reason gcc 3.3.5 wants
39  // this when using complex numbers (spetersen 04/20/06)
40  const unsigned int seven = 7;
41 
42  return AutoPtr<QBase>
43  (new QClough(dim,
44  static_cast<Order>
45  (std::max(static_cast<unsigned int>
46  (this->default_quadrature_order()), seven + extraorder))));
47  }
48 
49  return AutoPtr<QBase>
50  (new QGauss(dim, static_cast<Order>(this->default_quadrature_order()
51  + extraorder)));
52 }
53 
54 } // namespace libMesh

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

Hosted By:
SourceForge.net Logo