quadrature_grid_1D.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 // C++ includes
00021 
00022 // Local includes
00023 #include "libmesh/quadrature_grid.h"
00024 
00025 namespace libMesh
00026 {
00027 
00028 
00029 
00030 void QGrid::init_1D(const ElemType,
00031                     unsigned int)
00032 {
00033   //----------------------------------------------------------------------
00034   // 1D quadrature rules
00035 
00036   // We ignore p - the grid rule is just for experimentation
00037 
00038   _points.resize(_order + 1);
00039   _weights.resize(_order + 1);
00040   for (int i = 0; i != _order + 1; ++i)
00041     {
00042       _points[i](0) = 2.0 * (double)i / (double)_order - 1.0;
00043       _weights[i] = 2.0 / (double)(_order + 1);
00044     }
00045   return;
00046 }
00047 
00048 } // namespace libMesh
00049 
00050
Site Created By: libMesh Developers
Last modified: February 05 2013 19:54:48 UTC

Hosted By: