epetra/1Dfem/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"

#include "1DfemTransient.H"

class ExampleApplication1Dfem : public EpetraExt::ModelEvaluator {
public:

  // Constructor
  ExampleApplication1Dfem(Teuchos::RefCountPtr<Epetra_Comm> &epetra_comm_ptr, Teuchos::ParameterList &params);

  // Initialization
  void initialize(Teuchos::RefCountPtr<Epetra_Comm> &epetra_comm_ptr, Teuchos::ParameterList &params);


  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;

  Teuchos::RefCountPtr<Epetra_Vector> get_exact_solution( double t ) const;


private:

    // Epetra Comm:
    Teuchos::RefCountPtr<Epetra_Comm> epetra_comm_ptr_;
    // Epetra Map:
    Teuchos::RefCountPtr<const Epetra_Map> epetra_map_ptr_;
    
    // Global number of unknowns:
    int numElements_;

    Teuchos::RefCountPtr<Epetra_CrsGraph> W_graph_;

    // This ModelEvaluator object will own a TransientInterface object
    Teuchos::RefCountPtr<TransientInterface> problemInterfacePtr_;
};

#endif // EXAMPLE_APPLICATION_HPP

---