Midpoint.h

/*
 * CommunicationScheme.h
 *
 *  Created on: 14.07.2017
 *      Author: sauermann
 */
 
#pragma once
 
#include "Simulation.h"
#include "ZonalMethod.h"
 
class Midpoint: public ZonalMethod {
public:
    Midpoint() {
    }
    virtual ~~Midpoint() {
    }
 
    virtual std::vector<HaloRegion> getHaloImportForceExportRegions(HaloRegion& initialRegion, double cutoffRadius,
                                                                    bool coversWholeDomain[3],
                                                                    double cellLength[3]) override {
        // the midpoint traversal is cell based, so the halo region has to be cellLength wide.
        return getLeavingExportRegions(initialRegion, cellLength, coversWholeDomain);
    }
 
    virtual std::vector<HaloRegion> getHaloExportForceImportRegions(HaloRegion& initialRegion, double cutoffRadius,
                                                                    bool coversWholeDomain[3],
                                                                    double cellLength[3]) override {
        const std::function<bool(const int[3])> condition = [](const int[3])->bool {
            // no condition for leaving particles.
                return true;
            };
        // the midpoint traversal is cell based, so the halo region has to be cellLength wide.
        return getHaloRegionsConditionalInside(initialRegion, cellLength, coversWholeDomain, condition);
    }
};