2 * $Id: qplot.c,v 1.2 2003-11-24 17:31:28 pcg Exp $
14 # define GRX_SCALE 10000
17 #define DEFAULT_DATA_FILE "data"
20 axis_round(Real * min, Real * max, Real * grid_spacing)
24 logspace = (int)(log10((*max - *min) / 10.0) + 0.5);
25 *grid_spacing = pow(10, (double)logspace);
26 *min = (Real) ((int)(*min / (*grid_spacing))) * (*grid_spacing);
27 *max = (Real) ((int)(*max / (*grid_spacing)) + 1) * (*grid_spacing);
40 main(int argc, char **argv)
43 int Do_Start = 1, tmp;
44 int m, p, i, j, n, nchars, theight, twidth, xclick, yclick;
45 int downx = 1000, downy = 1000, upx, upy;
47 Real xmax, xmin, ymax, ymin, xdiff, ydiff, xgrid_spacing,
50 LineFunction linetype = StartLine;
51 char axis[100], line[256];
54 x = (Real *)malloc(1000000 * sizeof(Real));
55 y = (Real *)malloc(1000000 * sizeof(Real));
56 nls = (Real *)malloc(1000 * sizeof(Real));
57 if (x == NULL || y == NULL || nls == NULL) {
58 fprintf(stderr, "Can't allocate initial memory\n");
62 ymax = xmax = -HUGE_VAL;
63 ymin = xmin = HUGE_VAL;
65 for (i = 1; i < argc; i++) {
66 if (*argv[i] == '-') {
67 if (!strcmp(argv[i], "-nl"))
68 linetype = StartPoint;
69 else if (argv[i][1] == '\0') /* use stdin */
73 fprintf(stderr, "Usage:\n\t %s [options] [file]\n\n", argv[0]);
75 "where options include:\n"
76 " -pt plot with points instead of lines\n\n");
79 "file name `-' specifies stdin\n"
80 "if no file name is specified, "
81 "the default is \"%s\"\n\n", DEFAULT_DATA_FILE);
89 if (file && !strcmp(file, "-")) {
94 file = DEFAULT_DATA_FILE;
96 if ((fd = fopen(file, "r")) == NULL) {
97 fprintf(stderr, "%s: can't open file \"%s\"\n", argv[0], file);
103 while (fgets(line, sizeof(line), fd) != NULL) {
104 if (sscanf(line, "%f %f", &x[m], &y[m]) == 2) {
107 else if (x[m] < xmin)
111 else if (y[m] < ymin)
124 signal(SIGTERM, nice_end);
125 signal(SIGSTOP, nice_end);
126 signal(SIGTSTP, nice_end);
127 signal(SIGINT, nice_end);
128 signal(SIGQUIT, nice_end);
129 if (!InitializeGraphics(1))
134 axis_round(&xmin, &xmax, &xgrid_spacing);
135 axis_round(&ymin, &ymax, &ygrid_spacing);
137 id = CreateWin(0, 0, GRX_SCALE, GRX_SCALE);
139 fprintf(stderr, "Help id = 0\n");
142 /* Fill the window in black for real eye-catching graphics! */
145 FillArea(0, 0, GRX_SCALE, GRX_SCALE);
148 /* draw outline box in white */
151 /* Draw outline box */
160 /* Draw the data - either lines or dots */
161 xdiff = 8000 / (xmax - xmin);
162 ydiff = 8000 / (ymax - ymin);
164 for (i = j = 0; j < p; j++) {
167 ForeColor(j % 6 + 1);
168 while (((x[i] < xmin) || (x[i] > xmax) ||
169 (y[i] < ymin) || (y[i] > ymax)) && (i < nls[j]))
175 Extend(1000 + (x[i] - xmin) * xdiff,
176 9000 - (y[i] - ymin) * ydiff);
184 while ((i < nls[j]) &&
185 ((x[i] < xmin) || (x[i] > xmax) ||
186 (y[i] < ymin) || (y[i] > ymax)))
193 /* Do axis labels in black */
195 QueryWin(id, &twidth, &theight);
196 PlaceText(id, GRX_SCALE / 2, 0, HCENTER_TEXT | TOP_TEXT, file);
197 PlaceText(id, GRX_SCALE / 2, GRX_SCALE, HCENTER_TEXT | BOTTOM_TEXT,
199 PlaceText(id, 0, GRX_SCALE / 2, LEFT_TEXT | VCENTER_TEXT, "Y");
200 sprintf(axis, "%f", ymax);
201 nchars = 1000 / twidth;
203 PlaceText(id, GRX_SCALE / 10, GRX_SCALE / 10,
204 RIGHT_TEXT | TOP_TEXT, axis);
205 sprintf(axis, "%f", ymin);
207 PlaceText(id, GRX_SCALE / 10, 9 * GRX_SCALE / 10,
208 RIGHT_TEXT | BOTTOM_TEXT, axis);
209 sprintf(axis, "%f", xmax);
210 PlaceText(id, 9 * GRX_SCALE / 10, 9 * GRX_SCALE / 10,
211 HCENTER_TEXT | TOP_TEXT, axis);
212 sprintf(axis, "%f", xmin);
213 PlaceText(id, GRX_SCALE / 10, 9 * GRX_SCALE / 10,
214 HCENTER_TEXT | TOP_TEXT, axis);
218 n = WaitForCarriageReturn(&winclick, &xclick, &yclick);
237 xmin = (xmax - xmin) * (downx - 1000) / (8000) + xmin;
238 xmax = (xmax - xmin) * (upx - 1000) / (8000) + xmin;
239 ymax = ymax - (ymax - ymin) * (downy - 1000) / (8000);
240 ymin = ymax - (ymax - ymin) * (upy - 1000) / (8000);
243 } while (n && (n != 2));
245 nice_end(EXIT_SUCCESS);
249 /*----------------------- end-of-file (C source) -----------------------*/