FEApp_ConstantNodeBCStrategy.hpp

Go to the documentation of this file.

// $Id: FEApp_ConstantNodeBCStrategy.hpp,v 1.2.2.2 2007/08/14 00:19:05 etphipp Exp $ 
// $Source: /space/CVS/Trilinos/packages/sacado/example/FEApp/FEApp_ConstantNodeBCStrategy.hpp,v $ 
// @HEADER
// ***********************************************************************
// 
//                           Sacado Package
//                 Copyright (2006) Sandia Corporation
// 
// Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
// the U.S. Government retains certain rights in this software.
// 
// This library is free software; you can redistribute it and/or modify
// it under the terms of the GNU Lesser General Public License as
// published by the Free Software Foundation; either version 2.1 of the
// License, or (at your option) any later version.
//  
// This library is distributed in the hope that it will be useful, but
// WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
// Lesser General Public License for more details.
//  
// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
// USA
// Questions? Contact David M. Gay (dmgay@sandia.gov) or Eric T. Phipps
// (etphipp@sandia.gov).
// 
// ***********************************************************************
// @HEADER

#ifndef FEAPP_CONSTANTNODEBCSTRATEGY_HPP
#define FEAPP_CONSTANTNODEBCSTRATEGY_HPP

#include "FEApp_AbstractNodeBCStrategy.hpp"
#include "Sacado_ScalarParameterLibrary.hpp"

namespace FEApp {

  template <typename ScalarT>
  class ConstantNodeBCStrategy : 
    public FEApp::AbstractNodeBCStrategy<ScalarT> {
  public:

    ConstantNodeBCStrategy(
    unsigned int solution_index, 
    unsigned int residual_index,
    const ScalarT& value,
    unsigned int bc_id,
    const Teuchos::RCP<Sacado::ScalarParameterLibrary>& paramLib);

    virtual ~ConstantNodeBCStrategy();

    const std::vector<unsigned int>& getOffsets() const;

    virtual void evaluateResidual(const std::vector<ScalarT>* dot,
          const std::vector<ScalarT>& solution,
          std::vector<ScalarT>& residual) const;

    void setValue(const ScalarT& value, bool mark_constant);

    const ScalarT& getValue() const;

  private:
    
    ConstantNodeBCStrategy(const ConstantNodeBCStrategy&);

    ConstantNodeBCStrategy& operator=(const ConstantNodeBCStrategy&);

  protected:
    
    unsigned int sol_index;

    unsigned int res_index;

    ScalarT val;

    std::vector<unsigned int> offsets;

  };

  class ConstantNodeBCStrategy_TemplateBuilder {
  public:
    ConstantNodeBCStrategy_TemplateBuilder(
         unsigned int solution_index, 
         unsigned int residual_index,
         double value,
         unsigned int bc_id,
         const Teuchos::RCP<Sacado::ScalarParameterLibrary>& paramLib) :
      sol_index(solution_index), res_index(residual_index), val(value),
      bcid(bc_id), pl(paramLib) {}
    template <typename T>
    Teuchos::RCP<FEApp::AbstractNodeBCStrategy_NTBase> build() const {
      return Teuchos::rcp( new ConstantNodeBCStrategy<T>(sol_index, 
               res_index, 
               val,
               bcid,
               pl));
    }
  protected:
    unsigned int sol_index;
    unsigned int res_index;
    double val;
    unsigned int bcid;
    Teuchos::RCP<Sacado::ScalarParameterLibrary> pl;
  };

  template <typename ScalarT>
  class ConstantNodeBCParameter : 
    public Sacado::ScalarParameterEntry<ScalarT> {

  public:

    ConstantNodeBCParameter(
       const Teuchos::RCP< ConstantNodeBCStrategy<ScalarT> >& s) : 
      bc(s) {}

    virtual ~ConstantNodeBCParameter() {}

    virtual void setRealValue(double value) { 
      setValueAsConstant(ScalarT(value)); }

    virtual void setValueAsConstant(const ScalarT& value) { 
      bc->setValue(value, true); }
    
    virtual void setValueAsIndependent(const ScalarT& value) { 
      bc->setValue(value, false); }
    
    virtual const ScalarT& getValue() const { return bc->getValue(); }
    
  protected:  
    
    Teuchos::RCP< ConstantNodeBCStrategy<ScalarT> > bc;

  };

}

// Include implementation
#ifndef SACADO_ETI
#include "FEApp_ConstantNodeBCStrategyImpl.hpp"
#endif 

#endif // FEAPP_CONSTANTNODEBCSTRATEGY_HPP

---