RetinotopicOrientationConnectionPredicate.cpp

Go to the documentation of this file.

#include "RetinotopicOrientationConnectionPredicate.h"

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


RetinotopicOrientationConnectionPredicate::RetinotopicOrientationConnectionPredicate(
    double C,
    int max_input_syn,
    int src_shape_x, int src_shape_y,
    double src_scale,
    RandomDistribution &nsub,
    RandomDistribution &phi,
    RandomDistribution &f,
    double aspectratio_min,
    double aspectratio_max)
    : m_C(C)
    , m_max_input_syn(max_input_syn)
    , m_distance(false, 0,
                 src_shape_x, src_shape_y,
                 1, 1,
                 src_scale,
                 1.0,
                 Point2D<double>((src_shape_x-1)/2.0, (src_shape_y-1)/2.0),
                 Point2D<double>(0.0, 0.0))
    , m_nsub(nsub.clone())
    , m_phi(phi.clone())
    , m_f(f.clone())
    , m_aspectratio_min(aspectratio_min)
    , m_aspectratio_max(aspectratio_max)
    , m_unirnd(0,1)
{
}


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

    size_t nidx, sidx;
    double pi = acos(0.0) * 2;             // found no pi constant in cmath

    for(nidx=0; nidx < m_destinationPopulation->size(); nidx++)
    {
        const SimObjectAttributes & dst_attr = m_destinationPopulation->getAttributes(nidx);
        Point2D<double> d(m_destinationPopulation->getLocation(nidx));

        Point2D<double> dfp(dst_attr.get<double>(string("feature_pos_x")),
                            dst_attr.get<double>(string("feature_pos_y")));

        double theta = dst_attr.get<double>(string("orientation"));
        double preference = dst_attr.get<double>(string("preference"));
        double aspectratio = preference*(m_aspectratio_max-m_aspectratio_min) + m_aspectratio_min;

        double nsub = (*m_nsub)(*rnd);
        double phi = (*m_phi)(*rnd);
        double f = (*m_f)(*rnd);

        // see Troyer 98 for details
        double lambda_x = nsub/2.448/4.0/f;
        double lambda_y = aspectratio/2.448/4.0/f;

        vector<double> gabor(m_sourcePopulation->size());
        double z_min=0;

        for(sidx=0; sidx < m_sourcePopulation->size(); sidx++)
        {               
            Point3D<double> spos(m_sourcePopulation->getLocation(sidx));

            if(sidx==0)
                z_min=spos.z();
            z_min=min(z_min, spos.z());

            Point2D<double> s(spos);
            Point2D<double> vpos = m_distance.diff(s, dfp);

/*
            if (sidx==0 && nidx==0) {
                    cout<<"s:"<<s.x()<<","<<s.y()<< endl;       
                    cout<<"dfp:"<<dfp.x()<<","<<dfp.y()<< endl; 
                    cout<<"vpos:"<<vpos.x()<<","<<vpos.y()<< endl;
                    cout<<"theta:"<<theta<< endl;
                }
*/          

            double Phi = 2.0*pi*f * (vpos.x()*cos(theta) + vpos.y()*sin(theta)) + phi;

            vpos.rotate(theta);
            double edist2 = vpos.x()*vpos.x()/(lambda_x*lambda_x) + vpos.y()*vpos.y()/(lambda_y*lambda_y);
            gabor[sidx] = exp(-edist2/2.0) * cos(Phi);

/*            
            if (sidx==0 && nidx==0) {
                    cout<<"vpos_rotated:"<<vpos.x()<<","<<vpos.y()<< endl;
                    cout<<"edist2:"<<edist2<< endl;
                    cout<<"preference:"<<preference<< endl;
                    cout<<"aspectratio:"<<aspectratio<< endl;
                    cout<<"aspectratio_min:"<<m_aspectratio_min<< endl;
                    cout<<"aspectratio_max:"<<m_aspectratio_max<< endl;
                    cout<<"lambda_x:"<<lambda_x<<", lambda_y:"<<lambda_y<< endl;
                    cout<<"f:"<<f<<", phi:"<<phi<< endl;
                    cout<<"Phi:"<<Phi<< endl;
                    cout<<"gabor:"<<gabor[sidx]<< endl;
                }
*/            
        }

        vector<double> prob(m_sourcePopulation->size());
        double p_sum=0.0;
        for(sidx=0; sidx < m_sourcePopulation->size(); sidx++)
        {
            Point3D<double> spos(m_sourcePopulation->getLocation(sidx));
            if (spos.z() == z_min || spos.z() == z_min+2)                    // ON cell
                prob[sidx]= (gabor[sidx] > 0.0) ? gabor[sidx] : 0.0;
            else                                                         // OFF cell
                prob[sidx]= (gabor[sidx] < 0.0) ? fabs(gabor[sidx]) : 0.0;

            p_sum += prob[sidx];
/*            
            if (sidx==0 && nidx==0) {
                    cout<<"pos.z:"<<spos.z()<< endl;    
                    cout<<"prob:"<<prob[sidx]<< endl;   
                }
*/
        }
        
        for(sidx=0; sidx < m_sourcePopulation->size(); sidx++)
        {
/*
                double real_prob=m_C*prob[sidx]/p_sum*m_max_input_syn;
            if (sidx==0 && nidx==0) {
                cout<<"p.sum():"<<p_sum<<endl;
                    cout<<"real_prob:"<<real_prob<< endl;
                }
*/              
                
            if (m_unirnd(*rnd) < m_C*prob[sidx]/p_sum*m_max_input_syn)
                m_connections.insert(pair<size_t, size_t>(sidx, nidx));
        }
    }
}

---