interface_communicator_mpi.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Philipp Bucher, Jordi Cotela
//
// See Master-Thesis P.Bucher
// "Development and Implementation of a Parallel
//  Framework for Non-Matching Grid Mapping"
 
#pragma once
 
// System includes
 
// External includes
 
// Project includes
#include "custom_searching/interface_communicator.h"
 
namespace Kratos
{
 
 
 
class KRATOS_API(MAPPING_APPLICATION) InterfaceCommunicatorMPI : public InterfaceCommunicator
{
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(InterfaceCommunicatorMPI);
 
    typedef std::vector<std::vector<double>> BufferTypeDouble;
    typedef std::vector<std::vector<char>> BufferTypeChar;
 
 
    InterfaceCommunicatorMPI(ModelPart& rModelPartOrigin,
                             MapperLocalSystemPointerVector& rMapperLocalSystems,
                             Parameters SearchSettings);
 
    virtual ~InterfaceCommunicatorMPI()
    {
    }
 
 
    std::string Info() const override
    {
        std::stringstream buffer;
        buffer << "InterfaceCommunicatorMPI" ;
        return buffer.str();
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "InterfaceCommunicatorMPI";
    }
 
    void PrintData(std::ostream& rOStream) const override {}
 
 
protected:
 
    void InitializeSearch(const MapperInterfaceInfoUniquePointerType& rpRefInterfaceInfo) override;
 
    // This function constructs the InterfaceObjects on the Destination
    // In serial it only does it once, whereas in MPI this involves Data-Exchange!
    // Imagine a sliding interface, there the partitions might change!
    void InitializeSearchIteration(const MapperInterfaceInfoUniquePointerType& rpRefInterfaceInfo) override;
 
    void FinalizeSearchIteration(const MapperInterfaceInfoUniquePointerType& rpRefInterfaceInfo) override;
 
 
private:
 
    std::vector<double> mGlobalBoundingBoxes;
    // xmax, xmin,  ymax, ymin,  zmax, zmin
 
    int mCommRank;
    int mCommSize;
 
    std::vector<int> mSendSizes;
    std::vector<int> mRecvSizes;
 
    BufferTypeDouble mSendBufferDouble;
    BufferTypeDouble mRecvBufferDouble;
 
    BufferTypeChar mSendBufferChar;
    BufferTypeChar mRecvBufferChar;
 
 
    std::size_t GetBufferSizeEstimate() const
    {
        return mrMapperLocalSystems.size() / mCommSize;
    }
 
    void ComputeGlobalBoundingBoxes();
 
    template< typename TDataType >
    int ExchangeDataAsync(
        const std::vector<std::vector<TDataType>>& rSendBuffer,
        std::vector<std::vector<TDataType>>& rRecvBuffer);
 
 
}; // Class InterfaceCommunicatorMPI
 
 
 
}  // namespace Kratos.