FEApp_SGGaussQuadJacobianGlobalFill.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_SGGAUSSQUADJACOBIANGLOBALFILL_HPP
00033 #define FEAPP_SGGAUSSQUADJACOBIANGLOBALFILL_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 
00043 namespace FEApp {
00044 
00045   class SGGaussQuadJacobianGlobalFill : public GlobalFill<SGJacobianType> {
00046   public:
00047 
00049     typedef FEApp::EvaluationTraits::apply<SGJacobianType>::type ScalarT;
00050     
00052     SGGaussQuadJacobianGlobalFill(
00053       const Teuchos::RCP<const FEApp::Mesh>& elementMesh,
00054       const Teuchos::RCP<const FEApp::AbstractQuadrature>& quadRule,
00055       const Teuchos::RCP< FEApp::AbstractPDE<SGJacobianType> >& pdeEquations,
00056       const std::vector< Teuchos::RCP<FEApp::NodeBC> >& nodeBCs,
00057       bool is_transient,
00058       const Teuchos::RCP<const Stokhos::OrthogPolyBasis<int,double> >& sgBasis,
00059       const Teuchos::RCP<const Stokhos::Quadrature<int,double> >& sgQuad,
00060       const Teuchos::RCP< FEApp::AbstractPDE<JacobianType> >& jacPDEEquations,
00061       const ParamVec* p,
00062       double alpha, double beta);
00063   
00065     virtual ~SGGaussQuadJacobianGlobalFill();
00066 
00068     virtual void 
00069     computeGlobalFill(FEApp::AbstractInitPostOp<SGJacobianType>& initPostOp);
00070 
00071   private:
00072 
00074     SGGaussQuadJacobianGlobalFill(const SGGaussQuadJacobianGlobalFill&);
00075 
00077     SGGaussQuadJacobianGlobalFill& operator=(const SGGaussQuadJacobianGlobalFill&);
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<JacobianType> > jacPDE;
00088     const ParamVec* p;
00089     double alpha;
00090     double beta;
00091     const Teuchos::Array< Teuchos::Array<double> >& quad_points;
00092     const Teuchos::Array<double>& quad_weights;
00093     const Teuchos::Array< Teuchos::Array<double> >& quad_values;
00094     const Teuchos::Array<double>& norms;
00095     unsigned int sg_size;
00096     unsigned int nqp;
00097     std::vector<FadType> x;
00098     std::vector<FadType>* xdot;
00099     std::vector<FadType> f;
00100 
00101     std::vector<double> xqp;
00102     std::vector<double> xdotqp;
00103     std::vector<double> pqp;
00104     std::vector<double> fqp;
00105 
00106     std::vector<double> qv;
00107     std::vector<double> sqv;
00108 
00109     std::vector<double> sg_x;
00110     std::vector<double> sg_xdot;
00111     std::vector<double> sg_p;
00112     std::vector<double> sg_f;
00113 
00115     Teuchos::BLAS<int,double> blas;
00116 
00117   };
00118 
00119 }
00120 
00121 #endif // SG_ACTIVE
00122 
00123 #endif // SGGAUSSQUADRESIDUALGLOBALFILL_HPP

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