HDF5Import_decl.hpp

#ifndef HDF5IMPORT_DECL_hpp
#define HDF5IMPORT_DECL_hpp
 
#include <fstream>
#include "feddlib/core/FEDDCore.hpp"
#include "feddlib/core/LinearAlgebra/MultiVector.hpp"
// Trilinos
#include <Teuchos_Array.hpp>
 
#include <Epetra_ConfigDefs.h>
#ifdef HAVE_MPI
#include <mpi.h>
#include <Epetra_MpiComm.h>
#include <Epetra_SerialComm.h>
#else
#include <Epetra_SerialComm.h>
#endif
 
#include <Epetra_Map.h>
#include <Epetra_MultiVector.h>
 
#include <EpetraExt_HDF5.h>
#include <hdf5.h>

 
namespace FEDD {
template <class SC = default_sc, class LO = default_lo, class GO = default_go, class NO = default_no>
class HDF5Import {
public:
    typedef Teuchos::RCP<Epetra_Map> EpetraMapPtr_Type;
 
    typedef EpetraExt::HDF5 HDF5_Type;
    typedef Teuchos::RCP<HDF5_Type> HDF5Ptr_Type;
    
    typedef Teuchos::Comm<int> Comm_Type;
    typedef Teuchos::RCP<const Comm_Type> CommConstPtr_Type;
    
    typedef Map<LO,GO,NO> Map_Type;
    typedef Teuchos::RCP<const Map_Type> MapConstPtr_Type;
    
    typedef MultiVector<SC,LO,GO,NO> MultiVector_Type;
    typedef Teuchos::RCP<MultiVector_Type> MultiVectorPtr_Type;
 

    HDF5Import(MapConstPtr_Type readMap, std::string inputFilename);

    MultiVectorPtr_Type readVariablesHDF5(std::string varName);
 
    // Closing Importer
    void closeImporter();
 
  protected:
    
    HDF5Ptr_Type hdf5importer_; 
    CommConstPtr_Type comm_;
    Teuchos::RCP<Epetra_MpiComm> commEpetra_;
    
    // ------------------------
    // READ 
    // ------------------------
    std::string inputFilename_; 
    EpetraMapPtr_Type  readMap_;
    Epetra_MultiVector* u_import_Epetra_;
    MultiVectorPtr_Type u_import_Tpetra_;
 
};
 
}
 
#endif