quadrature_gm_3D.C

Go to the documentation of this file.
00001 // The libMesh Finite Element Library.
00002 // Copyright (C) 2002-2012 Benjamin S. Kirk, John W. Peterson, Roy H. Stogner
00003 
00004 // This library is free software; you can redistribute it and/or
00005 // modify it under the terms of the GNU Lesser General Public
00006 // License as published by the Free Software Foundation; either
00007 // version 2.1 of the License, or (at your option) any later version.
00008 
00009 // This library is distributed in the hope that it will be useful,
00010 // but WITHOUT ANY WARRANTY; without even the implied warranty of
00011 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
00012 // Lesser General Public License for more details.
00013 
00014 // You should have received a copy of the GNU Lesser General Public
00015 // License along with this library; if not, write to the Free Software
00016 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
00017 
00018 
00019 
00020 // Local includes
00021 #include "libmesh/quadrature_gm.h"
00022 
00023 namespace libMesh
00024 {
00025 
00026 
00027 
00028 void QGrundmann_Moller::init_3D(const ElemType type_in,
00029                                 unsigned int p)
00030 {
00031   // Nearly all GM rules contain negative weights, so if you are not
00032   // allowing rules with negative weights, we cannot continue!
00033   if (!allow_rules_with_negative_weights)
00034     {
00035       libMesh::err << "You requested a Grundmann-Moller rule but\n"
00036                     << "are not allowing rules with negative weights!\n"
00037                     << "Either select a different quadrature class or\n"
00038                     << "set allow_rules_with_negative_weights==true."
00039                     << std::endl;
00040 
00041       libmesh_error();
00042     }
00043 
00044   switch (type_in)
00045     {
00046     case TET4:
00047     case TET10:
00048       {
00049         // Untested above _order=23 but should work...
00050         gm_rule( (_order + 2*p)/2 );
00051         return;
00052 
00053       } // end case TET4, TET10
00054 
00055 
00056 
00057       //---------------------------------------------
00058       // Unsupported element type
00059     default:
00060       {
00061         libMesh::err << "ERROR: Unsupported element type: " << type_in << std::endl;
00062         libmesh_error();
00063       }
00064     } // end switch (type_in)
00065 
00066   // We must have returned or errored-out by this point.  If not,
00067   // throw an error now.
00068   libmesh_error();
00069   return;
00070 }
00071 
00072 } // namespace libMesh
Site Created By: libMesh Developers
Last modified: February 05 2013 19:54:48 UTC

Hosted By: