===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/engine/dalg.c,v
retrieving revision 1.11
retrieving revision 1.13
diff -u -p -r1.11 -r1.13
--- OpenXM_contrib2/asir2000/engine/dalg.c	2005/08/24 06:28:39	1.11
+++ OpenXM_contrib2/asir2000/engine/dalg.c	2006/01/05 00:21:20	1.13
@@ -1,5 +1,5 @@
 /*
- * $OpenXM: OpenXM_contrib2/asir2000/engine/dalg.c,v 1.10 2005/08/02 07:16:42 noro Exp $
+ * $OpenXM: OpenXM_contrib2/asir2000/engine/dalg.c,v 1.12 2005/10/12 14:43:36 noro Exp $
 */
 
 #include "ca.h"
@@ -651,6 +651,7 @@ NODE inv_or_split_dalg(DAlg a,DAlg *c)
 	struct order_spec *current_spec;
 	struct oEGT eg0,eg1;
 	extern struct oEGT eg_le;
+	extern int DP_Print;
 
 	if ( !(nf=current_numberfield) )
 		error("invdalg : current_numberfield is not set");
@@ -669,6 +670,7 @@ NODE inv_or_split_dalg(DAlg a,DAlg *c)
 	simp = (DAlg *)ALLOCA(dim*sizeof(DAlg));
 	current_spec = dp_current_spec; initd(nf->spec);
 	for ( i = 0; i < dim; i++ ) {
+		if ( DP_Print ) { fprintf(asir_out,"."); fflush(asir_out); }
 		m = mb[i];
 		for ( j = i-1; j >= 0; j-- )
 			if ( dp_redble(m,mb[j]) )
@@ -706,7 +708,7 @@ NODE inv_or_split_dalg(DAlg a,DAlg *c)
 		}
 	}
 	get_eg(&eg0);
-	rank = generic_gauss_elim_hensel(mobj,&sol,&dnsol,&rinfo,&cinfo);
+	rank = generic_gauss_elim_hensel_dalg(mobj,&sol,&dnsol,&rinfo,&cinfo);
 	get_eg(&eg1); add_eg(&eg_le,&eg0,&eg1);
 	if ( cinfo[0] == dim ) {
 		/* the input is invertible */
@@ -724,21 +726,12 @@ NODE inv_or_split_dalg(DAlg a,DAlg *c)
 		return 0;
 	} else {
 		/* the input is not invertible */
-		nparam = (dim+1)-rank;
-		/* the index 'dim' should not be in cinfo[] */
+		nparam = sol->col;
 		solmat = (Q **)BDY(sol);
-		for ( k = 0; k < nparam; k++ )
-			if ( cinfo[k] == dim )
-				error("invdalg : cannot happen");
 		nd0 = 0;
 		for ( k = 0; k < nparam; k++ ) {
-			m = mb[cinfo[k]];
-			for ( ndt = nd0; ndt; ndt = NEXT(ndt) ) {
-				if ( dp_redble(m,(DP)BDY(ndt)) ) break;
-			}
-			/* skip a redundunt basis element */
-			if ( ndt ) continue;
 			/* construct a new basis element */
+			m = mb[cinfo[k]];
 			mp0 = 0;
 			NEXTMP(mp0,mp);
 			chsgnq(dnsol,&dn1); mp->c = (P)dn1;