---

GnuPlotterTimeSeriesSharedBFactorObserver.cpp

Go to the documentation of this file.

// GnuPlotterTimeSeriesSharedBFactorObserver.cpp: implementation of the GnuPlotterTimeSeriesSharedBFactorObserver class.
//

#include "GnuPlotterTimeSeriesSharedBFactorObserver.h"

// Construction/Destruction

GnuPlotterTimeSeriesSharedBFactorObserver::GnuPlotterTimeSeriesSharedBFactorObserver()
{

}

GnuPlotterTimeSeriesSharedBFactorObserver::~GnuPlotterTimeSeriesSharedBFactorObserver()
{

}

// default behavior is to plot things in terms of (arbitrary) experimental
// units

// Evidently Gnuplot has a problem plotting more than 10 datasets and segfaults instead.  -JJW 11/01/03

void GnuPlotterTimeSeriesSharedBFactorObserver::Update(Subject *theChangedSubject)
{
        // static cast may be dangerous
        SharedBFactorNetworkMinimizable *nM = (SharedBFactorNetworkMinimizable *)theChangedSubject;

        // set the plot title
        char titlestring[100];
        char commandstring[1000];
        sprintf(titlestring,"set title \"%s\"\n",title);
        SendPlotCommand(titlestring);
        SendPlotCommand("set grid\n");
        SendPlotCommand("set xlabel \"time (min)\"\n");
        //SendPlotCommand("set ylabel \"[Chemical] (molecules)\"\n");
        SendPlotCommand("set ylabel \"Concentration/Activity (arb. units)\"\n");

        // write the experimental data to the temporary file
        freopen(buffername,"w",buffer);

        std::vector< std::vector<int> *> whichSimToPrint;
        // dump blocks of data/simulation to files, separated by double newlines
        // EXPERIMENT
        int nExp = nM->GetNExp();
        whichSimToPrint.reserve(nExp);

        std::vector<DataPoint *> nullVector;
        std::vector<int > nChem;
        std::vector<int > nTimeSeries;

        for(int expCounter=0; expCounter < nExp; expCounter++) {
                whichSimToPrint.push_back(new std::vector<int >);
                nChem.push_back(nM->GetExperiment(expCounter)->GetChemicalTimeSeriesData()->GetNChemicals());
                nTimeSeries.push_back(nM->GetExperiment(expCounter)->GetChemicalTimeSeriesData()->GetNTimeSeries());

                for(int i = 0; i < nChem[expCounter]; i++)
                {
                        if( (nM->GetExperiment(expCounter)->GetChemicalTimeSeriesData()->GetTimeSeries(i)) != nullVector )
                        {
                                whichSimToPrint[expCounter]->push_back(i);
                                for(int j = 0; j < (nM->GetExperiment(expCounter)->GetChemicalTimeSeriesData()->GetTimeSeries(i)).size(); j++)
                                {
                                        double time = nM->GetExperiment(expCounter)->GetChemicalTimeSeriesData()->GetTimeSeries(i)[j]->GetTime();
                                        double data = nM->GetExperiment(expCounter)->GetChemicalTimeSeriesData()->GetTimeSeries(i)[j]->GetDatum();
                                        double error = nM->GetExperiment(expCounter)->GetChemicalTimeSeriesData()->GetTimeSeries(i)[j]->GetError();
                                        fprintf(buffer,"%e\t%e\t%e\n",time,data,error);
                                }
                                fprintf(buffer,"\n\n");
                        }
                }

                // SIMULATION
                int maxTime = nM->GetCellObserver(expCounter)->GetNumberOfTimeSteps();
                for(int i = 0; i < whichSimToPrint[expCounter]->size(); i++)
                {
                        for(int j = 0; j < (maxTime-1); j++)
                        {
                                int target = (*whichSimToPrint[expCounter])[i];
                                double time = double (j);
                                double sim = nM->GetCellObserver(expCounter)->GetAverageConcentration()[target][j];
                                if(nM->GetConversionFactor(target)->IsFactorNeeded())
                                {
                                        double scale = nM->GetConversionFactor(target)->GetFactorValue();
                                        sim = sim*scale;
                                }
                                fprintf(buffer,"%e\t%e\n",time,sim);
                        }
                        fprintf(buffer,"\n\n");
                }
        }

        // now plot everything
        // first block (done this way to keep the commas straight)
        int index=0;
        char namestring[100];
        sprintf(namestring,nM->GetExperiment(0)->GetReactionNetwork()->GetChemical((*whichSimToPrint[0])[0])->GetName().c_str());
        sprintf(commandstring,"plot \"%s\" index %i t \"%s\" w errorbars lw 2 ps 2",buffername,index,namestring);
        index++;

        //Rest of experimental data for experiment 0
        for(int i = 1; i < nTimeSeries[0]; i++)
        {
                char tempbuffer[200];
                sprintf(namestring,nM->GetExperiment(0)->GetReactionNetwork()->GetChemical((*whichSimToPrint[0])[i])->GetName().c_str());
                sprintf(tempbuffer,",\"%s\" index %i t \"%s\" w errorbars lw 2 ps 2",buffername,index,namestring);
                strcat(commandstring,tempbuffer);
                index++;
        }

        //Simulation for experiment 0
        for(int i = 0; i < nTimeSeries[0]; i++)
        {
                char tempbuffer[200];
                sprintf(tempbuffer,",\"%s\" index %i notitle w lines lt %i lw 4",buffername,index,index-nTimeSeries[0]+1);
                strcat(commandstring,tempbuffer);
                index++;
        }

        for (int expCounter=1;expCounter < nExp; expCounter++) {
                // the rest - experiment first
                for(int i = 0; i < nTimeSeries[expCounter]; i++)
                {
                        char tempbuffer[200];
                        sprintf(namestring,nM->GetExperiment(expCounter)->GetReactionNetwork()->GetChemical((*whichSimToPrint[expCounter])[i])->GetName().c_str());
                        sprintf(tempbuffer,",\"%s\" index %i t \"%s\" w errorbars lw 2 ps 2",buffername,index,namestring);
                        strcat(commandstring,tempbuffer);
                        index++;
                }

                for(int i = 0; i < nTimeSeries[expCounter]; i++)
                {
                        char tempbuffer[200];
                        sprintf(tempbuffer,",\"%s\" index %i notitle w lines lt %i lw 4",buffername,index,index-nTimeSeries[expCounter]+1);
                        strcat(commandstring,tempbuffer);
                        index++;
                }

        }

        strcat(commandstring,"\n");
        // change single backslashes to forward slashes to keep windows happy
#ifdef _WIN32
        int iterator = 0;
        while(commandstring[iterator] != '\0')
        {
                if(commandstring[iterator] == '\\') {commandstring[iterator] = '/';}
                iterator++;
        }
#endif

        fflush(buffer);
        SendPlotCommand(commandstring);

        return;
}

---