/* The function plot_cfvr.c displays cardiac function and venous return
   curves by writing Gnuplot commands to the gpcommands file in the /tmp
   directory and then invoking gnuplot on the file.

   Function arguments:
         outputfilein - MATLAB array of the prefix name of the file containing             
                        the cardiac function/venous return numerics 
         covrformatin - MATLAB array indicating the desired plotting format
                       
   Function outputs: none */

#include <stdio.h>
#include <sys/types.h>
#include <dirent.h>
#include <unistd.h>
#include <string.h>
#include <matlab.h>

void mlfPlot_cfvr(const mxArray *outputfilein, const mxArray *covrformatin)
{

    /* Declaring variables. */
    char command[128], *outputfile;
    FILE *fp, *fpnew;
    double covrformat;

    /* Converting MATLAB types to C types. */
    outputfile=mxArrayToString(outputfilein);
    covrformat=mxGetScalar(covrformatin);

    /* If only the curve of the current simulation is desired for display
       or there is currently no /tmp/gpcommands file. */
    fp = fopen("/tmp/gpcommands","r");
    if (covrformat == 0 | covrformat == 2 | fp == 0)
    {
	
	fclose(fp);
	/* Delete existing gpcommands file. */
	system("rm -f /tmp/gpcommands");

	/* Creating new gpcommands file. */
	fpnew = fopen("/tmp/gpcommands","w");
	sprintf(command,"set tics out");
	fprintf(fpnew,"%s\n",command);
	sprintf(command,"set title ''");
	fprintf(fpnew,"%s\n",command);

	/* For cardiac function or venous return curve display. */
	if (covrformat == 0)
	{
	    sprintf(command,"set xlabel 'mPra [mmHg]'");
	}
	/* For cardiac output versus systemic arterial pressure curve display. */
	else
	{
	    sprintf(command,"set xlabel 'mPa [mmHg]'");
	}    
        fprintf(fpnew,"%s\n",command);

	/* For cardiac function or venous return curve display. */
	if (covrformat == 0)
	{
	    sprintf(command,"set ylabel 'mql,mqv [l/min]'");
	}
	/* For cardiac output versus systemic arterial pressure curve display. */
	else
	{
	    sprintf(command,"set ylabel 'mql [l/min]'");
	}    
	fprintf(fpnew,"%s\n",command);
	sprintf(command,"set xrange [:]");
	fprintf(fpnew,"%s\n",command);
	sprintf(command,"set yrange [:]");
	fprintf(fpnew,"%s\n",command);

	/* For cardiac function or venous return curve display. */
	if (covrformat == 0) 
	{
	    fprintf(fpnew,"plot '%s.txt' using 1:2 notitle with lines\n",outputfile);
	}
	/* For cardiac output versus systemic arterial pressure curve display. */
	else
	{
	    fprintf(fpnew,"plot '%s.txt' using 3:2 notitle with lines\n",outputfile);
	}
	fclose(fpnew);

	/* Copying gpcommands to a temporary file. */
	system("cp /tmp/gpcommands /tmp/gpcommands.orig");

	/* Adding a pause command to keep display window open until user
	   desires to close it to the gpcommands file. */
	fpnew = fopen("/tmp/gpcommands","a");
	sprintf(command,"pause -1 \"Press RETURN to dismiss the plot window: \"");
	fprintf(fpnew,"%s\n",command);
	fclose(fpnew);

	/* Invoking gnuplot. */
	system("gnuplot /tmp/gpcommands < /dev/tty");


	/* Copying temporary file -- without pause command -- to gpcommands. */
	system("mv -f /tmp/gpcommands.orig /tmp/gpcommands");

    }

    /* If multiple curves per a window are desired and a /tmp/gpcommands file 
       currently exists. */
    else
    {

	fclose(fp);

	/* Adding new Gnuplot commands to existing gpcommands file. */

	/* For cardiac function or venous return curve display. */
	if (covrformat == 1)
	{
	    sprintf(command,"replot \'%s.txt\' using 1:2 notitle with lines",outputfile);
	}
	/* For cardiac output versus systemic arterial pressure curve display. */
	else
	{
	    sprintf(command,"replot \'%s.txt\' using 3:2 notitle with lines",outputfile);
	}
	fpnew = fopen("/tmp/gpcommands","a");
	fprintf(fpnew,"%s\n",command);
	fclose(fpnew);

	/* Copying gpcommands to a temporary file. */
	system("cp /tmp/gpcommands /tmp/gpcommands.orig");

	/* Adding a pause command to keep display window open until user
	   desires to close it to the gpcommands file. */
	fpnew = fopen("/tmp/gpcommands","a");
	sprintf(command,"pause -1 \"Press RETURN to dismiss the plot window: ");
	fprintf(fpnew,"%s\n",command);
	fclose(fpnew);

	/* Invoking gnuplot. */
	system("gnuplot /tmp/gpcommands < /dev/tty");

	/* Copying temporary file -- without pause command -- to gpcommands. */
	system("mv -f /tmp/gpcommands.orig /tmp/gpcommands");

    }

    /* Freeing allocated dynamic memory and returning to rcvsim.m */
    mxFree(outputfile);
    return;

}