src/tetmesh/paralleltetmesh.h

Go to the documentation of this file.

//
// C++ Interface: paralleltetmesh
//
// Description: 
//
//
// Author: Dag Evensberget <deg.evensberget@psi.ch>, (C) 2005
//
// Copyright: See COPYING file that comes with this distribution
//
//
#ifndef ParallelTetMesh_h
#define ParallelTetMesh_h

#include <algorithm>
#include <map>
#include <vector>
#include <utility>
#include "mpi.h"
#include "entity.h"
#include "tetmesh.h"
#include "utility/parallel_tools.h"

namespace mesh {

class ParallelTetMesh : public TetMesh {

public:
    ParallelTetMesh(MPI_Comm _comm);

    virtual ~ParallelTetMesh();
    
    void initPoint(int nofPoint);
    
    void insertPoint(int id, double x, double y, double z);  
    
    void initTet(int nofTet);
    
    void insertTet(id_t global_id, id_t id0, id_t id1, id_t id2, id_t id3, id_t material);
    
    id_t lookupFace(id_t node0, id_t node1, id_t node2);

    void fix_edge_orientation();
    
    void finalize_mesh(int nofSym);
    
    id_t point_gid(id_t lid) const;
    
    id_t face_gid(id_t lid) const;
    
    id_t edge_gid(id_t lid) const;
    
    MPI_Comm get_comm();
    
    id_t get_dof_removed_offset(id_t my_offset);
    
    id_t get_dof_removed_offset();
    
    id_t get_num_global_points() const {
        return _num_global_points;
    }
    
    id_t get_num_global_edges() const {
        return _num_global_edges;    
    }
    
    id_t get_num_global_faces() const {
        return _num_global_faces;    
    }
    
        id_t get_num_global_tets();
 
    bool is_owned_point(id_t gid) const;
    
    bool is_owned_edge(id_t gid) const;

    bool is_owned_face(id_t gid) const;

    const std::vector<id_t>& get_edge_gids() const {
        return _edge_gids;
    }
    
    const std::vector<id_t>& get_node_gids() const {
        return _point_gids;
    }
    
private:

    void synchronize();

    void set_face_gids();
    
    void set_edge_gids();
    
    void set_point_ownership();
    
    std::vector<bool> _point_owned_by_me;

    utility::triple<id_t, id_t, id_t> get_face_corner_gids(const Face& face) const;
    
    std::pair<id_t, id_t> get_edge_corner_gids(const Edge& edge) const;
 
    MPI_Comm _comm;
    
    int _mpi_size;
    
    int _mpi_rank;
    
    int _tet_counter;

    int _point_counter;
    
    int _num_global_points;
    
    id_t _num_global_edges;
        
    id_t _num_global_faces;
        
    std::vector<id_t> _face_gids;
    
    std::vector<id_t> _edge_gids;
    
    std::vector<id_t> _point_gids;
    
    std::map<int,int> _point_local_id;
        
    id_t _my_edges[4];

    id_t _my_faces[4];
};

}
#endif

---