version 1.3, 2003/02/14 22:29:08 |
version 1.7, 2016/06/29 08:16:11 |
|
|
/* |
/* |
* $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.6 2015/08/14 13:51:54 fujimoto Exp $ |
*/ |
*/ |
#if defined(INTERVAL) |
#if defined(INTERVAL) |
#include <float.h> |
#include <float.h> |
Line 94 double NatToRealUp(N a, int *expo) |
|
Line 94 double NatToRealUp(N a, int *expo) |
|
#ifdef vax |
#ifdef vax |
t = m[0]; m[0] = m[1]; m[1] = t; itod(m); |
t = m[0]; m[0] = m[1]; m[1] = t; itod(m); |
#endif |
#endif |
#if defined(MIPSEL) || defined(TOWNS) || defined(linux) || defined(__alpha) || defined(__FreeBSD__) || defined(__NetBSD__) |
#if defined(MIPSEL) || defined(TOWNS) || defined(linux) || defined(__alpha) || defined(__FreeBSD__) || defined(__NetBSD__) || defined(ANDROID) |
t = m[0]; m[0] = m[1]; m[1] = t; |
t = m[0]; m[0] = m[1]; m[1] = t; |
#endif |
#endif |
return *((double *)m); |
return *((double *)m); |
Line 124 static double Q2doubleDown(Q a) |
|
Line 124 static double Q2doubleDown(Q a) |
|
#ifdef vax |
#ifdef vax |
s = p[0]; p[0] = p[1]; p[1] = s; itod(p); |
s = p[0]; p[0] = p[1]; p[1] = s; itod(p); |
#endif |
#endif |
#if defined(MIPSEL) || defined(TOWNS) || defined(linux) || defined(__alpha) || defined(__FreeBSD__) || defined(__NetBSD__) |
#if defined(MIPSEL) || defined(TOWNS) || defined(linux) || defined(__alpha) || defined(__FreeBSD__) || defined(__NetBSD__) | defined(ANDROID) |
s = p[0]; p[0] = p[1]; p[1] = s; |
s = p[0]; p[0] = p[1]; p[1] = s; |
#endif |
#endif |
FPMINUSINF |
FPMINUSINF |
Line 161 static double Q2doubleUp(Q a) |
|
Line 161 static double Q2doubleUp(Q a) |
|
#ifdef vax |
#ifdef vax |
s = p[0]; p[0] = p[1]; p[1] = s; itod(p); |
s = p[0]; p[0] = p[1]; p[1] = s; itod(p); |
#endif |
#endif |
#if defined(MIPSEL) || defined(TOWNS) || defined(linux) || defined(__alpha) || defined(__FreeBSD__) || defined(__NetBSD__) |
#if defined(MIPSEL) || defined(TOWNS) || defined(linux) || defined(__alpha) || defined(__FreeBSD__) || defined(__NetBSD__) || defined(ANDROID) |
s = p[0]; p[0] = p[1]; p[1] = s; |
s = p[0]; p[0] = p[1]; p[1] = s; |
#endif |
#endif |
#if 0 |
#if 0 |
Line 493 void pwritvd(Num a, Num e, IntervalDouble *c) |
|
Line 493 void pwritvd(Num a, Num e, IntervalDouble *c) |
|
#endif |
#endif |
} else if ( !INT(e) ) { |
} else if ( !INT(e) ) { |
#if defined(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 |