[BACK]Return to plotf.c CVS log [TXT][DIR] Up to [local] / OpenXM_contrib2 / asir2000 / plot

Diff for /OpenXM_contrib2/asir2000/plot/plotf.c between version 1.28 and 1.29

version 1.28, 2014/05/16 08:49:53 version 1.29, 2014/06/27 07:58:29
Line 45 
Line 45 
  * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,   * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
  * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.   * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
  *   *
  * $OpenXM: OpenXM_contrib2/asir2000/plot/plotf.c,v 1.27 2014/05/13 20:02:40 ohara Exp $   * $OpenXM: OpenXM_contrib2/asir2000/plot/plotf.c,v 1.28 2014/05/16 08:49:53 ohara Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "parse.h"  #include "parse.h"
Line 78  void PineqnandD(NODE,Obj *), PineqnandQ(NODE,Obj *), P
Line 78  void PineqnandD(NODE,Obj *), PineqnandQ(NODE,Obj *), P
 void PineqnorD(NODE,Obj *), PineqnorQ(NODE,Obj *), PineqnorB(NODE,Obj *);  void PineqnorD(NODE,Obj *), PineqnorQ(NODE,Obj *), PineqnorB(NODE,Obj *);
 void PineqnxorD(NODE,Obj *), PineqnxorQ(NODE,Obj *), PineqnxorB(NODE,Obj *);  void PineqnxorD(NODE,Obj *), PineqnxorQ(NODE,Obj *), PineqnxorB(NODE,Obj *);
 void PconplotD(NODE,Obj *),PconplotQ(NODE,Obj *),PconplotB(NODE,Obj *);  void PconplotD(NODE,Obj *),PconplotQ(NODE,Obj *),PconplotB(NODE,Obj *);
   void PpolarplotD(NODE,Obj *);
 void PplotoverD(NODE,Obj *),PplotoverQ(NODE,Obj *),PplotoverB(NODE,Obj *);  void PplotoverD(NODE,Obj *),PplotoverQ(NODE,Obj *),PplotoverB(NODE,Obj *);
 void ifplot_mainNG(NODE,char *,Obj *);  void ifplot_mainNG(NODE,char *,Obj *);
 void conplot_mainNG(NODE,char *,Obj *);  void conplot_mainNG(NODE,char *,Obj *);
Line 88  struct ftab plot_tab[]={
Line 89  struct ftab plot_tab[]={
 #if defined(INTERVAL)  #if defined(INTERVAL)
         {ITVIFPLOT,Pitvifplot,-8},          {ITVIFPLOT,Pitvifplot,-8},
 #endif  #endif
         {POLARPLOT,Ppolarplot,-6},  
         {OBJ_CP,Pobj_cp,4},          {OBJ_CP,Pobj_cp,4},
           {POLARPLOT,Ppolarplot,-6},
           {POLARPLOTD,PpolarplotD,-7},
         {IFPLOT,Pifplot,-7},          {IFPLOT,Pifplot,-7},
         {IFPLOTD,PifplotD,-8},          {IFPLOTD,PifplotD,-8},
         {IFPLOTQ,PifplotQ,-8},          {IFPLOTQ,PifplotQ,-8},
Line 118  struct ftab plot_tab[]={
Line 120  struct ftab plot_tab[]={
         {PLOTOVERD,PplotoverD,-4},          {PLOTOVERD,PplotoverD,-4},
         {PLOTOVERQ,PplotoverQ,-4},          {PLOTOVERQ,PplotoverQ,-4},
         {PLOTOVERB,PplotoverB,-4},          {PLOTOVERB,PplotoverB,-4},
         {MEMORY_IFPLOT,Pmemory_ifplot,-6},          {MEMORY_IFPLOT,Pmemory_ifplot,-7},
         {MEMORY_CONPLOT,Pmemory_conplot,-7},          {MEMORY_CONPLOT,Pmemory_conplot,-7},
         {MEMORY_PLOT,Pmemory_plot,-5},          {MEMORY_PLOT,Pmemory_plot,-7},
         {DRAWCIRCLE,Pdrawcircle,6},          {DRAWCIRCLE,Pdrawcircle,6},
         {OPEN_CANVAS,Popen_canvas,-3},          {OPEN_CANVAS,Popen_canvas,-3},
         {CLEAR_CANVAS,Pclear_canvas,2},          {CLEAR_CANVAS,Pclear_canvas,2},
Line 137  char *pfn[]={
Line 139  char *pfn[]={
         CONPLOTD,CONPLOTQ,CONPLOTB,ITVIFPLOT,          CONPLOTD,CONPLOTQ,CONPLOTB,ITVIFPLOT,
         PLOTOVERD,PLOTOVERQ,PLOTOVERB,          PLOTOVERD,PLOTOVERQ,PLOTOVERB,
         MEMORY_IFPLOT,MEMORY_CONPLOT,MEMORY_PLOT,ARRAYPLOT,OPEN_CANVAS,          MEMORY_IFPLOT,MEMORY_CONPLOT,MEMORY_PLOT,ARRAYPLOT,OPEN_CANVAS,
         DRAWCIRCLE,DRAW_OBJ,DRAW_STRING,OBJ_CP,CLEAR_CANVAS};          DRAWCIRCLE,DRAW_OBJ,DRAW_STRING,OBJ_CP,CLEAR_CANVAS,POLARPLOTD};
 /*  /*
         IFPLOT:0,CONPLOT:1,PLOT:2,INTERACTIVE:3,POLARPLOT:4,PLOTOVER:5,          IFPLOT:0,CONPLOT:1,PLOT:2,INTERACTIVE:3,POLARPLOT:4,PLOTOVER:5,
         IFPLOTD:6,IFPLOTQ:7,IFPLOTB:8,INEQND:9,INEQNQ:10,INEQNB:11,          IFPLOTD:6,IFPLOTQ:7,IFPLOTB:8,INEQND:9,INEQNQ:10,INEQNB:11,
Line 148  char *pfn[]={
Line 150  char *pfn[]={
         PLOTOVERD:25,PLOTOVERQ:26,PLOTOVERB:27,          PLOTOVERD:25,PLOTOVERQ:26,PLOTOVERB:27,
         MEMORY_IFPLOT:28,MEMORY_CONPLOT:29,MEMORY_PLOT:30,ARRAYPLOT:31,          MEMORY_IFPLOT:28,MEMORY_CONPLOT:29,MEMORY_PLOT:30,ARRAYPLOT:31,
         OPEN_CANVAS:32,DRAWCIRCLE:33,DRAW_OBJ:34,DRAW_STRING:35,OBJ_CP:36,          OPEN_CANVAS:32,DRAWCIRCLE:33,DRAW_OBJ:34,DRAW_STRING:35,OBJ_CP:36,
         CLEAR_CANVAS:37          CLEAR_CANVAS:37,POLARPLOTD:38
 */  */
 int modeNO(char *fn){  int modeNO(char *fn){
         int i;          int i;
Line 292  void ifplot_main(NODE arg,int is_memory,char *fn,Obj *
Line 294  void ifplot_main(NODE arg,int is_memory,char *fn,Obj *
                 Pox_pop_cmo(arg,rp);                  Pox_pop_cmo(arg,rp);
         } else {          } else {
                 MKSTR(fname,fn);                  MKSTR(fname,fn);
 /*  
 printf("%s\n",fn);  
 ListCheck("xrange",xrange);  
 ListCheck("yrange",yrange);  
 ListCheck("geom",geom);  
 */  
                 arg=mknode(8,s_id,fname,poly,xrange,yrange,zrange,geom,wname);                  arg=mknode(8,s_id,fname,poly,xrange,yrange,zrange,geom,wname);
                 Pox_rpc(arg,&t);                  Pox_rpc(arg,&t);
                 *rp=(Obj)s_id;                  *rp=(Obj)s_id;
         }          }
 }  }
   
 void Pobj_cp(NODE arg,Obj *rp){  
         //copy canvas  
         //1:and,3:copy,6:xor,7:or  
         Q sysid,index_A,index_B,op_code;  
         STRING fname;  
         Obj t;  
         sysid=(Q)ARG0(arg);  
         index_A=(Q)ARG1(arg);  
         index_B=(Q)ARG2(arg);  
         op_code=(Q)ARG3(arg);  
         MKSTR(fname,OBJ_CP);  
         arg=mknode(5,sysid,fname,index_A,index_B,op_code);  
         Pox_rpc(arg,&t);  
         *rp=(Obj)sysid;  
 }  
   
 void Pconplot(NODE arg,Obj *rp){conplot_main(arg, 0, rp);}  void Pconplot(NODE arg,Obj *rp){conplot_main(arg, 0, rp);}
 void Pmemory_conplot(NODE arg,Obj *rp){conplot_main(arg,1,rp);}  void Pmemory_conplot(NODE arg,Obj *rp){conplot_main(arg,1,rp);}
   
 void conplot_main(NODE arg,int is_memory,Obj *rp){  void conplot_main(NODE arg,int is_memory,Obj *rp){
         int stream,ri,i;          Q m2,p2,w300,s_id;
         Q m2,p2,w300,s_id,mxgc;          NODE defrange;
         NODE defrange,n,n0,n1,n2,n3;  
         LIST xrange,yrange,zrange,range[3],list,geom;          LIST xrange,yrange,zrange,range[3],list,geom;
         VL vl,vl0;          VL vl,vl0;
         V v[2],av[2];          V v[2],av[2];
         P poly,var;          int stream,ri,i;
           P poly;
           P var;
           NODE n,n0;
         STRING fname,wname;          STRING fname,wname;
         Obj t;          Obj t;
   
         STOQ(-2,m2);STOQ(2,p2);MKNODE(n,p2,0);MKNODE(defrange,m2,n);          STOQ(-2,m2); STOQ(2,p2);
         poly=0;vl=0;geom=0;wname=0;stream=-1;ri=0;v[0]=v[1]=0;          MKNODE(n,p2,0); MKNODE(defrange,m2,n);
         for(;arg;arg=NEXT(arg))          poly = 0; vl = 0; geom = 0; wname = 0; stream = -1; ri = 0;
                 if(!BDY(arg)) stream=0;          v[0] = v[1] = 0;
                 else switch (OID(BDY(arg))){          for ( ; arg; arg = NEXT(arg) )
                   if ( !BDY(arg) )
                           stream = 0;
                   else
                   switch ( OID(BDY(arg)) ) {
                         case O_P:                          case O_P:
                                 poly=(P)BDY(arg);                                  poly = (P)BDY(arg);
                                 get_vars_recursive((Obj)poly,&vl);                                  get_vars_recursive((Obj)poly,&vl);
                                 for(vl0=vl,i=0;vl0;vl0=NEXT(vl0)){                                  for ( vl0 = vl, i = 0; vl0; vl0 = NEXT(vl0) )
                                         if(vl0->v->attr==(pointer)V_IND){                                          if ( vl0->v->attr == (pointer)V_IND )
                                                 if(i>=2){                                                  if ( i >= 2 )
                                                         error("conplot : invalid argument");                                                          error("ifplot : invalid argument");
                                                 } else v[i++]=vl0->v;                                                  else
                                         }                                                          v[i++] = vl0->v;
                                 }  
                                 break;                                  break;
                         case O_LIST:                          case O_LIST:
                                 list=(LIST)BDY(arg);                                  list = (LIST)BDY(arg);
                                 if(OID(BDY(BDY(list)))==O_P)                                  if ( OID(BDY(BDY(list))) == O_P )
                                         if(ri>2) error("conplot : invalid argument");                                          if ( ri > 2 )
                                         else range[ri++]=list;                                                  error("ifplot : invalid argument");
                                 else geom=list;                                          else
                                                   range[ri++] = list;
                                   else
                                           geom = list;
                                 break;                                  break;
                         case O_N:                          case O_N:
                                 stream=QTOS((Q)BDY(arg));                                  stream = QTOS((Q)BDY(arg)); break;
                                 break;  
                         case O_STR:                          case O_STR:
                                 wname=(STRING)BDY(arg);                                  wname = (STRING)BDY(arg); break;
                                 break;  
                         default:                          default:
                                 error("conplot : invalid argument");                                  error("ifplot : invalid argument"); break;
                                 break;  
                 }                  }
         if(!poly) error("conplot : invalid argument");          if ( !poly )
         // list format var,num,num[,num]                  error("ifplot : invalid argument");
         switch (ri){          switch ( ri ) {
                 case 0:                  case 0:
                         if(!v[1]) error("conplot : please specify all variables");                          if ( !v[1] )
                         MKV(v[0],var);MKNODE(n,var,defrange);MKLIST(xrange,n);                                  error("ifplot : please specify all variables");
                         MKV(v[1],var);MKNODE(n,var,defrange);MKLIST(yrange,n);                          MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n);
                         STOQ(MAXGC,mxgc);MKNODE(n1,mxgc,0);MKNODE(n2,p2,n1);MKNODE(n3,p2,n2);                          MKV(v[1],var); MKNODE(n,var,defrange); MKLIST(yrange,n);
                         MKLIST(zrange,n3);                          MKNODE(n,0,defrange); MKLIST(zrange,n);
                         break;                          break;
                 case 1:                  case 1:
                         if(!v[1]) error("conplot : please specify all variables");                          if ( !v[1] )
                         av[0]=VR((P)BDY(BDY(range[0])));                                  error("ifplot : please specify all variables");
                         if(v[0]==av[0]){                          av[0] = VR((P)BDY(BDY(range[0])));
                                 xrange=range[0];                          if ( v[0] == av[0] ) {
                                 MKV(v[1],var);MKNODE(n,var,defrange);MKLIST(yrange,n);                                  xrange = range[0];
                                 STOQ(MAXGC,mxgc);MKNODE(n1,mxgc,0);MKNODE(n2,p2,n1);MKNODE(n3,p2,n2);                                  MKV(v[1],var); MKNODE(n,var,defrange); MKLIST(yrange,n);
                                 MKLIST(zrange,n3);                                  MKNODE(n,0,defrange); MKLIST(zrange,n);
                         } else if(v[1]==av[0]){                          } else if ( v[1] == av[0] ) {
                                 MKV(v[0],var);MKNODE(n,var,defrange);MKLIST(xrange,n);                                  MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n);
                                 yrange=range[0];                                  yrange = range[0];
                                 STOQ(MAXGC,mxgc);MKNODE(n1,mxgc,0);MKNODE(n2,p2,n1);MKNODE(n3,p2,n2);                                  MKNODE(n,0,defrange); MKLIST(zrange,n);
                                 MKLIST(zrange,n3);  
                         } else {                          } else {
                                 MKV(v[0],var);MKNODE(n,var,defrange);MKLIST(xrange,n);                                  MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(xrange,n);
                                 MKV(v[1],var);MKNODE(n,var,defrange);MKLIST(yrange,n);                                  MKV(v[1],var); MKNODE(n,var,defrange); MKLIST(yrange,n);
                                 n=NEXT((NODE)BDY(range[0]));                                  zrange = range[0];
                                 for(i=0,n1=n;n1;n1=NEXT(n1),i++);  
                                 if(i==3)MKLIST(zrange,n);  
                                 else if(i==2){  
                                         n1=(NODE)BDY(n);  
                                         n2=(NODE)(BDY(NEXT(n)));  
                                         STOQ(MAXGC,mxgc);MKNODE(n,mxgc,0);  
                                         MKNODE(n3,n1,n);  
                                         MKNODE(n1,n2,n3);  
                                         MKLIST(zrange,n1);  
                                 } else error("conplot : zrange error");  
                         }                          }
                         break;                          break;
                 case 2:                  case 2: case 3:
                 case 3:                          av[0] = VR((P)BDY(BDY(range[0])));
                         av[0]=VR((P)BDY(BDY(range[0])));                          av[1] = VR((P)BDY(BDY(range[1])));
                         av[1]=VR((P)BDY(BDY(range[1])));                          if ( ((v[0] == av[0]) && (!v[1] || v[1] == av[1])) ||
                         if(((v[0]==av[0])&&(!v[1]||v[1]==av[1]))||                                   ((v[0] == av[1]) && (!v[1] || v[1] == av[0])) ) {
                                 ((v[0]==av[1])&&(!v[1]||v[1]==av[0]))){                                  xrange = range[0]; yrange = range[1];
                                 xrange=range[0];                                  if ( ri == 3 )
                                 yrange=range[1];                                          zrange = range[2];
                                 if(ri==3) MKLIST(zrange,NEXT(BDY(range[2])));  
                                 else {                                  else {
                                         for(i=0,n=(NODE)BDY(range[0]);n;n=NEXT(n),i++);                                          MKNODE(n,0,defrange); MKLIST(zrange,n);
                                         MKNODE(n,0,defrange);  
                                         MKLIST(zrange,n);  
                                 }                                  }
                         } else error("conplot : invalid argument");                          } else
                                   error("ifplot : invalid argument");
                         break;                          break;
                 default:                  default:
                         error("conplot : cannot happen");                          error("ifplot : cannot happen"); break;
                         break;  
         }          }
         /* conplot in ox_plot requires          /* conplot in ox_plot requires
            [s_id (Q),             [s_id (Q),
                 formula (Obj),                  formula (Obj),
                 xrange=[x,xmin,xmax] (LIST),                  xrange=[x,xmin,xmax] (LIST),
                 yrange=[y,ymin,ymax] (LIST),                  yrange=[y,ymin,ymax] (LIST),
                 zrange=[z,zmin,zmax,nstep] (LIST),z:dummy var                  zrange=[z,zmin,zmax] (LIST),
                 geom=[xsize,ysize] (LIST),                  geom=[xsize,ysize] (LIST),
                 wname=name (STRING)]                  wname=name (STRING)]
         */          */
         stream=validate_ox_plot_stream(stream);  
           stream = validate_ox_plot_stream(stream);
         STOQ(stream,s_id);          STOQ(stream,s_id);
         if(!geom){          if ( !geom ) {
                 STOQ(300,w300);                  STOQ(300,w300);
                 MKNODE(n0,w300,0);                  MKNODE(n0,w300,0); MKNODE(n,w300,n0); MKLIST(geom,n);
                 MKNODE(n,w300,n0);  
                 MKLIST(geom,n);  
         }          }
         if(is_memory){          if ( is_memory ) {
                 MKSTR(fname,MEMORY_PLOT);                  MKSTR(fname,"memory_plot");
                 arg=mknode(7,s_id,fname,poly,xrange,yrange,zrange,geom);                  arg = mknode(7,s_id,fname,poly,xrange,yrange,zrange,geom);
                 Pox_rpc(arg,&t);                  Pox_rpc(arg,&t);
                 arg=mknode(1,s_id);                  arg = mknode(1,s_id);
                 Pox_pop_cmo(arg,rp);                  Pox_pop_cmo(arg,rp);
         } else {          } else {
                 //list check  
                 for(i=0,n=(NODE)BDY(xrange);n;i++,n=NEXT(n));  
                 if(i!=3)error("conplot : xrange error");  
                 for(i=0,n=(NODE)BDY(yrange);n;i++,n=NEXT(n));  
                 if(i!=3)error("conplot : yrange error");  
                 for(i=0,n=(NODE)BDY(zrange);n;i++,n=NEXT(n));  
                 if(i!=3)error("conplot : xrange error");  
                 MKSTR(fname,CONPLOT);                  MKSTR(fname,CONPLOT);
                 arg=mknode(8,s_id,fname,poly,xrange,yrange,zrange,geom,wname);                  arg = mknode(8,s_id,fname,poly,xrange,yrange,zrange,geom,wname);
                 Pox_rpc(arg,&t);                  Pox_rpc(arg,&t);
                 *rp=(Obj)s_id;                  *rp = (Obj)s_id;
         }          }
 }  }
   
Line 562  void plot_main(NODE arg,int is_memory,char *fn,Obj *rp
Line 525  void plot_main(NODE arg,int is_memory,char *fn,Obj *rp
 #define Pi 3.14159265358979323846264  #define Pi 3.14159265358979323846264
   
 void Ppolarplot(NODE arg,Obj *rp){  void Ppolarplot(NODE arg,Obj *rp){
         Q m2,p2,w300,defstep,s_id,color;          Q m2,p2,w300,s_id;
         NODE defrange,n,n0,n1,n2;          NODE defrange,n,n0;
         LIST range,geom,list[2];          LIST zrange,range[1],geom,list;
         VL vl,vl0;          VL vl,vl0;
         V v[1],av;          V v[1],av[1];
         int stream,ri,i,len,iNo,lNo,vNo,sNo,pfine,findG;          int stream,ri,i;
         P poly,var;          P poly,var;
         STRING fname,wname;          STRING fname,wname;
         Real pi2;          Real pi2;
         Obj t,frst,sec,thr;          Obj t;
         char ebuf[BUFSIZ];  
   
         iNo=lNo=sNo=findG=0;pfine=TRUE;          MKReal(2*Pi,pi2);MKNODE(n,pi2,0); MKNODE(defrange,0,n);
         poly=0;vl=0;geom=0;wname=0;color=0;stream=-1;ri=0;v[0]=0;          poly=0;vl=0;geom=0;wname=0;stream=-1;ri=0;v[0]=0;
         for(;arg;arg=NEXT(arg)){          for(;arg;arg=NEXT(arg)){
                 if(!BDY(arg)) iNo++;                  if(!BDY(arg)) stream=0;
                 else switch(OID(BDY(arg))){                  else switch(OID(BDY(arg))){
                 case O_P: case O_R://formular                  case O_P: case O_R://formular
                         poly=(P)BDY(arg);                          poly=(P)BDY(arg);
                         get_vars_recursive((Obj)poly,&vl);                          get_vars_recursive((Obj)poly,&vl);
                         for(vl0=vl,vNo=0;vl0;vl0=NEXT(vl0)){                          for(vl0=vl,i=0;vl0;vl0=NEXT(vl0))
                                 if(vl0->v->attr==(pointer)V_IND){                                  if(vl0->v->attr==(pointer)V_IND)
                                         if(vNo>=1){                                          if(i>=1)error("polarplot : invalid argument");
                                                 sprintf(ebuf,"%s : invalaid argument",POLARPLOT);                                  else v[i++]=vl0->v;
                                                 error(ebuf);                          if(i!=1)error("polarplot : invalid argument");
                                         } else v[vNo++]=vl0->v;  
                                 }  
                         }  
                         if(vNo!=1){  
                                 sprintf(ebuf,"%s : only uni-variate formular",POLARPLOT);  
                                 error(ebuf);  
                         }  
                         break;                          break;
                 case O_N://color,id,idx,division  
                         switch (iNo){  
                         case 0://color arg  
                                 color=(Q)BDY(arg);  
                                 iNo++;  
                                 break;  
                         case 1://stream arg  
                                 stream=QTOS((Q)BDY(arg));  
                                 iNo++;  
                                 break;  
                         default://error  
                                 sprintf(ebuf,"%s : invalid number arguments",POLARPLOT);  
                                 error(ebuf);  
                                 break;  
                         }  
                         break;  
                 case O_LIST://range,geomerty                  case O_LIST://range,geomerty
                         if(lNo<2)list[lNo++]=(LIST)BDY(arg);                          list=(LIST)BDY(arg);
                         else {                          if(OID(BDY(BDY(list)))==O_P)
                                 sprintf(ebuf,"%s : invalid list argument",POLARPLOT);                                  if(ri>0)error("polarplot : invalid argument");
                                 error(ebuf);                                  else range[ri++]=list;
                         }                          else geom=list;
                         break;                          break;
                   case O_N:
                           stream=QTOS((Q)BDY(arg));
                           break;
                 case O_STR://wname                  case O_STR://wname
                         wname=(STRING)BDY(arg);                          wname=(STRING)BDY(arg);
                         sNo++;  
                         break;                          break;
                 default:                  default:
                           error("polarplot : invalid argument");
                         break;                          break;
                 }                  }
         }          }
         //formular check          //formular check
         if(!poly){          if(!poly)error("polarplot : invalid argument");
                 sprintf(ebuf,"%s : invalid ploy argument",POLARPLOT);          switch (ri){
                 error(ebuf);          case 0:
                   MKV(v[0],var); MKNODE(n,var,defrange); MKLIST(zrange,n);
                   break;
           case 1:
                   av[0]=VR((P)BDY(BDY(range[0])));
                   if(v[0]==av[0]) zrange = range[0];
                   else error("polarplot : invalid argument");
                   break;
           default:
                   error("polarplot : cannot happen");
                   break;
         }          }
         //vars check  
         get_vars_recursive((Obj)poly,&vl);  
         for(vl0=vl,vNo=0;vl0;vNo++,vl0=NEXT(vl0)){  
                 if(vl0->v->attr==(pointer)V_IND){  
                         if(vNo>=2){  
                                 sprintf(ebuf,"%s : invalid ploy argument",POLARPLOT);  
                                 error(ebuf);  
                         } else v[vNo]=vl0->v;  
                 }  
         }  
         //list check  
         for(i=0;i<lNo;i++){  
                 if(OID(BDY(BDY(list[i])))!=O_P){  
                         // list first value is number (geometry)  
                         for(len=0,n=(NODE)BDY(list[i]);n;len++,n=NEXT(n)){  
                                 if(len==0) frst=BDY(n);  
                                 else if(len==1) sec=BDY(n);  
                                 else {  
                                         sprintf(ebuf,"%s : geometry list too long",POLARPLOT);  
                                         error(ebuf);  
                                 }  
                         }  
                         if(len!=2){  
                                 sprintf(ebuf,"%s : geometry requierd 2 numbers", POLARPLOT);  
                                 error(ebuf);  
                         } else geom=list[i];  
                 } else {  
                         //list first value is var (range)  
                         av=VR((P)BDY(BDY(list[i])));  
                         if(v[0]==av)range=list[i];  
                         else {  
                                 sprintf(ebuf,"%s : invalid list length",POLARPLOT);  
                                 error(ebuf);  
                         }  
                 }  
         }  
         // set default  
         if(!range){  
                 STOQ(DEFAULTPOLARSTEP,defstep);MKReal(2*Pi,pi2);MKV(v[0],var);  
                 MKNODE(n,defstep,0);MKNODE(n1,pi2,n);MKNODE(n2,0,n1);  
                 MKNODE(defrange,var,n2);MKLIST(range,defrange);  
         }  
         if(!geom){  
                 STOQ(300,w300);MKNODE(n0,w300,0);MKNODE(n,w300,n0);MKLIST(geom,n);  
         }  
         stream=validate_ox_plot_stream(stream);          stream=validate_ox_plot_stream(stream);
         STOQ(stream,s_id);          STOQ(stream,s_id);
           if(!geom){
                   STOQ(300,w300);
                   MKNODE(n0,w300,0); MKNODE(n,w300,n0); MKLIST(geom,n);
           }
         MKSTR(fname,POLARPLOT);          MKSTR(fname,POLARPLOT);
         arg=mknode(7,s_id,fname,poly,color,range,geom,wname);          arg=mknode(8,s_id,fname,poly,NULLP,NULLP,zrange,geom,wname);
         Pox_rpc(arg,&t);          Pox_rpc(arg,&t);
         *rp=(Obj)s_id;          *rp=s_id;
 }  }
   
 void Pplotover(NODE arg,Obj *rp){  void Pplotover(NODE arg,Obj *rp){
Line 826  void PplotoverD(NODE arg,Obj *rp){plotover_mainNG(arg,
Line 737  void PplotoverD(NODE arg,Obj *rp){plotover_mainNG(arg,
 void PplotoverQ(NODE arg,Obj *rp){plotover_mainNG(arg,PLOTOVERQ,rp);}  void PplotoverQ(NODE arg,Obj *rp){plotover_mainNG(arg,PLOTOVERQ,rp);}
 void PplotoverB(NODE arg,Obj *rp){plotover_mainNG(arg,PLOTOVERB,rp);}  void PplotoverB(NODE arg,Obj *rp){plotover_mainNG(arg,PLOTOVERB,rp);}
   
   void Pobj_cp(NODE arg,Obj *rp){
           //copy canvas
           //1:and,3:copy,6:xor,7:or
           Q sysid,index_A,index_B,op_code;
           STRING fname;
           Obj t;
           sysid=(Q)ARG0(arg);
           index_A=(Q)ARG1(arg);
           index_B=(Q)ARG2(arg);
           op_code=(Q)ARG3(arg);
           MKSTR(fname,OBJ_CP);
           arg=mknode(5,sysid,fname,index_A,index_B,op_code);
           Pox_rpc(arg,&t);
           *rp=(Obj)sysid;
   }
   
   void PpolarplotD(NODE arg,Obj *rp){
           Q m2,p2,w300,defstep,s_id,color;
           NODE defrange,n,n0,n1,n2;
           LIST range,geom,list[2];
           VL vl,vl0;
           V v[1],av;
           int stream,ri,i,len,iNo,lNo,vNo,sNo,pfine,findG;
           P poly,var;
           STRING fname,wname;
           Real pi2;
           Obj t,frst,sec,thr;
           char ebuf[BUFSIZ];
   
           iNo=lNo=sNo=findG=0;pfine=TRUE;
           poly=0;vl=0;geom=0;wname=0;color=0;stream=-1;ri=0;v[0]=0;
           for(;arg;arg=NEXT(arg)){
                   if(!BDY(arg)) iNo++;
                   else switch(OID(BDY(arg))){
                   case O_P: case O_R://formular
                           poly=(P)BDY(arg);
                           get_vars_recursive((Obj)poly,&vl);
                           for(vl0=vl,vNo=0;vl0;vl0=NEXT(vl0)){
                                   if(vl0->v->attr==(pointer)V_IND){
                                           if(vNo>=1){
                                                   sprintf(ebuf,"%s : invalaid argument",POLARPLOT);
                                                   error(ebuf);
                                           } else v[vNo++]=vl0->v;
                                   }
                           }
                           if(vNo!=1){
                                   sprintf(ebuf,"%s : only uni-variate formular",POLARPLOT);
                                   error(ebuf);
                           }
                           break;
                   case O_N://color,id,idx,division
                           switch (iNo){
                           case 0://color arg
                                   color=(Q)BDY(arg);
                                   iNo++;
                                   break;
                           case 1://stream arg
                                   stream=QTOS((Q)BDY(arg));
                                   iNo++;
                                   break;
                           default://error
                                   sprintf(ebuf,"%s : invalid number arguments",POLARPLOT);
                                   error(ebuf);
                                   break;
                           }
                           break;
                   case O_LIST://range,geomerty
                           if(lNo<2)list[lNo++]=(LIST)BDY(arg);
                           else {
                                   sprintf(ebuf,"%s : invalid list argument",POLARPLOT);
                                   error(ebuf);
                           }
                           break;
                   case O_STR://wname
                           wname=(STRING)BDY(arg);
                           sNo++;
                           break;
                   default:
                           break;
                   }
           }
           //formular check
           if(!poly){
                   sprintf(ebuf,"%s : invalid ploy argument",POLARPLOT);
                   error(ebuf);
           }
           //vars check
           get_vars_recursive((Obj)poly,&vl);
           for(vl0=vl,vNo=0;vl0;vNo++,vl0=NEXT(vl0)){
                   if(vl0->v->attr==(pointer)V_IND){
                           if(vNo>=2){
                                   sprintf(ebuf,"%s : invalid ploy argument",POLARPLOT);
                                   error(ebuf);
                           } else v[vNo]=vl0->v;
                   }
           }
           //list check
           for(i=0;i<lNo;i++){
                   if(OID(BDY(BDY(list[i])))!=O_P){
                           // list first value is number (geometry)
                           for(len=0,n=(NODE)BDY(list[i]);n;len++,n=NEXT(n)){
                                   if(len==0) frst=BDY(n);
                                   else if(len==1) sec=BDY(n);
                                   else {
                                           sprintf(ebuf,"%s : geometry list too long",POLARPLOT);
                                           error(ebuf);
                                   }
                           }
                           if(len!=2){
                                   sprintf(ebuf,"%s : geometry requierd 2 numbers", POLARPLOT);
                                   error(ebuf);
                           } else geom=list[i];
                   } else {
                           //list first value is var (range)
                           av=VR((P)BDY(BDY(list[i])));
                           if(v[0]==av)range=list[i];
                           else {
                                   sprintf(ebuf,"%s : invalid list length",POLARPLOT);
                                   error(ebuf);
                           }
                   }
           }
           // set default
           if(!range){
                   STOQ(DEFAULTPOLARSTEP,defstep);MKReal(2*Pi,pi2);MKV(v[0],var);
                   MKNODE(n,defstep,0);MKNODE(n1,pi2,n);MKNODE(n2,0,n1);
                   MKNODE(defrange,var,n2);MKLIST(range,defrange);
           }
           if(!geom){
                   STOQ(300,w300);MKNODE(n0,w300,0);MKNODE(n,w300,n0);MKLIST(geom,n);
           }
           stream=validate_ox_plot_stream(stream);
           STOQ(stream,s_id);
           MKSTR(fname,POLARPLOTD);
           arg=mknode(7,s_id,fname,poly,color,range,geom,wname);
           Pox_rpc(arg,&t);
           *rp=(Obj)s_id;
   }
   
 void ifplot_mainNG(NODE arg,char *fn,Obj *rp){  void ifplot_mainNG(NODE arg,char *fn,Obj *rp){
         Q m2,p2,w300,mxgc,s_id,color;          Q m2,p2,w300,mxgc,s_id,color;
         NODE defrange,n,n0,n1,n2;          NODE defrange,n,n0,n1,n2;
Line 968  void ifplot_mainNG(NODE arg,char *fn,Obj *rp){
Line 1018  void ifplot_mainNG(NODE arg,char *fn,Obj *rp){
         stream=validate_ox_plot_stream(stream);          stream=validate_ox_plot_stream(stream);
         STOQ(stream,s_id);          STOQ(stream,s_id);
         MKSTR(fname,fn);          MKSTR(fname,fn);
 /*  
 printf("%s\n",fn);  
 ListCheck("xrange",xrange);  
 ListCheck("yrange",yrange);  
 ListCheck("zrange",zrange);  
 ListCheck("geom",geom);  
 printf("idv %d\n",idv);  
 */  
         arg=mknode(9,s_id,fname,poly,color,xrange,yrange,zrange,geom,wname);          arg=mknode(9,s_id,fname,poly,color,xrange,yrange,zrange,geom,wname);
         Pox_rpc(arg,&t);          Pox_rpc(arg,&t);
         *rp=(Obj)s_id;          *rp=(Obj)s_id;
 }  }
   
 void plotover_mainNG(NODE arg,char *fn,Obj *rp){  void plotover_mainNG(NODE arg,char *fn,Obj *rp){
         //[s_id (Q), w_id (Q), formula (Obj), color (Q)]  
         Q s_id,w_id,color;          Q s_id,w_id,color;
         P poly;          P poly;
         STRING fname;          STRING fname;
Line 1030  void plotover_mainNG(NODE arg,char *fn,Obj *rp){
Line 1071  void plotover_mainNG(NODE arg,char *fn,Obj *rp){
                 }                  }
         }          }
         MKSTR(fname,fn);          MKSTR(fname,fn);
           //[s_id (Q), w_id (Q), formula (Obj), color (Q)]
         arg=mknode(5,s_id,fname,w_id,poly,color);          arg=mknode(5,s_id,fname,w_id,poly,color);
 //printf("fn:%s s_id:%d w_id:%d color:%d\n",BDY(fname),QTOS(s_id),QTOS(w_id),QTOS(color));  
         Pox_rpc(arg,&t);          Pox_rpc(arg,&t);
         *rp=(Obj)s_id;          *rp=(Obj)s_id;
 }  }
Legend:
Removed from v.1.28  
changed lines
  Added in v.1.29
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>