---

dlogger.cpp

Go to the documentation of this file.

#include "dlogger.h"
#include <iomanip>

using namespace std;

namespace punnets_common {

const int lev_shift = 0;

void tlogger::activate(tscheduler &scheduler, ntime_t current_time)
{ 
    out << setprecision(15) << current_time;

    int lev = 0;
    
    for( vector< neuentry >::iterator i = neus.begin(); i != neus.end(); i++ )
    {
        if( i->delay != 0.0 )
            out << "\t" << setprecision(15) << current_time - i->delay;
        out << "\t" << setprecision(15) << 
            (i->neuron->getLastFire() >= current_time - step ? 
                i->lastthr + 1.0 : 
                i->neuron->getCurrentSigLevel(current_time) ) + lev + i->offset;
        i->lastthr = i->neuron->getCurrentThrLevel(current_time);
        if( i->logopt & showthr )
            out << "\t" << setprecision(15) << i->lastthr + lev + i->offset;
        if( i->logopt & showext )
            out << "\t" << setprecision(15) << i->neuron->getCurrentExtInput(current_time) + lev + i->offset;
        if( i->logopt & showpart )
        {
            if( i->neuron->getLastSimulate() >= current_time - step )
            {
                switch( i->neuron->getLastSimulateType() )
                {
                    case 0:
                    default:
                        out << " -10 = = ="; break;
                    case 1:
                        out << " = -10 = ="; break;
                    case 2:
                        out << " = = -10 ="; break;
                    case 3:
                        out << " = = = -10"; break;
                }
            }
            else
                out << " = = = =";
        }
        lev += lev_shift;
    }
    for( vector<pair<tsynapse_base *, ntime_t> >::iterator i = syns.begin(); i != syns.end(); i++ )
    {
        if( i->second != 0.0 )
            out << "\t" << setprecision(15) << current_time - i->second;
        out << "\t" << setprecision(15) << i->first->getWeight();// + lev;
//      lev += lev_shift;
    }
    out << endl;
    if( current_time + step < until + delaymax )
        scheduler.scheduleEvent(current_time + step, *this);
};

void tlogger::gnuplot_def(ostream &os, string logfile)
{
//  cerr << "min/max: " << delaymin << "/" << delaymax << endl;
//  os << "set term X11" << endl;
    os << "set y2tics" << endl;
    os << "plot [" << from << ":" << until << "]";

    string delim = "";
    int column = 1;
    for( vector< neuentry >::iterator i = neus.begin(); i != neus.end(); i++ )
    {
        int xcol = (i->delay == 0.0 ? 1 : ++column);
        os << delim << " \\" << endl << "\t'" << logfile << "' " << 
            "using " << xcol << ":" << (++column) << " axes x1y1 " << 
            "title \"" << i->neuron->getClassName() << " '" << i->neuron->getName() << "'";
        if( i->delay != 0.0 )
            os << "(delay " << i->delay << ")";
        if( i->offset != 0.0 )
            os << "(offset " << i->offset << ")";
        os << " signal\" with lines";
        delim = ",";

        if( i->logopt & showthr )
        {
            os << delim << " \\" << endl << "\t'" << logfile << "' " << 
                "using " << xcol << ":" << (++column) << " axes x1y1 " << 
                "title \"" << i->neuron->getClassName() << " '" << i->neuron->getName() << "'";
            if( i->delay != 0.0 )
                os << "(delay " << i->delay << ")";
            if( i->offset != 0.0 )
                os << "(offset " << i->offset << ")";
            os << " threshold\" with dots";
            delim = ",";
        }
        if( i->logopt & showext )
        {
            os << delim << " \\" << endl << "\t'" << logfile << "' " << 
                "using " << xcol << ":" << (++column) << " axes x1y1 " << 
                "title \"" << i->neuron->getClassName() << " '" << i->neuron->getName() << "'";
            if( i->delay != 0.0 )
                os << "(delay " << i->delay << ")";
            if( i->offset != 0.0 )
                os << "(offset " << i->offset << ")";
            os << " external\" with dots";
            delim = ",";
        }
        if( i->logopt & showpart )
        {
            os << delim << " \\" << endl << "\t'" << logfile << "' " << 
                "using " << xcol << ":" << (++column) << " axes x1y1 " << 
                "title \"" << i->neuron->getClassName() << " '" << i->neuron->getName() << "'";
            if( i->delay != 0.0 )
                os << "(delay " << i->delay << ")";
            if( i->offset != 0.0 )
                os << "(offset " << i->offset << ")";
            os << " partition 0\" with impulses";
            delim = ",";

            os << delim << " \\" << endl << "\t'" << logfile << "' " << 
                "using " << xcol << ":" << (++column) << " axes x1y1 " << 
                "title \"" << i->neuron->getClassName() << " '" << i->neuron->getName() << "'";
            if( i->delay != 0.0 )
                os << "(delay " << i->delay << ")";
            if( i->offset != 0.0 )
                os << "(offset " << i->offset << ")";
            os << " partition 1\" with impulses";

            os << delim << " \\" << endl << "\t'" << logfile << "' " << 
                "using " << xcol << ":" << (++column) << " axes x1y1 " << 
                "title \"" << i->neuron->getClassName() << " '" << i->neuron->getName() << "'";
            if( i->delay != 0.0 )
                os << "(delay " << i->delay << ")";
            if( i->offset != 0.0 )
                os << "(offset " << i->offset << ")";
            os << " partition 2\" with impulses";

            os << delim << " \\" << endl << "\t'" << logfile << "' " << 
                "using " << xcol << ":" << (++column) << " axes x1y1 " << 
                "title \"" << i->neuron->getClassName() << " '" << i->neuron->getName() << "'";
            if( i->delay != 0.0 )
                os << "(delay " << i->delay << ")";
            if( i->offset != 0.0 )
                os << "(offset " << i->offset << ")";
            os << " partition delta_t\" with impulses";
        }
    }
    for( vector<pair<tsynapse_base *, ntime_t> >::iterator i = syns.begin(); i != syns.end(); i++ )
    {
        int xcol = (i->second == 0.0 ? 1 : ++column);
        os << delim << " \\" << endl << "\t'" << logfile << "' " << 
            "using " << xcol << ":" << (++column) << " axes x1y2 " <<
            "title \"" << i->first->getClassName() << " '" << i->first->getSrc().getName() << "' -> '" << i->first->getDest().getName() << "'";
        if( i->second != 0.0 )
            os << "(delay " << i->second << ")";
        os << "\" with lines";
        delim = ",";
    }

    os << endl;
    os << "pause -1" << endl;
}


} // namespace punnets_common

---