Subversion Repositories wimsdev

Rev

Rev 18325 | Rev 18350 | Go to most recent revision | Blame | Compare with Previous | Last modification | View Log | RSS feed

  1. /*    Copyright (C) 1998-2003 XIAO, Gang of Universite de Nice - Sophia Antipolis
  2.  *
  3.  *  This program is free software; you can redistribute it and/or modify
  4.  *  it under the terms of the GNU General Public License as published by
  5.  *  the Free Software Foundation; either version 2 of the License, or
  6.  *  (at your option) any later version.
  7.  *
  8.  *  This program is distributed in the hope that it will be useful,
  9.  *  but WITHOUT ANY WARRANTY; without even the implied warranty of
  10.  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  11.  *  GNU General Public License for more details.
  12.  *
  13.  *  You should have received a copy of the GNU General Public License
  14.  *  along with this program; if not, write to the Free Software
  15.  *  Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
  16.  */
  17.  
  18. #include <errno.h>
  19. #include "flydraw.h"
  20. void myGdImageArc(gdImagePtr, int, int, int, int, double, double, int);
  21. const double DEG=M_PI/180;
  22.  
  23. /* Tikz things */
  24.  
  25. static char *tikz_color (int the_color)
  26. {
  27.   static char buf[50];
  28.   int r=gdImageRed(image,the_color);
  29.   int g=gdImageGreen(image,the_color);
  30.   int b=gdImageBlue(image,the_color);
  31.   sprintf(buf, "{rgb,255:red,%i;green,%i;blue,%i}",r,g,b);
  32.   return buf;
  33. }
  34.  
  35. static char *tikz_options (int the_color, int the_fill)
  36. {
  37.   static char buf[100];
  38.   if (the_color == gdBrushed)
  39.     the_color = tikz_brushColor;
  40.   sprintf (buf, "draw=%s,cap=round",tikz_color(the_color));
  41.   if (the_color != tikz_brushColor) strcat(buf, ",thin");
  42.   if (the_fill == -1 || the_fill == 2) strcat(buf, ",dashed");
  43.   if (the_fill == 1 || the_fill == 2)
  44.     sprintf(buf+strlen(buf), ",fill=%s", tikz_color(the_color));
  45.   return buf;
  46. }
  47.  
  48. #define flip(y) (sizey-1-(y))
  49.  
  50. /* bug in gdImageFillToBorder */
  51.  
  52. void patchgdImageFillToBorder (gdImagePtr im, int x, int y, int border, int color)
  53. {
  54.   if(x>=im->sx || x<0 || y>=im->sy || y<0) return;
  55.   gdImageFillToBorder(im,x,y,border,color);
  56. }
  57.  
  58. void patchgdImageFill (gdImagePtr im, int x, int y, int color)
  59. {
  60.   if(x>=im->sx || x<0 || y>=im->sy || y<0) return;
  61.   gdImageFill(im,x,y,color);
  62. }
  63.  
  64. #define DASHPIXELS 8
  65.  
  66. void myDashedLine(gdImagePtr im, int x1, int y1, int x2, int y2, int color){
  67.   /**
  68.    * emulates gdImageDashedLine with gdImageSetStyle and gdImageLine
  69.    * GK: As gdImageDashedLine is slightly broken in latest libgd libraries,
  70.    * GK: I implemented an emulation named "myDashedLine"
  71.    **/
  72.  
  73.   int styleDashed[DASHPIXELS],i;
  74.   for (i=0; i< DASHPIXELS/2; i++) styleDashed[i]=color;
  75.   for (; i< DASHPIXELS; i++) styleDashed[i]=gdTransparent;
  76.   gdImageSetStyle(im, styleDashed, DASHPIXELS);
  77.   gdImageLine(im, x1, y1, x2, y2, gdStyled);
  78. }
  79.  
  80.  
  81. /* File opening: with security */
  82. FILE *open4read(char *n)
  83. {
  84.   char *p, *p1, *p2, namebuf[2048];
  85.   int t;
  86.   FILE *f;
  87.   n=find_word_start(n);
  88.   if(*n==0) return NULL;
  89.   p=getenv("flydraw_filebase");
  90.   p1=n+strlen(n)-4;if(p1<n || strcasecmp(p1,".gif")!=0) t=1; else t=0;
  91.   if(p!=NULL && *p!=0) {
  92.     char pbuf[MAX_LINELEN+1];
  93.     snprintf(pbuf,sizeof(pbuf),"%s",p);
  94.     p=find_word_start(pbuf); if(strstr(p,"..")!=NULL) return NULL;
  95.     if(*n=='/' || strstr(n,"..")!=NULL) return NULL;
  96.     /* prohibit unusual file/dir names */
  97.     for(p1=p;*p1;p1++)
  98.       if(!isalnum(*p1) && !isspace(*p1) && strchr("~_-/.",*p1)==NULL)
  99.         return NULL;
  100.     for(p1=n;*p1;p1++)
  101.       if(!isalnum(*p1) && !isspace(*p1) && strchr("~_-/.",*p1)==NULL)
  102.         return NULL;
  103.     f=NULL;
  104.     for(p1=p; *p1; p1=find_word_start(p2)) {
  105.       p2=find_word_end(p1);
  106.       if(*p2) *p2++=0;
  107.       snprintf(namebuf,sizeof(namebuf),"%s/%s",p1,n);
  108.       f=fopen(namebuf,"r"); if(f!=NULL) goto imgtype;
  109.     }
  110.     p1=getenv("w_wims_session");
  111.     if(p1!=NULL && strncmp(n,"insert",6)==0) {
  112.       snprintf(namebuf,sizeof(namebuf),"../s2/%s/%s",p1,n);
  113.       f=fopen(namebuf,"r");
  114.     }
  115.   }
  116.   else {
  117.     snprintf(namebuf,sizeof(namebuf),"%s",n);
  118.     f=fopen(namebuf,"r");
  119.   }
  120.   imgtype:
  121.   if(t && f!=NULL) {
  122.     char tbuf[1024],sbuf[4096];
  123.     fclose(f); f=NULL;
  124.     p1=getenv("TMPDIR"); if(p1==NULL || *p1==0) p1=".";
  125.     snprintf(tbuf,sizeof(tbuf),"%s/drawfile_.gif",p1);
  126.     snprintf(sbuf,sizeof(sbuf),"convert %s %s",namebuf,tbuf);
  127.     if (system(sbuf)) fly_error("system_failed");
  128.       f=fopen(tbuf,"r");
  129.   }
  130.   return f;
  131. }
  132.  
  133. /* Does nothing; just a comment. */
  134. void obj_comment(objparm *pm)
  135. {
  136.   return;
  137. }
  138.  
  139. /* define image size */
  140. void obj_size(objparm *pm)
  141. {
  142.   sizex=rint(pm->pd[0]); sizey=rint(pm->pd[1]);
  143.   if(sizex<0 || sizey<0 || sizex>MAX_SIZE || sizey>MAX_SIZE) {
  144.     fly_error("bad_size"); return;
  145.   }
  146.   if(image!=NULL) {
  147.     fly_error("size_already_defined"); return;
  148.   }
  149.   image=gdImageCreate(sizex,sizey);
  150.   if(tikz_file) fprintf(tikz_file,"\\clip (0,0) rectangle (%i, %i);\n", sizex, sizey);
  151.   if(image==NULL) fly_error("image_creation_failure");
  152.   else {
  153.     color_white=gdImageColorAllocate(image,255,255,255);
  154.     color_black=gdImageColorAllocate(image,0,0,0);
  155.     color_bounder=gdImageColorAllocate(image,1,2,3);
  156.     color_frame=gdImageColorAllocate(image,254,254,254);
  157.   }
  158. }
  159.  
  160. /* new image */
  161. void obj_new(objparm *pm)
  162. {
  163.   if(image) {
  164.     gdImageDestroy(image);image=NULL;
  165.   }
  166.   if(pm->pcnt>=2) { obj_size(pm); gdImageRectangle(image,0,0,sizex-1,sizey-1,color_frame);}
  167.   else sizex=sizey=0;
  168.   saved=0;
  169. }
  170.  
  171. /* new image */
  172. void obj_existing(objparm *pm)
  173. {
  174.   FILE *inf;
  175.   char *pp;
  176.  
  177.   if(image) {
  178.     gdImageDestroy(image);image=NULL;
  179.   }
  180.   pp=find_word_start(pm->str);*find_word_end(pp)=0;
  181.   inf=open4read(pp);
  182.   if(inf==NULL) {
  183.     fly_error("file_not_exist"); return;
  184.   }
  185.   image=gdImageCreateFromGif(inf); fclose(inf);
  186.   if(image==NULL) {
  187.     fly_error("bad_gif"); return;
  188.   }
  189.   sizex=image->sx; sizey=image->sy;
  190.   saved=0;
  191. }
  192. /* tikzfile */
  193. void obj_tikzfile(objparm *pm)
  194. {
  195.   char *p;
  196.   if(tikz_file){
  197.     fprintf(tikz_file,"\\end{tikzpicture}\n");
  198.     fclose(tikz_file);
  199.   }
  200.   p=find_word_start(pm->str); *find_word_end(p)=0;
  201.   snprintf(tikzfilename,sizeof(tikzfilename),"%s",p);
  202.   tikz_file=fopen(tikzfilename,"w");
  203.   if(!tikz_file)
  204.     fly_error("tikz_nopen");
  205.   fprintf(tikz_file,"\\begin{tikzpicture}\[x=0.02cm, y=0.02cm]\n");
  206. }
  207.  
  208. /* segment */
  209. void obj_line(objparm *pm)
  210. {
  211.   scale(pm->pd,pm->p,2);
  212.   if(pm->fill==-1 || pm->fill==2 )
  213.     myDashedLine(image,pm->p[0],pm->p[1],pm->p[2],pm->p[3],pm->color[0]);
  214.   else
  215.     gdImageLine(image,pm->p[0],pm->p[1],pm->p[2],pm->p[3],pm->color[0]);
  216.   if (tikz_file)
  217.     fprintf(tikz_file, "\\draw\[%s] (%i, %i) -- (%i, %i);\n",
  218.       tikz_options(pm->color[0],pm->fill),pm->p[0],flip(pm->p[1]),pm->p[2],flip(pm->p[3]));
  219.   if(vimg_enable) vimg_line(scale_buf[0],scale_buf[1],scale_buf[2],scale_buf[3]);
  220. }
  221.  
  222. void _obj_arrow(objparm *pm, int twoside)
  223. {
  224.   int l,xx,yy;
  225.   gdPoint ii[3];
  226.   double dx,dy,length,dd[6];
  227.   scale(pm->pd,pm->p,2);
  228.   xx=ii[0].x=pm->p[2];yy=ii[0].y=pm->p[3];
  229.   l=pm->pd[4];if(l<0) l=0; if(l>200) l=200;
  230.   scale2(pm->pd[0]-pm->pd[2],pm->pd[1]-pm->pd[3],&dx,&dy);
  231.   length=sqrt(dx*dx+dy*dy);
  232.   if(length<3 || l<5) goto stem;
  233.   dd[0]=l*dx/length; dd[1]=l*dy/length;
  234.   #define fat 0.27
  235.   dd[2]=dd[0]+dd[1]*fat; dd[3]=dd[1]-dd[0]*fat;
  236.   dd[4]=dd[0]-dd[1]*fat; dd[5]=dd[1]+dd[0]*fat;
  237.   ii[1].x=rint(dd[2])+ii[0].x; ii[1].y=rint(dd[3])+ii[0].y;
  238.   ii[2].x=rint(dd[4])+ii[0].x; ii[2].y=rint(dd[5])+ii[0].y;
  239.   gdImageFilledPolygon(image, ii,3,pm->color[0]);
  240.   if(tikz_file)
  241.     fprintf(tikz_file, "\\draw\[%s] (%i, %i) -- (%i, %i) -- (%i, %i) -- cycle;\n",
  242.       tikz_options(pm->color[0],1),ii[0].x,flip(ii[0].y),ii[1].x,flip(ii[1].y),ii[2].x,flip(ii[2].y));
  243.   if(vimg_enable) vimg_polyline(scale_buf,3,1);
  244.   xx=rint(dd[0])+ii[0].x;yy=rint(dd[1])+ii[0].y;
  245.   if(twoside) {
  246.     ii[0].x=pm->p[0]; ii[0].y=pm->p[1];
  247.     ii[1].x=-rint(dd[2])+ii[0].x; ii[1].y=-rint(dd[3])+ii[0].y;
  248.     ii[2].x=-rint(dd[4])+ii[0].x; ii[2].y=-rint(dd[5])+ii[0].y;
  249.     gdImageFilledPolygon(image, ii,3,pm->color[0]);
  250.   if(tikz_file)
  251.     fprintf(tikz_file, "\\draw\[%s] (%i, %i) -- (%i, %i) -- (%i, %i) -- cycle;\n",
  252.       tikz_options(pm->color[0],1),ii[0].x,flip(ii[0].y),ii[1].x,flip(ii[1].y),ii[2].x,flip(ii[2].y));
  253.   }
  254.   stem: if(pm->fill==-1 || pm->fill==2)
  255.     myDashedLine(image,pm->p[0],pm->p[1],xx,yy,pm->color[0]);
  256.   else
  257.     gdImageLine(image,pm->p[0],pm->p[1],xx,yy,pm->color[0]);
  258.   if(tikz_file)
  259.     fprintf(tikz_file, "\\draw\[%s] (%i, %i) -- (%i, %i);\n",
  260.       tikz_options(pm->color[0],pm->fill),pm->p[0],flip(pm->p[1]),xx,flip(yy));
  261.   if(vimg_enable) vimg_line(scale_buf[0],scale_buf[1],scale_buf[2],scale_buf[3]);
  262. }
  263.  
  264. /* Arrow */
  265. void obj_arrow(objparm *pm)
  266. {
  267.   _obj_arrow(pm,0);
  268. }
  269.  
  270. /* 2-sided arrow */
  271. void obj_arrow2(objparm *pm)
  272. {
  273.    _obj_arrow(pm,1);
  274. }
  275.  
  276. void _obj_arrows(objparm *pm, int twoside)
  277. {
  278.   #define fat 0.27
  279.   int i,l,xx,yy;
  280.   gdPoint ii[3];
  281.   double dx,dy,length,dd[6];
  282.   l=pm->pd[0];if(l<0) l=0; if(l>200) l=200;
  283.   scale(pm->pd+1,pm->p+1,pm->pcnt/2);
  284.   for (i=1;i<pm->pcnt;i+=4){
  285.     xx=ii[0].x=pm->p[i+2];yy=ii[0].y=pm->p[i+3];
  286.     scale2(pm->pd[i]-pm->pd[i+2],pm->pd[i+1]-pm->pd[i+3],&dx,&dy);
  287.     length=sqrt(dx*dx+dy*dy);
  288.     if(length<3 || l<5) goto stem;
  289.     dd[0]=l*dx/length; dd[1]=l*dy/length;
  290.     dd[2]=dd[0]+dd[1]*fat; dd[3]=dd[1]-dd[0]*fat;
  291.     dd[4]=dd[0]-dd[1]*fat; dd[5]=dd[1]+dd[0]*fat;
  292.     ii[1].x=rint(dd[2])+ii[0].x; ii[1].y=rint(dd[3])+ii[0].y;
  293.     ii[2].x=rint(dd[4])+ii[0].x; ii[2].y=rint(dd[5])+ii[0].y;
  294.     gdImageFilledPolygon(image, ii,3,pm->color[0]);
  295.     if(tikz_file)
  296.   fprintf(tikz_file, "\\draw\[%s] (%i, %i) -- (%i, %i) -- (%i, %i) -- cycle;\n",
  297.         tikz_options(pm->color[0],1),ii[0].x,flip(ii[0].y),ii[1].x,flip(ii[1].y),ii[2].x,flip(ii[2].y));
  298.       xx=rint(dd[0])+ii[0].x;yy=rint(dd[1])+ii[0].y;
  299.       if(twoside) {
  300.   ii[0].x=pm->p[i]; ii[0].y=pm->p[i+1];
  301.   ii[1].x=-rint(dd[2])+ii[0].x; ii[1].y=-rint(dd[3])+ii[0].y;
  302.   ii[2].x=-rint(dd[4])+ii[0].x; ii[2].y=-rint(dd[5])+ii[0].y;
  303.   gdImageFilledPolygon(image, ii,3,pm->color[0]);
  304.   if(tikz_file)
  305.     fprintf(tikz_file, "\\draw\[%s] (%i, %i) -- (%i, %i) -- (%i, %i) -- cycle;\n",
  306.         tikz_options(pm->color[0],1),ii[0].x,flip(ii[0].y),ii[1].x,flip(ii[1].y),ii[2].x,flip(ii[2].y));
  307.       }
  308.     stem: if(pm->fill==-1 || pm->fill==2)
  309.   myDashedLine(image,pm->p[i],pm->p[i+1],xx,yy,pm->color[0]);
  310.       else
  311.   gdImageLine(image,pm->p[i],pm->p[i+1],xx,yy,pm->color[0]);
  312.       if(tikz_file)
  313.   fprintf(tikz_file, "\\draw\[%s] (%i, %i) -- (%i, %i);\n",
  314.     tikz_options(pm->color[0],pm->fill),pm->p[i],flip(pm->p[i+1]),xx,flip(yy));
  315.     }
  316. }
  317.  
  318. /* Arrows */
  319. void obj_arrows(objparm *pm)
  320. {
  321.   _obj_arrows(pm,0);
  322. }
  323.  
  324. /* 2-sided arrows */
  325. void obj_arrows2(objparm *pm)
  326. {
  327.    _obj_arrows(pm,1);
  328. }
  329.  
  330. /* horizontal line */
  331. void obj_hline(objparm *pm)
  332. {
  333.   scale(pm->pd,pm->p,1);
  334.   if(pm->fill==-1 || pm->fill==2)
  335.     myDashedLine(image,0,pm->p[1],sizex,pm->p[1],pm->color[0]);
  336.   else
  337.     gdImageLine(image,0,pm->p[1],sizex,pm->p[1],pm->color[0]);
  338.   if (tikz_file)
  339.     fprintf(tikz_file, "\\draw\[%s] (%i, %i) -- (%i, %i);\n",
  340.       tikz_options(pm->color[0],pm->fill),0,flip(pm->p[1]),sizex,flip(pm->p[1]));
  341. }
  342.  
  343. static void line_extend(int x1, int y1, int x2, int y2, int* xend, int *yend)
  344. {
  345.   double t1, t2, t3;
  346.   if (x2 != x1){t1 = (x2 > x1) ? sizex : 0; t1 = (t1 - x1)/(x2 - x1);}
  347.   else t1 = INFINITY;
  348.   if (y2 != y1){t2 = (y2 > y1) ? sizey : 0; t2 = (t2 - y1)/(y2 - y1);}
  349.   else t2 = INFINITY;
  350.   t3 = (t1 < t2) ? t1 : t2;
  351.   *xend = x1+t3*(x2-x1);
  352.   *yend = y1+t3*(y2-y1);
  353. }
  354.  
  355. /* halfline */
  356. void obj_halfline(objparm *pm)
  357. {
  358.   int xend,yend;
  359.   scale(pm->pd,pm->p,2);
  360.   line_extend(pm->p[0],pm->p[1],pm->p[2],pm->p[3],&xend,&yend);
  361.   if(pm->fill==-1 || pm->fill==2)
  362.     myDashedLine(image,pm->p[0],pm->p[1],xend,yend,pm->color[0]);
  363.   else
  364.     gdImageLine(image,pm->p[0],pm->p[1],xend,yend,pm->color[0]);
  365.   if (tikz_file)
  366.     fprintf(tikz_file, "\\draw\[%s] (%i, %i) -- (%i, %i);\n",
  367.       tikz_options(pm->color[0],pm->fill),
  368.       pm->p[0],flip(pm->p[1]),xend,flip(yend));
  369. }
  370.  
  371. void obj_fullline(objparm *pm)
  372. {
  373.   int xstart,ystart,xend,yend;
  374.   scale(pm->pd,pm->p,2);
  375.   line_extend(pm->p[0],pm->p[1],pm->p[2],pm->p[3],&xend,&yend);
  376.   line_extend(pm->p[2],pm->p[3],pm->p[0],pm->p[1],&xstart,&ystart);
  377.   if(pm->fill==-1 || pm->fill==2)
  378.     myDashedLine(image,xstart,ystart,xend,yend,pm->color[0]);
  379.   else
  380.     gdImageLine(image,xstart,ystart,xend,yend,pm->color[0]);
  381.   if (tikz_file)
  382.     fprintf(tikz_file, "\\draw\[%s] (%i, %i) -- (%i, %i);\n",
  383.       tikz_options(pm->color[0],pm->fill),xstart,flip(ystart),xend,flip(yend));
  384. }
  385.  
  386. void obj_fulllines(objparm *pm)
  387. {
  388.   int i,xstart,ystart,xend,yend;
  389.   scale(pm->pd,pm->p,pm->pcnt/2);
  390.   if (tikz_file)
  391.       fprintf(tikz_file, "\\draw\[%s]",tikz_options(pm->color[0],pm->fill));
  392.   for (i=0;i<pm->pcnt;i+=4){
  393.     line_extend(pm->p[i],pm->p[i+1],pm->p[i+2],pm->p[i+3],&xend,&yend);
  394.     line_extend(pm->p[i+2],pm->p[i+3],pm->p[i],pm->p[i+1],&xstart,&ystart);
  395.     if(pm->fill==-1 || pm->fill==2)
  396.       myDashedLine(image,xstart,ystart,xend,yend,pm->color[0]);
  397.     else
  398.       gdImageLine(image,xstart,ystart,xend,yend,pm->color[0]);
  399.     if (tikz_file)
  400.       fprintf(tikz_file, "(%i,%i)--(%i, %i)",xstart,flip(ystart),xend,flip(yend));
  401.   }
  402.   if (tikz_file)
  403.       fprintf(tikz_file, ";\n");
  404. }
  405.  
  406. /* vertical line */
  407. void obj_vline(objparm *pm)
  408. {
  409.   scale(pm->pd,pm->p,1);
  410.   if(pm->fill==-1 || pm->fill==2)
  411.     myDashedLine(image,pm->p[0],0,pm->p[0],sizey,pm->color[0]);
  412.   else
  413.     gdImageLine(image,pm->p[0],0,pm->p[0],sizey,pm->color[0]);
  414.   if (tikz_file)
  415.     fprintf(tikz_file, "\\draw\[%s] (%i, %i) -- (%i, %i);\n",
  416.       tikz_options(pm->color[0],pm->fill),pm->p[0],0,pm->p[0],sizey);
  417. }
  418.  
  419. /* dashed line */
  420. void obj_dline(objparm *pm)
  421. {
  422.   scale(pm->pd,pm->p,2);
  423.   myDashedLine(image,pm->p[0],pm->p[1],pm->p[2],pm->p[3], pm->color[0]);
  424.   if (tikz_file)
  425.     fprintf(tikz_file, "\\draw\[%s] (%i, %i) -- (%i, %i);\n",
  426.       tikz_options(pm->color[0],-1),
  427.       pm->p[0],flip(pm->p[1]),pm->p[2],flip(pm->p[3]));
  428. }
  429.  
  430. /* parallel lines.
  431.  * x1,y1,x2,y2,xv,yv,n,color */
  432. void obj_parallel(objparm *pm)
  433. {
  434.   int i, n, xi,yi;
  435.   double xv,yv;
  436.   n=pm->pd[6]; if(n<0) return; if(n>256) n=256;
  437.   scale(pm->pd,pm->p,3);
  438.   scale2(pm->pd[4],pm->pd[5],&xv,&yv);
  439.   if (tikz_file)
  440.     fprintf(tikz_file, "\\draw\[%s]", tikz_options(pm->color[0],pm->fill));
  441.   for(i=0;i<n;i++) {
  442.     xi=rint(i*xv); yi=rint(i*yv);
  443.     if(pm->fill==-1 || pm->fill==2)
  444.       myDashedLine(image,pm->p[0]+xi,pm->p[1]+yi,pm->p[2]+xi,pm->p[3]+yi,
  445.         pm->color[0]);
  446.     else
  447.       gdImageLine(image,pm->p[0]+xi,pm->p[1]+yi,pm->p[2]+xi,pm->p[3]+yi,
  448.         pm->color[0]);
  449.     if (tikz_file)
  450.       fprintf(tikz_file, "(%i,%i)--(%i, %i)",
  451.         pm->p[0]+xi,flip(pm->p[1]+yi),pm->p[2]+xi,flip(pm->p[3]+yi));
  452.     if(vimg_enable) vimg_line(scale_buf[0]+i*(scale_buf[4]-transx),
  453.       scale_buf[1]+i*(scale_buf[5]-transy),
  454.       scale_buf[2]+i*(scale_buf[4]-transx),
  455.       scale_buf[3]+i*(scale_buf[5]-transy));
  456.   }
  457.  if(tikz_file) fprintf(tikz_file,";\n");
  458. }
  459.  
  460. /* rectangle */
  461. void obj_rect(objparm *pm)
  462. {
  463.   int x1,y1,x2,y2;
  464.   scale(pm->pd,pm->p,2);
  465.   x1=min(pm->p[0],pm->p[2]); x2=max(pm->p[0],pm->p[2]);
  466.   y1=min(pm->p[1],pm->p[3]); y2=max(pm->p[1],pm->p[3]);
  467.   if(pm->fill!=0){
  468.     if(pm->fill==-1 || pm->fill==2){
  469.       myDashedLine(image,x1,y1,x1,y2,pm->color[0]);
  470.       myDashedLine(image,x1,y2,x2,y2,pm->color[0]);
  471.       myDashedLine(image,x2,y2,x2,y1,pm->color[0]);
  472.       myDashedLine(image,x2,y1,x1,y1,pm->color[0]);
  473.     }
  474.     if(pm->fill==1 || pm->fill==2)
  475.       gdImageFilledRectangle(image,x1,y1,x2,y2,pm->color[0]);
  476.   }
  477.   else
  478.     gdImageRectangle(image,x1,y1,x2,y2,pm->color[0]);
  479.   if(tikz_file)
  480.     fprintf(tikz_file,
  481.       "\\draw\[%s] (%i,%i) rectangle (%i,%i);\n",
  482.       tikz_options(pm->color[0],pm->fill),x1,flip(y1),x2,flip(y2));
  483.   if(vimg_enable) vimg_rect(scale_buf[0],scale_buf[1],scale_buf[2],scale_buf[3]);
  484. }
  485.  
  486. /* square */
  487. void obj_square(objparm *pm)
  488. {
  489.   int w,h;
  490.   scale(pm->pd,pm->p,1);
  491.   w=rint(pm->pd[2]); h=rint(pm->pd[2]);
  492.   if(pm->fill!=0){
  493.     if(pm->fill==-1 || pm->fill==2){
  494.       myDashedLine(image,pm->p[0],pm->p[1],pm->p[0]+w,pm->p[1],pm->color[0]);
  495.       myDashedLine(image,pm->p[0]+w,pm->p[1],pm->p[0]+w,pm->p[1]+h,pm->color[0]);
  496.       myDashedLine(image,pm->p[0]+w,pm->p[1]+h,pm->p[0],pm->p[1]+h,pm->color[0]);
  497.       myDashedLine(image,pm->p[0],pm->p[1]+h,pm->p[0],pm->p[1],pm->color[0]);
  498.     }
  499.     if(pm->fill==1 || pm->fill==2)
  500.     gdImageFilledRectangle(image,pm->p[0],pm->p[1],
  501.         pm->p[0]+w,pm->p[1]+h,pm->color[0]);
  502.   }
  503.   else
  504.     gdImageRectangle(image,pm->p[0],pm->p[1],
  505.                pm->p[0]+w,pm->p[1]+h,pm->color[0]);
  506.   if(tikz_file){
  507.     fprintf(tikz_file,
  508.       "\\draw\[%s] (%i,%i) rectangle (%i,%i);\n",
  509.       tikz_options(pm->color[0],pm->fill),pm->p[0],flip(pm->p[1]),
  510.       pm->p[0]+w,flip(pm->p[1]+h));
  511.   }
  512.   if(vimg_enable) vimg_rect(scale_buf[0],scale_buf[1],
  513.                       scale_buf[0]+pm->pd[2],scale_buf[1]+pm->pd[2]);
  514. }
  515.  
  516. /* triangle */
  517. void obj_triangle(objparm *pm)
  518. {
  519.   scale(pm->pd,pm->p,3);
  520.   int i,n=2;
  521.   if(pm->fill!=0){
  522.     if(pm->fill==-1 || pm->fill==2){
  523.       for(i=0;i<n;i+=1)
  524.         myDashedLine(image,pm->p[2*i],pm->p[2*i+1],pm->p[2*i+2],pm->p[2*i+3],pm->color[0]);
  525.       myDashedLine(image,pm->p[2*n],pm->p[2*n+1],pm->p[0],pm->p[1],pm->color[0]);
  526.     }
  527.     if(pm->fill==1 || pm->fill==2)
  528.       gdImageFilledPolygon(image,(gdPointPtr) pm->p,3,pm->color[0]);
  529.   }
  530.   else
  531.     gdImagePolygon(image,(gdPointPtr) pm->p,3,pm->color[0]);
  532.   if (tikz_file)
  533.     fprintf(tikz_file, "\\draw\[%s] (%i, %i) -- (%i, %i) -- (%i, %i) -- cycle;\n",
  534.       tikz_options(pm->color[0],pm->fill),pm->p[0],flip(pm->p[1]),pm->p[2],flip(pm->p[3]),pm->p[4],flip(pm->p[5]));
  535.   if(vimg_enable) vimg_polyline(scale_buf,3,1);
  536. }
  537.  
  538. /* polygon */
  539. void obj_poly(objparm *pm)
  540. {
  541.   int cnt,i;
  542.   cnt=(pm->pcnt)/2;
  543.   scale(pm->pd,pm->p,cnt);
  544.   if(pm->fill==1 || pm->fill==2)
  545.     gdImageFilledPolygon(image,(gdPointPtr) pm->p,cnt,pm->color[0]);
  546.   else
  547.     gdImagePolygon(image,(gdPointPtr) pm->p,cnt,pm->color[0]);
  548.   if(tikz_file){
  549.     fprintf(tikz_file,"\\draw\[%s] (%i,%i) --",tikz_options(pm->color[0],pm->fill),pm->p[0],flip(pm->p[1]));
  550.     for (i= 1; i<cnt; i++)
  551.         fprintf(tikz_file,"(%i,%i)--", pm->p[2*i],flip(pm->p[2*i+1]));
  552.     fprintf(tikz_file," cycle;\n");
  553.   }
  554.   if(vimg_enable) vimg_polyline(scale_buf,cnt,1);
  555. }
  556.  
  557. /* rays */
  558. void obj_rays(objparm *pm)
  559. {
  560.   int i, n;
  561.   n=(pm->pcnt)/2;
  562.   scale(pm->pd,pm->p,n);
  563.   if (tikz_file)
  564.     fprintf(tikz_file, "\\draw\[%s]",tikz_options(pm->color[0],pm->fill));
  565.   for(i=2;i<2*n;i+=2) {
  566.     if(pm->fill==-1 || pm->fill==2)
  567.       myDashedLine(image,pm->p[0],pm->p[1],pm->p[i],pm->p[i+1],pm->color[0]);
  568.     else
  569.       gdImageLine(image,pm->p[0],pm->p[1],pm->p[i],pm->p[i+1],pm->color[0]);
  570.     if (tikz_file)
  571.       fprintf(tikz_file, "(%i,%i) -- (%i,%i)",
  572.         pm->p[0],flip(pm->p[1]),pm->p[i],flip(pm->p[i+1]));
  573.     if(vimg_enable) vimg_line(scale_buf[0],scale_buf[1],
  574.                         scale_buf[i],scale_buf[i+1]);
  575.   }
  576.   if (tikz_file) fprintf(tikz_file,";\n");
  577. }
  578. /* crosshair */
  579. void obj_crosshair(objparm *pm)
  580. {
  581.   scale(pm->pd,pm->p,2);
  582.   gdImageLine(image,pm->p[0]+width2,pm->p[1]+width2,pm->p[0]-width2,pm->p[1]-width2,pm->color[0]);
  583.   gdImageLine(image,pm->p[0]-width2,pm->p[1]+width2,pm->p[0]+width2,pm->p[1]-width2,pm->color[0]);
  584.   if (tikz_file){
  585.     fprintf(tikz_file, "\\draw\[%s] (%i, %i) -- (%i, %i);\n",
  586.       tikz_options(pm->color[0],pm->fill),pm->p[0]+width2,flip(pm->p[1]+width2),pm->p[0]-width2,flip(pm->p[1]-width2));
  587.     fprintf(tikz_file, "\\draw\[%s] (%i, %i) -- (%i, %i);\n",
  588.       tikz_options(pm->color[0],pm->fill),pm->p[0]-width2,flip(pm->p[1]+width2),pm->p[0]+width2,flip(pm->p[1]-width2));
  589.   }
  590. }
  591. /* crosshairs */
  592.  
  593. void obj_crosshairs(objparm *pm)
  594. {
  595.   int i, n;
  596.   n=(pm->pcnt)/2;
  597.   scale(pm->pd,pm->p,n);
  598.   if (tikz_file)
  599.       fprintf(tikz_file, "\\draw\[%s]",tikz_options(pm->color[0],pm->fill));
  600.   for(i=0;i<2*n;i+=2) {
  601.     gdImageLine(image,pm->p[i]+width2,pm->p[i+1]+width2,pm->p[i]-width2,pm->p[i+1]-width2,pm->color[0]);
  602.     gdImageLine(image,pm->p[i]-width2,pm->p[i+1]+width2,pm->p[i]+width2,pm->p[i+1]-width2,pm->color[0]);
  603.     if (tikz_file){
  604.       fprintf(tikz_file, "(%i, %i) -- (%i, %i) ",
  605.         pm->p[i]+width2,flip(pm->p[i+1]+width2),pm->p[i]-width2,flip(pm->p[i+1]-width2));
  606.       fprintf(tikz_file, "(%i, %i) -- (%i, %i)",
  607.         pm->p[i]-width2,flip(pm->p[i+1]+width2),pm->p[i]+width2,flip(pm->p[i+1]-width2));
  608.     }
  609.   }
  610.   if (tikz_file) fprintf(tikz_file,";\n");
  611. }
  612.  
  613.  
  614. /* segments */
  615. void obj_lines(objparm *pm)
  616. {
  617.  int i, n;
  618.  n=(pm->pcnt)/2;
  619.  scale(pm->pd,pm->p,n);
  620.  if (tikz_file){
  621.    fprintf(tikz_file,"\\draw\[%s] (%i,%i)",
  622.      tikz_options(pm->color[0],pm->fill),pm->p[0],flip(pm->p[1]));
  623.    for(i=2;i<2*n;i+=2)
  624.      fprintf(tikz_file, "--(%i,%i)",pm->p[i],flip(pm->p[i+1]));
  625.    fprintf(tikz_file, ";\n");
  626.  }
  627.  for(i=2;i<2*n;i+=2){
  628.   if(pm->fill==-1 || pm->fill==2 )
  629.     myDashedLine(image,pm->p[i-2],pm->p[i-1],pm->p[i],pm->p[i+1],pm->color[0]);
  630.   else
  631.     gdImageLine(image,pm->p[i-2],pm->p[i-1],pm->p[i],pm->p[i+1],pm->color[0]);
  632.   if(vimg_enable) vimg_polyline(scale_buf,n,0);
  633.   }
  634. }
  635. /*segments from x1,y1 to x2,y2, x3,y3 to x4,y4, etc */
  636. void obj_segments(objparm *pm)
  637. {
  638.   int i, n;
  639.   n=(pm->pcnt)/2;
  640.   scale(pm->pd,pm->p,n);
  641.   if (tikz_file)
  642.     fprintf(tikz_file, "\\draw\[%s]",tikz_options(pm->color[0],pm->fill));
  643.   for(i=0;i<2*n;i+=4){
  644.     if(pm->fill==-1 || pm->fill==2)
  645.       myDashedLine(image,pm->p[i],pm->p[i+1],pm->p[i+2],pm->p[i+3],pm->color[0]);
  646.     else
  647.       gdImageLine(image,pm->p[i],pm->p[i+1],pm->p[i+2],pm->p[i+3],pm->color[0]);
  648.     if (tikz_file)
  649.       fprintf(tikz_file, "(%i,%i)--(%i,%i)",
  650.        pm->p[i],flip(pm->p[i+1]),pm->p[i+2],flip(pm->p[i+3]));
  651.   }
  652.   if(tikz_file) fprintf(tikz_file, ";\n");
  653. }
  654. /* for geodesics, return data for obj_hypgeods + a point on the arc
  655.   to be used for filling hyperbolic triangle*/
  656. int hypgeodaux(double *q, double* res)
  657. {
  658.   double r,cx,cy,a1,a2,a3,tmp,
  659.     nx = -q[0]*q[2]*q[2]+(q[0]*q[0]+q[1]*q[1]+1)*q[2]-q[0]*q[3]*q[3]-q[0],
  660.     ny = -q[1]*q[2]*q[2]-q[1]*q[3]*q[3]+(q[0]*q[0]+q[1]*q[1]+1)*q[3]-q[1],
  661.     dy = -2*q[1]*q[2]+2*q[0]*q[3];
  662.   if (dy*dy*1e4 < nx*nx+ny*ny){
  663.     res[5]=(q[0]+q[2])/2;
  664.     res[6]=(q[1]+q[3])/2;
  665.     return 0;}
  666.   cx = ny/dy; cy=-nx/dy;
  667.   a1 = atan2(q[1]-cy, q[0]-cx);
  668.   a2 = atan2(q[3]-cy, q[2]-cx);
  669.   if (fabs(a2-a1)>M_PI){if(a1<a2) a1+=2*M_PI; else a2+=2*M_PI;};
  670.   a3 = (a1+a2)/2;
  671.   if(a1>a2) {tmp=a1; a1=a2; a2=tmp;}
  672.   r = sqrt(cx*cx+cy*cy-1);
  673.   res[0]=cx;
  674.   res[1]=cy;
  675.   res[2]=r;
  676.   res[3]=a1;
  677.   res[4]=a2;
  678.   res[5]=cx+r*cos(a3);
  679.   res[6]=cy+r*sin(a3);
  680.   return 1;
  681. }
  682. /* hyperbolic geodesics from x1,y1 to x2,y2, x3,y3 to x4,y4 in Poincaré disk */
  683.  
  684. void obj_hypgeods(objparm *pm)
  685. {
  686.   int i, *qq=pm->p;
  687.   double rx,ry,res[7];
  688.   for (i = 0; i < pm->pcnt; i+=4)
  689.     if (hypgeodaux(pm->pd+i,res)){
  690.       scale2(res[2],res[2],&rx,&ry);
  691.       scale(res,qq,1);
  692.       myGdImageArc(image,qq[0],qq[1],rx,ry,res[3]/DEG,res[4]/DEG,pm->color[0]);
  693.     }
  694.     else {
  695.       scale(pm->pd+i,qq,2);
  696.       gdImageLine(image,qq[0],qq[1],qq[2],qq[3],pm->color[0]);
  697.     }
  698. }
  699. /* hyperbolique triangle, can be filled */
  700. void obj_hyptriangle(objparm *pm)
  701. {
  702.   double rx,ry,data[8],res[7];
  703.   int i, *qq=pm->p;
  704.   for(i=0; i<8; i++) data[i]=pm->pd[i%6];
  705.   for(i=0; i<3; i++)
  706.     if(hypgeodaux(data+2*i,res)) {
  707.       scale2(res[2],res[2],&rx,&ry);
  708.       scale(res,qq,1);
  709.       myGdImageArc(image,qq[0],qq[1],rx,ry,res[3]/DEG,res[4]/DEG,pm->color[0]);
  710.     }
  711.     else {
  712.       scale(data+2*i,qq,2);
  713.       gdImageLine(image,qq[0],qq[1],qq[2],qq[3],pm->color[0]);
  714.     }
  715.   if(pm->fill){
  716.     data[0]=res[5];
  717.     data[1]=res[6];
  718.     hypgeodaux(data,res);
  719.     scale(res+5,qq,1);
  720.     patchgdImageFill(image,qq[0],qq[1],pm->color[0]);
  721.   }
  722. }
  723. /* hyperbolic geodesic through two points */
  724. void obj_hypline(objparm *pm)
  725. {
  726.   double rx,ry,res[7];
  727.   if (hypgeodaux(pm->pd,res)){
  728.       double alpha=atan(1/res[2]), beta = atan2(res[0], res[1]);
  729.       scale(pm->pd, pm->p, 1);
  730.       scale2(res[2],res[2],&rx,&ry);
  731.       myGdImageArc(image, pm->p[0], pm->p[1], rx, ry, 180+(beta-alpha)/DEG, 180+(beta+alpha)/DEG,pm->color[0]);
  732.   }
  733.   else {
  734.     double gamma = atan2(pm->pd[0],pm->pd[1]);
  735.     res[0]=cos(gamma);
  736.     res[1]=sin(gamma);
  737.     res[2]=-cos(gamma);
  738.     res[3]=-sin(gamma);
  739.     scale(res,pm->p,2);
  740.     gdImageLine(image,pm->p[0],pm->p[1],pm->p[2],pm->p[3],pm->color[0]);
  741.   }
  742. }
  743.  
  744. void obj_hyplines(objparm *pm)
  745. {
  746.   double rx,ry,res[7];
  747.   int i;
  748.   for (i=0; i < pm->pcnt; i+=4)
  749.     if (hypgeodaux(pm->pd+i,res)) {
  750.       double alpha=atan(1/res[2]), beta = atan2(res[0],res[1]);
  751.       scale(pm->pd+i,pm->p,1);
  752.       scale2(res[2],res[2],&rx,&ry);
  753.       myGdImageArc(image,pm->p[0],pm->p[1],rx,ry,180+(beta-alpha)/DEG,180+(beta+alpha)/DEG,pm->color[0]);
  754.     }
  755.     else {
  756.       double gamma = atan2(pm->pd[i],pm->pd[i+1]);
  757.       res[0]=cos(gamma);
  758.       res[1]=sin(gamma);
  759.       res[2]=-cos(gamma);
  760.       res[3]=-sin(gamma);
  761.       scale(res,pm->p,2);
  762.       gdImageLine(image,pm->p[0],pm->p[1],pm->p[2],pm->p[3],pm->color[0]);
  763.     }
  764. }
  765. /* segments */
  766. void obj_dlines(objparm *pm)
  767. {
  768.  int i, n;
  769.  n=(pm->pcnt)/2;
  770.  scale(pm->pd,pm->p,n);
  771.  if (tikz_file){
  772.    fprintf(tikz_file,"\\draw\[%s] (%i,%i)",
  773.      tikz_options(pm->color[0],pm->fill),pm->p[0],flip(pm->p[1]));
  774.    for(i=2;i<2*n;i+=2)
  775.      fprintf(tikz_file, "--(%i,%i)",pm->p[i],flip(pm->p[i+1]));
  776.    fprintf(tikz_file, ";\n");
  777.  }
  778.  for(i=2;i<2*n;i+=2)
  779.    myDashedLine(image,pm->p[i-2],pm->p[i-1],pm->p[i],pm->p[i+1],pm->color[0]);
  780.  if(vimg_enable) vimg_polyline(scale_buf,n,0);
  781. }
  782.  
  783. /* points */
  784. void obj_points(objparm *pm)
  785. {
  786.  int i, n;
  787.  n=(pm->pcnt)/2;
  788.  scale(pm->pd,pm->p,n);
  789.  for(i=0;i<2*n;i+=2) gdImageSetPixel(image,pm->p[i],pm->p[i+1],pm->color[0]);
  790.  if(tikz_file) {
  791.    fprintf(tikz_file, "\\draw\[%s]", tikz_options(pm->color[0],pm->fill));
  792.    for(i=0;i<2*n;i+=2)
  793.      fprintf(tikz_file,"(%i,%i)circle(1pt)",pm->p[i],flip(pm->p[i+1]));
  794.    fprintf(tikz_file,";\n");
  795.  }
  796. }
  797.  
  798. /* lattice.
  799.  * x0,y0,xv1,yv1,xv2,yv2,n1,n2,color */
  800. void obj_lattice(objparm *pm)
  801. {
  802.  int n1,n2,i1,i2,xi1,yi1,xi2,yi2;
  803.  double xv1,xv2,yv1,yv2;
  804.  n1=pm->pd[6];n2=pm->pd[7]; if(n1<0 || n2<0) return;
  805.  if(n1>256) n1=256;
  806.  if(n2>256) n2=256;
  807.  scale(pm->pd,pm->p,1);
  808.  scale2(pm->pd[2],pm->pd[3],&xv1,&yv1);
  809.  scale2(pm->pd[4],pm->pd[5],&xv2,&yv2);
  810.  if(tikz_file)
  811.    fprintf(tikz_file,"\\draw[%s]", tikz_options(pm->color[0],0));
  812.  for(i1=0;i1<n1;i1++) {
  813.    xi1=rint(i1*xv1)+pm->p[0]; yi1=rint(i1*yv1)+pm->p[1];
  814.    for(i2=0;i2<n2;i2++) {
  815.      xi2=i2*xv2+xi1;yi2=i2*yv2+yi1;
  816.      gdImageSetPixel(image,xi2,yi2,pm->color[0]);
  817.      if(tikz_file) fprintf(tikz_file,"(%i,%i)circle(1pt)",xi2,flip(yi2));
  818.    }
  819.  }
  820.  if(tikz_file) fprintf(tikz_file,";\n");
  821. }
  822.  
  823. /* arc */
  824. /*trouble with gdImageArc as the angles are of type int */
  825. void myGdImageArc(gdImagePtr im, int cx, int cy, int rx, int ry, double t1, double t2, int color)
  826. {
  827.   int i,nb;
  828.   double dt, r=rx>ry?rx:ry;
  829.   t2-=t1;
  830.   t2-=360*floor(t2/360);
  831.   if(t2==0) t2=360;
  832.   t1*=DEG;t2*=DEG;
  833.   nb = r*t2/3;
  834.   dt = t2/nb;
  835.   for (i=0; i<nb; ++i,t1+=dt)
  836.     gdImageLine(im,cx+rx*cos(t1),cy+ry*sin(t1),
  837.     cx+rx*cos(t1+dt),cy+ry*sin(t1+dt),color);
  838. }
  839.  
  840. void obj_arc(objparm *pm)
  841. {
  842.   scale(pm->pd,pm->p,1);
  843.   int rx=pm->pd[2]*xscale/2, ry=pm->pd[3]*yscale/2;
  844.   myGdImageArc(image, pm->p[0], pm->p[1], rx, ry, pm->pd[4],pm->pd[5],pm->color[0]);
  845.   if(tikz_file){
  846.     pm->pd[5]-=pm->pd[4];
  847.     pm->pd[5]-=360*floor(pm->pd[5]/360);
  848.     pm->pd[5]+=pm->pd[4];
  849.     fprintf(tikz_file,
  850. /*    "\\draw\[%s, domain=%f:%f] plot ({%i+%f*cos(\\x)}, {%i+%f*sin(\\x)});\n",*/
  851.       "\\draw \[%s] (%i,%i) arc (%f:%f:%i and %i);\n",
  852.         tikz_options(pm->color[0],pm->fill),
  853.         (int)rint(pm->p[0]+rx*cos(DEG*pm->pd[4])),
  854.         flip((int)rint(pm->p[1]+ry*sin(DEG*pm->pd[4]))),
  855.         pm->pd[4],pm->pd[5],rx,-ry);
  856.   }
  857.   /*FIXME echelle mauvaise*/
  858.   if(vimg_enable) vimg_arc(scale_buf[0],scale_buf[1],
  859.                      0.5*pm->pd[2],0.5*pm->pd[3],pm->pd[4],pm->pd[5]);
  860. }
  861.  
  862. void obj_angle(objparm *pm)
  863. {
  864. /* dans la doc: mettre les angles dans l'ordre ???? demander qu'ils soient entre 0 et 360 */
  865.   double dpts[6];
  866.   int pts[6], wx=rint(2*pm->pd[2]*xscale), wy=rint(2*pm->pd[2]*yscale);
  867.   dpts[0]=pm->pd[0];
  868.   dpts[1]=pm->pd[1];
  869.   dpts[2]=pm->pd[0]+pm->pd[2]*cos(DEG*pm->pd[3]);
  870.   dpts[3]=pm->pd[1]+pm->pd[2]*sin(DEG*pm->pd[3]);
  871.   dpts[4]=pm->pd[0]+pm->pd[2]*cos(DEG*pm->pd[4]);
  872.   dpts[5]=pm->pd[1]+pm->pd[2]*sin(DEG*pm->pd[4]);
  873.   scale(dpts, pts, 3);
  874.     gdImageLine(image,pts[0],pts[1],pts[2],pts[3],pm->color[0]);
  875.     gdImageLine(image,pts[0],pts[1],pts[4],pts[5],pm->color[0]);
  876.   gdImageArc(image,pts[0],pts[1],wx,wy,pm->pd[3],pm->pd[4],pm->color[0]);
  877.   if(tikz_file) {
  878.     fprintf(tikz_file, "\\draw\[%s, domain=%f:%f] plot ({%i+%f*cos(\\x)}, {%i+%f*sin(\\x)});\n",
  879.       tikz_options(pm->color[0],pm->fill),pm->pd[3],pm->pd[4],pts[0],wx*0.5,flip(pts[1]),wy*(-0.5));
  880.     fprintf(tikz_file, "\\draw\[%s] (%i, %i) -- (%i, %i) (%i, %i) -- (%i, %i);\n",
  881.         tikz_options(pm->color[0],pm->fill),pts[0],flip(pts[1]),pts[2],flip(pts[3]),pts[0],flip(pts[1]),pts[4],flip(pts[5]));
  882.   }
  883. }
  884. /* Ellipse: centre 0,1, width 2, hight 3, color 4,5,6 */
  885. void obj_ellipse(objparm *pm)
  886. {
  887.   scale(pm->pd,pm->p,1);
  888.   pm->p[2]=pm->pd[2]*xscale; pm->p[3]=pm->pd[3]*yscale;
  889.   if(pm->fill==1 || pm->fill==2) {
  890.     gdImageArc(image,pm->p[0],pm->p[1],pm->p[2],pm->p[3],0,360,
  891.              color_bounder);
  892.     patchgdImageFillToBorder(image,pm->p[0],pm->p[1],
  893.                     color_bounder,pm->color[0]);
  894.   }
  895.   gdImageArc(image,pm->p[0],pm->p[1],rint(pm->p[2]),rint(pm->p[3]),0,360,pm->color[0]);
  896.   if(tikz_file) fprintf(tikz_file,"\\draw\[%s] (%i,%i) ellipse (%i and %i);\n\n",
  897.     tikz_options(pm->color[0],pm->fill),pm->p[0],flip(pm->p[1]),(int) rint(0.5*pm->p[2]),(int) rint(0.5*pm->p[3]));
  898.   if(vimg_enable) vimg_ellipse(scale_buf[0],scale_buf[1],0.5*pm->pd[2],0.5*pm->pd[3]);
  899. }
  900.  
  901. /* Circle radius in pixels*/
  902. void obj_circle(objparm *pm)
  903. {
  904.   scale(pm->pd,pm->p,1);
  905.   pm->p[2]=rint(pm->pd[2]); pm->p[3]=rint(pm->pd[2]);
  906.   if(pm->fill==1 || pm->fill==2) {
  907.     gdImageArc(image,pm->p[0],pm->p[1],pm->p[2],pm->p[3],0,360,
  908.              color_bounder);
  909.     patchgdImageFillToBorder(image,pm->p[0],pm->p[1],
  910.                     color_bounder,pm->color[0]);
  911.   }
  912.   gdImageArc(image,pm->p[0],pm->p[1],pm->p[2],pm->p[3],0,360,pm->color[0]);
  913.   if(tikz_file) fprintf(tikz_file, "\\draw\[%s] (%i, %i) circle (%i);\n",
  914.       tikz_options(pm->color[0],pm->fill),pm->p[0],flip(pm->p[1]),(int) rint(pm->pd[2]/2.));
  915. }
  916. /* Circle, radius in xrange */
  917. void obj_circles(objparm *pm)
  918. {
  919.   int i, n;
  920.   n=(pm->pcnt)/3;
  921.   for (i=0; i<n; ++i) scale(pm->pd+3*i, pm->p+3*i, 1);
  922.   if (tikz_file) fprintf(tikz_file,"\\draw\[%s]",tikz_options(pm->color[0],pm->fill));
  923.   for(i=0;i<3*n;i+=3){
  924.     if(pm->fill==1 || pm->fill==2) {
  925.       gdImageArc(image,pm->p[i],pm->p[i+1],
  926.         rint(2*pm->pd[i+2]*xscale),rint(2*pm->pd[i+2]*xscale),0,360,
  927.           color_bounder);
  928.       patchgdImageFillToBorder(image,pm->p[i],pm->p[i+1],
  929.                     color_bounder,pm->color[0]);
  930.     }
  931.     gdImageArc(image,pm->p[i],pm->p[i+1],
  932.         rint(2*pm->pd[i+2]*xscale),rint(2*pm->pd[i+2]*xscale),0,360,pm->color[0]);
  933.     if (tikz_file){
  934.       fprintf(tikz_file, "(%i, %i) circle (%i and %i)",
  935.         pm->p[i],flip(pm->p[i+1]),(int) rint(pm->pd[i+2]*xscale),(int) rint(pm->pd[i+2]*xscale));
  936.     }
  937.   }
  938.   if (tikz_file) fprintf(tikz_file,";\n");
  939. }
  940. /* ellipse, width in xrange,height in yrange */
  941. void obj_ellipses(objparm *pm)
  942. {
  943.   int i, n;
  944.   n=(pm->pcnt)/4;
  945.   for (i=0; i<n; ++i) scale(pm->pd+4*i, pm->p+4*i, 1);
  946.   if (tikz_file) fprintf(tikz_file,"\\draw\[%s]",tikz_options(pm->color[0],pm->fill));
  947.   for(i=0;i<4*n;i+=4){
  948.     if(pm->fill==1 || pm->fill==2){
  949.       gdImageArc(image,pm->p[i],pm->p[i+1],
  950.         rint(pm->pd[i+2]*xscale),rint(pm->pd[i+2]*xscale),0,360,
  951.         color_bounder);
  952.       patchgdImageFillToBorder(image,pm->p[i],pm->p[i+1],
  953.         color_bounder,pm->color[0]);
  954.     } else
  955.       gdImageArc(image,pm->p[i],pm->p[i+1],rint(pm->pd[i+2]*xscale),rint(pm->pd[i+3]*yscale),0,360,pm->color[0]);
  956.     if (tikz_file){
  957.       fprintf(tikz_file, "(%i, %i) circle (%i and %i)",
  958.         pm->p[i],flip(pm->p[i+1]),(int) rint(pm->pd[i+2]*xscale/2),(int) rint(pm->pd[i+3]*yscale/2));
  959.     }
  960.   }
  961.   if (tikz_file) fprintf(tikz_file,";\n");
  962. }
  963. typedef enum tikz_fill_options {grid, dot, hatch, diamond} tfo;
  964. int compar(const void *a, const void *b) {return *(int *)b - *(int *)a;}
  965.  
  966. static void tikz_fill(int x, int y, int nx, int ny, int wall, int the_color, tfo opt)
  967. {
  968.   int numpix=sizex*sizey, top = 0, a, c, cy, lx=0, rx=0, seed, nt,ct=0;
  969.   int *stack = xmalloc(numpix*sizeof(int));
  970.   int *stack2 = xmalloc(numpix*sizeof(int));
  971.   if (nx==0) nx=1;
  972.   if (ny==0) ny=1;
  973.   nt = abs(nx*ny);
  974.   seed = gdImageGetPixel(image,x,y);
  975.   char *check = xmalloc(numpix);
  976.   while (numpix--) check[numpix]=0;
  977.   a=x+y*sizex;
  978.   check[a]=1;
  979.   stack[top++]=a;
  980.   while(top)
  981.     {
  982.       a = stack[--top];
  983.       x = a%sizex; y = a/sizex;
  984.       c = gdImageGetPixel(image,x,y);
  985.       if ((wall && (c == wall))||(!wall && c!=seed))
  986.         continue;
  987.       if (x>0 && !check[a-1]) {check[a-1]=1; stack[top++]=a-1;}
  988.       if (x+1<sizex && !check[a+1]) {check[a+1]=1; stack[top++]=a+1;}
  989.       if (y>0 && !check[a-sizex]) {check[a-sizex]=1; stack[top++]=a-sizex;}
  990.       if (y+1<sizey && !check[a+sizex]) {check[a+sizex]=1; stack[top++]=a+sizex;}
  991.       switch(opt)
  992.   {
  993.   case grid: if(((x%nx)!=(nx/2))&&((y%ny)!=(ny/2))) continue; break;
  994.   case dot: if(x%nx!=nx/2||y%ny!=ny/2) continue; break;
  995.   case hatch: if((ny*x-nx*y+nt)%nt) continue; break;
  996.   case diamond: if((ny*x+nx*y)%nt && (ny*x-nx*y+nt)%nt)continue;
  997.   default: break;
  998.   }
  999.       stack2[ct++]=a;
  1000.     }
  1001.   free(stack);
  1002.   free(check);
  1003.   if (ct==0) return;
  1004.   fprintf(tikz_file,"\\draw\[%s]", tikz_options(the_color,1));
  1005.   qsort(stack2, ct, sizeof(int), compar);
  1006.   cy=-1;
  1007.   while (ct--)
  1008.     {
  1009.       a = stack2[ct];
  1010.       x = a%sizex; y = a/sizex;
  1011.       if (y != cy || x != rx+1)
  1012.       {
  1013.   if (cy >= 0)
  1014.     fprintf(tikz_file,
  1015.       "(%.2f,%.2f)rectangle(%.2f,%.2f)",
  1016.       lx-0.0,flip(cy-0.0),rx+0.0,flip(cy+0.0));
  1017.   cy = y;
  1018.   lx = rx = x;
  1019.       }
  1020.       else
  1021.   rx++;
  1022.     }
  1023.   fprintf(tikz_file,
  1024.     "(%.2f,%.2f)rectangle(%.2f,%.2f);\n",
  1025.      lx-0.0,flip(cy-0.0),rx+0.0,flip(cy+0.0));
  1026.   free(stack2);
  1027. }
  1028. void obj_dashed(objparm *pm){ dashed=1;}
  1029. void obj_filled(objparm *pm){ filled=1;}
  1030. void obj_noreset(objparm *pm){ noreset=1;}
  1031. void obj_reset(objparm *pm){ noreset=0;dashed=0;filled=0;}
  1032.  
  1033. /* flood fill */
  1034. void obj_fill(objparm *pm)
  1035. {
  1036.   scale(pm->pd,pm->p,1);
  1037.   if (tikz_file) tikz_fill(pm->p[0],pm->p[1],1,1,0,pm->color[0],grid);
  1038.   patchgdImageFill(image,pm->p[0],pm->p[1],pm->color[0]);
  1039. }
  1040.  
  1041. /* flood fill to border*/
  1042. void obj_fillb(objparm *pm)
  1043. {
  1044.   scale(pm->pd,pm->p,1);
  1045.   if (tikz_file) tikz_fill(pm->p[0],pm->p[1],1,1,pm->color[0],pm->color[1],grid);
  1046.   patchgdImageFillToBorder(image,pm->p[0],pm->p[1],pm->color[0],pm->color[1]);
  1047. }
  1048.  
  1049. gdImagePtr himg;
  1050.  
  1051. int makehatchimage(int x, int y, int px, int py, int col)
  1052. {
  1053.   int c1,c2,r,g,b;
  1054.   gdImagePtr saveimg;
  1055.   himg=gdImageCreate(x,y);
  1056.   c1=gdImageGetPixel(image,px,py);
  1057.   r=gdImageRed(image,c1); g=gdImageGreen(image,c1); b=gdImageBlue(image,c1);
  1058.   if(r>=255) r--; else r++; if(g>=255) g--; else g++; if(b>=255) b--; else b++;
  1059.   c1=gdImageColorAllocate(himg,r,g,b);
  1060.   r=gdImageRed(image,col); g=gdImageGreen(image,col); b=gdImageBlue(image,col);
  1061.   c2=gdImageColorAllocate(himg,r,g,b);
  1062.   if(width>1) {
  1063.     savew=-1; saveimg=image;
  1064.     image=himg; c2=widthcolor(width,c2); image=saveimg;
  1065.     c2=gdBrushed; savew=-1;
  1066.   }
  1067.   return c2;
  1068. }
  1069.  
  1070. /* flood fill with hatching */
  1071. void obj_hatchfill(objparm *pm)
  1072. {
  1073.   int nx,ny,ax,ay, dir, c;
  1074.   scale(pm->pd,pm->p,1);
  1075.   nx=pm->pd[2]; ny=pm->pd[3]; ax=abs(nx); ay=abs(ny);
  1076.   if(nx==0 && ny==0) {fly_error("bad displacement vector"); return;}
  1077.   if((nx>0 && ny>0) || (nx<0 && ny<0)) dir=1; else dir=-1;
  1078.   if(ax==0) {ax=100; dir=2;}
  1079.   if(ay==0) {ay=100; dir=3;}
  1080.   c=makehatchimage(ax,ay,pm->p[0],pm->p[1],pm->color[0]);
  1081.   switch(dir) {
  1082.     case -1: {
  1083.       gdImageLine(himg,0,ay-1,ax-1,0,c);
  1084.       if(width>1) {
  1085.         gdImageLine(himg,-ax,ay-1,-1,0,c);
  1086.         gdImageLine(himg,ax,ay-1,2*ax-1,0,c);
  1087.         gdImageLine(himg,0,-1,ax-1,-ay,c);
  1088.         gdImageLine(himg,0,2*ay-1,ax-1,ay,c);
  1089.       }
  1090.       break;
  1091.     }
  1092.     case 1: {
  1093.       gdImageLine(himg,0,0,ax-1,ay-1,c);
  1094.       if(width>1) {
  1095.         gdImageLine(himg,-ax,0,-1,ay-1,c);
  1096.         gdImageLine(himg,ax,0,2*ax-1,ay-1,c);
  1097.         gdImageLine(himg,0,-ay,ax-1,-1,c);
  1098.         gdImageLine(himg,0,ay,ax-1,2*ay-1,c);
  1099.       }
  1100.       break;
  1101.     }
  1102.     case 2: gdImageLine(himg,0,ay/2,ax-1,ay/2,c); break;
  1103.     case 3: gdImageLine(himg,ax/2,0,ax/2,ay-1,c); break;
  1104.   }
  1105.   if (tikz_file) tikz_fill(pm->p[0],pm->p[1],nx,ny,0,pm->color[0],hatch);
  1106.   gdImageSetTile(image,himg);
  1107.   patchgdImageFill(image,pm->p[0],pm->p[1],gdTiled);
  1108.   gdImageDestroy(himg);
  1109.   if(tiled) gdImageSetTile(image,tileimg);
  1110. }
  1111.  
  1112. /* flood fill with grid */
  1113. void obj_gridfill(objparm *pm)
  1114. {
  1115.   int nx,ny, c;
  1116.   scale(pm->pd,pm->p,1);
  1117.   nx=pm->pd[2]; ny=pm->pd[3]; nx=abs(nx); ny=abs(ny);
  1118.   if(nx==0 && ny==0) {fly_error("bad grid size"); return;}
  1119.   c=makehatchimage(nx,ny,pm->p[0],pm->p[1],pm->color[0]);
  1120.   if (tikz_file) tikz_fill(pm->p[0],pm->p[1],nx,ny,0,pm->color[0],grid);
  1121.   gdImageLine(himg,0,ny/2,nx-1,ny/2,c); gdImageLine(himg,nx/2,0,nx/2,ny-1,c);
  1122.   gdImageSetTile(image,himg);
  1123.   patchgdImageFill(image,pm->p[0],pm->p[1],gdTiled);
  1124.   gdImageDestroy(himg);
  1125.   if(tiled) gdImageSetTile(image,tileimg);
  1126. }
  1127.  
  1128. /* flood fill with double hatching */
  1129. void obj_diafill(objparm *pm)
  1130. {
  1131.   int nx,ny, c;
  1132.   scale(pm->pd,pm->p,1);
  1133.   nx=pm->pd[2]; ny=pm->pd[3]; nx=abs(nx); ny=abs(ny);
  1134.   if(nx==0 && ny==0) {fly_error("bad grid size"); return;}
  1135.   c=makehatchimage(nx,ny,pm->p[0],pm->p[1],pm->color[0]);
  1136.   if (tikz_file) tikz_fill(pm->p[0],pm->p[1],nx,ny,0,pm->color[0],diamond);
  1137.   gdImageLine(himg,0,0,nx-1,ny-1,c); gdImageLine(himg,0,ny-1,nx-1,0,c);
  1138.   gdImageSetTile(image,himg);
  1139.   patchgdImageFill(image,pm->p[0],pm->p[1],gdTiled);
  1140.   gdImageDestroy(himg);
  1141.   if(tiled) gdImageSetTile(image,tileimg);
  1142. }
  1143.  
  1144. /* flood fill with dots */
  1145. void obj_dotfill(objparm *pm)
  1146. {
  1147.   int nx,ny, c;
  1148.   scale(pm->pd,pm->p,1);
  1149.   nx=pm->pd[2]; ny=pm->pd[3]; nx=abs(nx); ny=abs(ny);
  1150.   if(nx==0 && ny==0) {fly_error("bad grid size"); return;}
  1151.   c=makehatchimage(nx,ny,pm->p[0],pm->p[1],pm->color[0]);
  1152.   if (tikz_file) tikz_fill(pm->p[0],pm->p[1],nx,ny,0,pm->color[0],dot);
  1153.   gdImageSetPixel(himg,nx/2,ny/2,c);
  1154.   gdImageSetTile(image,himg);
  1155.   patchgdImageFill(image,pm->p[0],pm->p[1],gdTiled);
  1156.   gdImageDestroy(himg);
  1157.   if(tiled) gdImageSetTile(image,tileimg);
  1158. }
  1159.  
  1160. struct {
  1161.   char *name;
  1162.   gdFontPtr *fpt;
  1163.   char *ftikz;
  1164. } fonttab[]={
  1165.   {"tiny",  &gdFontTiny,"0.3"},
  1166.   {"small", &gdFontSmall,"0.5"},
  1167.   {"medium",&gdFontMediumBold,"0.7"},
  1168.   {"large", &gdFontLarge,"0.9"},
  1169.   {"giant", &gdFontGiant,"1"},
  1170.   {"huge",  &gdFontGiant,"1"}
  1171. };
  1172.  
  1173. #define fonttab_no (sizeof(fonttab)/sizeof(fonttab[0]))
  1174.  
  1175. /* string */
  1176. void obj_string(objparm *pm)
  1177. {
  1178.   char *pp, *pe, *p2;
  1179.   int i;
  1180.   pp=pm->str; pe=strchr(pp,','); if(pe==NULL) {
  1181.     fly_error("too_few_parms"); return;
  1182.   }
  1183.   *pe++=0; pp=find_word_start(pp); *find_word_end(pp)=0;
  1184.   pe=find_word_start(pe); strip_trailing_spaces(pe);
  1185.   if(*pp) {
  1186.     for(i=0;i<fonttab_no && strcmp(pp,fonttab[i].name)!=0; i++);
  1187.     if(i>=fonttab_no) i=1;
  1188.   }
  1189.   else i=1;
  1190.   scale(pm->pd,pm->p,1);
  1191.   if(*pe=='"') {
  1192.     p2=strchr(pe+1,'"');
  1193.     if(p2 && *(p2+1)==0) {*p2=0; pe++;}
  1194.   }
  1195.   if(pm->fill){
  1196.     gdImageStringUp(image,*(fonttab[i].fpt),pm->p[0],pm->p[1], (unsigned char*) pe,
  1197.         pm->color[0]);
  1198.     if(tikz_file) fprintf(tikz_file,"\\draw (%i,%i) node[scale=%s,rotate=90,color=%s,anchor=north west,inner sep=0pt] {\\(%s\\)};\n",
  1199.       pm->p[0],flip(pm->p[1]),fonttab[i].ftikz,tikz_color(pm->color[0]),(unsigned char*) pe);
  1200.   }
  1201.   else {
  1202.     gdImageString(image,*(fonttab[i].fpt),pm->p[0],pm->p[1], (unsigned char*) pe,
  1203.         pm->color[0]);
  1204.     if(tikz_file) fprintf(tikz_file,"\\draw (%i,%i) node[scale=%s,color=%s,anchor=north west,inner sep=0pt] {\\(%s\\)};\n",
  1205.       pm->p[0],flip(pm->p[1]),fonttab[i].ftikz,tikz_color(pm->color[0]),(unsigned char*) pe);
  1206.   }
  1207. }
  1208. /*FIXME; giant does not exist in tikz and ... samething as huge */
  1209. /* point */
  1210. void obj_point(objparm *pm)
  1211. {
  1212.   scale(pm->pd,pm->p,1);
  1213.   gdImageSetPixel(image,pm->p[0],pm->p[1],pm->color[0]);
  1214.   if(tikz_file) fprintf(tikz_file,"\\draw[%s] (%i,%i) circle (1pt);\n",
  1215.     tikz_options(pm->color[0],0),pm->p[0],flip(pm->p[1]));
  1216. }
  1217.  
  1218. /* copy an image file */
  1219. void obj_copy(objparm *pm)
  1220. {
  1221.   char *pp;
  1222.   FILE *inf;
  1223.   gdImagePtr insimg;
  1224.  
  1225.   pp=find_word_start(pm->str);*find_word_end(pp)=0;
  1226.   inf=open4read(pp);
  1227.   if(inf==NULL) {
  1228.     fly_error("file_does_not_exist"); return;
  1229.   }
  1230.   insimg=gdImageCreateFromGif(inf); fclose(inf);
  1231.   if(insimg==NULL) {
  1232.     fly_error("bad_gif"); return;
  1233.   }
  1234.   scale(pm->pd,pm->p,1);
  1235.   if(pm->pd[2]<0 && pm->pd[3]<0 && pm->pd[4]<0 && pm->pd[5]<0)
  1236.     gdImageCopy(image,insimg,pm->p[0],pm->p[1],0,0,
  1237.             insimg->sx,insimg->sy);
  1238.   else
  1239.     gdImageCopy(image,insimg,pm->p[0],pm->p[1],pm->pd[2],pm->pd[3],
  1240.             pm->pd[4]-pm->pd[2],pm->pd[5]-pm->pd[3]);
  1241.   if(tikz_file) fprintf(tikz_file,
  1242.     "\\node[anchor=north west,inner sep=0pt] at (%i,%i) {\\includegraphics\[width=.05\\linewidth]{%s}};\n",
  1243.       pm->p[0],flip(pm->p[1]),pm->str);
  1244.     gdImageDestroy(insimg);
  1245. /*FIXME position non correcte; on ne peut pas utiliser une image gif */
  1246. }
  1247.  
  1248. /* copy an image file, with resizing */
  1249. void obj_copyresize(objparm *pm)
  1250. {
  1251.   char *pp;
  1252.   FILE *inf;
  1253.   gdImagePtr insimg;
  1254.  
  1255.   pp=find_word_start(pm->str);*find_word_end(pp)=0;
  1256.   inf=open4read(pp);
  1257.   if(inf==NULL) {
  1258.     fly_error("file_not_found"); return;
  1259.   }
  1260.   insimg=gdImageCreateFromGif(inf); fclose(inf);
  1261.   if(insimg==NULL) {
  1262.     fly_error("bad_gif"); return;
  1263.   }
  1264.   scale(pm->pd+4,pm->p+4,2);
  1265.   if(pm->pd[0]<0 && pm->pd[1]<0 && pm->pd[2]<0 && pm->pd[3]<0)
  1266.     gdImageCopyResized(image,insimg,pm->p[4],pm->p[5],0,0,
  1267.                  pm->p[6]-pm->p[4]+1,pm->p[7]-pm->p[5]+1,
  1268.                  insimg->sx,insimg->sy);
  1269.   else
  1270.     gdImageCopyResized(image,insimg,pm->p[4],pm->p[5],pm->pd[0],pm->pd[1],
  1271.                  pm->p[6]-pm->p[4]+1,pm->p[7]-pm->p[5]+1,
  1272.                  pm->pd[2]-pm->pd[0]+1,pm->pd[3]-pm->pd[1]+1);
  1273.   gdImageDestroy(insimg);
  1274. }
  1275.  
  1276. /* set brush or tile */
  1277. void obj_setbrush(objparm *pm)
  1278. {
  1279.   char *pp;
  1280.   FILE *inf;
  1281.   gdImagePtr insimg;
  1282.  
  1283.   pp=find_word_start(pm->str); *find_word_end(pp)=0;
  1284.   inf=open4read(pp); if(inf==NULL) {
  1285.     fly_error("file_not_found"); return;
  1286.   }
  1287.   insimg=gdImageCreateFromGif(inf); fclose(inf);
  1288.   if(insimg==NULL) {
  1289.     fly_error("bad_gif"); return;
  1290.   }
  1291.   if(pm->fill) {
  1292.     gdImageSetTile(image,insimg); tiled=1; tileimg=insimg;
  1293.   }
  1294.   else {
  1295.     gdImageSetBrush(image,insimg);brushed=1; brushimg=insimg;
  1296.   }
  1297. }
  1298.  
  1299. /* kill brush */
  1300. void obj_killbrush(objparm *pm)
  1301. {
  1302.   if(brushimg) gdImageDestroy(brushimg);
  1303.   brushed=0; brushimg=NULL;
  1304. }
  1305.  
  1306. /* kill tile */
  1307. void obj_killtile(objparm *pm)
  1308. {
  1309.   if(tileimg) gdImageDestroy(tileimg);
  1310.   tiled=0; tileimg=NULL;
  1311. }
  1312.  
  1313. /* set style */
  1314. void obj_setstyle(objparm *pm)
  1315. {
  1316.   int i,t;
  1317.   t=pm->pcnt/3; if(t<=0) {
  1318.     fly_error("too_few_parms"); return;
  1319.   }
  1320.   for(i=0;i<t;i++) {
  1321.     if(pm->pd[3*i]<0 || pm->pd[3*i+1]<0 || pm->pd[3*i+2]<0)
  1322.       pm->p[i]=gdTransparent;
  1323.     else
  1324.       pm->p[i]=getcolor(pm->pd[3*i],pm->pd[3*i+1],pm->pd[3*i+2]);
  1325.   }
  1326.   gdImageSetStyle(image,pm->p,t); styled=1;
  1327. }
  1328.  
  1329. /* kill style */
  1330. void obj_killstyle(objparm *pm)
  1331. {
  1332.   styled=0;
  1333. }
  1334.  
  1335. /* set transparent */
  1336. void obj_transp(objparm *pm)
  1337. {
  1338.   gdImageColorTransparent(image,pm->color[0]);
  1339. }
  1340.  
  1341. /* set interlace */
  1342. void obj_interlace(objparm *pm)
  1343. {
  1344.   gdImageInterlace(image,1);
  1345. }
  1346.  
  1347. /* set linewidth */
  1348. void obj_linewidth(objparm *pm)
  1349. {
  1350.   if(pm->pd[0]<1 || pm->pd[0]>255) fly_error("bad_parms");
  1351.   width=pm->pd[0];
  1352.   if(tikz_file) fprintf(tikz_file,"\\pgfsetlinewidth{%fpt}\n",width*0.7);
  1353. }
  1354.  
  1355. /* set crosshairsize */
  1356. void obj_crosshairsize(objparm *pm)
  1357. {
  1358.   if(pm->pd[0]<1 || pm->pd[0]>255) fly_error("bad_parms");
  1359.   else width2=pm->pd[0];
  1360. }
  1361.  
  1362. /* set x range */
  1363. void obj_xrange(objparm *pm)
  1364. {
  1365.   double dd;
  1366.   dd=pm->pd[1]-pm->pd[0];
  1367.   xstart=pm->pd[0]; xscale=sizex/dd;
  1368. }
  1369.  
  1370. /* set y range */
  1371. void obj_yrange(objparm *pm)
  1372. {
  1373.   double dd;
  1374.   dd=pm->pd[1]-pm->pd[0];
  1375.   ystart=pm->pd[1]; yscale=-sizey/dd;
  1376. }
  1377.  
  1378. /* set x et y range */
  1379. void obj_range(objparm *pm)
  1380. {
  1381.   double d1,d2;
  1382.   d1=pm->pd[1]-pm->pd[0]; d2=pm->pd[3]-pm->pd[2];
  1383.   xstart=pm->pd[0]; xscale=(sizex-1)/d1;
  1384.   ystart=pm->pd[3]; yscale=-(sizey-1)/d2;
  1385. }
  1386.  
  1387. /* set t range */
  1388. void obj_trange(objparm *pm)
  1389. {
  1390.   /*double dd;
  1391.   dd=pm->pd[1]-pm->pd[0];*/
  1392.   tstart=pm->pd[0]; tend=pm->pd[1];
  1393.   tranged=1;
  1394. }
  1395.  
  1396. /* set tstep (plotting step) */
  1397. void obj_tstep(objparm *pm)
  1398. {
  1399.   int dd;
  1400.   dd=pm->pd[0];
  1401.   if(dd<3) {
  1402.     fly_error("bad_step"); return;
  1403.   }
  1404.   if(dd>2000) dd=2000;
  1405.   tstep=dd;
  1406. }
  1407.  
  1408. /* set plotjump (plot jump break threashold) */
  1409. void obj_plotjump(objparm *pm)
  1410. {
  1411.   int dd;
  1412.   dd=pm->pd[0];
  1413.   if(dd<3) dd=3;
  1414.   if(dd>MAX_SIZE) dd=MAX_SIZE;
  1415.   plotjump=dd;
  1416. }
  1417.  
  1418. /* plot a curve, either parametric or explicit */
  1419. void _obj_plot(objparm *pm,int dash)
  1420. {
  1421.   int i,j,n,dist,xx,yy,varpos;
  1422.   char p1[MAX_LINELEN+1], p2[MAX_LINELEN+1];
  1423.   char *varn, *pp;
  1424.   double dc[2],t,v;
  1425.   int ic[2],oc[2];
  1426.  
  1427.   n=itemnum(pm->str);
  1428.   if(n<1) {
  1429.     fly_error("bad_parms"); return;
  1430.   }
  1431.   fnd_item(pm->str,1,p1); fnd_item(pm->str,2,p2);
  1432.   if(n==1 && !tranged) v=sizex/xscale/tstep;
  1433.   else v=(tend-tstart)/tstep;
  1434.   if(n==1) varn="x"; else varn="t";
  1435.   for(pp=varchr(p1,varn); pp; pp=varchr(pp,varn)) {
  1436.     string_modify(p1,pp,pp+strlen(varn),EV_T);
  1437.     pp+=strlen(EV_T);
  1438.   }
  1439.   for(pp=varchr(p2,varn); pp; pp=varchr(pp,varn)) {
  1440.     string_modify(p2,pp,pp+strlen(varn),EV_T);
  1441.       pp+=strlen(EV_T);
  1442.   }
  1443.   varpos=eval_getpos(EV_T);
  1444.   if(varpos<0) return; /* unknown error */
  1445.   evalue_compile(p1); evalue_compile(p2);
  1446.   if(vimg_enable) vimg_plotstart();
  1447.   if (tikz_file) fprintf(tikz_file,"\%\%begin plot\n\\draw\[%s]",
  1448.     tikz_options(pm->color[0],-dash));
  1449.   for(i=j=0;i<=tstep;i++) {
  1450.     if(n==1) {
  1451.       if(tranged) t=tstart+i*v; else t=xstart+i*v;
  1452.       eval_setval(varpos,t); dc[0]=t; dc[1]=strevalue(p1);
  1453.     }
  1454.     else {
  1455.       t=tstart+i*v; eval_setval(varpos,t);
  1456.       dc[0]=strevalue(p1); dc[1]=strevalue(p2);
  1457.     }
  1458.     if(!isfinite(dc[0]) || !isfinite(dc[1])) ic[0]=ic[1]=-BOUND;
  1459.     else scale(dc,ic,1);
  1460.     if(vimg_enable) vimg_plot1 (scale_buf[0],scale_buf[1]);
  1461.     if(j==0) {
  1462.       gdImageSetPixel(image,ic[0],ic[1],pm->color[0]); j++;
  1463.       if (tikz_file)
  1464.         fprintf(tikz_file,"(%i,%i)rectangle(%i,%i)",
  1465.           ic[0],flip(ic[1]),ic[0]+1,flip(ic[1]+1));
  1466.     }
  1467.     else {
  1468.       xx=ic[0]-oc[0]; yy=ic[1]-oc[1];
  1469.       dist=sqrt(xx*xx+yy*yy);
  1470.       if(dist>0){
  1471.         if(dist>plotjump || dist==1) {
  1472.           gdImageSetPixel(image,ic[0],ic[1],pm->color[0]);
  1473.           if (tikz_file)
  1474.             fprintf(tikz_file,"(%i,%i)rectangle(%i,%i)",
  1475.               ic[0],flip(ic[1]),ic[0]+1,flip(ic[1]+1));
  1476.         }
  1477.         else {
  1478.           if (dash==1)
  1479.             myDashedLine(image,oc[0],oc[1],ic[0],ic[1],pm->color[0]);
  1480.           else
  1481.             gdImageLine(image,oc[0],oc[1],ic[0],ic[1],pm->color[0]);
  1482.           if (tikz_file)
  1483.             fprintf(tikz_file, "(%i,%i)--(%i,%i)",
  1484.                 oc[0],flip(oc[1]),ic[0],flip(ic[1]));
  1485.         }
  1486.       }
  1487.     }
  1488.     memmove(oc,ic,sizeof(oc));
  1489.   }
  1490.   if (tikz_file) fprintf(tikz_file,";\n\%\%end plot\n");
  1491.   if(vimg_enable) vimg_plotend();
  1492. }
  1493.  
  1494. void obj_plot(objparm *pm) { _obj_plot( pm,0) ;}
  1495. void obj_dplot(objparm *pm) { _obj_plot( pm,1) ; }
  1496.  
  1497. /* set levelcurve granularity */
  1498. void obj_levelstep(objparm *pm)
  1499. {
  1500.   int dd;
  1501.   dd=pm->pd[0];
  1502.   if(dd<1) return;
  1503.   if(dd>16) dd=16;
  1504.   lstep=dd;
  1505. }
  1506.  
  1507. /* level curve */
  1508. void obj_levelcurve(objparm *pm)
  1509. {
  1510.   char fn[MAX_LINELEN+1],tc[MAX_LINELEN+1];
  1511.   int n,i;
  1512.   double d;
  1513.   leveldata *ld;
  1514.  
  1515.   ld=xmalloc(sizeof(leveldata)+16);
  1516.   ld->xname="x"; ld->yname="y";
  1517.   ld->xsize=sizex; ld->ysize=sizey; ld->datacnt=0;
  1518.   ld->xrange[0]=xstart; ld->xrange[1]=sizex/xscale+xstart;
  1519.   ld->yrange[0]=sizey/yscale+ystart; ld->yrange[1]=ystart;
  1520.   ld->grain=lstep;
  1521.   fnd_item(pm->str,1,fn); ld->fn=fn;
  1522.   n=itemnum(pm->str);
  1523.   if(n<=1) {ld->levelcnt=1; ld->levels[0]=0;}
  1524.   else {
  1525.     if(n>LEVEL_LIM+1) n=LEVEL_LIM+1;
  1526.     for(i=0;i<n-1;i++) {
  1527.       fnd_item(pm->str,i+2,tc);
  1528.       d=strevalue(tc);
  1529.       if(isfinite(d)) ld->levels[i]=d; else ld->levels[i]=0;
  1530.     }
  1531.     ld->levelcnt=n-1;
  1532.   }
  1533.   levelcurve(ld);
  1534.   if (tikz_file)
  1535.     fprintf(tikz_file,"\%\%begin levelcurve\n\\draw\[%s]"
  1536.       ,tikz_options(pm->color[0],1));
  1537.   for(i=0;i<ld->datacnt;i++) {
  1538.     gdImageSetPixel(image,ld->xdata[i],ld->ydata[i],pm->color[0]);
  1539.     if (tikz_file)
  1540.       fprintf(tikz_file,
  1541.         "(%i,%i)rectangle(%i,%i)",
  1542.         ld->xdata[i],flip(ld->ydata[i]),ld->xdata[i]+1,flip(ld->ydata[i]+1));
  1543.   }
  1544.   if (tikz_file) fprintf(tikz_file,";\n\%\%end levelcurve\n");
  1545.   free(ld);
  1546. }
  1547.  
  1548. /* set animation step */
  1549. void obj_animstep(objparm *pm)
  1550. {
  1551.   animstep=pm->pd[0];
  1552.   setvar("animstep",pm->pd[0]);
  1553. }
  1554.  
  1555. /* Starts a line counting */
  1556. void obj_linecount(objparm *pm)
  1557. {
  1558.   linecnt=pm->pd[0];
  1559. }
  1560.  
  1561. /* marks end of execution */
  1562. void obj_end(objparm *pm)
  1563. {
  1564.   fly_error("successful_end_of_execution");
  1565. }
  1566.  
  1567. /* output */
  1568. void obj_output(objparm *pm)
  1569. {
  1570.   char *p, namebuf[1024];
  1571.   p=find_word_start(pm->str); *find_word_end(p)=0;
  1572.   snprintf(namebuf,sizeof(namebuf),"%s",imagefilename);
  1573.   snprintf(imagefilename,sizeof(imagefilename),"%s",p);
  1574.   output();
  1575.   snprintf(imagefilename,sizeof(imagefilename),"%s",namebuf);
  1576. }
  1577.  
  1578. /* vimgfile */
  1579. void obj_vimgfile(objparm *pm)
  1580. {
  1581.   char *p;
  1582.   p=find_word_start(pm->str); *find_word_end(p)=0;
  1583.   snprintf(vimgfilename,sizeof(vimgfilename),"%s",p);
  1584.   if(vimg_ready) vimg_close();
  1585. }
  1586.  
  1587. /* vimg enable/disable */
  1588. void obj_vimg(objparm *pm)
  1589. {
  1590.   vimg_enable=pm->pd[0];
  1591.   if(vimg_enable>0 && vimg_ready==0) vimg_init();
  1592. }
  1593.  
  1594. /* Set affine transformation */
  1595. void obj_affine(objparm *pm)
  1596. {
  1597.   int i;
  1598.   for(i=0;i<4;i++) matrix[i]=pm->pd[i];
  1599.   transx=pm->pd[4]; transy=pm->pd[5];
  1600.   transform=1;
  1601. }
  1602.  
  1603. /* Set affine transformation */
  1604. void obj_rotation(objparm *pm)
  1605. {
  1606.   double r;
  1607.   r=M_PI*pm->pd[0]/180;
  1608.   matrix[0]=matrix[3]=cos(r);
  1609.   matrix[1]=-sin(r); matrix[2]=sin(r);
  1610.   transform=1;
  1611. }
  1612.  
  1613. /* Set linear transformation */
  1614. void obj_linear(objparm *pm)
  1615. {
  1616.   int i;
  1617.   for(i=0;i<4;i++) matrix[i]=pm->pd[i];
  1618.   transform=1;
  1619. }
  1620.  
  1621. /* Set translation transformation */
  1622. void obj_translation(objparm *pm)
  1623. {
  1624.   transx=pm->pd[0]; transy=pm->pd[1];
  1625. }
  1626.  
  1627. /* kill affine transformation */
  1628. void obj_killaffine(objparm *pm)
  1629. {
  1630.   matrix[0]=matrix[3]=1;
  1631.   matrix[1]=matrix[2]=transx=transy=0;
  1632.   transform=0;
  1633. }
  1634.  
  1635. /* kill affine transformation */
  1636. void obj_killlinear(objparm *pm)
  1637. {
  1638.   matrix[0]=matrix[3]=1;
  1639.   matrix[1]=matrix[2]=0;
  1640.   transform=0;
  1641. }
  1642.  
  1643. /* kill affine transformation */
  1644. void obj_killtranslation(objparm *pm)
  1645. {
  1646.   transx=transy=0;
  1647. }
  1648.  
  1649. /***** Les modifs de Jean-Christophe Leger Fev 2006 *****/
  1650.  
  1651. /* arguments: un numero de matrice entre 1 et JC_NB_MATRICES, une liste de 4 nombres reels pour la matrice */
  1652. void obj_setmatrix(objparm *pm)
  1653. {
  1654.   int nummatrix = (int) (pm->pd[0]);
  1655.   if((nummatrix < 1) ||(nummatrix>JC_NB_MATRICES)) {
  1656.     fly_error("bad_matrix_number");
  1657.     return;
  1658.   }
  1659.   matrices_pavage[nummatrix][0] = pm->pd[1];
  1660.   matrices_pavage[nummatrix][1] = pm->pd[2];
  1661.   matrices_pavage[nummatrix][2] = pm->pd[3];
  1662.   matrices_pavage[nummatrix][3] = pm->pd[4];
  1663. }
  1664.  
  1665. /* arguments: un numero de matrice entre 1 et JC_NB_MATRICES */
  1666. void obj_resetmatrix(objparm *pm)
  1667. {
  1668.   int nummatrix = (int) (pm->pd[0]);
  1669.   if((nummatrix < 1) ||(nummatrix>JC_NB_MATRICES)) {
  1670.     fly_error("bad_matrix_number");
  1671.     return;
  1672.   }
  1673.   matrices_pavage[nummatrix][0] = 1;
  1674.   matrices_pavage[nummatrix][1] = 0;
  1675.   matrices_pavage[nummatrix][2] = 0;
  1676.   matrices_pavage[nummatrix][3] = 1;
  1677.  
  1678. }
  1679. /* arguments: un numero de vecteur entre 1 et JC_NB_MATRICES, une liste de 2 nombres reels pour le vecteur */
  1680. void obj_setvector(objparm *pm)
  1681. {
  1682.   int nummatrix = (int) (pm->pd[0]);
  1683.   if((nummatrix < 1) ||(nummatrix>JC_NB_MATRICES)) {
  1684.     fly_error("bad_vector_number");
  1685.     return;
  1686.   }
  1687.   vecteurs_pavage[nummatrix][0] = pm->pd[1];
  1688.   vecteurs_pavage[nummatrix][1] = pm->pd[2];
  1689. }
  1690.  
  1691. /* arguments: un numero de matrice entre 1 et JC_NB_MATRICES */
  1692. void obj_resetvector(objparm *pm)
  1693. {
  1694.   int nummatrix = (int) (pm->pd[0]);
  1695.   if((nummatrix < 1) ||(nummatrix>JC_NB_MATRICES)) {
  1696.     fly_error("bad_vector_number");
  1697.     return;
  1698.   }
  1699.   vecteurs_pavage[nummatrix][0] = 0;
  1700.   vecteurs_pavage[nummatrix][1] = 0;
  1701. }
  1702.  
  1703. /* arguments: un numero de vecteur entre 1 et JC_NB_MATRICES, une liste de 6 nombres reels pour la matrice et le vecteur */
  1704. void obj_settransform(objparm *pm)
  1705. {
  1706.   int nummatrix = (int) (pm->pd[0]);
  1707.   if((nummatrix < 1) ||(nummatrix>JC_NB_MATRICES)) {
  1708.     fly_error("bad_vector_number");
  1709.     return;
  1710.   }
  1711.   matrices_pavage[nummatrix][0] = pm->pd[1];
  1712.   matrices_pavage[nummatrix][1] = pm->pd[2];
  1713.   matrices_pavage[nummatrix][2] = pm->pd[3];
  1714.   matrices_pavage[nummatrix][3] = pm->pd[4];
  1715.   vecteurs_pavage[nummatrix][0] = pm->pd[5];
  1716.   vecteurs_pavage[nummatrix][1] = pm->pd[6];
  1717. }
  1718.  
  1719.  
  1720. /* arguments: un numero de matrice entre 1 et JC_NB_MATRICES */
  1721. void obj_resettransform(objparm *pm)
  1722. {
  1723.   int nummatrix = (int) (pm->pd[0]);
  1724.   if((nummatrix < 1) ||(nummatrix>JC_NB_MATRICES)) {
  1725.     fly_error("bad_vector_number");
  1726.     return;
  1727.   }
  1728.   vecteurs_pavage[nummatrix][0] = 0;
  1729.   vecteurs_pavage[nummatrix][1] = 0;
  1730.   matrices_pavage[nummatrix][0] = 1;
  1731.   matrices_pavage[nummatrix][1] = 0;
  1732.   matrices_pavage[nummatrix][2] = 0;
  1733.   matrices_pavage[nummatrix][3] = 1;
  1734. }
  1735.  
  1736. /* arguments: une liste de 6 nombres reels pour les coordonnees de l'origine et des vecteurs de base */
  1737. void obj_setparallelogram(objparm *pm)
  1738. {
  1739.   int i;
  1740.   for(i=0;i<6;i++)  parallogram_fonda[i]=pm->pd[i];
  1741. }
  1742.  
  1743. /* arguments :aucun */
  1744. void obj_resetparallelogram(objparm *pm)
  1745. {
  1746.   parallogram_fonda[0]=0;
  1747.   parallogram_fonda[1]=0;
  1748.   parallogram_fonda[2]=100;
  1749.   parallogram_fonda[3]=0;
  1750.   parallogram_fonda[4]=0;
  1751.   parallogram_fonda[5]=100;
  1752. }
  1753.  
  1754. /* argument : la liste des numeros des matrices a faire agir, le nom du fichier de l'image a inclure */
  1755. void obj_multicopy(objparm *pm)
  1756. {
  1757.   char *pp,*pe,*cptr;
  1758.   FILE *inf;
  1759.   gdImagePtr insimg;
  1760.   int i,j;
  1761.   int imax,jmax;
  1762.   int c; /* la couleur reperee */
  1763.   int mat;
  1764.   int nummatrices[JC_NB_MATRICES];
  1765.   double pt[2]; /* les coordonnees math. du point repere dans le parallelogramme */
  1766.   double ptr[2]; /* les coordonnees math. du point transforme */
  1767.   int pti[2]; /* les coordonnees img du point a placer sur le canevas */
  1768.   int t;
  1769.   /* gestion palette de couleur de l'image a inserer */
  1770.   int colorMap[gdMaxColors];
  1771.   int comptmat=0; /* combien de matrices en tout ? */
  1772.  
  1773.   for (i=0; (i<gdMaxColors); i++) {
  1774.     colorMap[i] = (-1);
  1775.   }
  1776.  
  1777.   /* Gestion des parametres la liste est censee finir avec le nom du fichier */
  1778.   t=pm->pcnt-1; /* il faut enlever le nom de fichier ! c'est le nombre de parametres en plus */
  1779.   pp=find_word_start(pm->str);
  1780.   /* visiblement find_word_start n'arrive pas a trouver le dernier mot dans un contexte ou le nombre de parameters est variable
  1781.      * on cherche donc l'emplacement de la derniere virgule:
  1782.      * il y en a une car t>0, on retrouve ensuite le debut de mot
  1783.      */
  1784.   for(pe=pp;*pe!=0;pe++);/* trouve la fin de la chaine */
  1785.   if(t>0){
  1786.     for(cptr = pe; cptr > pp && *cptr != ','; cptr--);
  1787.     pp=find_word_start(cptr+1);
  1788.   }
  1789.   pe=find_word_end(pp);*pe=0; /* strip junk final */
  1790.   //      printf("pp contient %s\n",pp);
  1791.  
  1792.  
  1793.   inf=open4read(pp);
  1794.   if(inf==NULL) {
  1795.     fly_error("file_not_exist"); return;
  1796.   }
  1797.   insimg=gdImageCreateFromGif(inf); fclose(inf);
  1798.   if(insimg==NULL) {
  1799.     fly_error("bad_gif"); return;
  1800.   }
  1801.  
  1802.   /* On recupere les numeros des matrices/vecteurs a faire agir,
  1803.     * s'il n'y en a pas, on les fait toutes agir
  1804.   */
  1805.   for(i=0;i<t && i< JC_NB_MATRICES;i++) {
  1806.     if(pm->pd[i]>=1 && pm->pd[i]<=JC_NB_MATRICES){
  1807.       nummatrices[comptmat] = pm->pd[i];
  1808.       comptmat++;
  1809.     }
  1810.   }
  1811.   if(t<=0){
  1812.     for(i=0;i<JC_NB_MATRICES;i++) {
  1813.       nummatrices[i] = i+1;
  1814.     }
  1815.     comptmat=JC_NB_MATRICES;
  1816.   }
  1817.  
  1818.  
  1819.   imax = gdImageSX(insimg);
  1820.   jmax = gdImageSY(insimg);
  1821.  
  1822.   for(i=0;i<imax;i++){
  1823.     for(j=0;j<jmax;j++){
  1824.       int nc;
  1825.       c=gdImageGetPixel(insimg,i,j);
  1826.       /* Le code suivant est une copie du code dans gdImageCopy  Added 7/24/95: support transparent copies */
  1827.       if (gdImageGetTransparent(insimg) != c) {
  1828.         /* Have we established a mapping for this color? */
  1829.         if (colorMap[c] == (-1)) {
  1830.           /* If it's the same image, mapping is trivial, dans notre cas ca n'arrive jamais... */
  1831.           if (image == insimg) {
  1832.             nc = c;
  1833.           } else {
  1834.             /* First look for an exact match */
  1835.             nc = gdImageColorExact(image,
  1836.                    insimg->red[c], insimg->green[c],
  1837.                            insimg->blue[c]);
  1838.           }
  1839.           if (nc == (-1)) {
  1840.           /* No, so try to allocate it */
  1841.             nc = gdImageColorAllocate(image,
  1842.                              insimg->red[c], insimg->green[c],
  1843.                              insimg->blue[c]);
  1844.              /* If we're out of colors, go for the closest color */
  1845.             if (nc == (-1)) {
  1846.               nc = gdImageColorClosest(image,
  1847.                               insimg->red[c], insimg->green[c],
  1848.                               insimg->blue[c]);
  1849.             }
  1850.           }
  1851.           colorMap[c] = nc;
  1852.         }
  1853.  
  1854.         pt[0]= i*(parallogram_fonda[2])/(imax-1)+(jmax-j)*(parallogram_fonda[4])/(jmax-1)+parallogram_fonda[0];
  1855.         pt[1]= i*(parallogram_fonda[3])/(imax-1)+(jmax-j)*(parallogram_fonda[5])/(jmax-1)+parallogram_fonda[1];
  1856.         for(mat=0;mat<comptmat;mat++){
  1857.           double *matricecourante = matrices_pavage[nummatrices[mat]];
  1858.           double *vecteurcourant = vecteurs_pavage[nummatrices[mat]];
  1859.           ptr[0] = matricecourante[0]*pt[0]+matricecourante[1]*pt[1]+vecteurcourant[0];
  1860.           ptr[1] = matricecourante[2]*pt[0]+matricecourante[3]*pt[1]+vecteurcourant[1];
  1861.           scale(ptr,pti,1);
  1862.           gdImageSetPixel(image,pti[0],pti[1],colorMap[c]);
  1863.         }
  1864.       }
  1865.     }
  1866.   }
  1867.   gdImageDestroy(insimg);
  1868. }
  1869.  
  1870. /**** Fin modifs JC Fev 06 ******/
  1871.  
  1872. /* get color from value */
  1873. int calc_color(char *p, struct objtab *o)
  1874. {
  1875.   int k, cc[3];
  1876.   char buf[MAX_LINELEN+1];
  1877.   for(k=0;k<3;k++) {
  1878.     fnd_item(p,k+1,buf);
  1879.     cc[k]=strevalue(buf);
  1880.   }
  1881.   collapse_item(p,3);
  1882.   if(cc[0]==-1 && cc[1]==-1 && cc[2]==-255) {
  1883.  
  1884.     if(brushed && o->subst&1) return gdBrushed;
  1885.     else return color_black;
  1886.   }
  1887.   if(cc[0]==-1 && cc[1]==-255 && cc[2]==-1) {
  1888.     if(styled && o->subst&2)  return gdStyled;
  1889.     else return color_black;
  1890.   }
  1891.   if(cc[0]==-255 && cc[1]==-1 && cc[2]==-1) {
  1892.     if(tiled && o->fill_tag==1)  return gdTiled;
  1893.     else return color_black;
  1894.   }
  1895.   return getcolor(cc[0],cc[1],cc[2]);
  1896. }
  1897.  
  1898. /* Routine to parse parameters */
  1899. int parse_parms(char *p,objparm *pm,struct objtab *o)
  1900. {
  1901.   char buf[MAX_LINELEN+1];
  1902.   char *p1, *p2;
  1903.   int j,t,c,c1,c2;
  1904.  
  1905.   c=o->color_pos;c1=c2=0;
  1906.   pm->color[0]=pm->color[1]=0;
  1907.   if(c>0) c1=c;
  1908.   if(c<0) c2=-c;
  1909.   c=c1+c2;
  1910.   t=itemnum(p);if(t<o->required_parms+3*c) return -1;
  1911.   if(c1>0 && t>o->required_parms+3*c) t=o->required_parms+3*c;
  1912.   pm->pcnt=t-3*c;
  1913.   if(pm->pcnt>MAX_PARMS) pm->pcnt=MAX_PARMS;
  1914.   if(c2>0) {
  1915.     for(j=0;j<2 && j<c2; j++) pm->color[j]=calc_color(p,o);
  1916.   }
  1917.   snprintf(buf,sizeof(buf),"%s",p);
  1918.   for(j=0, p1=buf; j<pm->pcnt; j++, p1=p2) {
  1919.     p2=find_item_end(p1); if(*p2) *p2++=0;
  1920.     p1=find_word_start(p1);
  1921.     if(*p1) pm->pd[j]=strevalue(p1); else pm->pd[j]=0;
  1922.     if(!isfinite(pm->pd[j])) {
  1923.         if(j<o->required_parms) return -1;
  1924.         else {pm->pd[j]=0;break;}
  1925.     }
  1926.   }
  1927.   collapse_item(p,o->required_parms);
  1928.   if(c1>0) {
  1929.     for(j=0;j<c1 && j<2; j++) pm->color[j]=calc_color(p,o);
  1930.   }
  1931.   if(width>1 && o->subst&1 && pm->color[0]!=gdBrushed
  1932.       && pm->color[0]!=gdStyled && pm->color[0]!=gdTiled) {
  1933.     pm->color[1]=pm->color[0];
  1934.     pm->color[0]=widthcolor(width,pm->color[0]);
  1935.   }
  1936.   pm->fill=o->fill_tag;
  1937.   if(o->required_parms!=0){
  1938.     if(dashed) {
  1939.       if(pm->fill==1) pm->fill=2;
  1940.       if(pm->fill==0) pm->fill=-1;
  1941.        if(noreset==0) dashed=0;
  1942.     }
  1943.     if(filled) {
  1944.       if(pm->fill==-1) pm->fill=2;
  1945.       if(pm->fill==0) pm->fill=1;
  1946.       if(noreset==0) filled=0;
  1947.     }
  1948.   }
  1949.   ovlstrcpy(pm->str,p); return 0;
  1950. }
  1951.  
  1952. /* Create the struct objparm pm corresponding to the objparm p
  1953.  * Execute a command. Returns 0 if OK.
  1954.  */
  1955. int obj_main(char *p)
  1956. {
  1957.   int i;
  1958.   char *pp,*name,*pt;
  1959.   char tbuf2[MAX_LINELEN+1];
  1960.   struct objparm pm;
  1961.  
  1962.   p=find_word_start(p);
  1963.   if(*p==exec_prefix_char || *p==0) return 0; /* comment */
  1964.   name=p;
  1965.   pp=find_name_end(p);
  1966.   pt=find_word_start(pp); if(*pt=='=') *pt=' ';
  1967.   if(*pt==':' && *(pt+1)=='=') *pt=*(pt+1)=' ';
  1968.   pp=find_word_end(p);
  1969.   if(*pp!=0) {
  1970.     *(pp++)=0; pp=find_word_start(pp);
  1971.   }
  1972.   if(strlen(p)==1 || (strlen(p)==2 && isdigit(*(p+1)))) {
  1973.     setvar(p,strevalue(pp)); return 0;
  1974.   }
  1975.   /* search the number of the object */
  1976.   i=search_list(objtab,obj_no,sizeof(objtab[0]),name);
  1977.   if(i<0) {
  1978.     fly_error("bad_cmd"); return 1;
  1979.   }
  1980.   if(image==NULL && (objtab[i].color_pos || objtab[i].required_parms>2)) {
  1981.     fly_error("image_not_defined"); return 1;
  1982.   }
  1983.   ovlstrcpy(tbuf2,pp);
  1984.   if(objtab[i].color_pos || objtab[i].routine==obj_setstyle) {
  1985.     substit(tbuf2);
  1986.   }
  1987.   if(parse_parms(tbuf2,&pm,objtab+i)!=0) fly_error("bad_parms");
  1988.   else objtab[i].routine(&pm);
  1989.   return 0;
  1990. }
  1991.