Annotation of OpenXM_contrib/pari-2.2/src/graph/plotgnuplot.c, Revision 1.1
1.1 ! noro 1: /* $Id: plotgnuplot.c,v 1.8 2001/03/13 17:56:53 karim Exp $
! 2:
! 3: Copyright (C) 2000 The PARI group.
! 4:
! 5: This file is part of the PARI/GP package.
! 6:
! 7: PARI/GP is free software; you can redistribute it and/or modify it under the
! 8: terms of the GNU General Public License as published by the Free Software
! 9: Foundation. It is distributed in the hope that it will be useful, but WITHOUT
! 10: ANY WARRANTY WHATSOEVER.
! 11:
! 12: Check the License for details. You should have received a copy of it, along
! 13: with the package; see the file 'COPYING'. If not, write to the Free Software
! 14: Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. */
! 15:
! 16: /*******************************************************************/
! 17: /* */
! 18: /* HI-RES PLOT. GNUPLOT INTERFACE */
! 19: /* written by Ilya Zakharevich */
! 20: /* */
! 21: /*******************************************************************/
! 22: # include "pari.h"
! 23: #include "rect.h"
! 24: #define croak(str) err(talker,str)
! 25: #define SET_OPTIONS_FROM_STRING
! 26: #define GNUPLOT_OUTLINE_STDOUT
! 27: #define DONT_POLLUTE_INIT
! 28: #include "Gnuplot.h"
! 29:
! 30: #ifdef __EMX__
! 31: # define DEF_TERM "pm"
! 32: #else
! 33: # define DEF_TERM (getenv("DISPLAY") ? "X11" : "dumb")
! 34: #endif
! 35:
! 36: void
! 37: rectdraw0(long *w, long *x, long *y, long lw, long do_free)
! 38: {
! 39: double *ptx,*pty;
! 40: long i,j,x0,y0, hjust, vjust, hgap, vgap, hgapsize, vgapsize;
! 41: long good, seen_graph = 0;
! 42: int point_type = -1, line_type = 0;
! 43: PariRect *e;
! 44: RectObj *p1;
! 45: int strdir = RoSTdirLEFT, can_justify = 1, shift = 0, xstart, xend;
! 46:
! 47: (void)do_free;
! 48: PARI_get_plot(0);
! 49:
! 50: hgapsize = h_unit; vgapsize = v_unit;
! 51: if (hgapsize == 1)
! 52: hgapsize = 2; /* Vertical direction is subjectively different! */
! 53: /* Find the info about the *actual* x and y-coords of the
! 54: rectangles. Use the first rectangle with has_graph attribute. */
! 55:
! 56: for(i=0;i<lw;i++) {
! 57: e=rectgraph[w[i]];
! 58: if (RHasGraph(e)) {
! 59: set_mouse_feedback_rectangle(
! 60: x[i], x[i] + RXsize(e) - 1,
! 61: w_height - 1 - y[i] - (RYsize(e) - 1), w_height - 1 - y[i],
! 62: (0 - RXshift(e))/RXscale(e),
! 63: (RXsize(e) - 1 - RXshift(e))/RXscale(e),
! 64: (RYsize(e) - 1 - RYshift(e))/RYscale(e),
! 65: (0 - RYshift(e))/RYscale(e)
! 66: );
! 67: seen_graph = 1;
! 68: break;
! 69: }
! 70: }
! 71: if (!seen_graph) /* Put some reasonable values */
! 72: set_mouse_feedback_rectangle( 0, w_width - 1, 0, w_height - 1,
! 73: 0, 0, 0, 0);
! 74:
! 75: #if 0
! 76: graphics(); /* Switch on terminal. */
! 77: #else
! 78: term_start_plot(); /* Switch on terminal. */
! 79: #endif
! 80: linetype(line_type); /* X does not work otherwise. */
! 81: setpointsize(pointsize);
! 82:
! 83: for(i=0;i<lw;i++)
! 84: {
! 85: e=rectgraph[w[i]]; p1=RHead(e); x0=x[i]; y0=y[i];
! 86: while(p1)
! 87: {
! 88: switch(RoType(p1))
! 89: {
! 90: case ROt_PT:
! 91: point(DTOL(RoPTx(p1)+x0), w_height - 1 - DTOL(RoPTy(p1) + y0),
! 92: point_type);
! 93: break;
! 94: case ROt_LN:
! 95: move(DTOL(RoLNx1(p1)+x0), w_height - 1 - DTOL(RoLNy1(p1) + y0));
! 96: vector(DTOL(RoLNx2(p1)+x0), w_height - 1 - DTOL(RoLNy2(p1) + y0));
! 97: break;
! 98: case ROt_BX:
! 99: move(DTOL(RoBXx1(p1)+x0), w_height - 1 - DTOL(RoBXy1(p1) + y0));
! 100: vector(DTOL(RoBXx2(p1)+x0), w_height - 1 - DTOL(RoBXy1(p1) + y0));
! 101: vector(DTOL(RoBXx2(p1)+x0), w_height - 1 - DTOL(RoBXy2(p1) + y0));
! 102: vector(DTOL(RoBXx1(p1)+x0), w_height - 1 - DTOL(RoBXy2(p1) + y0));
! 103: vector(DTOL(RoBXx1(p1)+x0), w_height - 1 - DTOL(RoBXy1(p1) + y0));
! 104: break;
! 105: case ROt_MP:
! 106: ptx=RoMPxs(p1);
! 107: pty=RoMPys(p1);
! 108: for(j=0;j<RoMPcnt(p1);j++)
! 109: {
! 110: point(DTOL(ptx[j] + x0), w_height - 1 - DTOL(pty[j] + y0),
! 111: point_type);
! 112: }
! 113: break;
! 114: case ROt_ML:
! 115: ptx=RoMLxs(p1);
! 116: pty=RoMLys(p1);
! 117: j = 0;
! 118: if (DTOL(ptx[j]+x0) < 0 || DTOL(ptx[j]+x0) >= w_width
! 119: || DTOL(pty[j] + y0) < 0 || DTOL(pty[j] + y0) >= w_height) {
! 120: good = 0;
! 121: } else {
! 122: move(DTOL(ptx[j]+x0), w_height - 1 - DTOL(pty[j] + y0));
! 123: good = 1;
! 124: }
! 125: for(j=1;j<RoMLcnt(p1);j++)
! 126: {
! 127: if (good) {
! 128: if (DTOL(ptx[j] + x0) < 0 || DTOL(ptx[j]+x0) >= w_width
! 129: || DTOL(pty[j] + y0) < 0 || DTOL(pty[j] + y0) >= w_height) {
! 130: good = 0;
! 131: } else {
! 132: vector(DTOL(ptx[j] + x0), w_height - 1 - DTOL(pty[j] + y0));
! 133: }
! 134: } else {
! 135: if (DTOL(ptx[j] + x0) < 0 || DTOL(ptx[j] + x0) >= w_width
! 136: || DTOL(pty[j] + y0) < 0 || DTOL(pty[j] + y0) >= w_height) {
! 137: } else {
! 138: move(DTOL(ptx[j]+x0), w_height - 1 - DTOL(pty[j] + y0));
! 139: good = 1;
! 140: }
! 141: }
! 142: }
! 143: break;
! 144: case ROt_ST:
! 145: hjust = RoSTdir(p1) & RoSTdirHPOS_mask;
! 146: vjust = RoSTdir(p1) & RoSTdirVPOS_mask;
! 147: hgap = RoSTdir(p1) & RoSTdirHGAP;
! 148: if (hgap)
! 149: hgap = (hjust == RoSTdirLEFT) ? hgapsize : -hgapsize;
! 150: vgap = RoSTdir(p1) & RoSTdirVGAP;
! 151: if (vgap)
! 152: vgap = (vjust == RoSTdirBOTTOM) ? vgapsize : -vgapsize;
! 153: if (vjust != RoSTdirVCENTER)
! 154: vgap += ((vjust == RoSTdirTOP) ? -1 : 1) * (f_height - 1)/2;
! 155: if (strdir != hjust) {
! 156: shift = (hjust == RoSTdirLEFT ? 0 :
! 157: (hjust == RoSTdirRIGHT ? 2 : 1));
! 158: can_justify = justify_text(shift); /* 1 for LEFT */
! 159: strdir = RoSTdir(p1);
! 160: }
! 161: xstart = DTOL(RoSTx(p1) + x0) + hgap
! 162: - (can_justify ? 0
! 163: : ((RoSTl(p1) * pari_plot.fwidth - 1) * shift / 2));
! 164: xend = xstart + (can_justify ? 0 : RoSTl(p1) * pari_plot.fwidth - 1);
! 165: if (xstart < 0 || xend >= w_width
! 166: || DTOL(RoSTy(p1) + y0) < 0
! 167: || DTOL(RoSTy(p1) + y0) >= w_height) {
! 168: } else {
! 169: put_text(xstart,
! 170: w_height - 1 - DTOL(RoSTy(p1) + y0) + vgap,
! 171: RoSTs(p1));
! 172: }
! 173: break;
! 174: case ROt_PTT:
! 175: point_type = RoPTTpen(p1);
! 176: break;
! 177: case ROt_PTS:
! 178: pointsize = RoPTSsize(p1);
! 179: setpointsize(pointsize);
! 180: break;
! 181: case ROt_LNT:
! 182: linetype(RoLNTpen(p1));
! 183: break;
! 184: default: break;
! 185: }
! 186: p1=RoNext(p1);
! 187: }
! 188: }
! 189: #if 0
! 190: text(); /* Reset terminal */
! 191: #else
! 192: term_end_plot(); /* Reset terminal. */
! 193: #endif
! 194: }
! 195:
! 196: void
! 197: PARI_get_plot(long fatal)
! 198: {
! 199: (void)fatal;
! 200: if (pari_plot.init) {
! 201: return;
! 202: }
! 203: term_set( DEF_TERM );
! 204: }
! 205:
! 206:
! 207: long
! 208: term_set(char *s)
! 209: {
! 210: char *t, *size = NULL;
! 211: double x, y;
! 212: static int had_error;
! 213:
! 214: setup_gpshim();
! 215: if (*s == 0)
! 216: s = pari_plot.name;
! 217: t = s;
! 218: if (t[1] == '\0' && t[0] == '?') {
! 219: list_terms();
! 220: return 1;
! 221: }
! 222: while (*t && !(*t == ' ' || *t == '\t' || *t == '\n' || *t == '='))
! 223: t++;
! 224: if ((t-s) > PLOT_NAME_LEN)
! 225: err(talker,"name \"%s\" for terminal too long", s);
! 226: if (*pari_plot.name && !had_error
! 227: && (strlen(pari_plot.name) != t - s /* As strcmp() without \0 at end */
! 228: || (strncmp(pari_plot.name, s, t-s) != 0)) )
! 229: reset();
! 230: strncpy(pari_plot.name,s,t-s);
! 231: pari_plot.name[t-s] = '\0';
! 232:
! 233: had_error = 1;
! 234: if (!termset( pari_plot.name ))
! 235: err(talker,"error setting terminal \"%s\"", pari_plot.name);
! 236: had_error = 0;
! 237:
! 238: if (*t == '=') {
! 239: size = ++t;
! 240: x = atof(size);
! 241: while (*t && !(*t == ' ' || *t == '\t' || *t == '\n' || *t == ','))
! 242: t++;
! 243: if (*t != ',')
! 244: err(talker, "Terminal size directive without ','");
! 245: y = atof(++t);
! 246: while (*t && !(*t == ' ' || *t == '\t' || *t == '\n'))
! 247: t++;
! 248: plotsizes_scale(x*(1 + 1e-6)/termprop(xmax),
! 249: y*(1 + 1e-6)/termprop(ymax)); /* Later - truncated! */
! 250: } else {
! 251: plotsizes_scale(1,1);
! 252: }
! 253:
! 254: /* *Needed*, say, by gif output: */
! 255: set_options_from(t);
! 256:
! 257: #if 0
! 258: gptable_init(); /* Init terminal. */
! 259: #else
! 260: term_init();
! 261: #endif
! 262:
! 263: setpointsize(pointsize);
! 264:
! 265: w_width = scaled_xmax();
! 266: w_height = scaled_ymax();
! 267: f_height = termprop(v_char);
! 268: f_width = termprop(h_char);
! 269: h_unit = termprop(h_tic);
! 270: v_unit = termprop(v_tic);
! 271: pari_plot.init = 1;
! 272:
! 273: return 1;
! 274: }
! 275:
! 276: long
! 277: plot_outfile_set(char *s) {
! 278: int normal = (strcmp(s,"-") == 0);
! 279:
! 280: setup_gpshim();
! 281: /* Delegate all the hard work to term_set_output() */
! 282:
! 283: if (normal)
! 284: term_set_output(NULL);
! 285: else { /* term_set_output() needs
! 286: a malloced string */
! 287: char *s1 = (char*) malloc(strlen(s) + 1);
! 288:
! 289: strcpy(s1,s);
! 290: term_set_output(s1);
! 291: }
! 292: return 1;
! 293: }
! 294:
! 295: void
! 296: set_pointsize(double d)
! 297: {
! 298: pointsize = d;
! 299: if (pari_plot.init)
! 300: setpointsize(d);
! 301: }
! 302:
! 303: #ifdef DYNAMIC_PLOTTING_RUNTIME_LINK
! 304: #include <dlfcn.h>
! 305:
! 306: get_term_ftable_t *
! 307: get_term_ftable_get(void) /* Establish runtime link with gnuplot engine */
! 308: {
! 309: char *s = getenv("GNUPLOT_DRAW_DLL"), buf[4096];
! 310: void *h, *f;
! 311: int mode = RTLD_LAZY;
! 312:
! 313: #ifdef RTLD_GLOBAL
! 314: mode |= RTLD_GLOBAL;
! 315: #endif
! 316:
! 317: if (!s)
! 318: s = DYNAMIC_PLOTTING_RUNTIME_LINK;
! 319: h = dlopen(s, mode);
! 320: if (!h) {
! 321: sprintf(buf,"Can't load Gnuplot drawing engine from '%s': %s", s, dlerror());
! 322: croak(buf);
! 323: return 0;
! 324: }
! 325: f = dlsym(h, "get_term_ftable");
! 326: if (f)
! 327: return (get_term_ftable_t *)f;
! 328: sprintf(buf, "Can't resolve 'get_term_ftable' function from Gnuplot drawing engine '%s': %s", s, dlerror());
! 329: croak(buf);
! 330: return 0;
! 331: }
! 332: #endif
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to plotgnuplot.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / pari-2.2 / src / graph