FEApp_SGGaussQuadResidualGlobalFill.hpp

Go to the documentation of this file.
00001 // $Id$ 
00002 // $Source$ 
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 #ifndef FEAPP_SGGAUSSQUADRESIDUALGLOBALFILL_HPP
00033 #define FEAPP_SGGAUSSQUADRESIDUALGLOBALFILL_HPP
00034 
00035 #include "FEApp_TemplateTypes.hpp"
00036 #if SG_ACTIVE
00037 
00038 #include "FEApp_GlobalFill.hpp"
00039 #include "Stokhos_OrthogPolyBasis.hpp"
00040 #include "Stokhos_Quadrature.hpp"
00041 #include "Sacado_ScalarParameterVector.hpp"
00042 #include "Teuchos_BLAS.hpp"
00043 
00044 namespace FEApp {
00045 
00046   class SGGaussQuadResidualGlobalFill : public GlobalFill<SGResidualType> {
00047   public:
00048 
00050     typedef FEApp::EvaluationTraits::apply<SGResidualType>::type ScalarT;
00051     
00053     SGGaussQuadResidualGlobalFill(
00054       const Teuchos::RCP<const FEApp::Mesh>& elementMesh,
00055       const Teuchos::RCP<const FEApp::AbstractQuadrature>& quadRule,
00056       const Teuchos::RCP< FEApp::AbstractPDE<SGResidualType> >& pdeEquations,
00057       const std::vector< Teuchos::RCP<FEApp::NodeBC> >& nodeBCs,
00058       bool is_transient,
00059       const Teuchos::RCP<const Stokhos::OrthogPolyBasis<int,double> >& sgBasis,
00060       const Teuchos::RCP<const Stokhos::Quadrature<int,double> >& sgQuad,
00061       const Teuchos::RCP< FEApp::AbstractPDE<ResidualType> >& resPDEEquations,
00062       const ParamVec* pvec);
00063   
00065     virtual ~SGGaussQuadResidualGlobalFill();
00066 
00068     virtual void 
00069     computeGlobalFill(FEApp::AbstractInitPostOp<SGResidualType>& initPostOp);
00070 
00071   private:
00072 
00074     SGGaussQuadResidualGlobalFill(const SGGaussQuadResidualGlobalFill&);
00075 
00077     SGGaussQuadResidualGlobalFill& operator=(const SGGaussQuadResidualGlobalFill&);
00078 
00079   protected:
00080     
00082     Teuchos::RCP<const Stokhos::OrthogPolyBasis<int,double> > sg_basis;
00083 
00085     Teuchos::RCP<const Stokhos::Quadrature<int,double> > sg_quad;
00086 
00087     Teuchos::RCP< FEApp::AbstractPDE<ResidualType> > residPDE;
00088     const ParamVec* p;
00089     const Teuchos::Array< Teuchos::Array<double> >& quad_points;
00090     const Teuchos::Array<double>& quad_weights;
00091     const Teuchos::Array< Teuchos::Array<double> >& quad_values;
00092     const Teuchos::Array<double>& norms;
00093     unsigned int sg_size;
00094     unsigned int nqp;
00095     std::vector<double> x;
00096     std::vector<double>* xdot;
00097     std::vector<double> f;
00098 
00099     std::vector<double> xqp;
00100     std::vector<double> xdotqp;
00101     std::vector<double> pqp;
00102     std::vector<double> fqp;
00103 
00104     std::vector<double> qv;
00105     std::vector<double> sqv;
00106 
00107     std::vector<double> sg_x;
00108     std::vector<double> sg_xdot;
00109     std::vector<double> sg_p;
00110     std::vector<double> sg_f;
00111 
00113     Teuchos::BLAS<int,double> blas;
00114 
00115   };
00116 
00117 }
00118 
00119 #endif // SG_ACTIVE
00120 
00121 #endif // SGGAUSSQUADRESIDUALGLOBALFILL_HPP

Generated on Wed May 12 21:39:33 2010 for Sacado Package Browser (Single Doxygen Collection) by  doxygen 1.4.7