quadrature_gm_3D.C
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 // Local includes
21 #include "libmesh/quadrature_gm.h"
22 
23 namespace libMesh
24 {
25 
26 
27 
29  unsigned int p)
30 {
31  // Nearly all GM rules contain negative weights, so if you are not
32  // allowing rules with negative weights, we cannot continue!
34  {
35  libMesh::err << "You requested a Grundmann-Moller rule but\n"
36  << "are not allowing rules with negative weights!\n"
37  << "Either select a different quadrature class or\n"
38  << "set allow_rules_with_negative_weights==true."
39  << std::endl;
40 
41  libmesh_error();
42  }
43 
44  switch (type_in)
45  {
46  case TET4:
47  case TET10:
48  {
49  // Untested above _order=23 but should work...
50  gm_rule( (_order + 2*p)/2 );
51  return;
52 
53  } // end case TET4, TET10
54 
55 
56 
57  //---------------------------------------------
58  // Unsupported element type
59  default:
60  {
61  libMesh::err << "ERROR: Unsupported element type: " << type_in << std::endl;
62  libmesh_error();
63  }
64  } // end switch (type_in)
65 
66  // We must have returned or errored-out by this point. If not,
67  // throw an error now.
68  libmesh_error();
69  return;
70 }
71 
72 } // namespace libMesh

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

Hosted By:
SourceForge.net Logo