HalfShell.h

/*
 * CommunicationScheme.h
 *
 *  Created on: 14.07.2017
 *      Author: sauermann
 */
 
#pragma once
 
#include "ZonalMethod.h"
 
class HalfShell: public ZonalMethod {
public:
    HalfShell() = default;
    ~~HalfShell() override = default;
 
    std::vector<HaloRegion> getHaloImportForceExportRegions(HaloRegion& initialRegion, double cutoffRadius,
                                                            bool coversWholeDomain[3], double cellLength[3]) override {
        auto condition = [](const int d[3])->bool {
            // Here, we cannot directly apply the stencil from the HS traversal, as for multiple cells also the
            // other directions are needed. Visualize: Apply the stencil in the z-direction for the lowest cell.
            return d[2] >= 0;
        };
        return getHaloRegionsConditional(initialRegion, cutoffRadius, coversWholeDomain, condition);
    }
 
    std::vector<HaloRegion> getHaloExportForceImportRegions(HaloRegion& initialRegion, double cutoffRadius,
                                                            bool coversWholeDomain[3], double cellLength[3]) override {
        auto condition = [](const int d[3])->bool {
            // Here, we cannot directly apply the stencil from the HS traversal, as for multiple cells also the
            // other directions are needed. Visualize: Apply the stencil in the z-direction for the lowest cell.
            return d[2] <= 0;
        };
        return getHaloRegionsConditionalInside(initialRegion, cutoffRadius, coversWholeDomain, condition);
    }
};