par_geom.h

Go to the documentation of this file.

#ifndef PAR_GEOM_FILE
#define PAR_GEOM_FILE
#include "geom.h"
#include "vdop.h"
#include <array>
#include <functional>             // function; C++11
#include <iostream>
#include <map>
#include <memory>                 // shared_ptr
#include <mpi.h>                  // MPI
#include <string>
#include <vector>
 
class ParMesh: public Mesh
{
public:
    explicit ParMesh(int ndim, int nvert_e = 0, int ndof_e = 0, int nedge_e = 0, MPI_Comm const &icomm = MPI_COMM_WORLD);
 
    ParMesh(ParMesh const &) = default;
 
    ParMesh &operator=(ParMesh const &) = delete;
 
 
    virtual ~ParMesh();
 
    explicit ParMesh(std::string const &sname, MPI_Comm const &icomm = MPI_COMM_WORLD);
 
    void VecAccu(std::vector<double> &w) const;
 
    double dscapr(std::vector<double> const &x, std::vector<double> const &y) const
    {
        return par_scalar(x, y, _icomm);
    }
 
    void Visualize(std::vector<double> const &v) const override;
 
private:
    std::vector<int> ReadElementSubdomains(std::string const &dname);
 
 
    void Transform_Local2Global_Vertex(int myrank, std::vector<int> const &t2d);
 
 
    void Generate_VectorAdd();
 
    bool CheckInterfaceExchange_InPlace() const;
    bool CheckInterfaceExchange() const;
    bool CheckInterfaceAdd_InPlace() const;
    bool CheckInterfaceAdd() const;
 
 
public:
    int MyRank() const
    {
        return _myrank;
    }
 
    int NumProcs() const
    {
        return _numprocs;
    }
 
    MPI_Comm GetCommunicator() const
    {
        return _icomm;
    }
 
private:
    // Don't use  &_icomm  ==> Error
    MPI_Comm const _icomm;                     
    int _numprocs;                             
    int _myrank;                               
    std::vector<int> _v_l2g;                   
    std::vector<int> _t_l2g;                   
    std::map<int, int> _v_g2l;                 
    std::map<int, int> _t_g2l;                 
 
    //std::vector<int> e_l2g;                    //!< edges:     local to global mapping
 
    std::vector<int> _valence;                 
    //  MPI_Alltoallv needs:
    mutable std::vector<double> _sendbuf;      
    std::vector<int>    _sendcounts;           
    std::vector<int>    _sdispls;              
    //
    //  We need to map the interface vertices onto the sendbuffer:
    std::vector<int> _loc_itf;                 
    std::vector<int> _gloc_itf;                
    std::vector<int> _buf2loc;                 
 
 
};
 
 
#endif