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

Diff for /OpenXM_contrib2/asir2000/builtin/strobj.c between version 1.4 and 1.14

version 1.4, 2000/02/07 06:04:56 version 1.14, 2004/03/03 09:25:30
Line 1 
Line 1 
 /* $OpenXM: OpenXM_contrib2/asir2000/builtin/strobj.c,v 1.3 2000/02/07 05:21:32 noro Exp $ */  /*
    * Copyright (c) 1994-2000 FUJITSU LABORATORIES LIMITED
    * All rights reserved.
    *
    * FUJITSU LABORATORIES LIMITED ("FLL") hereby grants you a limited,
    * non-exclusive and royalty-free license to use, copy, modify and
    * redistribute, solely for non-commercial and non-profit purposes, the
    * computer program, "Risa/Asir" ("SOFTWARE"), subject to the terms and
    * conditions of this Agreement. For the avoidance of doubt, you acquire
    * only a limited right to use the SOFTWARE hereunder, and FLL or any
    * third party developer retains all rights, including but not limited to
    * copyrights, in and to the SOFTWARE.
    *
    * (1) FLL does not grant you a license in any way for commercial
    * purposes. You may use the SOFTWARE only for non-commercial and
    * non-profit purposes only, such as academic, research and internal
    * business use.
    * (2) The SOFTWARE is protected by the Copyright Law of Japan and
    * international copyright treaties. If you make copies of the SOFTWARE,
    * with or without modification, as permitted hereunder, you shall affix
    * to all such copies of the SOFTWARE the above copyright notice.
    * (3) An explicit reference to this SOFTWARE and its copyright owner
    * shall be made on your publication or presentation in any form of the
    * results obtained by use of the SOFTWARE.
    * (4) In the event that you modify the SOFTWARE, you shall notify FLL by
    * e-mail at risa-admin@sec.flab.fujitsu.co.jp of the detailed specification
    * for such modification or the source code of the modified part of the
    * SOFTWARE.
    *
    * THE SOFTWARE IS PROVIDED AS IS WITHOUT ANY WARRANTY OF ANY KIND. FLL
    * MAKES ABSOLUTELY NO WARRANTIES, EXPRESSED, IMPLIED OR STATUTORY, AND
    * EXPRESSLY DISCLAIMS ANY IMPLIED WARRANTY OF MERCHANTABILITY, FITNESS
    * FOR A PARTICULAR PURPOSE OR NONINFRINGEMENT OF THIRD PARTIES'
    * RIGHTS. NO FLL DEALER, AGENT, EMPLOYEES IS AUTHORIZED TO MAKE ANY
    * MODIFICATIONS, EXTENSIONS, OR ADDITIONS TO THIS WARRANTY.
    * UNDER NO CIRCUMSTANCES AND UNDER NO LEGAL THEORY, TORT, CONTRACT,
    * OR OTHERWISE, SHALL FLL BE LIABLE TO YOU OR ANY OTHER PERSON FOR ANY
    * DIRECT, INDIRECT, SPECIAL, INCIDENTAL, PUNITIVE OR CONSEQUENTIAL
    * DAMAGES OF ANY CHARACTER, INCLUDING, WITHOUT LIMITATION, DAMAGES
    * ARISING OUT OF OR RELATING TO THE SOFTWARE OR THIS AGREEMENT, DAMAGES
    * FOR LOSS OF GOODWILL, WORK STOPPAGE, OR LOSS OF DATA, OR FOR ANY
    * DAMAGES, EVEN IF FLL SHALL HAVE BEEN INFORMED OF THE POSSIBILITY OF
    * SUCH DAMAGES, OR FOR ANY CLAIM BY ANY OTHER PARTY. EVEN IF A PART
    * OF THE SOFTWARE HAS BEEN DEVELOPED BY A THIRD PARTY, THE THIRD PARTY
    * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
    * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
    *
    * $OpenXM: OpenXM_contrib2/asir2000/builtin/strobj.c,v 1.13 2004/02/26 10:07:55 noro Exp $
   */
 #include "ca.h"  #include "ca.h"
 #include "parse.h"  #include "parse.h"
 #include "ctype.h"  #include "ctype.h"
 #if PARI  #if defined(PARI)
 #include "genpari.h"  #include "genpari.h"
   #  if !(PARI_VERSION_CODE > 131588)
 extern jmp_buf environnement;  extern jmp_buf environnement;
   #  endif
 #endif  #endif
   #include <string.h>
   
 extern char *parse_strp;  extern char *parse_strp;
   
 void Prtostr(), Pstrtov(), Peval_str();  void Prtostr(), Pstrtov(), Peval_str();
 void Pstrtoascii(), Pasciitostr();  void Pstrtoascii(), Pasciitostr();
   void Pstr_len(), Pstr_chr(), Psub_str();
   void Pwrite_to_tb();
   void Ptb_to_string();
   void Pclear_tb();
   void Pstring_to_tb();
   void Pquotetotex_tb();
   void Pquotetotex();
   void fnodetotex_tb(FNODE f,TB tb);
   char *symbol_name(char *name);
   void tb_to_string(TB tb,STRING *rp);
   void fnodenodetotex_tb(NODE n,TB tb);
   void fargstotex_tb(char *opname,FNODE f,TB tb);
   
 struct ftab str_tab[] = {  struct ftab str_tab[] = {
         {"rtostr",Prtostr,1},          {"rtostr",Prtostr,1},
Line 17  struct ftab str_tab[] = {
Line 81  struct ftab str_tab[] = {
         {"eval_str",Peval_str,1},          {"eval_str",Peval_str,1},
         {"strtoascii",Pstrtoascii,1},          {"strtoascii",Pstrtoascii,1},
         {"asciitostr",Pasciitostr,1},          {"asciitostr",Pasciitostr,1},
           {"str_len",Pstr_len,1},
           {"str_chr",Pstr_chr,3},
           {"sub_str",Psub_str,3},
           {"write_to_tb",Pwrite_to_tb,2},
           {"clear_tb",Pclear_tb,1},
           {"tb_to_string",Ptb_to_string,1},
           {"string_to_tb",Pstring_to_tb,1},
           {"quotetotex_tb",Pquotetotex_tb,2},
           {"quotetotex",Pquotetotex,1},
         {0,0,0},          {0,0,0},
 };  };
   
   void write_tb(char *s,TB tb)
   {
           if ( tb->next == tb->size ) {
                   tb->size *= 2;
                   tb->body = (char **)REALLOC(tb->body,tb->size*sizeof(char *));
           }
           tb->body[tb->next] = s;
           tb->next++;
   }
   
   void Pwrite_to_tb(NODE arg,Q *rp)
   {
           int i;
   
           asir_assert(ARG1(arg),O_TB,"write_to_tb");
           write_tb(BDY((STRING)ARG0(arg)),ARG1(arg));
           *rp = 0;
   }
   
   void Pquotetotex(NODE arg,STRING *rp)
   {
           TB tb;
   
           NEWTB(tb);
           fnodetotex_tb(BDY((QUOTE)ARG0(arg)),tb);
           tb_to_string(tb,rp);
   }
   
   void Pquotetotex_tb(NODE arg,Q *rp)
   {
           int i;
           TB tb;
   
           asir_assert(ARG1(arg),O_TB,"quotetotex_tb");
           fnodetotex_tb(BDY((QUOTE)ARG0(arg)),ARG1(arg));
           *rp = 0;
   }
   
   void Pstring_to_tb(NODE arg,TB *rp)
   {
           TB tb;
   
           asir_assert(ARG0(arg),O_STR,"string_to_tb");
           NEWTB(tb);
           tb->body[0] = BDY((STRING)ARG0(arg));
           tb->next++;
           *rp = tb;
   }
   
   void Ptb_to_string(NODE arg,STRING *rp)
   {
           TB tb;
   
           asir_assert(ARG0(arg),O_TB,"tb_to_string");
           tb = (TB)ARG0(arg);
           tb_to_string(tb,rp);
   }
   
   void tb_to_string(TB tb,STRING *rp)
   {
           int j,len;
           char *all,*p,*q;
   
           for ( j = 0, len = 0; j < tb->next; j++ )
                   len += strlen(tb->body[j]);
           all = (char *)MALLOC_ATOMIC((len+1)*sizeof(char));
           for ( j = 0, p = all; j < tb->next; j++ )
                   for ( q = tb->body[j]; *q; *p++ = *q++ );
           *p = 0;
           MKSTR(*rp,all);
   }
   
   void Pclear_tb(NODE arg,Q *rp)
   {
           TB tb;
           int j;
   
           asir_assert(ARG0(arg),O_TB,"clear_tb");
           tb = (TB)ARG0(arg);
           for ( j = 0; j < tb->next; j++ )
                   tb->body[j] = 0;
           tb->next = 0;
           *rp = 0;
   }
   
   void Pstr_len(arg,rp)
   NODE arg;
   Q *rp;
   {
           STRING str;
           int r;
   
           str = (STRING)ARG0(arg);
           asir_assert(str,O_STR,"str_chr");
           r = strlen(BDY(str));
           STOQ(r,*rp);
   }
   
   void Pstr_chr(arg,rp)
   NODE arg;
   Q *rp;
   {
           STRING str,terminator;
           Q start;
           char *p,*ind;
           int chr,spos,r;
   
           str = (STRING)ARG0(arg);
           start = (Q)ARG1(arg);
           terminator = (STRING)ARG2(arg);
           asir_assert(str,O_STR,"str_chr");
           asir_assert(start,O_N,"str_chr");
           asir_assert(terminator,O_STR,"str_chr");
           p = BDY(str);
           spos = QTOS(start);
           chr = BDY(terminator)[0];
           if ( spos > (int)strlen(p) )
                   r = -1;
           else {
                   ind = strchr(p+spos,chr);
                   if ( ind )
                           r = ind-p;
                   else
                           r = -1;
           }
           STOQ(r,*rp);
   }
   
   void Psub_str(arg,rp)
   NODE arg;
   STRING *rp;
   {
           STRING str;
           Q head,tail;
           char *p,*r;
           int spos,epos,len;
   
           str = (STRING)ARG0(arg);
           head = (Q)ARG1(arg);
           tail = (Q)ARG2(arg);
           asir_assert(str,O_STR,"sub_str");
           asir_assert(head,O_N,"sub_str");
           asir_assert(tail,O_N,"sub_str");
           p = BDY(str);
           spos = QTOS(head);
           epos = QTOS(tail);
           len = strlen(p);
           if ( (spos >= len) || (epos < spos) ) {
                   *rp = 0; return;
           }
           if ( epos >= len )
                   epos = len-1;
           len = epos-spos+1;
           r = (char *)MALLOC(len+1);
           strncpy(r,p+spos,len);
           r[len] = 0;
           MKSTR(*rp,r);
   }
   
 void Pstrtoascii(arg,rp)  void Pstrtoascii(arg,rp)
 NODE arg;  NODE arg;
 LIST *rp;  LIST *rp;
Line 75  Obj *rp;
Line 307  Obj *rp;
 {  {
         FNODE fnode;          FNODE fnode;
         char *cmd;          char *cmd;
 #if PARI  #if defined(PARI)
           void recover(int);
   
         recover(0);          recover(0);
   #  if !(PARI_VERSION_CODE > 131588)
         if ( setjmp(environnement) ) {          if ( setjmp(environnement) ) {
                 avma = top; recover(1);                  avma = top; recover(1);
                 resetenv("");                  resetenv("");
         }          }
   #  endif
 #endif  #endif
         cmd = BDY((STRING)ARG0(arg));          cmd = BDY((STRING)ARG0(arg));
         exprparse(0,cmd,&fnode);          exprparse_create_var(0,cmd,&fnode);
         *rp = eval(fnode);          *rp = eval(fnode);
 }  }
   
Line 95  STRING *rp;
Line 331  STRING *rp;
         int len;          int len;
   
         len = estimate_length(CO,ARG0(arg));          len = estimate_length(CO,ARG0(arg));
         b = (char *)MALLOC(len+1);          b = (char *)MALLOC_ATOMIC(len+1);
         soutput_init(b);          soutput_init(b);
         sprintexpr(CO,ARG0(arg));          sprintexpr(CO,ARG0(arg));
         MKSTR(*rp,b);          MKSTR(*rp,b);
Line 105  void Pstrtov(arg,rp)
Line 341  void Pstrtov(arg,rp)
 NODE arg;  NODE arg;
 P *rp;  P *rp;
 {  {
         char *p,*t;          char *p;
   
         p = BDY((STRING)ARG0(arg));          p = BDY((STRING)ARG0(arg));
 #if 0  #if 0
Line 121  P *rp;
Line 357  P *rp;
 #else  #else
         makevar(p,rp);          makevar(p,rp);
 #endif  #endif
   }
   
   char *symbol_name(char *name)
   {
           /* XXX */
           return name;
   }
   
   void fnodetotex_tb(FNODE f,TB tb)
   {
           NODE n,t,t0;
           char vname[BUFSIZ];
           char *opname;
           Obj obj;
           int i,len;
   
           write_tb(" ",tb);
           if ( !f ) {
                   write_tb("0",tb);
                   return;
           }
           switch ( f->id ) {
                   /* unary operators */
                   case I_NOT: case I_PAREN: case I_MINUS:
                           switch ( f->id ) {
                                   case I_NOT:
                                           write_tb("\\neg (",tb);
                                           fnodetotex_tb((FNODE)FA0(f),tb);
                                           write_tb(")",tb);
                                           break;
                                   case I_PAREN:
                                           write_tb("(",tb);
                                           fnodetotex_tb((FNODE)FA0(f),tb);
                                           write_tb(")",tb);
                                           break;
                                   case I_MINUS:
                                           write_tb("-",tb);
                                           fnodetotex_tb((FNODE)FA0(f),tb);
                                           break;
                           }
                           break;
   
                   /* binary operators */
                   case I_BOP: case I_COP: case I_LOP: case I_AND: case I_OR:
                           /* arg list */
                           /* I_AND, I_OR => FA0(f), FA1(f) */
                           /* otherwise   => FA1(f), FA2(f) */
   
                           /* op */
                           switch ( f->id ) {
                                   case I_BOP:
                                           opname = ((ARF)FA0(f))->name;
                                           if ( !strcmp(opname,"+") ) {
                                                   fnodetotex_tb((FNODE)FA1(f),tb);
                                                   write_tb(opname,tb);
                                                   fnodetotex_tb((FNODE)FA2(f),tb);
                                           } else if ( !strcmp(opname,"-") ) {
                                                   if ( FA1(f) ) fnodetotex_tb((FNODE)FA1(f),tb);
                                                   write_tb(opname,tb);
                                                   fnodetotex_tb((FNODE)FA2(f),tb);
                                           } else if ( !strcmp(opname,"*") ) {
                                                   fnodetotex_tb((FNODE)FA1(f),tb);
                                                   write_tb(" ",tb);
                                                   fnodetotex_tb((FNODE)FA2(f),tb);
                                           } else if ( !strcmp(opname,"/") ) {
                                                   write_tb("\\frac{",tb);
                                                   fnodetotex_tb((FNODE)FA1(f),tb);
                                                   write_tb("} {",tb);
                                                   fnodetotex_tb((FNODE)FA2(f),tb);
                                                   write_tb("}",tb);
                                           } else if ( !strcmp(opname,"^") ) {
                                                   fnodetotex_tb((FNODE)FA1(f),tb);
                                                   write_tb("^{",tb);
                                                   fnodetotex_tb((FNODE)FA2(f),tb);
                                                   write_tb("} ",tb);
                                           } else if ( !strcmp(opname,"%") ) {
                                                   fnodetotex_tb((FNODE)FA1(f),tb);
                                                   write_tb(" {\\rm mod}\\, ",tb);
                                                   fnodetotex_tb((FNODE)FA2(f),tb);
                                           } else
                                                   error("invalid binary operator");
   
                                   case I_COP:
                                           switch( (cid)FA0(f) ) {
                                                   case C_EQ:
                                                           fnodetotex_tb((FNODE)FA1(f),tb);
                                                           write_tb(" = ",tb);
                                                           fnodetotex_tb((FNODE)FA2(f),tb);
                                                           break;
                                                   case C_NE:
                                                           fnodetotex_tb((FNODE)FA1(f),tb);
                                                           write_tb(" \\neq ",tb);
                                                           fnodetotex_tb((FNODE)FA2(f),tb);
                                                           break;
                                                   case C_GT:
                                                           fnodetotex_tb((FNODE)FA1(f),tb);
                                                           write_tb(" \\gt ",tb);
                                                           fnodetotex_tb((FNODE)FA2(f),tb);
                                                           break;
                                                   case C_LT:
                                                           fnodetotex_tb((FNODE)FA1(f),tb);
                                                           write_tb(" \\lt ",tb);
                                                           fnodetotex_tb((FNODE)FA2(f),tb);
                                                           break;
                                                   case C_GE:
                                                           fnodetotex_tb((FNODE)FA1(f),tb);
                                                           write_tb(" \\geq ",tb);
                                                           fnodetotex_tb((FNODE)FA2(f),tb);
                                                           break;
                                                   case C_LE:
                                                           fnodetotex_tb((FNODE)FA1(f),tb);
                                                           write_tb(" \\leq ",tb);
                                                           fnodetotex_tb((FNODE)FA2(f),tb);
                                                           break;
                                           }
                                           break;
   
                                   case I_LOP:
                                           switch( (lid)FA0(f) ) {
                                                   case L_EQ:
                                                           fnodetotex_tb((FNODE)FA1(f),tb);
                                                           write_tb(" = ",tb);
                                                           fnodetotex_tb((FNODE)FA2(f),tb);
                                                           break;
                                                   case L_NE:
                                                           fnodetotex_tb((FNODE)FA1(f),tb);
                                                           write_tb(" \\neq ",tb);
                                                           fnodetotex_tb((FNODE)FA2(f),tb);
                                                           break;
                                                   case L_GT:
                                                           fnodetotex_tb((FNODE)FA1(f),tb);
                                                           write_tb(" \\gt ",tb);
                                                           fnodetotex_tb((FNODE)FA2(f),tb);
                                                           break;
                                                   case L_LT:
                                                           fnodetotex_tb((FNODE)FA1(f),tb);
                                                           write_tb(" \\lt ",tb);
                                                           fnodetotex_tb((FNODE)FA2(f),tb);
                                                           break;
                                                   case L_GE:
                                                           fnodetotex_tb((FNODE)FA1(f),tb);
                                                           write_tb(" \\geq ",tb);
                                                           fnodetotex_tb((FNODE)FA2(f),tb);
                                                           break;
                                                   case L_LE:
                                                           fnodetotex_tb((FNODE)FA1(f),tb);
                                                           write_tb(" \\leq ",tb);
                                                           fnodetotex_tb((FNODE)FA2(f),tb);
                                                           break;
                                                   case L_AND:
                                                           fnodetotex_tb((FNODE)FA1(f),tb);
                                                           write_tb(" {\\rm \\ and\\ } ",tb);
                                                           fnodetotex_tb((FNODE)FA2(f),tb);
                                                           break;
                                                   case L_OR:
                                                           fnodetotex_tb((FNODE)FA1(f),tb);
                                                           write_tb(" {\\rm \\ or\\ } ",tb);
                                                           fnodetotex_tb((FNODE)FA2(f),tb);
                                                           break;
                                                   case L_NOT:
                                                           /* XXX : L_NOT is a unary operator */
                                                           write_tb("\\neg (",tb);
                                                           fnodetotex_tb((FNODE)FA1(f),tb);
                                                           write_tb(")",tb);
                                                           return;
                                           }
                                           break;
   
                                   case I_AND:
                                           fnodetotex_tb((FNODE)FA0(f),tb);
                                           write_tb(" {\\rm \\ and\\ } ",tb);
                                           fnodetotex_tb((FNODE)FA1(f),tb);
                                           break;
   
                                   case I_OR:
                                           fnodetotex_tb((FNODE)FA0(f),tb);
                                           write_tb(" {\\rm \\ or\\ } ",tb);
                                           fnodetotex_tb((FNODE)FA1(f),tb);
                                           break;
                           }
                           break;
   
                   /* ternary operators */
                   case I_CE:
                           error("fnodetotex_tb : not implemented yet");
                           break;
   
                   /* lists */
                   case I_LIST:
                           write_tb(" [ ",tb);
                           n = (NODE)FA0(f);
                           fnodenodetotex_tb(n,tb);
                           write_tb("]",tb);
                           break;
   
                   /* function */
                   case I_FUNC: case I_CAR: case I_CDR: case I_EV:
                           switch ( f->id ) {
                                   case I_FUNC:
                                           opname = symbol_name(((FUNC)FA0(f))->name);
                                           write_tb(opname,tb);
                                           write_tb("(",tb);
                                           fargstotex_tb(opname,FA1(f),tb);
                                           write_tb(")",tb);
                                           break;
                                   case I_CAR:
                                           opname = symbol_name("car");
                                           write_tb(opname,tb);
                                           write_tb("(",tb);
                                           fargstotex_tb(opname,FA0(f),tb);
                                           write_tb(")",tb);
                                           break;
                                   case I_CDR:
                                           opname = symbol_name("cdr");
                                           write_tb(opname,tb);
                                           write_tb("(",tb);
                                           fargstotex_tb(opname,FA0(f),tb);
                                           write_tb(")",tb);
                                           break;
                                   case I_EV:
                                           n = (NODE)FA0(f);
                                           for ( t0 = 0, i = 0; n; n = NEXT(n), i++ ) {
                                                   sprintf(vname,"x_{%d}^{",i);
                                                   write_tb(vname,tb);
                                                   fnodetotex_tb((FNODE)BDY(n),tb);
                                                   write_tb("} ",tb);
                                           }
                                           break;
                           }
                           break;
   
                   case I_STR:
                           write_tb((char *)FA0(f),tb);
                           break;
   
                   case I_FORMULA:
                           obj = (Obj)FA0(f);
                           if ( obj && OID(obj) == O_P ) {
                                   opname = symbol_name(VR((P)obj)->name);
                           } else {
                                   len = estimate_length(CO,obj);
                                   opname = (char *)MALLOC_ATOMIC(len+1);
                                   soutput_init(opname);
                                   sprintexpr(CO,obj);
                           }
                           write_tb(opname,tb);
                           break;
   
                   case I_PVAR:
                           if ( FA1(f) )
                                   error("fnodetotex_tb : not implemented yet");
                           GETPVNAME(FA0(f),opname);
                           write_tb(opname,tb);
                           break;
   
                   default:
                           error("fnodetotex_tb : not implemented yet");
           }
   }
   
   void fnodenodetotex_tb(NODE n,TB tb)
   {
           for ( ; n; n = NEXT(n) ) {
                   fnodetotex_tb((FNODE)BDY(n),tb);
                   if ( NEXT(n) ) write_tb(", ",tb);
           }
   }
   
   void fargstotex_tb(char *name,FNODE f,TB tb)
   {
           NODE n;
   
           if ( !strcmp(name,"matrix") ) {
                   error("fargstotex_tb : not implemented yet");
           } else if ( !strcmp(name,"vector") ) {
                   error("fargstotex_tb : not implemented yet");
           } else {
                   if ( f->id == I_LIST ) {
                           n = (NODE)FA0(f);
                           fnodenodetotex_tb(n,tb);
                   } else
                           fnodetotex_tb(f,tb);
           }
 }  }

Legend:
Removed from v.1.4  
changed lines
  Added in v.1.14

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>