FEApp_SGDakotaResidualGlobalFill.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_SGDAKOTARESIDUALGLOBALFILL_HPP
00033 #define FEAPP_SGDAKOTARESIDUALGLOBALFILL_HPP
00034 
00035 #include "FEApp_TemplateTypes.hpp"
00036 #if SG_ACTIVE
00037 
00038 #include "Stokhos_ConfigDefs.h"
00039 #ifdef HAVE_DAKOTA
00040 
00041 #include "FEApp_GlobalFill.hpp"
00042 #include "Stokhos_OrthogPolyBasis.hpp"
00043 #include "Sacado_ScalarParameterVector.hpp"
00044 
00045 // Dakota includes
00046 #include "system_defs.h"
00047 #include "ParallelLibrary.H"
00048 #include "ProblemDescDB.H"
00049 #include "DakotaStrategy.H"
00050 #include "DakotaModel.H"
00051 #include "FEApp_DakotaElementResidualInterface.hpp"
00052 
00053 namespace FEApp {
00054 
00055   class SGDakotaResidualGlobalFill : public GlobalFill<SGResidualType> {
00056   public:
00057 
00059     typedef FEApp::EvaluationTraits::apply<SGResidualType>::type ScalarT;
00060     
00062     SGDakotaResidualGlobalFill(
00063       const Teuchos::RCP<const FEApp::Mesh>& elementMesh,
00064       const Teuchos::RCP<const FEApp::AbstractQuadrature>& quadRule,
00065       const Teuchos::RCP< FEApp::AbstractPDE<SGResidualType> >& pdeEquations,
00066       const std::vector< Teuchos::RCP<FEApp::NodeBC> >& nodeBCs,
00067       bool is_transient,
00068       const Teuchos::RCP<const Stokhos::OrthogPolyBasis<int,double> >& sgBasis,
00069       const Teuchos::RCP< FEApp::AbstractPDE<ResidualType> >& resPDEEquations,
00070       const ParamVec* pvec);
00071   
00073     virtual ~SGDakotaResidualGlobalFill();
00074 
00076     virtual void 
00077     computeGlobalFill(FEApp::AbstractInitPostOp<SGResidualType>& initPostOp);
00078 
00079   private:
00080 
00082     SGDakotaResidualGlobalFill(const SGDakotaResidualGlobalFill&);
00083 
00085     SGDakotaResidualGlobalFill& operator=(const SGDakotaResidualGlobalFill&);
00086 
00087   protected:
00088     
00090     Teuchos::RCP<const Stokhos::OrthogPolyBasis<int,double> > sg_basis;
00091     Teuchos::RCP< FEApp::AbstractPDE<ResidualType> > residPDE;
00092     const ParamVec* p;
00093     unsigned int sg_size;
00094     Dakota::ParallelLibrary parallel_lib;
00095     Dakota::ProblemDescDB problem_db;
00096     Dakota::Strategy selected_strategy;
00097     FEApp::DakotaElementResidualInterface *elemInterface;
00098   };
00099 
00100 }
00101 
00102 #endif // HAVE_DAKOTA
00103 
00104 #endif // SG_ACTIVE
00105 
00106 #endif // SGDAKOTARESIDUALGLOBALFILL_HPP

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