WebSVN – wimsdev – Diff – /trunk/wims/src/Flydraw/objects.c

  *  Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
  */
 #include <errno.h>
 #include "flydraw.h"
-void myGdImageArc(gdImagePtr, double, double, double, double, double, double, int);
 const double DEG=M_PI/180;
 /* Tikz things */
 static char *tikz_color (int the_color)
       fprintf(tikz_file, "(%i,%i)--(%i,%i)",
        pm->p[i],flip(pm->p[i+1]),pm->p[i+2],flip(pm->p[i+3]));
+  }
   if(tikz_file) fprintf(tikz_file, ";\n");
+}
-/* for geodesics, return data for obj_hypgeods + a point on the arc
-  to be used for filling hyperbolic triangle*/
-int hypgeodaux(double *q, double* res)
-{
-  double r,cx,cy,a1,a2,a3,tmp,
-    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],
-    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],
-    dy = -2*q[1]*q[2]+2*q[0]*q[3];
-  if (dy*dy*1e4 < nx*nx+ny*ny){
-    res[5]=(q[0]+q[2])/2;
-    res[6]=(q[1]+q[3])/2;
-    return 0;}
-  cx = ny/dy; cy=-nx/dy;
-  a1 = atan2(q[1]-cy, q[0]-cx);
-  a2 = atan2(q[3]-cy, q[2]-cx);
-  if (fabs(a2-a1)>M_PI){if(a1<a2) a1+=2*M_PI; else a2+=2*M_PI;};
-  a3 = (a1+a2)/2;
-  if(a1>a2) {tmp=a1; a1=a2; a2=tmp;}
-  r = sqrt(cx*cx+cy*cy-1);
-  res[0]=cx;
-  res[1]=cy;
-  res[2]=r;
-  res[3]=a1;
-  res[4]=a2;
-  res[5]=cx+r*cos(a3);
-  res[6]=cy+r*sin(a3);
-  return 1;
-}
-/* hyperbolic geodesics from x1,y1 to x2,y2, x3,y3 to x4,y4 in Poincaré disk */
-void obj_hypgeods(objparm *pm)
-{
-  int i, *qq=pm->p;
-  double rx,ry,res[7];
-  for (i = 0; i < pm->pcnt; i+=4)
-    if (hypgeodaux(pm->pd+i,res)){
-      scale2(res[2],res[2],&rx,&ry);
-      scale(res,qq,1);
-      myGdImageArc(image,qq[0],qq[1],rx,ry,res[3]/DEG,res[4]/DEG,pm->color[0]);
-    }
-    else {
-      scale(pm->pd+i,qq,2);
-      gdImageLine(image,qq[0],qq[1],qq[2],qq[3],pm->color[0]);
-    }
-}
-/* hyperbolique triangle, can be filled */
-void obj_hyptriangle(objparm *pm)
-{
-  double rx,ry,data[8],res[7];
-  int i, *qq=pm->p;
-  for(i=0; i<8; i++) data[i]=pm->pd[i%6];
-  for(i=0; i<3; i++)
-    if(hypgeodaux(data+2*i,res)) {
-      scale2(res[2],res[2],&rx,&ry);
-      scale(res,qq,1);
-      myGdImageArc(image,qq[0],qq[1],rx,ry,res[3]/DEG,res[4]/DEG,pm->color[0]);
-    }
-    else {
-      scale(data+2*i,qq,2);
-      gdImageLine(image,qq[0],qq[1],qq[2],qq[3],pm->color[0]);
-    }
-  if(pm->fill){
-    data[0]=res[5];
-    data[1]=res[6];
-    hypgeodaux(data,res);
-    scale(res+5,qq,1);
-    patchgdImageFill(image,qq[0],qq[1],pm->color[0]);
-  }
-}
-/* complete hyperbolic geodesic through two points */
-void obj_hyplines(objparm *pm)
-{
-  double rx,ry,res[7],*pd = pm->pd;
-  int i;
-  for (i=0; i < pm->pcnt; i+=4,pd+=4)
-    if (hypgeodaux(pd,res)) {
-      double alpha=atan(1/res[2]), beta = atan2(res[1],res[0]);
-      scale(res,pm->p,1);
-      scale2(res[2],res[2],&rx,&ry);
-      myGdImageArc(image,pm->p[0],pm->p[1],rx,ry,180+(beta-alpha)/DEG,180+(beta+alpha)/DEG,pm->color[0]);
-    }
-    else {
-      double gamma;
-      if (pd[1]*pd[1]+pd[0]*pd[0] > pd[2]*pd[2]+pd[3]*pd[3])
-        gamma = atan2(pd[1],pd[0]);
-      else
-        gamma = atan2(pd[3],pd[2]);
-      res[0]=cos(gamma);
-      res[1]=sin(gamma);
-      res[2]=-cos(gamma);
-      res[3]=-sin(gamma);
-      scale(res,pm->p,2);
-      gdImageLine(image,pm->p[0],pm->p[1],pm->p[2],pm->p[3],pm->color[0]);
-    }
-}
 /* segments */
 void obj_dlines(objparm *pm)
+{
  int i, n;
  n=(pm->pcnt)/2;
   if(tikz_file){
     pm->pd[5]-=pm->pd[4];
     pm->pd[5]-=360*floor(pm->pd[5]/360);
     pm->pd[5]+=pm->pd[4];
     fprintf(tikz_file,
-/*    "\\draw\[%s, domain=%f:%f] plot ({%i+%f*cos(\\x)}, {%i+%f*sin(\\x)});\n",*/
       "\\draw \[%s] (%i,%i) arc (%f:%f:%i and %i);\n",
         tikz_options(pm->color[0],pm->fill),
         (int)rint(pm->p[0]+rx*cos(DEG*pm->pd[4])),
         flip((int)rint(pm->p[1]+ry*sin(DEG*pm->pd[4]))),
         pm->pd[4],pm->pd[5],rx,-ry);
   {"giant", &gdFontGiant,"1"},
   {"huge",  &gdFontGiant,"1"}
 };
 #define fonttab_no (sizeof(fonttab)/sizeof(fonttab[0]))
 /* string */
 void obj_string(objparm *pm)
+{
   char *pp, *pe, *p2;
   int i;
   if(pm->fill) {
     gdImageSetTile(image,insimg); tiled=1; tileimg=insimg;
+  }
   else {
     gdImageSetBrush(image,insimg);brushed=1; brushimg=insimg;
   }
 }
 /* kill brush */
 void obj_killbrush(objparm *pm)
+{
   if(brushimg) gdImageDestroy(brushimg);
   styled=0;
+}
 /* set transparent */
 void obj_transp(objparm *pm)
 {
   gdImageColorTransparent(image,pm->color[0]);
+}
 /* set interlace */
 void obj_interlace(objparm *pm)
   if (tikz_file) fprintf(tikz_file,";\n\%\%end plot\n");
   if(vimg_enable) vimg_plotend();
+}
 void obj_plot(objparm *pm) { _obj_plot( pm,0) ;}
-void obj_dplot(objparm *pm) { _obj_plot( pm,1) ; }
+void obj_dplot(objparm *pm) { _obj_plot( pm,1) ;}
 /* set levelcurve granularity */
 void obj_levelstep(objparm *pm)
+{
   int dd;
 /***** Les modifs de Jean-Christophe Leger Fev 2006 *****/
 /* arguments: un numero de matrice entre 1 et JC_NB_MATRICES, une liste de 4 nombres reels pour la matrice */
 void obj_setmatrix(objparm *pm)
 {
   int nummatrix = (int) (pm->pd[0]);
   if((nummatrix < 1) ||(nummatrix>JC_NB_MATRICES)) {
     fly_error("bad_matrix_number");
     return;
   }
   matrices_pavage[nummatrix][0] = pm->pd[1];
   matrices_pavage[nummatrix][1] = pm->pd[2];
   matrices_pavage[nummatrix][2] = pm->pd[3];
   matrices_pavage[nummatrix][3] = pm->pd[4];
 }
 /* arguments: un numero de matrice entre 1 et JC_NB_MATRICES */
 void obj_resetmatrix(objparm *pm)
+{
   int nummatrix = (int) (pm->pd[0]);
   if((nummatrix < 1) ||(nummatrix>JC_NB_MATRICES)) {
     substit(tbuf2);
+  }
   if(parse_parms(tbuf2,&pm,objtab+i)!=0) fly_error("bad_parms");
   else objtab[i].routine(&pm);
   return 0;
+}
+/* for hyperbolic geodesics, return data for obj_hypgeods + a point on the arc
+  to be used for filling hyperbolic triangle*/
+int hypgeodaux(double *q, double* res)
+{
+  double r,cx,cy,a1,a2,a3,tmp,
+    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],
+    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],
+    dy = -2*q[1]*q[2]+2*q[0]*q[3];
+  if (dy*dy*1e4 < nx*nx+ny*ny){
+    res[5]=(q[0]+q[2])/2;
+    res[6]=(q[1]+q[3])/2;
+    return 0;}
+  cx = ny/dy; cy=-nx/dy;
+  a1 = atan2(q[1]-cy, q[0]-cx);
+  a2 = atan2(q[3]-cy, q[2]-cx);
+  if (fabs(a2-a1)>M_PI){if(a1<a2) a1+=2*M_PI; else a2+=2*M_PI;};
+  a3 = (a1+a2)/2;
+  if(a1>a2) {tmp=a1; a1=a2; a2=tmp;}
+  r = sqrt(cx*cx+cy*cy-1);
+  res[0]=cx;
+  res[1]=cy;
+  res[2]=r;
+  res[3]=a1;
+  res[4]=a2;
+  res[5]=cx+r*cos(a3);
+  res[6]=cy+r*sin(a3);
+  return 1;
+}
+/* hyperbolic geodesics from x1,y1 to x2,y2, x3,y3 to x4,y4 in Poincaré disk */
+void obj_hypgeods(objparm *pm)
+{
+  int i, *qq=pm->p;
+  double rx,ry,res[7];
+  if (tikz_file) fprintf(tikz_file, "\\draw[%s]",tikz_options(pm->color[0],pm->fill));
+  for (i = 0; i < pm->pcnt; i+=4)
+    if (hypgeodaux(pm->pd+i,res)){
+      scale2(res[2],res[2],&rx,&ry);
+      scale(res,qq,1);
+      myGdImageArc(image,qq[0],qq[1],rx,ry,res[3]/DEG,res[4]/DEG,pm->color[0]);
+      if(tikz_file) fprintf(tikz_file,"(%i,%i)arc(%f:%f:%i and %i)",
+        qq[0],flip(qq[1]),res[3]/DEG,res[4]/DEG,(int)rint(rx),-(int)rint(ry));
+    }
+    else {
+      scale(pm->pd+i,qq,2);
+      gdImageLine(image,qq[0],qq[1],qq[2],qq[3],pm->color[0]);
+      if(tikz_file) fprintf(tikz_file,"(%i,%i)--(%i,%i)",
+          qq[0],flip(qq[1]),qq[2],flip(qq[3]));
+    }
+  if(tikz_file) fprintf(tikz_file,";\n");
+}
+/* hyperbolic triangle, can be filled */
+void obj_hyptriangle(objparm *pm)
+{
+  double rx,ry,data[8],res[7];
+  int i, *qq=pm->p;
+  for(i=0; i<8; i++) data[i]=pm->pd[i%6];
+  if (tikz_file) fprintf(tikz_file, "\\draw[%s]",tikz_options(pm->color[0],pm->fill));
+  for(i=0; i<3; i++)
+    if(hypgeodaux(data+2*i,res)){
+      scale2(res[2],res[2],&rx,&ry);
+      scale(res,qq,1);
+      myGdImageArc(image,qq[0],qq[1],rx,ry,res[3]/DEG,res[4]/DEG,pm->color[0]);
+      if(tikz_file)
+        fprintf(tikz_file, "(%i,%i) arc (%f:%f:%i and %i)",
+          (int)rint(qq[0]+rx*cos(res[3])),
+          flip((int)rint(qq[1]+ry*sin(res[3]))),
+          res[3]/DEG,res[4]/DEG,(int)rint(rx),-(int)rint(ry));
+    }
+    else {
+      scale(data+2*i,qq,2);
+      gdImageLine(image,qq[0],qq[1],qq[2],qq[3],pm->color[0]);
+      if(tikz_file) fprintf(tikz_file,"(%i,%i)--(%i,%i)",
+            qq[0],flip(qq[1]),qq[2],flip(qq[3]));
+    }
+    if(tikz_file) fprintf(tikz_file,";\n");
+    if(pm->fill){
+      data[0]=res[5];
+      data[1]=res[6];
+      hypgeodaux(data,res);
+      scale(res+5,qq,1);
+      //patchgdImageFill(image,qq[0],qq[1],pm->color[0]);
+      patchgdImageFillToBorder(image,qq[0],qq[1],pm->color[0],pm->color[0]);
+      if(tikz_file) tikz_fill(qq[0],qq[1],1,1,0,pm->color[0],grid);
+    }
+}
+/* complete hyperbolic geodesic through the two first points, then the following two points, etc */
+void obj_hyplines(objparm *pm)
+{
+  double rx,ry,res[7],*pd = pm->pd;
+  int i;
+  if (tikz_file) fprintf(tikz_file, "\\draw[%s]",tikz_options(pm->color[0],0));
+  for (i=0; i < pm->pcnt; i+=4,pd+=4)
+    if (hypgeodaux(pd,res)) {
+      double alpha=atan(1/res[2]), beta = atan2(res[1],res[0]);
+      scale(res,pm->p,1);
+      scale2(res[2],res[2],&rx,&ry);
+      myGdImageArc(image,pm->p[0],pm->p[1],rx,ry,180+(beta-alpha)/DEG,180+(beta+alpha)/DEG,pm->color[0]);
+      if (tikz_file) fprintf(tikz_file, "(%i,%i) arc (%f:%f:%i and %i);\n",
+        (int)rint(pm->p[0]-rx*cos(beta-alpha)), flip((int)rint(pm->p[1]-ry*sin(beta-alpha))),
+          (beta-alpha)/DEG,(beta+alpha)/DEG, (int)rint(rx),-(int)rint(ry));
+    }
+    else {
+      double gamma;
+      if (pd[1]*pd[1]+pd[0]*pd[0] > pd[2]*pd[2]+pd[3]*pd[3])
+        gamma = atan2(pd[1],pd[0]);
+      else
+        gamma = atan2(pd[3],pd[2]);
+        res[0]=cos(gamma);
+        res[1]=sin(gamma);
+        res[2]=-cos(gamma);
+        res[3]=-sin(gamma);
+        scale(res,pm->p,2);
+        gdImageLine(image,pm->p[0],pm->p[1],pm->p[2],pm->p[3],pm->color[0]);
+        if(tikz_file) fprintf(tikz_file, "(%i,%i)--(%i,%i)",
+          pm->p[0],flip(pm->p[1]),pm->p[2],flip(pm->p[3]));
+    }
+    if (tikz_file) fprintf(tikz_file, ";\n");
+}

Subversion Repositories wimsdev

(root)/trunk/wims/src/Flydraw/objects.c – Rev 18351 → 18352