=================================================================== RCS file: /home/cvs/OpenXM_contrib2/asir2000/engine/d-itv.c,v retrieving revision 1.3 retrieving revision 1.5 diff -u -p -r1.3 -r1.5 --- OpenXM_contrib2/asir2000/engine/d-itv.c 2003/02/14 22:29:08 1.3 +++ OpenXM_contrib2/asir2000/engine/d-itv.c 2015/08/08 14:19:41 1.5 @@ -1,5 +1,5 @@ /* - * $OpenXM: OpenXM_contrib2/asir2000/engine/d-itv.c,v 1.2 2002/01/08 04:14:37 kondoh Exp $ + * $OpenXM: OpenXM_contrib2/asir2000/engine/d-itv.c,v 1.4 2009/03/27 14:42:29 ohara Exp $ */ #if defined(INTERVAL) #include @@ -28,6 +28,9 @@ void printbinint(int d) } } fprintf(stderr,"\n"); +#if defined(__MINGW32__) || defined(__MINGW64__) + fflush(stderr); +#endif } #endif @@ -48,6 +51,9 @@ double NatToRealUp(N a, int *expo) #if defined(ITVDEBUG) fprintf(stderr," %d : tail = %d\n", j, tail); +#if defined(__MINGW32__) || defined(__MINGW64__) + fflush(stderr); +#endif printbinint(p[j]); #endif kk = (1<< (BSH - tail)) - 1; @@ -493,12 +499,7 @@ void pwritvd(Num a, Num e, IntervalDouble *c) #endif } else if ( !INT(e) ) { #if defined(PARI) && 0 - GEN pa,pe,z; - 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); + gpui_ri((Obj)a,(Obj)c,(Obj *)c); #else error("pwritvd : can't calculate a fractional power"); #endif