OpenXM_contrib2/asir2000/engine/t-itv.c - diff

Return to t-itv.c CVS log

Up to [local] / OpenXM_contrib2 / asir2000 / engine

version 1.2, 2002/01/08 04:14:37

version 1.5, 2009/03/27 14:42:29

Line 1

* $OpenXM: OpenXM_contrib2/asir2000/engine/t-itv.c,v 1.1 2000/12/22 10:03:28 saito Exp $

* $OpenXM: OpenXM_contrib2/asir2000/engine/t-itv.c,v 1.4 2003/10/20 07:18:42 saito Exp $

#if defined(INTERVAL)

#include "ca.h"

#include "base.h"

#if PARI

#if defined(PARI)

#include "genpari.h"

#endif

Line 227 void pwritvp(Itv a, Num e, Itv *c)

else if ( NID(a) <= N_B )

pwrnum(0,a,e,c);

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

error("pwritv : can't calculate a fractional power");

#endif

Line 330 void miditvp(Itv a, Num *b)

Line 325 void miditvp(Itv a, Num *b)

{

Num ai,as;

Num t;

Q TWO;

if ( ! a ) *b = 0;

else if ( (NID(a) <= N_B) )