===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/io/bload.c,v
retrieving revision 1.12
retrieving revision 1.17
diff -u -p -r1.12 -r1.17
--- OpenXM_contrib2/asir2000/io/bload.c	2003/02/14 22:29:15	1.12
+++ OpenXM_contrib2/asir2000/io/bload.c	2015/08/04 06:20:45	1.17
@@ -44,26 +44,26 @@
  * OF THE SOFTWARE HAS BEEN DEVELOPED BY A THIRD PARTY, THE THIRD PARTY
  * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
  * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
- * $OpenXM: OpenXM_contrib2/asir2000/io/bload.c,v 1.11 2002/01/08 04:14:38 kondoh Exp $ 
+ * $OpenXM: OpenXM_contrib2/asir2000/io/bload.c,v 1.16 2009/03/16 16:43:03 ohara Exp $ 
 */
 #include "ca.h"
 #include "parse.h"
 #include "com.h"
-#if defined(PARI)
-#include "genpari.h"
-int get_lg(GEN);
-#endif
 
 extern VL file_vl;
 
+void loadnbp(FILE *s,NBP *p);
+
 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)
 void loaditv();
 void loaditvd();
-void (*nloadf[])() = { loadq, loadreal, 0, loadbf, loaditv, loaditvd, 0, loaditv, loadcplx, loadmi, loadlm, loadgf2n, loadgfpn, loadgfs, loadgfsn };
+void (*nloadf[])() = { loadq, loadreal, 0, loadbf, loaditv, loaditvd, 0, loaditv, loadcplx, loadmi, loadlm, loadgf2n, loadgfpn, loadgfs, loadgfsn, loaddalg };
 #else
-void (*nloadf[])() = { loadq, loadreal, 0, loadbf, loadcplx, loadmi, loadlm, loadgf2n, loadgfpn, loadgfs, loadgfsn };
+void (*nloadf[])() = { loadq, loadreal, 0, loadbf, loadcplx, loadmi, loadlm, loadgf2n, loadgfpn, loadgfs, loadgfsn, loaddalg };
 #endif
 
 void loadobj(FILE *s,Obj *p)
@@ -120,45 +120,27 @@ void loadreal(FILE *s,Real *p)
 
 void loadbf(FILE *s,BF *p)
 {
-#if defined(PARI)
-	GEN z;
-	unsigned int uexpo,lexpo;
-	char dmy;
-	int sign;
-	unsigned int len;
-	BF q;
+	BF r;
+  char dmy;
+	int sgn,prec;
+  UL exp;
 
+	int len;
 	read_char(s,&dmy);
-	read_int(s,&sign);
-	read_int(s,&uexpo);
-	read_int(s,&lexpo);
-
-#if defined(LONG_IS_32BIT)
-	if ( uexpo )
-		error("loadbf : exponent too large");
+	NEWBF(r);
+	read_int(s,&sgn);
+	read_int(s,&prec);
+	read_int64(s,&exp);
 	read_int(s,&len);
-	NEWBF(q,len+2);
-	z = (GEN)BDY(q);
-	settyp(z,t_REAL);
-	setlg(z,len+2);
-	setsigne(z,(long)sign);
-	setexpo(z,(long)lexpo);
-	read_intarray(s,(int *)(z+2),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);
+  mpfr_init2(r->body,prec);
+  MPFR_SIGN(r->body) = sgn;
+	MPFR_EXP(r->body) = (int)exp;
+#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
-	*p = q;
-#else
-	error("loadbf : PARI is not combined");
-#endif
+	*p = r;
 }
 
 #if defined(INTERVAL)
@@ -265,6 +247,17 @@ void loadgfsn(FILE *s,GFSN *p)
 	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;
@@ -452,4 +445,22 @@ void loadgfmmat(FILE *s,GFMMAT *p) 
 	for ( i = 0; i < row; i++ )
 		read_intarray(s,a[i],col);
 	*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);
 }