/* $Id: plstripc.c,v 1.11 2005/04/27 06:44:25 rlaboiss Exp $
*
* Plots a simple stripchart.
*
* Copyright (C) 2004 Alan W. Irwin
*
* This file is part of PLplot.
*
* PLplot is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Library Public License as published
* by the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* PLplot is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU Library General Public License for more details.
*
* You should have received a copy of the GNU Library General Public License
* along with PLplot; if not, write to the Free Software
* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
*
*
* ToDo: better way of clearing plot. search for `plvsta'.
*/
#include "plplotP.h"
/* Data declarations for stripcharts. */
#define PEN 4
typedef struct {
PLFLT xmin, xmax, ymin, ymax, xjump, xlen;
PLFLT wxmin, wxmax, wymin, wymax; /* FIXME - some redundancy might exist */
char *xspec, *yspec, *labx, *laby, *labtop;
PLINT y_ascl, acc, colbox, collab;
PLFLT xlpos, ylpos;
PLFLT *x[PEN], *y[PEN];
PLINT npts[PEN], nptsmax[PEN];
PLINT colline[PEN], styline[PEN];
char *legline[PEN];
} PLStrip;
static int sid; /* strip id number */
#define MAX_STRIPC 1000 /* Max allowed */
static PLStrip *strip[MAX_STRIPC]; /* Array of pointers */
static PLStrip *stripc; /* current strip chart */
/* Generates a complete stripchart plot. */
static void
plstrip_gen(PLStrip *strip);
/* draw legend */
static void
plstrip_legend(PLStrip *strip, int flag);
/*--------------------------------------------------------------------------*\
* plstripc
*
* Create 1d stripchart.
\*--------------------------------------------------------------------------*/
void
c_plstripc( PLINT *id, char *xspec, char *yspec,
PLFLT xmin, PLFLT xmax, PLFLT xjump, PLFLT ymin, PLFLT ymax,
PLFLT xlpos, PLFLT ylpos,
PLINT y_ascl, PLINT acc,
PLINT colbox, PLINT collab,
PLINT *colline, PLINT *styline, char *legline[],
char *labx, char *laby, char *labtop )
{
int i;
if (!plsc->termin) {
plabort("Stripchart plot requires a terminal output device");
return;
}
/* Get a free strip id and allocate it */
for (i = 0; i < MAX_STRIPC; i++)
if (strip[i] == NULL)
break;
if (i == MAX_STRIPC) {
plabort("plstripc: Cannot create new strip chart");
*id = -1;
return;
}
else {
sid = *id = i;
strip[sid] = (PLStrip *) calloc(1, (size_t) sizeof(PLStrip));
if (strip[sid] == NULL) {
plabort("plstripc: Out of memory.");
*id = -1;
return;
}
}
/* Fill up the struct with all relevant info */
stripc = strip[sid];
for (i=0; i<PEN; i++) {
stripc->npts[i] = 0;
stripc->nptsmax[i] = 100;
stripc->colline[i] = colline[i];
stripc->styline[i] = styline[i];
stripc->legline[i] = plstrdup(legline[i]);
stripc->x[i] = (PLFLT *) malloc((size_t) sizeof(PLFLT) * stripc->nptsmax[i]);;
stripc->y[i] = (PLFLT *) malloc((size_t) sizeof(PLFLT) * stripc->nptsmax[i]);;
if (stripc->x[i] == NULL || stripc->y[i] == NULL) {
plabort("plstripc: Out of memory.");
plstripd(sid);
*id = -1;
return;
}
}
stripc->xlpos = xlpos; /* legend position [0..1] */
stripc->ylpos = ylpos;
stripc->xmin = xmin; /* initial bounding box */
stripc->xmax = xmax;
stripc->ymin = ymin;
stripc->ymax = ymax;
stripc->xjump = xjump; /* jump x step(%) when x attains xmax (xmax is then set to xmax+xjump) */
stripc->xlen = xmax - xmin; /* length of x scale */
stripc->y_ascl = y_ascl; /* autoscale y between x jump scale */
stripc->acc = acc; /* accumulate plot (not really stripchart) */
stripc->xspec = plstrdup(xspec); /* x axis specification */
stripc->yspec = plstrdup(yspec);
stripc->labx = plstrdup(labx); /* x label */
stripc->laby = plstrdup(laby);
stripc->labtop = plstrdup(labtop); /* title */
stripc->colbox = colbox; /* box color */
stripc->collab = collab; /* label color */
/* Generate the plot */
plstrip_gen(stripc);
plstrip_legend(stripc,1);
}
static void plstrip_legend(PLStrip *stripc, int first)
{
int i;
PLFLT sc, dy;
/* draw legend */
plgchr(&sc, &dy);
sc = dy = dy/100;
plwind(-0.01, 1.01, -0.01, 1.01);
for (i=0; i<PEN; i++) {
if (stripc->npts[i] || first) {
plcol(stripc->colline[i]); pllsty(stripc->styline[i]);
pljoin(stripc->xlpos, stripc->ylpos - sc, stripc->xlpos + 0.1, stripc->ylpos - sc);
plcol(stripc->collab);
plptex(stripc->xlpos + 0.11, stripc->ylpos - sc, 0., 0., 0, stripc->legline[i]);sc += dy;
}
}
plwind(stripc->xmin, stripc->xmax, stripc->ymin, stripc->ymax);
plflush();
}
/*--------------------------------------------------------------------------*\
* plstrip_gen
*
* Generates a complete stripchart plot. Used either initially or
* during rescaling.
\*--------------------------------------------------------------------------*/
PLFLT oxm,oxM, oym,oyM;
static void plstrip_gen( PLStrip *strip )
{
int i;
PLFLT x[]={0.,1.,1.,0.}, y[]={0.,0.,1.,1.};
/* Set up window */
plvpor(0,1,0,1);
plwind(0,1,0,1);
plcol(0);plpsty(0);
plfill(4, &x[0], &y[0]);
plvsta();
/* Draw box and same window dimensions */
strip->wxmin=strip->xmin; strip->wxmax=strip->xmax;
strip->wymin=strip->ymin; strip->wymax=strip->ymax; /* FIXME - can exist some redundancy here */
plwind(strip->xmin, strip->xmax, strip->ymin, strip->ymax);
pllsty(1);
plcol(strip->colbox);
plbox(strip->xspec, 0.0, 0, strip->yspec, 0.0, 0);
plcol(strip->collab);
pllab(strip->labx, strip->laby, strip->labtop);
for (i=0; i<PEN; i++) {
if (strip->npts[i] > 0) {
plcol(strip->colline[i]);pllsty(strip->styline[i]);
plline(strip->npts[i], strip->x[i], strip->y[i]);
}
}
plstrip_legend(strip,0);
}
/*--------------------------------------------------------------------------*\
* plstripa
*
* Add a point to a stripchart.
* Allocates memory and rescales as necessary.
\*--------------------------------------------------------------------------*/
void c_plstripa( PLINT id, PLINT p, PLFLT x, PLFLT y )
{
int j, yasc=0, istart;
if (p >= PEN) {
plabort("Non existent pen");
return;
}
if ((id < 0) || (id >= MAX_STRIPC) ||
((stripc = strip[id]) == NULL)) {
plabort("Non existent stripchart");
return;
}
/* Add new point, allocating memory if necessary */
if (++stripc->npts[p] > stripc->nptsmax[p]) {
stripc->nptsmax[p] += 32;
stripc->x[p] = (PLFLT *) realloc((void *) stripc->x[p], sizeof(PLFLT)*stripc->nptsmax[p]);
stripc->y[p] = (PLFLT *) realloc((void *) stripc->y[p], sizeof(PLFLT)*stripc->nptsmax[p]);
if (stripc->x[p] == NULL || stripc->y[p] == NULL) {
plabort("plstripc: Out of memory.");
plstripd(id);
return;
}
}
stripc->x[p][stripc->npts[p]-1] = x;
stripc->y[p][stripc->npts[p]-1] = y;
stripc->xmax = x;
if (stripc->y_ascl == 1 && (y > stripc->ymax || y < stripc->ymin))
yasc=1;
if (y > stripc->ymax)
stripc->ymax = stripc->ymin + 1.1*(y - stripc->ymin);
if (y < stripc->ymin)
stripc->ymin = stripc->ymax - 1.1*(stripc->ymax - y);
/* Now either plot new point or regenerate plot */
if (stripc->xmax - stripc->xmin < stripc->xlen) {
if( yasc == 0) {
/* If user has changed subwindow, make shure we have the correct one */
plvsta();
plwind(stripc->wxmin, stripc->wxmax, stripc->wymin, stripc->wymax); /* FIXME - can exist some redundancy here */
plcol(stripc->colline[p]); pllsty(stripc->styline[p]);
if ((stripc->npts[p]-2) < 0)
plP_movwor(stripc->x[p][stripc->npts[p]-1], stripc->y[p][stripc->npts[p]-1]);
else
plP_movwor(stripc->x[p][stripc->npts[p]-2], stripc->y[p][stripc->npts[p]-2]);
plP_drawor(stripc->x[p][stripc->npts[p]-1], stripc->y[p][stripc->npts[p]-1]);
plflush();
}
else {
stripc->xmax = stripc->xmin + stripc->xlen;
plstrip_gen(stripc);
}
}
else {
/* Regenerating plot */
if (stripc->acc == 0) {
for (j=0; j<PEN; j++) {
if (stripc->npts[j] > 0) {
istart = 0;
while (stripc->x[j][istart] < stripc->xmin + stripc->xlen*stripc->xjump)
istart++;
stripc->npts[j] = stripc->npts[j] - istart;
memcpy( &stripc->x[j][0], &stripc->x[j][istart], (stripc->npts[j])*sizeof(PLFLT));
memcpy( &stripc->y[j][0], &stripc->y[j][istart], (stripc->npts[j])*sizeof(PLFLT));
}
}
} else
stripc->xlen = stripc->xlen * (1 + stripc->xjump);
stripc->xmin = stripc->x[p][0];
stripc->xmax = stripc->xmax + stripc->xlen*stripc->xjump;
plstrip_gen(stripc);
}
}
/*--------------------------------------------------------------------------*\
* plstripd
*
* Deletes and releases memory used by a stripchart.
\*--------------------------------------------------------------------------*/
void c_plstripd( PLINT id )
{
int i;
if ((id < 0) || (id >= MAX_STRIPC) ||
((stripc = strip[id]) == NULL)) {
plabort("Non existent stripchart");
return;
}
for (i=0; i<PEN; i++) {
if (stripc->npts[i]) {
free((void *) stripc->x[i]);
free((void *) stripc->y[i]);
free(stripc->legline[i]);
}
}
free(stripc->xspec);
free(stripc->yspec);
free(stripc->labx);
free(stripc->laby);
free(stripc->labtop);
free((void *) stripc);
strip[id] = NULL;
}
syntax highlighted by Code2HTML, v. 0.9.1