Demo of multiple stream/window capability

 #pragma _ipath "./"
 #include "plcdemos.h"
 
 static PLFLT x[101], y[101];
 static PLFLT xscale, yscale, xoff, yoff, xs[6], ys[6];
 static PLINT space0 = 0, mark0 = 0, space1 = 1500, mark1 = 1500;
 
 void plot1();
 void plot2();
 void plot3();
 void plot4();
 void plot5();
 
 
 int
 main(int argc, char *argv[])
 {
     int i, digmax;
     int xleng0 = 400, yleng0 = 300, xoff0 = 200, yoff0 = 200;
     int xleng1 = 400, yleng1 = 300, xoff1 = 500, yoff1 = 500;
 
 
     char geometry_master[] = "500x410+100+200";
     char geometry_slave[]  = "500x410+650+200";
 
     char driver[80];
 
 
     (void) plParseOpts(&argc, argv, PL_PARSE_FULL);
 
     plgdev(driver);
 
     printf("Demo of multiple output streams via the %s driver.\n", driver);
     printf("Running with the second stream as slave to the first.\n");
     printf("\n");
 
 
     plsetopt("geometry", geometry_master);
 
     plsdev(driver);
     plssub(2, 2);
     plinit();
 
 
     plsstrm(1);
 
 
     plsetopt("geometry", geometry_slave);
     plspause(0);
     plsdev(driver);
     plinit();
 
 
     plsstrm(0);
 
     xscale = 6.;
     yscale = 1.;
     xoff = 0.;
     yoff = 0.;
     plot1();
 
 
     xscale = 1.;
     yscale = 1.e+6;
     plot1();
 
 
     xscale = 1.;
     yscale = 1.e-6;
     digmax = 2;
     plsyax(digmax, 0);
     plot1();
 
 
     xscale = 1.;
     yscale = 0.0014;
     yoff = 0.0185;
     digmax = 5;
     plsyax(digmax, 0);
     plot1();
 
 
     plsstrm(1);
     plot4();
     pleop();
 
 
     plsstrm(0);
     plot2();
     plot3();
 
 
     plsstrm(1);
     plot5();
     pleop();
 
 
     plsstrm(0);
     pleop();
 
 
     plend();
     exit(0);
 }
 
 
 void
 plot1(void)
 {
     int i;
     PLFLT xmin, xmax, ymin, ymax;
 
     for (i = 0; i < 60; i++) {
 	x[i] = xoff + xscale * (i + 1) / 60.0;
 	y[i] = yoff + yscale * pow(x[i], 2.);
     }
 
     xmin = x[0];
     xmax = x[59];
     ymin = y[0];
     ymax = y[59];
 
     for (i = 0; i < 6; i++) {
 	xs[i] = x[i * 10 + 3];
 	ys[i] = y[i * 10 + 3];
     }
 
 
     plcol0(1);
     plenv(xmin, xmax, ymin, ymax, 0, 0);
     plcol0(6);
     pllab("(x)", "(y)", "#frPLplot Example 1 - y=x#u2");
 
 
     plcol0(9);
     plpoin(6, xs, ys, 9);
 
 
     plcol0(4);
     plline(60, x, y);
     plflush();
 }
 
 
 
 void
 plot2(void)
 {
     int i;
 
 
     plcol0(1);
     plenv(-2.0, 10.0, -0.4, 1.2, 0, 1);
     plcol0(2);
     pllab("(x)", "sin(x)/x", "#frPLplot Example 1 - Sinc Function");
 
 
     for (i = 0; i < 100; i++) {
 	x[i] = (i - 19.0) / 6.0;
 	y[i] = 1.0;
 	if (x[i] != 0.0)
 	    y[i] = sin(x[i]) / x[i];
     }
 
 
     plcol0(3);
     plline(100, x, y);
     plflush();
 }
 
 
 void
 plot3(void)
 {
     int i;
 
 
     pladv(0);
 
 
     plvsta();
     plwind(0.0, 360.0, -1.2, 1.2);
 
 
     plcol0(1);
     plbox("bcnst", 60.0, 2, "bcnstv", 0.2, 2);
 
 
     plstyl(1, &mark1, &space1);
     plcol0(2);
     plbox("g", 30.0, 0, "g", 0.2, 0);
     plstyl(0, &mark0, &space0);
 
     plcol0(3);
     pllab("Angle (degrees)", "sine", "#frPLplot Example 1 - Sine function");
 
     for (i = 0; i < 101; i++) {
 	x[i] = 3.6 * i;
 	y[i] = sin(x[i] * PI / 180.0);
     }
 
     plcol0(4);
     plline(101, x, y);
     plflush();
 }
 
 
 void
 plot4(void)
 {
     int i, j;
     PLFLT dtr, theta, dx, dy, r;
     char text[3];
     PLFLT x0[361], y0[361];
     PLFLT x[361], y[361];
 
     dtr = PI / 180.0;
     for (i = 0; i <= 360; i++) {
 	x0[i] = cos(dtr * i);
 	y0[i] = sin(dtr * i);
     }
 
 
     plenv(-1.3, 1.3, -1.3, 1.3, 1, -2);
     for (i = 1; i <= 10; i++) {
 	for (j = 0; j <= 360; j++) {
 	    x[j] = 0.1 * i * x0[j];
 	    y[j] = 0.1 * i * y0[j];
 	}
 
 
 	plline(361, x, y);
     }
 
     plcol0(2);
     for (i = 0; i <= 11; i++) {
 	theta = 30.0 * i;
 	dx = cos(dtr * theta);
 	dy = sin(dtr * theta);
 
 
 	pljoin(0.0, 0.0, dx, dy);
 	sprintf(text, "%d", ROUND(theta));
 
 
 	if (dx >= -0.00001)
 	    plptex(dx, dy, dx, dy, -0.15, text);
 	else
 	    plptex(dx, dy, -dx, -dy, 1.15, text);
     }
 
 
     for (i = 0; i <= 360; i++) {
 	r = sin(dtr * (5 * i));
 	x[i] = x0[i] * r;
 	y[i] = y0[i] * r;
     }
     plcol0(3);
     plline(361, x, y);
 
     plcol0(4);
     plmtex("t", 2.0, 0.5, 0.5,
 	   "#frPLplot Example 3 - r(#gh)=sin 5#gh");
     plflush();
 }
 
 
 
 #define XPTS      35
 #define YPTS      46
 #define XSPA      2./(XPTS-1)
 #define YSPA      2./(YPTS-1)
 
 PLFLT tr[6] =
 {XSPA, 0.0, -1.0, 0.0, YSPA, -1.0};
 
 void
 mypltr(PLFLT x, PLFLT y, PLFLT *tx, PLFLT *ty, void *pltr_data)
 {
     *tx = tr[0] * x + tr[1] * y + tr[2];
     *ty = tr[3] * x + tr[4] * y + tr[5];
 }
 
 static PLFLT clevel[11] =
 {-1., -.8, -.6, -.4, -.2, 0, .2, .4, .6, .8, 1.};
 
 void
 plot5(void)
 {
     int i, j;
     PLFLT xx, yy;
     PLFLT **z, **w;
     static PLINT mark = 1500, space = 1500;
 
 
     plAlloc2dGrid(&z, XPTS, YPTS);
     plAlloc2dGrid(&w, XPTS, YPTS);
 
     for (i = 0; i < XPTS; i++) {
 	xx = (double) (i - (XPTS / 2)) / (double) (XPTS / 2);
 	for (j = 0; j < YPTS; j++) {
 	    yy = (double) (j - (YPTS / 2)) / (double) (YPTS / 2) - 1.0;
 	    z[i][j] = xx * xx - yy * yy;
 	    w[i][j] = 2 * xx * yy;
 	}
     }
 
     plenv(-1.0, 1.0, -1.0, 1.0, 0, 0);
     plcol0(2);
     plcont(z, XPTS, YPTS, 1, XPTS, 1, YPTS, clevel, 11, mypltr, NULL);
     plstyl(1, &mark, &space);
     plcol0(3);
     plcont(w, XPTS, YPTS, 1, XPTS, 1, YPTS, clevel, 11, mypltr, NULL);
     plcol0(1);
     pllab("X Coordinate", "Y Coordinate", "Streamlines of flow");
     plflush();
 
     plFree2dGrid(z, XPTS, YPTS);
     plFree2dGrid(w, XPTS, YPTS);
 
 }