ZonalMethod.h

/*
 * ZonalMethod.hpp
 *
 *  Created on: Feb 28, 2017
 *      Author: seckler
 */
 
#pragma once
 
#include <functional>
#include <vector>
#include "parallel/HaloRegion.h"
 
class ZonalMethod {
public:
    ZonalMethod();
 
    virtual ~~ZonalMethod();
 
    virtual std::vector<HaloRegion> getHaloImportForceExportRegions(HaloRegion& initialRegion, double cutoffRadius,
                                                                    bool coversWholeDomain[3],
                                                                    double cellLength[3]) = 0;
 
    virtual std::vector<HaloRegion> getHaloExportForceImportRegions(HaloRegion& initialRegion, double cutoffRadius,
                                                                    bool coversWholeDomain[3],
                                                                    double cellLength[3]) = 0;
 
    virtual std::vector<HaloRegion> getLeavingExportRegions(HaloRegion& initialRegion, double cutoffRadius,
                                                            bool coversWholeDomain[3]);
 
    virtual std::vector<HaloRegion> getLeavingExportRegions(HaloRegion& initialRegion, double cutoffRadius[3],
                                                            bool coversWholeDomain[3]);
 
protected:
    std::vector<HaloRegion> getHaloRegionsConditional(HaloRegion& initialRegion, double cutoffRadius,
                                                      bool coversWholeDomain[3],
                                                      const std::function<bool(const int[3])>& condition);
 
    std::vector<HaloRegion> getHaloRegionsConditional(HaloRegion& initialRegion, const double cutoffRadius[3],
                                                      bool coversWholeDomain[3],
                                                      const std::function<bool(const int[3])>& condition);
 
    std::vector<HaloRegion> getHaloRegionsConditionalInside(HaloRegion& initialRegion, double cutoffRadius,
                                                            bool coversWholeDomain[3],
                                                            const std::function<bool(const int[3])>& condition);
 
    std::vector<HaloRegion> getHaloRegionsConditionalInside(HaloRegion& initialRegion, const double cutoffRadius[3],
                                                            bool coversWholeDomain[3],
                                                            const std::function<bool(const int[3])>& condition);
};