version 1.211, 2013/09/26 08:55:11 |
version 1.213, 2013/09/27 02:45:17 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.210 2013/09/26 00:38:47 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.212 2013/09/27 02:35:15 noro Exp $ */ |
|
|
#include "nd.h" |
#include "nd.h" |
|
|
Line 6262 NODE nd_f4_red_q_main(ND_pairs sp0,int nsp,int trace,U |
|
Line 6262 NODE nd_f4_red_q_main(ND_pairs sp0,int nsp,int trace,U |
|
rank = nd_gauss_elim_q(spmat,spsugar,sprow,spcol,colstat); |
rank = nd_gauss_elim_q(spmat,spsugar,sprow,spcol,colstat); |
w = (pointer *)ALLOCA(rank*sizeof(pointer)); |
w = (pointer *)ALLOCA(rank*sizeof(pointer)); |
for ( i = 0; i < rank; i++ ) { |
for ( i = 0; i < rank; i++ ) { |
|
#if 0 |
w[rank-i-1] = (pointer)vect_to_ndv_q(spmat[i],spcol,col,rhead,s0vect); |
w[rank-i-1] = (pointer)vect_to_ndv_q(spmat[i],spcol,col,rhead,s0vect); |
SG((NDV)w[rank-i-1]) = spsugar[i]; |
SG((NDV)w[rank-i-1]) = spsugar[i]; |
|
#else |
|
w[i] = (pointer)vect_to_ndv_q(spmat[i],spcol,col,rhead,s0vect); |
|
SG((NDV)w[i]) = spsugar[i]; |
|
#endif |
/* GCFREE(spmat[i]); */ |
/* GCFREE(spmat[i]); */ |
} |
} |
#if 0 |
#if 0 |