version 1.15, 2006/08/09 02:40:47 |
version 1.17, 2015/08/04 06:20:45 |
|
|
* OF THE SOFTWARE HAS BEEN DEVELOPED BY A THIRD PARTY, THE THIRD PARTY |
* OF THE SOFTWARE HAS BEEN DEVELOPED BY A THIRD PARTY, THE THIRD PARTY |
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* $OpenXM: OpenXM_contrib2/asir2000/io/bload.c,v 1.14 2004/12/10 07:36:35 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/io/bload.c,v 1.16 2009/03/16 16:43:03 ohara Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "parse.h" |
#include "parse.h" |
#include "com.h" |
#include "com.h" |
#if defined(PARI) |
|
#include "genpari.h" |
|
int get_lg(GEN); |
|
#endif |
|
|
|
extern VL file_vl; |
extern VL file_vl; |
|
|
Line 124 void loadreal(FILE *s,Real *p) |
|
Line 120 void loadreal(FILE *s,Real *p) |
|
|
|
void loadbf(FILE *s,BF *p) |
void loadbf(FILE *s,BF *p) |
{ |
{ |
#if defined(PARI) |
BF r; |
GEN z; |
char dmy; |
unsigned int uexpo,lexpo; |
int sgn,prec; |
char dmy; |
UL exp; |
int sign; |
|
unsigned int len; |
|
unsigned long expo; |
|
BF q; |
|
|
|
|
int len; |
read_char(s,&dmy); |
read_char(s,&dmy); |
read_int(s,&sign); |
NEWBF(r); |
read_int(s,&uexpo); |
read_int(s,&sgn); |
read_int(s,&lexpo); |
read_int(s,&prec); |
|
read_int64(s,&exp); |
#if defined(LONG_IS_32BIT) |
|
if ( uexpo ) |
|
error("loadbf : exponent too large"); |
|
read_int(s,&len); |
read_int(s,&len); |
NEWBF(q,len+2); |
mpfr_init2(r->body,prec); |
z = (GEN)BDY(q); |
MPFR_SIGN(r->body) = sgn; |
settyp(z,t_REAL); |
MPFR_EXP(r->body) = (int)exp; |
setlg(z,len+2); |
#if SIZEOF_LONG == 4 |
setsigne(z,(long)sign); |
read_intarray(s,(int *)r->body->_mpfr_d,len); |
setexpo(z,(long)lexpo); |
#else /* SIZEOF_LONG == 8 */ |
read_intarray(s,(int *)(z+2),len); |
read_longarray(s,(long *)r->body->_mpfr_d,len); |
#elif defined(LONG_IS_64BIT) |
|
expo = (((UL)uexpo)<<32)|((UL)lexpo); |
|
read_int(s,&len); |
|
NEWBF(q,(len+5)/2); /* 2+(len+1)/2 */ |
|
z = (GEN)BDY(q); |
|
settyp(z,t_REAL); |
|
setlg(z,(len+5)/2); |
|
setsigne(z,(long)sign); |
|
setexpo(z,(long)expo); |
|
read_longarray(s,z+2,len); |
|
#endif |
#endif |
*p = q; |
*p = r; |
#else |
|
error("loadbf : PARI is not combined"); |
|
#endif |
|
} |
} |
|
|
#if defined(INTERVAL) |
#if defined(INTERVAL) |