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

Diff for /OpenXM_contrib2/asir2000/builtin/parif.c between version 1.6 and 1.7

version 1.6, 2000/12/05 01:24:51 version 1.7, 2000/12/22 09:58:32
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/parif.c,v 1.5 2000/11/15 00:07:20 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/builtin/parif.c,v 1.6 2000/12/05 01:24:51 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "parse.h"  #include "parse.h"
Line 63  void ritopa_i(N,int,GEN *);
Line 63  void ritopa_i(N,int,GEN *);
 void Peval(),Psetprec(),p_pi(),p_e(),p_mul(),p_gcd();  void Peval(),Psetprec(),p_pi(),p_e(),p_mul(),p_gcd();
 void asir_cgiv(GEN);  void asir_cgiv(GEN);
   
   #if defined(INTERVAL) || 1
   void Psetprecword();
   #endif
   
 struct ftab pari_tab[] = {  struct ftab pari_tab[] = {
         {"eval",Peval,-2}, {"setprec",Psetprec,-1}, {0,0,0},          {"eval",Peval,-2},
           {"setprec",Psetprec,-1},
   #if defined(INTERVAL) || 1
           {"setprecword",Psetprecword,-1},
   #endif
           {0,0,0},
 };  };
   
 #define MKPREC(a,i,b) (argc(a)==(i)?mkprec(QTOS((Q)(b))):prec)  #define MKPREC(a,i,b) (argc(a)==(i)?mkprec(QTOS((Q)(b))):prec)
Line 83  void f(NODE,Obj *);\
Line 92  void f(NODE,Obj *);\
 void f(ar,rp) NODE ar; Obj *rp;\  void f(ar,rp) NODE ar; Obj *rp;\
 { GEN _pt1_,_pt2_,_pt3_; CALLPARI2(pf,ARG0(ar),ARG1(ar),MKPREC(ar,3,ARG2(ar)),rp); }  { GEN _pt1_,_pt2_,_pt3_; CALLPARI2(pf,ARG0(ar),ARG1(ar),MKPREC(ar,3,ARG2(ar)),rp); }
   
   #if defined(INTERVAL)
   #define PREC_CONV       pariK1
   #else
 #if defined(LONG_IS_32BIT)  #if defined(LONG_IS_32BIT)
 #define PREC_CONV               0.103810253  #define PREC_CONV               0.103810253
 #endif  #endif
 #if defined(LONG_IS_64BIT)  #if defined(LONG_IS_64BIT)
 #define PREC_CONV               0.051905126  #define PREC_CONV               0.051905126
 #endif  #endif
   #endif
   
 /* XXX : we should be more careful when we free PARI pointers. */  /* XXX : we should be more careful when we free PARI pointers. */
   
Line 123  Obj *rp;
Line 136  Obj *rp;
         int p;          int p;
         Q q;          Q q;
   
   #if defined(INTERVAL) || 1
           p = (int)((prec-2)/PREC_CONV); STOQ(p,q); *rp = (Obj)q;
           if ( arg ) {
                   asir_assert(ARG0(arg),O_N,"setprec");
                   p = QTOS((Q)ARG0(arg));
                   if ( p > 0 )
                           prec = (long)(p*PREC_CONV+3);
           }
   #else
         p = (int)((prec-3)/PREC_CONV); STOQ(p,q); *rp = (Obj)q;          p = (int)((prec-3)/PREC_CONV); STOQ(p,q); *rp = (Obj)q;
         if ( arg ) {          if ( arg ) {
                 asir_assert(ARG0(arg),O_N,"setprec");                  asir_assert(ARG0(arg),O_N,"setprec");
                 prec = mkprec(QTOS((Q)ARG0(arg)));                  prec = mkprec(QTOS((Q)ARG0(arg)));
         }          }
   #endif
 }  }
   
   #if defined(INTERVAL) || 1
   void Psetprecword(arg,rp)
   NODE arg;
   Obj *rp;
   {
           int p;
           Q q;
   
           p = (int)((prec-2)); STOQ(p,q); *rp = (Obj)q;
           if ( arg ) {
                   asir_assert(ARG0(arg),O_N,"setprecword");
                   p = QTOS((Q)ARG0(arg));
                   if ( p > 0 ) {
                           prec = p + 2;
                   }
           }
   }
   #endif
   
 void p_pi(arg,rp)  void p_pi(arg,rp)
 NODE arg;  NODE arg;
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>