[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.10 and 1.16

version 1.10, 2003/02/14 22:29:07 version 1.16, 2004/03/04 03:31:28
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/builtin/strobj.c,v 1.9 2002/12/09 00:42:13 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/builtin/strobj.c,v 1.15 2004/03/04 01:41:32 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "parse.h"  #include "parse.h"
 #include "ctype.h"  #include "ctype.h"
 #if defined(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>  #include <string.h>
   
Line 61  extern char *parse_strp;
Line 63  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 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 Ptrim_tb();
   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);
   void dp_trim_tb(TB tb);
   
 struct ftab str_tab[] = {  struct ftab str_tab[] = {
         {"rtostr",Prtostr,1},          {"rtostr",Prtostr,1},
Line 71  struct ftab str_tab[] = {
Line 86  struct ftab str_tab[] = {
         {"str_len",Pstr_len,1},          {"str_len",Pstr_len,1},
         {"str_chr",Pstr_chr,3},          {"str_chr",Pstr_chr,3},
         {"sub_str",Psub_str,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},
           {"trim_tb",Ptrim_tb,2},
           {"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 dp_trim_tb(TB tb)
   {
           int i,j,l,onei,onej;
           char *p;
   
           if ( !tb->next ) return;
           /* number    1 =>  number */
           onei = -1; onej = -1;
           for ( i = tb->next-1; i >= 0; i-- ) {
                   p = tb->body[i];
                   l = strlen(p);
                   if ( onei < 0 ) {
                           for ( j = l-1; j >= 0 && isspace(p[j]); j-- );
                           if ( j < 0 ) continue;
                           if ( p[j] != '1' ) return;
                           /* we found 1 */
                           onei = i; onej = j; j--;
                   } else
                           j = l-1;
                   /* we search the previous non-space character */
                   for ( ; j >= 0 && isspace(p[j]); j-- );
                   if ( j < 0 ) continue;
                   if ( p[j] == '+' || p[j] == '-' ) return;
                   else break;
           }
           if ( i < 0 ) return;
           /* the previous character is a number */
           l = strlen(tb->body[onei]);
           p = (char *)MALLOC_ATOMIC(l+1);
           strcpy(p,tb->body[onei]);
           p[onej] = ' ';
           tb->body[onei] = p;
   }
   
   void Ptrim_tb(NODE arg,Q *rp)
   {
           asir_assert(ARG0(arg),O_TB,"trim_tb");
           if ( !INT(ARG1(arg)) ) error("trim_tb : invalid argument");
           switch ( QTOS((Q)ARG1(arg)) ) {
                   case O_DP:
                           dp_trim_tb((TB)ARG0(arg));
                           break;
                   default:
                           break;
           }
           *rp = 0;
   }
   
   void Pwrite_to_tb(NODE arg,Q *rp)
   {
           int i;
           Obj obj;
           TB tb;
   
           asir_assert(ARG1(arg),O_TB,"write_to_tb");
           obj = ARG0(arg);
           if ( !obj )
                   write_tb("",ARG1(arg));
           else if ( OID(obj) == O_STR )
                   write_tb(BDY((STRING)obj),ARG1(arg));
           else if ( OID(obj) == O_TB ) {
                   tb = (TB)obj;
                   for ( i = 0; i < tb->next; i++ )
                           write_tb(tb->body[i],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)  void Pstr_len(arg,rp)
 NODE arg;  NODE arg;
 Q *rp;  Q *rp;
 {  {
         STRING str;          Obj obj;
         int r;          TB tb;
           int r,i;
   
         str = (STRING)ARG0(arg);          obj = (Obj)ARG0(arg);
         asir_assert(str,O_STR,"str_chr");          if ( !obj || (OID(obj) != O_STR && OID(obj) != O_TB) )
         r = strlen(BDY(str));                  error("str_len : invalid argument");
           if ( OID(obj) == O_STR)
                   r = strlen(BDY((STRING)obj));
           else if ( OID(obj) == O_TB ) {
                   tb = (TB)obj;
                   for ( r = i = 0; i < tb->next; i++ )
                           r += strlen(tb->body[i]);
           }
         STOQ(r,*rp);          STOQ(r,*rp);
 }  }
   
Line 207  Obj *rp;
Line 381  Obj *rp;
         void recover(int);          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_create_var(0,cmd,&fnode);          exprparse_create_var(0,cmd,&fnode);
Line 225  STRING *rp;
Line 401  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 251  P *rp;
Line 427  P *rp;
 #else  #else
         makevar(p,rp);          makevar(p,rp);
 #endif  #endif
   }
   
   struct TeXSymbol {
           char *text;
           char *symbol;
   };
   
   static struct TeXSymbol texsymbol[] = {
    {"sin","\\sin"},
    {"cos","\\cos"},
    {"tan","\\tan"},
    {"sinh","\\sinh"},
    {"cosh","\\cosh"},
    {"tanh","\\tanh"},
    {"exp","\\exp"},
    {"log","\\log"},
   
   /* Greek Letters (lower case) */
    {"alpha","\\alpha"},
    {"beta","\\beta"},
    {"gamma","\\gamma"},
    {"delta","\\delta"},
    {"epsilon","\\epsilon"},
    {"varepsilon","\\varepsilon"},
    {"zeta","\\zeta"},
    {"eta","\\eta"},
    {"theta","\\theta"},
    {"vartheta","\\vartheta"},
    {"iota","\\iota"},
    {"kappa","\\kappa"},
    {"lambda","\\lambda"},
    {"mu","\\mu"},
    {"nu","\\nu"},
    {"xi","\\xi"},
    {"pi","\\pi"},
    {"varpi","\\varpi"},
    {"rho","\\rho"},
    {"sigma","\\sigma"},
    {"varsigma","\\varsigma"},
    {"tau","\\tau"},
    {"upsilon","\\upsilon"},
    {"phi","\\phi"},
    {"varphi","\\varphi"},
    {"chi","\\chi"},
    {"omega","\\omega"},
   
   /* Greek Letters, (upper case) */
    {"ggamma","\\Gamma"},
    {"ddelta","\\Delta"},
    {"ttheta","\\Theta"},
    {"llambda","\\Lambda"},
    {"xxi","\\Xi"},
    {"ppi","\\Pi"},
    {"ssigma","\\Sigma"},
    {"uupsilon","\\Upsilon"},
    {"pphi","\\Phi"},
    {"ppsi","\\Psi"},
    {"oomega","\\Omega"},
   
    /* Our own mathematical functions */
    {"algebra_tensor","\\otimes"},
    {"base_where","{\\rm \\ where \\ }"},
    /* Mathematical constants */
    {"c_pi","\\pi"},
    {"c_i","\\sqrt{-1}"},
   
    /* Temporary  */
    {"dx","\\partial"},
    {0,0}
   };
   
   char *symbol_name(char *name)
   {
           int i;
   
           for ( i = 0; texsymbol[i].text; i++ )
                   if ( !strcmp(texsymbol[i].text,name) )
                           return texsymbol[i].symbol;
           return name;
   }
   
   void fnodetotex_tb(FNODE f,TB tb)
   {
           NODE n,t,t0;
           char vname[BUFSIZ];
           char *opname;
           Obj obj;
           int i,len,allzero;
           FNODE fi;
   
           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);
                                           allzero = 1;
                                           for ( t0 = 0, i = 0; n; n = NEXT(n), i++ ) {
                                                   fi = (FNODE)BDY(n);
                                                   if ( fi->id == I_FORMULA && !FA0(fi) ) continue;
                                                   allzero = 0;
                                                   if ( fi->id == I_FORMULA && UNIQ(FA0(fi)) ) {
                                                           sprintf(vname,"x_{%d}",i);
                                                           len = strlen(vname);
                                                           opname = MALLOC_ATOMIC(len+1);
                                                           strcpy(opname,vname);
                                                           write_tb(opname,tb);
                                                   } else {
                                                           sprintf(vname,"x_{%d}^{",i);
                                                           len = strlen(vname);
                                                           opname = MALLOC_ATOMIC(len+1);
                                                           strcpy(opname,vname);
                                                           write_tb(opname,tb);
                                                           fnodetotex_tb((FNODE)BDY(n),tb);
                                                           write_tb("} ",tb);
                                                   }
                                           }
                                           /* XXX */
                                           if ( allzero )
                                                   write_tb(" 1 ",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.10  
changed lines
  Added in v.1.16
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>