basicExample/ExampleApplication.hpp

Go to the documentation of this file.

//@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., 59 Temple Place, Suite 330, Boston, MA 02111-1307
// USA
// Questions? Contact Todd S. Coffey (tscoffe@sandia.gov)
//
// ***********************************************************************
//@HEADER

#ifndef EXAMPLE_APPLICATION_HPP
#define EXAMPLE_APPLICATION_HPP

#include "EpetraExt_ModelEvaluator.h"
#include "Teuchos_ParameterList.hpp"
#include "Rythmos_ConfigDefs.h"
#include "Epetra_Map.h"
#include "Epetra_Vector.h"
#include "Epetra_Comm.h"
#include "Epetra_CrsGraph.h"

class ExampleApplication : public EpetraExt::ModelEvaluator {
public:

  // Constructor
  ExampleApplication(Teuchos::RefCountPtr<Epetra_Comm> &epetra_comm_ptr, Teuchos::ParameterList &params);
  // return ODE decay coefficient
  Teuchos::RefCountPtr<const Epetra_Vector> get_coeff() const;
  //
  Teuchos::RefCountPtr<const Epetra_Vector> get_exact_solution(double t) const;


  Teuchos::RefCountPtr<const Epetra_Map> get_x_map() const;
  Teuchos::RefCountPtr<const Epetra_Map> get_f_map() const;
  Teuchos::RefCountPtr<const Epetra_Vector> get_x_init() const;
  Teuchos::RefCountPtr<const Epetra_Vector> get_x_dot_init() const;
  Teuchos::RefCountPtr<Epetra_Operator> create_W() const;
  InArgs createInArgs() const;
  OutArgs createOutArgs() const;
  void evalModel( const InArgs& inArgs, const OutArgs& outArgs ) const;
  double evalR(const double& t, const double& lambda, const double& s) const;


private:

    // Epetra Comm:
    Teuchos::RefCountPtr<Epetra_Comm> epetra_comm_ptr_;
    // Epetra Map:
    Teuchos::RefCountPtr<Epetra_Map> epetra_map_ptr_;
    
    // Mode
    bool implicit_;
    // Global number of unknowns:
    int numElements_;
    // Coefficients for ODE
    double lambda_min_;
    double lambda_max_;
    double coeff_s_; // Coefficient for forcing term
    std::string lambda_fit_;
    Teuchos::RefCountPtr<Epetra_Vector> lambda_ptr_;
    // Constant initial condition for the problem:
    double x0_;

    Teuchos::RefCountPtr<Epetra_CrsGraph> W_graph_;

};

#endif // EXAMPLE_APPLICATION_HPP

---