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

Diff for /OpenXM_contrib2/asir2000/engine/pari.c between version 1.4 and 1.12

version 1.4, 2000/12/05 01:24:52 version 1.12, 2014/07/31 08:01: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/engine/pari.c,v 1.3 2000/08/22 05:04:06 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/engine/pari.c,v 1.11 2011/12/21 19:38:19 ohara Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #if PARI  
   #if defined(PARI)
 #include "base.h"  #include "base.h"
 #include <math.h>  #include <math.h>
 #include "genpari.h"  #include "genpari.h"
Line 58  void patori_i(GEN,N *);
Line 59  void patori_i(GEN,N *);
 void ritopa(Obj,GEN *);  void ritopa(Obj,GEN *);
 void ritopa_i(N,int,GEN *);  void ritopa_i(N,int,GEN *);
   
   // PARI_VERSION(2,2,12) == 131596
   #if PARI_VERSION_CODE >= 131596
   #define prec precreal
   #endif
 extern long prec;  extern long prec;
 extern int paristack;  extern int paristack;
   
 void risa_pari_init() {  long get_pariprec() {
         char buf[BUFSIZ];          return prec;
         int i;  }
   void set_pariprec(long p) {
           prec = p;
   }
   
   void risa_pari_init() {
         pari_init(paristack,2);          pari_init(paristack,2);
         prec = 4;          set_pariprec(4);
 }  }
   
 void create_pari_variable(index)  void create_pari_variable(index)
Line 79  int index;
Line 88  int index;
         if ( index > max_varn ) {          if ( index > max_varn ) {
                 for ( i = max_varn+1; i <= index; i++ ) {                  for ( i = max_varn+1; i <= index; i++ ) {
                         sprintf(name,"x%d",i);                          sprintf(name,"x%d",i);
   #if (PARI_VERSION_CODE < 131594)
                         fetch_named_var(name,0);                          fetch_named_var(name,0);
   #else
                           fetch_named_var(name);
   #endif
                 }                  }
                 max_varn = index;                  max_varn = index;
         }          }
Line 91  GEN a;
Line 104  GEN a;
         return lg(a);          return lg(a);
 }  }
   
   void gpui_ri(Obj a, Obj e, Obj *c)
   {
           GEN pa,pe,z;
           long ltop,lbot;
   
           ltop = avma; ritopa(a,&pa); ritopa(e,&pe); lbot = avma;
           z = gerepile(ltop,lbot,gpui(pa,pe,get_pariprec()));
           patori(z,c); cgiv(z);
   }
   
 void ritopa(a,rp)  void ritopa(a,rp)
 Obj a;  Obj a;
 GEN *rp;  GEN *rp;
 {  {
         long ltop;          long ltop;
         GEN pnm,z,w;          GEN pnm,z,w,u;
         DCP dc;          DCP dc;
         int i,j,l,row,col;          int i,j,l,row,col;
         VL vl;          VL vl;
Line 114  GEN *rp;
Line 137  GEN *rp;
                                                 *rp = pnm;                                                  *rp = pnm;
                                         else {                                          else {
                                                 *rp = z = cgetg(3,4); z[1] = (long)pnm;                                                  *rp = z = cgetg(3,4); z[1] = (long)pnm;
                                                 ritopa_i(DN((Q)a),1,(GEN *)&z[2]);                                                  ritopa_i(DN((Q)a),1,&u); z[2] = u;
                                         }                                          }
                                         break;                                          break;
                                 case N_R:                                  case N_R:
Line 123  GEN *rp;
Line 146  GEN *rp;
                                         *rp = gcopy((GEN)BDY(((BF)a))); break;                                          *rp = gcopy((GEN)BDY(((BF)a))); break;
                                 case N_C:                                  case N_C:
                                         z = cgetg(3,6);                                          z = cgetg(3,6);
                                         ritopa((Obj)((C)a)->r,(GEN *)&z[1]); ritopa((Obj)((C)a)->i,(GEN *)&z[2]);                                          ritopa((Obj)((C)a)->r,&u); z[1] = u;
                                           ritopa((Obj)((C)a)->i,&u); z[2] = u;
                                         *rp = z;                                          *rp = z;
                                         break;                                          break;
                                 default:                                  default:
Line 140  GEN *rp;
Line 164  GEN *rp;
                         setlgef(z,l+3);                          setlgef(z,l+3);
                         for ( i = l+2; i >= 2; i-- )                          for ( i = l+2; i >= 2; i-- )
                                 z[i] = (long)gzero;                                  z[i] = (long)gzero;
                         for ( dc = DC((P)a); dc; dc = NEXT(dc) )                          for ( dc = DC((P)a); dc; dc = NEXT(dc) ) {
                                 ritopa((Obj)COEF(dc),(GEN *)&z[QTOS(DEG(dc))+2]);                                  ritopa((Obj)COEF(dc),&u); z[QTOS(DEG(dc))+2] = u;
                           }
                         break;                          break;
                 case O_VECT:                  case O_VECT:
                         l = ((VECT)a)->len; z = cgetg(l+1,17);                          l = ((VECT)a)->len; z = cgetg(l+1,17);
                         for ( i = 0; i < l; i++ )                          for ( i = 0; i < l; i++ ) {
                                 ritopa((Obj)BDY((VECT)a)[i],(GEN *)&z[i+1]);                                  ritopa((Obj)BDY((VECT)a)[i],&u); z[i+1] = u;
                           }
                         *rp = z;                          *rp = z;
                         break;                          break;
                 case O_MAT:                  case O_MAT:
                         row = ((MAT)a)->row; col = ((MAT)a)->col; z = cgetg(col+1,19);                          row = ((MAT)a)->row; col = ((MAT)a)->col; z = cgetg(col+1,19);
                         for ( j = 0; j < col; j++ ) {                          for ( j = 0; j < col; j++ ) {
                                 w = cgetg(row+1,18);                                  w = cgetg(row+1,18);
                                 for ( i = 0; i < row; i++ )                                  for ( i = 0; i < row; i++ ) {
                                         ritopa((Obj)BDY((MAT)a)[i][j],(GEN *)&w[i+1]);                                          ritopa((Obj)BDY((MAT)a)[i][j],&u); w[i+1] = u;
                                   }
                                 z[j+1] = (long)w;                                  z[j+1] = (long)w;
                         }                          }
                         *rp = z;                          *rp = z;
Line 176  Obj *rp;
Line 203  Obj *rp;
         N n,nm,dn;          N n,nm,dn;
         DCP dc0,dc;          DCP dc0,dc;
         P t;          P t;
         int s,i,j,l,row,col;          int s,i,j,l,row,col,ty;
         GEN b;          GEN b;
         VL vl;          VL vl;
   
         if ( gcmp0(a) )          if ( gcmp0(a) )
                 *rp = 0;                  *rp = 0;
         else {          else {
                 switch ( typ(a) ) {                  switch ( ty = typ(a) ) {
                         case 1: /* integer */                          case 1: /* integer */
                                 patori_i(a,&n); NTOQ(n,(char)signe(a),q); *rp = (Obj)q;                                  patori_i(a,&n); NTOQ(n,(char)signe(a),q); *rp = (Obj)q;
                                 break;                                  break;
Line 191  Obj *rp;
Line 218  Obj *rp;
                                 NEWBF(r,lg(a)); bcopy((char *)a,(char *)BDY(r),lg(a)*sizeof(long));                                  NEWBF(r,lg(a)); bcopy((char *)a,(char *)BDY(r),lg(a)*sizeof(long));
                                 *rp = (Obj)r;                                  *rp = (Obj)r;
                                 break;                                  break;
                           case 3: /* integermod */
                                   MKVECT(v,2);
                                   patori((GEN)a[1],(Obj *)&BDY(v)[0]);
                                   patori((GEN)a[2],(Obj *)&BDY(v)[1]);
                                   *rp = (Obj)v;
                                   break;
                         case 4: /* rational; reduced */                          case 4: /* rational; reduced */
                                 patori_i((GEN)a[1],&nm); patori_i((GEN)a[2],&dn);                                  patori_i((GEN)a[1],&nm); patori_i((GEN)a[2],&dn);
                                 s = signe(a[1])*signe(a[2]);                                  s = signe(a[1])*signe(a[2]);
Line 249  Obj *rp;
Line 282  Obj *rp;
         }          }
 }  }
   
 #if defined(LONG_IS_32BIT)  #if SIZEOF_LONG == 4
 void ritopa_i(a,s,rp)  void ritopa_i(a,s,rp)
 N a;  N a;
 int s;  int s;
Line 284  N *rp;
Line 317  N *rp;
         for ( j = 0, b = (unsigned int *)BD(z); j < l; j++ )          for ( j = 0, b = (unsigned int *)BD(z); j < l; j++ )
                 b[l-j-1] = ((unsigned int *)g)[j+2];                  b[l-j-1] = ((unsigned int *)g)[j+2];
 }  }
 #endif  
   
 #if defined(LONG_IS_64BIT)  #elif SIZEOF_LONG == 8
 void ritopa_i(a,s,rp)  void ritopa_i(a,s,rp)
 N a;  N a;
 int s;  int s;
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.12
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>