| version 1.1, 2018/09/19 05:45:07 |
version 1.2, 2018/09/28 08:20:28 |
|
|
| /* |
/* |
| * $OpenXM$ |
* $OpenXM: OpenXM_contrib2/asir2018/engine/d-itv.c,v 1.1 2018/09/19 05:45:07 noro Exp $ |
| */ |
*/ |
| #if defined(INTERVAL) |
#if defined(INTERVAL) |
| #include <float.h> |
#include <float.h> |
| Line 498 void pwritvd(Num a, Num e, IntervalDouble *c) |
|
| Line 498 void pwritvd(Num a, Num e, IntervalDouble *c) |
|
| error("pwritvd : can't calculate a fractional power"); |
error("pwritvd : can't calculate a fractional power"); |
| #endif |
#endif |
| } else { |
} else { |
| ei = QTOS((Q)e); |
ei = ZTOS((Q)e); |
| pwritv0d((IntervalDouble)a,ei,&t); |
pwritv0d((IntervalDouble)a,ei,&t); |
| if ( SGN((Q)e) < 0 ) |
if ( SGN((Q)e) < 0 ) |
| divnum(0,(Num)ONE,(Num)t,(Num *)c); |
divnum(0,(Num)ONE,(Num)t,(Num *)c); |