LateralLongRangeOrientationConnectionPredicate.cpp

Go to the documentation of this file.

#include "LateralLongRangeOrientationConnectionPredicate.h"

#include <iostream>
using std::cout;
using std::cerr;
using std::endl;

LateralLongRangeOrientationConnectionPredicate::LateralLongRangeOrientationConnectionPredicate(
    double C,
    double sigma21,
    double kappa2,
    double mu2,
    bool toroid,
    double toroid_off,
    int src_shape_x, int src_shape_y,
    int dst_shape_x, int dst_shape_y,
    double src_scale,
    double dst_scale,
    double src_center_x, double src_center_y,
    double dst_center_x, double dst_center_y
    )
    : m_C(C)
    , m_sigma21(sigma21)
    , m_kappa2(kappa2)
    , m_mu2(mu2)
    , m_distance(toroid, toroid_off,
                    src_shape_x, src_shape_y,
                    dst_shape_x, dst_shape_y,
                    src_scale, dst_scale,
                    Point2D<double>(src_center_x, src_center_y),
                    Point2D<double>(dst_center_x, dst_center_y))
    , m_unirnd(0.0, 1.0)
{
}

void LateralLongRangeOrientationConnectionPredicate::init(SimObjectPopulation const& src, SimObjectPopulation const& dst, RandomEngine *rnd)
{
    AugmentedConnectionDecisionPredicate::init(src, dst, rnd);
}

bool LateralLongRangeOrientationConnectionPredicate::decide(size_t src, size_t dst, RandomEngine *rnd)
{
    if (m_sourcePopulation->getID(src) == m_destinationPopulation->getID(dst))
        return false;

    Point2D<double> s(m_sourcePopulation->getLocation(src));
    Point2D<double> d(m_destinationPopulation->getLocation(dst));

    const SimObjectAttributes & src_attr = m_sourcePopulation->getAttributes(src);
    const SimObjectAttributes & dst_attr = m_destinationPopulation->getAttributes(dst);
    double sori=src_attr.get<double>(string("orientation"));
    double dori=dst_attr.get<double>(string("orientation"));

    double vdist2 = m_distance.sqr_distance(s, d);

    //mises distribution
    double phi = sori - dori - m_mu2;
    double V2 = exp(m_kappa2 * cos(2*phi))/exp(m_kappa2);

    //Gaussian part of ori-component
    double G21 = exp(-vdist2/(2.0 * m_sigma21*m_sigma21));

    return (m_unirnd(*rnd) < (V2 * G21 * m_C));
}

---