src/qplot.c

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
  
#include <stdio.h>
  #include <stdlib.h>
  #include <string.h>
  #include <signal.h>

  FILE * gnuplot = NULL;

  void cleanup (int number)
  {
    pclose (gnuplot);
    int ret = system ("rm -f .qplot??????");
    exit (1 + ret);
  }

  int main (int argc, char * argv[])
  {
    float fn;
    char command[80] = "gnuplot", * files[80], * commands[80];
    int i, nf = 0, nc = 0;
    for (i = 1; i < argc; i++)
      if (!strcmp (argv[i], "-e")) {
        if (i + 1 < argc)
  	commands[nc++] = argv[++i];
      }
      else if (argv[i][0] == '-') {
        strcat (command, " "); strcat (command, argv[i]);
      }
      else
        files[nf++] = argv[i];
    gnuplot = popen (command, "w");
    if (!gnuplot) {
      perror ("qplot: could not open pipe: ");
      return 1;
    }
    fprintf (gnuplot, "set terminal wxt noraise\n");
    for (i = 0; i < nf; i++)
      fprintf (gnuplot, "load '%s'\n", files[i]);
    for (i = 0; i < nc; i++)
      fprintf (gnuplot, "%s\n", commands[i]);
    signal (SIGINT, cleanup);
    while (fread (&fn, sizeof(float), 1, stdin) == 1) {
      int n = fn;
      float * v = malloc (sizeof (float)*(n + 1));
      if (v == NULL)
        break;
      char fname[] = ".qplotXXXXXX";
      int fd = mkstemp (fname), nr;
      if (fd < 0) {
        perror ("qplot: could not create temporary file: ");
        return 1;
      }
      FILE * fp = fdopen (fd, "w");
      fwrite (&fn, sizeof(float), 1, fp);
      if ((nr = fread (v, sizeof (float), n, stdin)) != n) {
        fprintf (stderr, "qplot: expecting %d y-values only got %d\n", n, nr);
        remove (fname);
        return 1;
      }
      fwrite (v, sizeof (float), n, fp);
      int i;
      for (i = 0; i < n; i++) {
        if ((nr = fread (v, sizeof (float), n + 1, stdin)) != n + 1) {
  	fprintf (stderr, "qplot: expecting %d x-z-values only got %d\n", n + 1, nr);
  	remove (fname);
  	return 1;
        }
        fwrite (v, sizeof (float), n + 1, fp);
      }
      free (v);
      fclose (fp);
      fprintf (gnuplot, "splot '%s' binary u 2:1:3 t ''\n", fname);
      fprintf (gnuplot, "! rm -f %s\n", fname);
      fflush (gnuplot);
    }
    pclose (gnuplot);
    int ret = system ("rm -f .qplot??????");
    return 0*ret;
  }