DomainDecomposition.h

#ifndef DOMAINDECOMPOSITION_H_
#define DOMAINDECOMPOSITION_H_
 
#include "DomainDecompMPIBase.h"
 
 
class DomainDecomposition : public DomainDecompMPIBase {
public:
    DomainDecomposition();
 
    // documentation see father class (DomainDecompBase.h)
    ~~DomainDecomposition();
 
    void readXML(XMLfileUnits& xmlconfig) override;
 
    // documentation see father class (DomainDecompBase.h)
    double getBoundingBoxMin(int dimension, Domain* domain) override;
 
    // documentation see father class (DomainDecompBase.h)
    double getBoundingBoxMax(int dimension, Domain* domain) override;
 
    void balanceAndExchange(double lastTraversalTime, bool forceRebalancing, ParticleContainer* moleculeContainer, Domain* domain) override;
 
    void initCommunicationPartners(double cutoffRadius, Domain * domain,ParticleContainer* moleculeContainer);
 
    //returns a vector of the neighbour ranks in x y and z direction (only neighbours connected by an area to local area)
    std::vector<int> getNeighbourRanks() override;
    
    //returns a vector of all 26 neighbour ranks in x y and z direction
    std::vector<int> getNeighbourRanksFullShell() override;
 
    //returns the ranks of the neighbours
    std::vector<std::vector<std::vector<int>>> getAllRanks();
 
 
    // documentation in base class
    void prepareNonBlockingStage(bool forceRebalancing,
            ParticleContainer* moleculeContainer, Domain* domain,
            unsigned int stageNumber) override;
 
    // documentation in base class
    void finishNonBlockingStage(bool forceRebalancing,
            ParticleContainer* moleculeContainer, Domain* domain,
            unsigned int stageNumber) override;
 
    // documentation in base class
    bool queryBalanceAndExchangeNonBlocking(bool forceRebalancing, ParticleContainer* moleculeContainer, Domain* domain, double etime) override;
 
    std::vector<CommunicationPartner> getNeighboursFromHaloRegion(Domain* domain, const HaloRegion& haloRegion, double cutoff) override;
 
private:
    void initMPIGridDims();
 
    int _gridSize[DIMgeom]; 
    int _coords[DIMgeom]; 
};
 
#endif /* DOMAINDECOMPOSITION_H_ */