pointcloud.cpp

/*---------------------------------------------------------------------------*\
 *
 *  bitpit
 *
 *  Copyright (C) 2015-2021 OPTIMAD engineering Srl
 *
 *  -------------------------------------------------------------------------
 *  License
 *  This file is part of bitpit.
 *
 *  bitpit is free software: you can redistribute it and/or modify it
 *  under the terms of the GNU Lesser General Public License v3 (LGPL)
 *  as published by the Free Software Foundation.
 *
 *  bitpit is distributed in the hope that it will be useful, but WITHOUT
 *  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
 *  FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
 *  License for more details.
 *
 *  You should have received a copy of the GNU Lesser General Public License
 *  along with bitpit. If not, see <http://www.gnu.org/licenses/>.
 *
\*---------------------------------------------------------------------------*/
 
#include "bitpit_common.hpp"
 
#include "pointcloud.hpp"
 
namespace bitpit {

 
#if BITPIT_ENABLE_MPI==1
PointCloud::PointCloud(MPI_Comm communicator)
    : PointKernel(communicator, ADAPTION_MANUAL, PARTITIONING_DISABLED)
#else
PointCloud::PointCloud()
    : PointKernel(ADAPTION_MANUAL)
#endif
{
}
 
#if BITPIT_ENABLE_MPI==1
PointCloud::PointCloud(int dimension, MPI_Comm communicator)
    : PointKernel(PatchManager::AUTOMATIC_ID, dimension, communicator, ADAPTION_MANUAL, PARTITIONING_DISABLED)
#else
PointCloud::PointCloud(int dimension)
    : PointKernel(PatchManager::AUTOMATIC_ID, dimension, ADAPTION_MANUAL)
#endif
{
}
 
#if BITPIT_ENABLE_MPI==1
PointCloud::PointCloud(int id, int dimension, MPI_Comm communicator)
    : PointKernel(id, dimension, communicator, ADAPTION_MANUAL, PARTITIONING_DISABLED)
#else
PointCloud::PointCloud(int id, int dimension)
    : PointKernel(id, dimension, ADAPTION_MANUAL)
#endif
{
}
 
#if BITPIT_ENABLE_MPI==1
PointCloud::PointCloud(std::istream &stream, MPI_Comm communicator)
    : PointKernel(communicator, ADAPTION_MANUAL, PARTITIONING_DISABLED)
#else
PointCloud::PointCloud(std::istream &stream)
    : PointKernel(ADAPTION_MANUAL)
#endif
{
    // Restore the patch
    restore(stream);
}

std::unique_ptr<PatchKernel> PointCloud::clone() const
{
    return std::unique_ptr<PointCloud>(new PointCloud(*this));
}

int PointCloud::_getDumpVersion() const
{
    const int DUMP_VERSION = 1;
 
    return DUMP_VERSION;
}

void PointCloud::_dump(std::ostream &stream) const
{
#if BITPIT_ENABLE_MPI==1
    // Dump works only for serial calculations
    if (getProcessorCount() != 1) {
        throw std::runtime_error ("Dump of lineunstructured is implemented only for serial calculations.");
    }
#endif
 
    // Save the vertices
    dumpVertices(stream);
 
    // Save the cells
    dumpCells(stream);
 
    // Save the interfaces
    dumpInterfaces(stream);
}

void PointCloud::_restore(std::istream &stream)
{
#if BITPIT_ENABLE_MPI==1
    // Restore works only for serial calculations
    if (getProcessorCount() != 1) {
        throw std::runtime_error ("Restore of lineunstructured is implemented only for serial calculations.");
    }
#endif
 
    // Restore the vertices
    restoreVertices(stream);
 
    // Restore the cells
    restoreCells(stream);
 
    // Restore the interfaces
    restoreInterfaces(stream);
}

long PointCloud::locatePoint(const std::array<double, 3> &point) const
{
    BITPIT_UNUSED(point);
 
    throw std::runtime_error ("The function 'locatePoint' is not implemented yet");
 
    return false;
}
 
}