00001 /* @HEADER@ */ 00002 // ************************************************************************ 00003 // 00004 // Sundance 00005 // Copyright (2005) Sandia Corporation 00006 // 00007 // Copyright (year first published) Sandia Corporation. Under the terms 00008 // of Contract DE-AC04-94AL85000 with Sandia Corporation, the U.S. Government 00009 // retains certain rights in this software. 00010 // 00011 // This library is free software; you can redistribute it and/or modify 00012 // it under the terms of the GNU Lesser General Public License as 00013 // published by the Free Software Foundation; either version 2.1 of the 00014 // License, or (at your option) any later version. 00015 // 00016 // This library is distributed in the hope that it will be useful, but 00017 // WITHOUT ANY WARRANTY; without even the implied warranty of 00018 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00019 // Lesser General Public License for more details. 00020 // 00021 // You should have received a copy of the GNU Lesser General Public 00022 // License along with this library; if not, write to the Free Software 00023 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 00024 // USA 00025 // Questions? Contact Kevin Long (krlong@sandia.gov), 00026 // Sandia National Laboratories, Livermore, California, USA 00027 // 00028 // ************************************************************************ 00029 /* @HEADER@ */ 00030 00031 #ifndef SUNDANCE_FOURIER_H 00032 #define SUNDANCE_FOURIER_H 00033 00034 #include "SundanceDefs.hpp" 00035 #include "Teuchos_RefCountPtr.hpp" 00036 #include "SundanceBasisFamilyBase.hpp" 00037 00038 namespace Sundance 00039 { 00040 /** 00041 * Fourier basis 00042 */ 00043 class Fourier : public ScalarBasis 00044 { 00045 public: 00046 /** Construct a Fourier basis with the specified maximum frequency */ 00047 Fourier(int N); 00048 00049 /** 00050 * \brief Inform caller as to whether a given cell type is supported 00051 */ 00052 bool supportsCellTypePair( 00053 const CellType& maximalCellType, 00054 const CellType& cellType 00055 ) const ; 00056 00057 /** */ 00058 void print(std::ostream& os) const ; 00059 00060 /** */ 00061 int order() const {return 2*N_+1;} 00062 00063 /** return the number of nodes for this basis on the given cell type */ 00064 int nReferenceDOFsWithoutFacets( 00065 const CellType& maximalCellType, 00066 const CellType& cellType 00067 ) const ; 00068 00069 /** */ 00070 void getReferenceDOFs( 00071 const CellType& maximalCellType, 00072 const CellType& cellType, 00073 Array<Array<Array<int> > >& dofs) const ; 00074 00075 /** */ 00076 void refEval( 00077 const CellType& cellType, 00078 const Array<Point>& pts, 00079 const SpatialDerivSpecifier& deriv, 00080 Array<Array<Array<double> > >& result, 00081 int verbosity=0) const ; 00082 00083 /* Handleable boilerplate */ 00084 GET_RCP(BasisFamilyBase); 00085 00086 private: 00087 static Array<int> makeRange(int low, int high); 00088 00089 /** evaluate on a line cell */ 00090 void evalOnLine(const Point& pt, 00091 const MultiIndex& deriv, 00092 Array<double>& result) const ; 00093 00094 /** the order of the basis*/ 00095 int N_; 00096 00097 00098 }; 00099 } 00100 00101 #endif