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

Diff for /OpenXM_contrib2/asir2000/engine/d-itv.c between version 1.2 and 1.5

version 1.2, 2002/01/08 04:14:37 version 1.5, 2015/08/08 14:19:41
Line 1 
Line 1 
 /*  /*
  * $OpenXM: OpenXM_contrib2/asir2000/engine/d-itv.c,v 1.1 2000/12/22 10:03:28 saito Exp $   * $OpenXM: OpenXM_contrib2/asir2000/engine/d-itv.c,v 1.4 2009/03/27 14:42:29 ohara Exp $
 */  */
 #if defined(INTERVAL)  #if defined(INTERVAL)
 #include <float.h>  #include <float.h>
 #include "ca.h"  #include "ca.h"
 #include "base.h"  #include "base.h"
 #if PARI  #if defined(PARI)
 #include "genpari.h"  #include "genpari.h"
 #endif  #endif
   
Line 28  void printbinint(int d)
Line 28  void printbinint(int d)
                 }                  }
         }          }
         fprintf(stderr,"\n");          fprintf(stderr,"\n");
   #if defined(__MINGW32__) || defined(__MINGW64__)
           fflush(stderr);
   #endif
 }  }
 #endif  #endif
   
Line 48  double NatToRealUp(N a, int *expo)
Line 51  double NatToRealUp(N a, int *expo)
   
 #if defined(ITVDEBUG)  #if defined(ITVDEBUG)
         fprintf(stderr," %d : tail = %d\n", j, tail);          fprintf(stderr," %d : tail = %d\n", j, tail);
   #if defined(__MINGW32__) || defined(__MINGW64__)
           fflush(stderr);
   #endif
         printbinint(p[j]);          printbinint(p[j]);
 #endif  #endif
                 kk = (1<< (BSH - tail)) - 1;                  kk = (1<< (BSH - tail)) - 1;
Line 492  void pwritvd(Num a, Num e, IntervalDouble *c)
Line 498  void pwritvd(Num a, Num e, IntervalDouble *c)
                 pwrnum(0,a,e,c);                  pwrnum(0,a,e,c);
 #endif  #endif
         } else if ( !INT(e) ) {          } else if ( !INT(e) ) {
 #if PARI && 0  #if defined(PARI) && 0
                 GEN pa,pe,z;                  gpui_ri((Obj)a,(Obj)c,(Obj *)c);
                 int ltop,lbot;  
   
                 ltop = avma; ritopa(a,&pa); ritopa(e,&pe); lbot = avma;  
                 z = gerepile(ltop,lbot,gpui(pa,pe,prec));  
                 patori(z,c); cgiv(z);  
 #else  #else
                 error("pwritvd : can't calculate a fractional power");                  error("pwritvd : can't calculate a fractional power");
 #endif  #endif
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.5
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>