version 1.6, 2000/12/22 10:03:30 |
version 1.18, 2017/08/31 09:11:04 |
|
|
* 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.5 2000/12/05 01:24:54 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/io/bload.c,v 1.17 2015/08/04 06:20:45 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "parse.h" |
#include "parse.h" |
#include "com.h" |
#include "com.h" |
#if PARI |
|
#include "genpari.h" |
|
int get_lg(GEN); |
|
#endif |
|
|
|
void loaderror(FILE *,ERR *); |
|
void loadui(FILE *,USINT *); |
|
void loaddp(FILE *,DP *); |
|
void loadstr(FILE *,char **); |
|
void loadstring(FILE *,STRING *); |
|
void loadmat(FILE *,MAT *); |
|
void loadvect(FILE *,VECT *); |
|
void loadlist(FILE *,LIST *); |
|
void loadr(FILE *,R *); |
|
void loadp(FILE *,P *); |
|
void loadgf2n(FILE *,GF2N *); |
|
void loadgfpn(FILE *,GFPN *); |
|
void loadlm(FILE *,LM *); |
|
void loadmi(FILE *,MQ *); |
|
void loadcplx(FILE *,C *); |
|
void loadbf(FILE *,BF *); |
|
void loadreal(FILE *,Real *); |
|
void loadq(FILE *,Q *); |
|
void loadnum(FILE *,Num *); |
|
void loadgfmmat(FILE *,GFMMAT *); |
|
void loadbytearray(FILE *,BYTEARRAY *); |
|
|
|
V loadpfins(FILE *); |
|
|
|
extern VL file_vl; |
extern VL file_vl; |
|
|
|
void loadnbp(FILE *s,NBP *p); |
|
|
void (*loadf[])() = { 0, loadnum, loadp, loadr, loadlist, loadvect, loadmat, |
void (*loadf[])() = { 0, loadnum, loadp, loadr, loadlist, loadvect, loadmat, |
loadstring, 0, loaddp, loadui, loaderror,0,0,0,loadgfmmat, loadbytearray }; |
loadstring, 0, loaddp, loadui, loaderror,0,0,0,loadgfmmat, |
|
loadbytearray, 0, 0, 0, 0, 0, 0, 0, 0, loadnbp }; |
|
|
#if defined(INTERVAL) |
#if defined(INTERVAL) |
int loaditv(); |
void loaditv(); |
int loaditvd(); |
void loaditvd(); |
void (*nloadf[])() = { loadq, loadreal, 0, loadbf, loaditv, loaditvd, 0, loaditv, loadcplx, loadmi, loadlm, loadgf2n, loadgfpn }; |
void (*nloadf[])() = { loadq, loadreal, 0, loadbf, loaditv, loaditvd, 0, loaditv, loadcplx, loadmi, loadlm, loadgf2n, loadgfpn, loadgfs, loadgfsn, loaddalg }; |
#else |
#else |
void (*nloadf[])() = { loadq, loadreal, 0, loadbf, loadcplx, loadmi, loadlm, loadgf2n, loadgfpn }; |
void (*nloadf[])() = { loadq, loadreal, 0, loadbf, loadcplx, loadmi, loadlm, loadgf2n, loadgfpn, loadgfs, loadgfsn, loaddalg }; |
#endif |
#endif |
|
|
void loadobj(s,p) |
void loadobj(FILE *s,Obj *p) |
FILE *s; |
|
Obj *p; |
|
{ |
{ |
short id; |
short id; |
|
|
|
|
(*loadf[id])(s,p); |
(*loadf[id])(s,p); |
} |
} |
|
|
void loadnum(s,p) |
void loadnum(FILE *s,Num *p) |
FILE *s; |
|
Num *p; |
|
{ |
{ |
char nid; |
char nid; |
|
|
|
|
(*nloadf[nid])(s,p); |
(*nloadf[nid])(s,p); |
} |
} |
|
|
void loadq(s,p) |
void loadq(FILE *s,Q *p) |
FILE *s; |
|
Q *p; |
|
{ |
{ |
int size[2]; |
int size[2]; |
char sgn; |
char sgn; |
|
|
NDTOQ(nm,dn,sgn,*p); |
NDTOQ(nm,dn,sgn,*p); |
} |
} |
|
|
void loadreal(s,p) |
void loadreal(FILE *s,Real *p) |
FILE *s; |
|
Real *p; |
|
{ |
{ |
Real q; |
Real q; |
char dmy; |
char dmy; |
|
|
*p = q; |
*p = q; |
} |
} |
|
|
void loadbf(s,p) |
void loadbf(FILE *s,BF *p) |
FILE *s; |
|
BF *p; |
|
{ |
{ |
#if PARI |
BF r; |
GEN z; |
char dmy; |
unsigned int uexpo,lexpo; |
int sgn,prec; |
UL expo; |
UL exp; |
char dmy; |
|
int sign; |
|
unsigned int len; |
|
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 defined(VISUAL) |
setsigne(z,(long)sign); |
#if !defined(_WIN64) |
setexpo(z,(long)lexpo); |
read_intarray(s,(int *)r->body->_mpfr_d,len); |
read_intarray(s,(int *)(z+2),len); |
#else |
#elif defined(LONG_IS_64BIT) |
read_longarray(s,(long long*)r->body->_mpfr_d,len); |
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; |
|
#else |
#else |
error("loadbf : PARI is not combined"); |
#if SIZEOF_LONG == 4 |
|
read_intarray(s,(int *)r->body->_mpfr_d,len); |
|
#else /* SIZEOF_LONG == 8 */ |
|
read_longarray(s,(long *)r->body->_mpfr_d,len); |
#endif |
#endif |
|
#endif |
|
*p = r; |
} |
} |
|
|
#if defined(INTERVAL) |
#if defined(INTERVAL) |
loaditv(s,p) |
void loaditv(FILE *s,Itv *p) |
FILE *s; |
|
Itv *p; |
|
{ |
{ |
Itv q; |
Itv q; |
char dmy; |
char dmy; |
|
|
*p = q; |
*p = q; |
} |
} |
|
|
loaditvd(s,p) |
void loaditvd(FILE *s,IntervalDouble *p) |
FILE *s; |
|
ItvD *p; |
|
{ |
{ |
ItvD q; |
IntervalDouble q; |
char dmy; |
char dmy; |
|
|
read_char(s,&dmy); |
read_char(s,&dmy); |
NEWItvD(q); |
NEWIntervalDouble(q); |
read_double(s,&INF(q)); |
read_double(s,&INF(q)); |
read_double(s,&SUP(q)); |
read_double(s,&SUP(q)); |
*p = q; |
*p = q; |
} |
} |
#endif |
#endif |
|
|
void loadcplx(s,p) |
void loadcplx(FILE *s,C *p) |
FILE *s; |
|
C *p; |
|
{ |
{ |
C q; |
C q; |
char dmy; |
char dmy; |
|
|
*p = q; |
*p = q; |
} |
} |
|
|
void loadmi(s,p) |
void loadmi(FILE *s,MQ *p) |
FILE *s; |
|
MQ *p; |
|
{ |
{ |
MQ q; |
MQ q; |
char dmy; |
char dmy; |
|
|
*p = q; |
*p = q; |
} |
} |
|
|
void loadlm(s,p) |
void loadlm(FILE *s,LM *p) |
FILE *s; |
|
LM *p; |
|
{ |
{ |
int size; |
int size; |
char dmy; |
char dmy; |
|
|
MKLM(body,*p); |
MKLM(body,*p); |
} |
} |
|
|
void loadgf2n(s,p) |
void loadgf2n(FILE *s,GF2N *p) |
FILE *s; |
|
GF2N *p; |
|
{ |
{ |
char dmy; |
char dmy; |
int len; |
int len; |
|
|
MKGF2N(body,*p); |
MKGF2N(body,*p); |
} |
} |
|
|
void loadgfpn(s,p) |
void loadgfpn(FILE *s,GFPN *p) |
FILE *s; |
|
GFPN *p; |
|
{ |
{ |
char dmy; |
char dmy; |
int d,i; |
int d,i; |
|
|
MKGFPN(body,*p); |
MKGFPN(body,*p); |
} |
} |
|
|
void loadp(s,p) |
void loadgfs(FILE *s,GFS *p) |
FILE *s; |
|
P *p; |
|
{ |
{ |
|
GFS q; |
|
char dmy; |
|
|
|
read_char(s,&dmy); |
|
NEWGFS(q); read_int(s,(int *)&CONT(q)); |
|
*p = q; |
|
} |
|
|
|
void loadgfsn(FILE *s,GFSN *p) |
|
{ |
|
char dmy; |
|
int d; |
|
UM body; |
|
|
|
read_char(s,&dmy); read_int(s,&d); |
|
body = UMALLOC(d); DEG(body) = d; |
|
read_intarray(s,COEF(body),d+1); |
|
MKGFSN(body,*p); |
|
} |
|
|
|
void loaddalg(FILE *s,DAlg *p) |
|
{ |
|
char dmy; |
|
Obj nm,dn; |
|
|
|
read_char(s,&dmy); |
|
loadobj(s,&nm); |
|
loadobj(s,&dn); |
|
MKDAlg((DP)nm,(Q)dn,*p); |
|
} |
|
|
|
void loadp(FILE *s,P *p) |
|
{ |
V v; |
V v; |
int n,vindex; |
int n,vindex; |
DCP dc,dc0; |
DCP dc,dc0; |
|
|
|
|
/* |name(str)|argc(int)|darray(intarray)|args| */ |
/* |name(str)|argc(int)|darray(intarray)|args| */ |
|
|
V loadpfins(s) |
V loadpfins(FILE *s) |
FILE *s; |
|
{ |
{ |
char *name; |
char *name; |
FUNC fp; |
FUNC fp; |
|
|
return v; |
return v; |
} |
} |
|
|
void loadr(s,p) |
void loadr(FILE *s,R *p) |
FILE *s; |
|
R *p; |
|
{ |
{ |
R r; |
R r; |
|
|
|
|
loadobj(s,(Obj *)&NM(r)); loadobj(s,(Obj *)&DN(r)); *p = r; |
loadobj(s,(Obj *)&NM(r)); loadobj(s,(Obj *)&DN(r)); *p = r; |
} |
} |
|
|
void loadlist(s,p) |
void loadlist(FILE *s,LIST *p) |
FILE *s; |
|
LIST *p; |
|
{ |
{ |
int n; |
int n; |
NODE tn,tn0; |
NODE tn,tn0; |
|
|
MKLIST(*p,tn0); |
MKLIST(*p,tn0); |
} |
} |
|
|
void loadvect(s,p) |
void loadvect(FILE *s,VECT *p) |
FILE *s; |
|
VECT *p; |
|
{ |
{ |
int i,len; |
int i,len; |
VECT vect; |
VECT vect; |
|
|
*p = vect; |
*p = vect; |
} |
} |
|
|
void loadmat(s,p) |
void loadmat(FILE *s,MAT *p) |
FILE *s; |
|
MAT *p; |
|
{ |
{ |
int row,col,i,j; |
int row,col,i,j; |
MAT mat; |
MAT mat; |
|
|
*p = mat; |
*p = mat; |
} |
} |
|
|
void loadstring(s,p) |
void loadstring(FILE *s,STRING *p) |
FILE *s; |
|
STRING *p; |
|
{ |
{ |
char *t; |
char *t; |
|
|
loadstr(s,&t); MKSTR(*p,t); |
loadstr(s,&t); MKSTR(*p,t); |
} |
} |
|
|
void loadstr(s,p) |
void loadstr(FILE *s,char **p) |
FILE *s; |
|
char **p; |
|
{ |
{ |
int len; |
int len; |
char *t; |
char *t; |
|
|
*p = t; |
*p = t; |
} |
} |
|
|
void loadbytearray(s,p) |
void loadbytearray(FILE *s,BYTEARRAY *p) |
FILE *s; |
|
BYTEARRAY *p; |
|
{ |
{ |
int len; |
int len; |
BYTEARRAY array; |
BYTEARRAY array; |
|
|
*p = array; |
*p = array; |
} |
} |
|
|
void loaddp(s,p) |
void loaddp(FILE *s,DP *p) |
FILE *s; |
|
DP *p; |
|
{ |
{ |
int nv,n,i,sugar; |
int nv,n,i,sugar; |
DP dp; |
DP dp; |
|
|
NEXT(m) = 0; MKDP(nv,m0,dp); dp->sugar = sugar; *p = dp; |
NEXT(m) = 0; MKDP(nv,m0,dp); dp->sugar = sugar; *p = dp; |
} |
} |
|
|
void loadui(s,u) |
void loadui(FILE *s,USINT *u) |
FILE *s; |
|
USINT *u; |
|
{ |
{ |
unsigned int b; |
unsigned int b; |
|
|
read_int(s,&b); MKUSINT(*u,b); |
read_int(s,&b); MKUSINT(*u,b); |
} |
} |
|
|
void loaderror(s,e) |
void loaderror(FILE *s,ERR *e) |
FILE *s; |
|
ERR *e; |
|
{ |
{ |
Obj b; |
Obj b; |
|
|
|
|
} |
} |
|
|
|
|
void loadgfmmat(s,p) |
void loadgfmmat(FILE *s,GFMMAT *p) |
FILE *s; |
|
GFMMAT *p; |
|
{ |
{ |
int i,j,row,col; |
int i,row,col; |
unsigned int **a; |
unsigned int **a; |
GFMMAT mat; |
GFMMAT mat; |
|
|
|
|
for ( i = 0; i < row; i++ ) |
for ( i = 0; i < row; i++ ) |
read_intarray(s,a[i],col); |
read_intarray(s,a[i],col); |
*p = mat; |
*p = mat; |
|
} |
|
|
|
void loadnbp(FILE *s,NBP *p) |
|
{ |
|
int n,i; |
|
NBM m; |
|
NODE r0,r; |
|
|
|
read_int(s,&n); |
|
for ( i = 0, r0 = 0; i < n; i++ ) { |
|
NEWNBM(m); |
|
loadobj(s,(Obj *)&m->c); |
|
read_int(s,&m->d); |
|
NEWNBMBDY(m,m->d); read_intarray(s,m->b,(m->d+31)/32); |
|
NEXTNODE(r0,r); BDY(r) = (pointer)m; |
|
} |
|
if ( r0 ) NEXT(r) = 0; |
|
MKNBP(*p,r0); |
} |
} |