Sacado Package Browser (Single Doxygen Collection): FEApp_SourceFunctionFactoryImpl.hpp Source File

00001 // $Id: FEApp_SourceFunctionFactoryImpl.hpp,v 1.5 2008/06/16 17:05:40 etphipp Exp $ 
00002 // $Source: /space/CVS/Trilinos/packages/sacado/example/FEApp/FEApp_SourceFunctionFactoryImpl.hpp,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 #include "Teuchos_TestForException.hpp"
00032 #include "FEApp_QuadraticSourceFunction.hpp"
00033 #include "FEApp_CubicSourceFunction.hpp"
00034 #include "FEApp_ExponentialSourceFunction.hpp"
00035 
00036 template <typename ScalarT>
00037 FEApp::SourceFunctionFactory<ScalarT>::SourceFunctionFactory(
00038       const Teuchos::RCP<Teuchos::ParameterList>& funcParams_,
00039       const Teuchos::RCP<Sacado::ScalarParameterLibrary>& paramLib_) :
00040   funcParams(funcParams_), paramLib(paramLib_)
00041 {
00042 }
00043 
00044 template <typename ScalarT>
00045 Teuchos::RCP< FEApp::AbstractSourceFunction<ScalarT> >
00046 FEApp::SourceFunctionFactory<ScalarT>::create()
00047 {
00048   Teuchos::RCP< FEApp::AbstractSourceFunction<ScalarT> > strategy;
00049 
00050   std::string& method = funcParams->get("Name", "Quadratic");
00051   if (method == "Quadratic") {
00052     double factor = funcParams->get("Nonlinear Factor", 1.0);
00053     strategy = 
00054       Teuchos::rcp(new FEApp::QuadraticSourceFunction<ScalarT>(factor,
00055                      paramLib));
00056   }
00057   else if (method == "Cubic") {
00058     double factor = funcParams->get("Nonlinear Factor", 1.0);
00059     strategy = 
00060       Teuchos::rcp(new FEApp::CubicSourceFunction<ScalarT>(factor,
00061                  paramLib));
00062   }
00063   else if (method == "Exponential") {
00064     double factor = funcParams->get("Nonlinear Factor", 1.0);
00065     strategy = 
00066       Teuchos::rcp(new FEApp::ExponentialSourceFunction<ScalarT>(factor,
00067                  paramLib));
00068   }
00069   else {
00070     TEST_FOR_EXCEPTION(true, Teuchos::Exceptions::InvalidParameter,
00071            std::endl << 
00072            "Error!  Unknown source function " << method << 
00073            "!" << std::endl << "Supplied parameter list is " << 
00074            std::endl << *funcParams);
00075   }
00076 
00077   return strategy;
00078 }