Rythmos_IntegratorBuilder_decl.hpp

//@HEADER
// ***********************************************************************
//
//                           Rythmos Package
//                 Copyright (2006) Sandia Corporation
//
// Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
// license for use of this work by or on behalf of the U.S. Government.
//
// 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., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301
// USA
// Questions? Contact Todd S. Coffey (tscoffe@sandia.gov)
//
// ***********************************************************************
//@HEADER
 
 
#ifndef Rythmos_INTEGRATOR_BUILDER_DECL_H
#define Rythmos_INTEGRATOR_BUILDER_DECL_H
 
// Rythmos classes:
#include "Rythmos_Types.hpp"
#include "Rythmos_IntegratorBase.hpp"
#include "Rythmos_IntegrationControlStrategyBase.hpp"
#include "Rythmos_StepperBuilder.hpp"
#include "Rythmos_StepControlStrategyBase.hpp"
#include "Rythmos_InterpolationBufferBase.hpp"
#include "Rythmos_InterpolationBufferAppenderBase.hpp"
#include "Rythmos_ErrWtVecCalcBase.hpp"
#include "Rythmos_InterpolatorBase.hpp"
#include "Rythmos_RKButcherTableauBuilder.hpp"
 
// Thyra classes:
#include "Thyra_LinearOpWithSolveFactoryBase.hpp"
#include "Thyra_NonlinearSolverBase.hpp"
#include "Thyra_ModelEvaluator.hpp"
 
// Teuchos:
#include "Teuchos_ObjectBuilder.hpp"
#include "Teuchos_ParameterListAcceptor.hpp"
 
namespace Rythmos {
 
 
template<class Scalar>
  class IntegratorBuilder : virtual public Teuchos::ParameterListAcceptor
{
public:
 
  IntegratorBuilder();
 
  virtual ~IntegratorBuilder();
 
  void setIntegratorFactory(
    const RCP<const AbstractFactory<IntegratorBase<Scalar> > > &integratorFactory,
    const std::string &integratorFactoryName
    );
 
  void setIntegrationControlFactory(
    const RCP<const AbstractFactory<IntegrationControlStrategyBase<Scalar> > > &integrationControlFactory,
    const std::string &integrationControlName
    );
 
  void setStepperBuilder(
    const RCP<StepperBuilder<Scalar> > &stepperBuilder
    );
 
  RCP<StepperBuilder<Scalar> > getStepperBuilder();
 
  void setRKButcherTableauBuilder(
      const RCP<RKButcherTableauBuilder<Scalar> > & rkbtBuilder
      );
 
  void setStepControlFactory(
    const RCP<const AbstractFactory<StepControlStrategyBase<Scalar> > > &stepControlStrategyFactory,
    const std::string &stepControlName
    );
 
  void setInterpolationBufferFactory(
    const RCP<const AbstractFactory<InterpolationBufferBase<Scalar> > > &interpolationBufferFactory,
    const std::string &interpolationBufferName
    );
 
  void setInterpolationBufferAppenderFactory(
    const RCP<const AbstractFactory<InterpolationBufferAppenderBase<Scalar> > > &interpolationBufferAppenderFactory,
    const std::string &interpolationBufferAppenderName
    );
 
  void setErrWtVecCalcFactory(
    const RCP<const AbstractFactory<ErrWtVecCalcBase<Scalar> > > &errWtVecCalcFactory,
    const std::string &errWtVecCalcFactoryName
    );
  
  void setInterpolatorFactory(
    const RCP<const AbstractFactory<InterpolatorBase<Scalar> > > &interpolatorFactory,
    const std::string &interpolatorFactoryName
    );
  
  void setWFactoryObject(
    const RCP<Thyra::LinearOpWithSolveFactoryBase<Scalar> > &wFactoryObject
    );
  
  RCP<IntegratorBase<Scalar> > create(
    const RCP<const Thyra::ModelEvaluator<Scalar> >& model,
    const Thyra::ModelEvaluatorBase::InArgs<Scalar>& initialCondition,
    const RCP<Thyra::NonlinearSolverBase<Scalar> >& nlSolver
    ) const;
  
 
  void setParameterList(const RCP<Teuchos::ParameterList> & paramList);
  
  RCP<const Teuchos::ParameterList> getValidParameters() const;
 
  RCP<ParameterList> getNonconstParameterList();
 
  RCP<ParameterList> unsetParameterList();
  
  RCP<const ParameterList> getParameterList() const;
 
 
private:
 
  // //////////////////////////////////////
  // Private data members
 
  RCP<Teuchos::ObjectBuilder<IntegratorBase<Scalar> > > integratorBuilder_;
  RCP<Teuchos::ObjectBuilder<IntegrationControlStrategyBase<Scalar> > > integrationControlBuilder_;
  RCP<StepperBuilder<Scalar> > stepperBuilder_;
  RCP<RKButcherTableauBuilder<Scalar> > rkbtBuilder_;
  RCP<Teuchos::ObjectBuilder<StepControlStrategyBase<Scalar> > > stepControlBuilder_;
  RCP<Teuchos::ObjectBuilder<InterpolationBufferBase<Scalar> > > interpolationBufferBuilder_;
  RCP<Teuchos::ObjectBuilder<InterpolationBufferAppenderBase<Scalar> > > interpolationBufferAppenderBuilder_;
  RCP<Teuchos::ObjectBuilder<ErrWtVecCalcBase<Scalar> > > errWtVecCalcBuilder_;
  RCP<Teuchos::ObjectBuilder<InterpolatorBase<Scalar> > > interpolatorBuilder_;
 
  RCP<Thyra::LinearOpWithSolveFactoryBase<Scalar> > wFactoryObject_;
 
  RCP<ParameterList> paramList_;
  mutable RCP<ParameterList> validPL_;
 
  // //////////////////////////////////////
  // Private member functions
 
  void initializeDefaults_();
 
};
 
 
template<class Scalar>
RCP<IntegratorBuilder<Scalar> > integratorBuilder();
 
 
template<class Scalar>
RCP<IntegratorBuilder<Scalar> >
integratorBuilder(const RCP<ParameterList> &paramList);
 
template<class Scalar>
RCP<IntegratorBase<Scalar> > createForwardSensitivityIntegrator(
    const RCP<const Thyra::ModelEvaluator<Scalar> >& model,
    const int& p_index,
    const Thyra::ModelEvaluatorBase::InArgs<Scalar>& model_ic,
    const RCP<Thyra::NonlinearSolverBase<Scalar> >& nlSolver,
    const RCP<ParameterList>& integratorBuilderPL
    );
 
} // namespace Rythmos
 
#endif //Rythmos_INTEGRATOR_BUILDER_DECL_H