00001 // $Id: FEApp_GaussianQuadrature2.cpp,v 1.2 2007/07/24 00:08:04 etphipp Exp $ 00002 // $Source: /space/CVS/Trilinos/packages/sacado/example/FEApp/FEApp_GaussianQuadrature2.cpp,v $ 00003 // @HEADER 00004 // *********************************************************************** 00005 // 00006 // Sacado Package 00007 // Copyright (2006) Sandia Corporation 00008 // 00009 // Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation, 00010 // the U.S. Government retains certain rights in this software. 00011 // 00012 // This library is free software; you can redistribute it and/or modify 00013 // it under the terms of the GNU Lesser General Public License as 00014 // published by the Free Software Foundation; either version 2.1 of the 00015 // License, or (at your option) any later version. 00016 // 00017 // This library is distributed in the hope that it will be useful, but 00018 // WITHOUT ANY WARRANTY; without even the implied warranty of 00019 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00020 // Lesser General Public License for more details. 00021 // 00022 // You should have received a copy of the GNU Lesser General Public 00023 // License along with this library; if not, write to the Free Software 00024 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 00025 // USA 00026 // Questions? Contact David M. Gay (dmgay@sandia.gov) or Eric T. Phipps 00027 // (etphipp@sandia.gov). 00028 // 00029 // *********************************************************************** 00030 // @HEADER 00031 00032 #include "FEApp_GaussianQuadrature2.hpp" 00033 #include <cmath> // for sqrt 00034 00035 FEApp::GaussianQuadrature2::GaussianQuadrature2() : 00036 qp(2), 00037 w(2) 00038 { 00039 qp[0] = -1.0 / std::sqrt(3.0); 00040 qp[1] = -qp[0]; 00041 00042 w[0] = 1.0; 00043 w[1] = 1.0; 00044 } 00045 00046 FEApp::GaussianQuadrature2::~GaussianQuadrature2() 00047 { 00048 } 00049 00050 unsigned int 00051 FEApp::GaussianQuadrature2::numPoints() const 00052 { 00053 return 2; 00054 } 00055 00056 const std::vector<double>& 00057 FEApp::GaussianQuadrature2::quadPoints() const 00058 { 00059 return qp; 00060 } 00061 00062 const std::vector<double>& 00063 FEApp::GaussianQuadrature2::weights() const 00064 { 00065 return w; 00066 }
1.4.7