GeometryNode.hpp

/* Copyright (c) 2013       David Sichau <mail"at"sichau"dot"eu>
 *               2013-2015  Simon Tanaka <tanakas"at"gmx"dot"ch>
 *
 * Permission is hereby granted, free of charge, to any person obtaining a copy
 * of this software and associated documentation files (the "Software"), to deal
 * in the Software without restriction, including without limitation the rights
 * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
 * copies of the Software, and to permit persons to whom the Software is
 * furnished to do so, subject to the following conditions:
 *
 * The above copyright notice and this permission notice shall be included in
 * all copies or substantial portions of the Software.
 *
 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND EXPRESS OR
 * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
 * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
 * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
 * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
 * THE SOFTWARE.
 */
#ifndef POLYGONPOINT_HPP
#define POLYGONPOINT_HPP

#include <LbmLib/include/nodes/LagrangianPoint.hpp>
#include <LbmLib/include/geometry/Connection.hpp>
#include <LbmLib/include/GlobalSimulationParameters.hpp>
#include <UtilLib/include/Log.hpp>
#include <boost/geometry/geometries/point.hpp>
#include <array>
#include <iostream>
#include <ostream>
#include <string>

namespace LbmLib {
namespace geometry {
class Connection;
}
}

namespace LbmLib {
namespace nodes {
class PhysicalNode;

class GeometryNode : public LagrangianPoint {
 public:
    GeometryNode(
            double x,
            double y,
            unsigned int id);
    ~GeometryNode();


    void writeNode(std::ostream* stream);

    void move();

    void addPhysicalNode(
            PhysicalNode* const physicalNode,
            unsigned int pos);

    virtual std::string getType() const;

    virtual void setPos(
            double x,
            double y);

    void collectVelocity();

    void distributeForce();

    bool isUpdateNeeded() const;

    void setForce(Field<double> f);

    void addForce(Field<double> f);

    unsigned int getId() const;

    template<std::size_t K>
    void setConnection(const std::shared_ptr<LbmLib::geometry::Connection> newconnection)
    {
        assert(K==0 || K==1);
        if (this->connections_[K] != nullptr) {
            if (newconnection != nullptr) {
                std::stringstream message;
                message << std::setprecision(12) <<
                           "GeometryNode::setConnection<" << K << ">() on point ID= "<< this->getId() <<
                           " (" << this->getXPos() << "," << this->getYPos() << "): overwrite non-nullptr entry (this should not happen)";
                lbm_fail(message.str().c_str());
            }
        }
        this->connections_[K] = newconnection;
    }

    template<std::size_t K>
    const std::shared_ptr<LbmLib::geometry::Connection>  getConnection() const
    {
        if (this->connections_[K] == nullptr) {
            std::stringstream message;
            message << std::setprecision(12) <<
                       "GeometryNode::getConnection<" << K << ">() on point ID= "<< this->getId() <<
                       " (" << this->getXPos() << "," << this->getYPos() << ")" <<
                       " at time " << Parameters.getCurrentIteration() <<
                       ": nullptr found (this should not happen)";
            LOG(UtilLib::logINFO) << message.str().c_str();
        }
        assert(K==0 || K==1);
        assert(this->connections_[K] != nullptr);
        return this->connections_[K];
    }

    unsigned int getDomainIdOfAdjacentConnections(void) const;



 private:
    Field<double> oldPosition_;
    Field<double> force_;
    void perturbatePosition();

    const unsigned int id_;

    std::array<PhysicalNode*,
            16> neighbourPhysicalNodes_;

    std::array<std::shared_ptr<LbmLib::geometry::Connection>,
            2> connections_;


};
}   // end namespace
}  // end namespace

#endif  // POLYGONPOINT_HPP