src/HDF5ParallelReader.h

Go to the documentation of this file.

#ifndef HDF5ParallelReader_H
#define HDF5ParallelReader_H

#include <sys/times.h>
#include <unistd.h>
#include <string>

#include <iostream>
#include <fstream>
#include <stdexcept>
#include "assert.h"
#include "hdf5.h"
#include "parmetis.h"
#include "entity.h"
#include "tetmeshbuilder.h"
#include "pbedefs.h"

#define REALMEMTYPE H5T_NATIVE_DOUBLE
#define INTMEMTYPE H5T_NATIVE_INT

using namespace std;

class HDF5ParallelReader {

public: 

    ~HDF5ParallelReader();
    
    HDF5ParallelReader(const char* inputFile, MeshBuilder* builder, MPI_Comm comm);
    
    void read();
    
    void export_dist(char* file);
    
    int main (int argc, char **argv);
    
private:
    
    double seconds(); 
    
    void print_tet_global_id(); 
    
    void read_tets(); 

    void read_mats(); 
    
    void partition(); 
    
    void sort_buffer(int* sendcnts);
    
    void distribute(); 
    
    void generate_point_global_id(int* slab_offsets, int* slab_lengths, int slab_count); 
    
    void read_points(); 
    
    void read_boundary(); 
    
    void open_file(); 
    
    void close_file(); 


private: // Private attributes
    
    const char* _inputFile;
    
    MeshBuilder* _builder;
    
    MPI_Comm _comm;
    
    MPI_Info _info;
    
    int _mpi_size;
    
    int _mpi_rank;
    
    double* _point_coords;
    
    int _numpoints_global;
    
    int _numpoints_local;
    
    int* _tetids;
    
    int _numtets_global;
    
    int _numtets_local;

    int* _matids;
    
    int* _faces;
    
    int _numbfaces_global;
    
    int _numbfaces_local;

    int* _tet_partition;

    int* _tet_global_id; 
    
    int* _point_global_id;

    hid_t _file_id;
    
    hid_t _plist_id;

};

#endif

---