version 1.65, 2013/12/20 02:02:23 |
version 1.67, 2015/08/08 14:19:41 |
|
|
* 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/builtin/array.c,v 1.64 2013/11/05 02:55:02 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/builtin/array.c,v 1.66 2015/03/04 08:44:26 ohara Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "base.h" |
#include "base.h" |
Line 231 int lu_gauss(int *ul,ent **u,int *ll,ent **l,int n,int |
|
Line 231 int lu_gauss(int *ul,ent **u,int *ll,ent **l,int n,int |
|
for ( i = 0; i < n; i++ ) ll2[i] = 0; |
for ( i = 0; i < n; i++ ) ll2[i] = 0; |
for ( i = 0; i < n; i++ ) { |
for ( i = 0; i < n; i++ ) { |
fprintf(stderr,"i=%d\n",i); |
fprintf(stderr,"i=%d\n",i); |
|
#if defined(__MINGW32__) || defined(__MINGW64__) |
|
fflush(stderr); |
|
#endif |
inv = invm(u[i][0].e,mod); |
inv = invm(u[i][0].e,mod); |
for ( k = i+1; k < n; k++ ) |
for ( k = i+1; k < n; k++ ) |
if ( u[k][0].j == n-i ) { |
if ( u[k][0].j == n-i ) { |
Line 895 void Pvtol(NODE arg,LIST *rp) |
|
Line 898 void Pvtol(NODE arg,LIST *rp) |
|
pointer *a; |
pointer *a; |
int len,i; |
int len,i; |
|
|
|
if ( OID(ARG0(arg)) == O_LIST ) { |
|
*rp = ARG0(arg); |
|
return; |
|
} |
asir_assert(ARG0(arg),O_VECT,"vtol"); |
asir_assert(ARG0(arg),O_VECT,"vtol"); |
v = (VECT)ARG0(arg); len = v->len; a = BDY(v); |
v = (VECT)ARG0(arg); len = v->len; a = BDY(v); |
for ( i = len - 1, n = 0; i >= 0; i-- ) { |
for ( i = len - 1, n = 0; i >= 0; i-- ) { |
Line 906 void Pvtol(NODE arg,LIST *rp) |
|
Line 913 void Pvtol(NODE arg,LIST *rp) |
|
void Pltov(NODE arg,VECT *rp) |
void Pltov(NODE arg,VECT *rp) |
{ |
{ |
NODE n; |
NODE n; |
VECT v; |
VECT v,v0; |
int len,i; |
int len,i; |
|
|
|
if ( OID(ARG0(arg)) == O_VECT ) { |
|
v0 = (VECT)ARG0(arg); len = v0->len; |
|
MKVECT(v,len); |
|
for ( i = 0; i < len; i++ ) { |
|
BDY(v)[i] = BDY(v0)[i]; |
|
} |
|
*rp = v; |
|
return; |
|
} |
asir_assert(ARG0(arg),O_LIST,"ltov"); |
asir_assert(ARG0(arg),O_LIST,"ltov"); |
n = (NODE)BDY((LIST)ARG0(arg)); |
n = (NODE)BDY((LIST)ARG0(arg)); |
len = length(n); |
len = length(n); |
Line 1762 int generic_gauss_elim_hensel(MAT mat,MAT *nmmat,Q *dn |
|
Line 1778 int generic_gauss_elim_hensel(MAT mat,MAT *nmmat,Q *dn |
|
} |
} |
} |
} |
} |
} |
|
#if defined(__MINGW32__) || defined(__MINGW64__) |
|
fflush(stderr); |
|
#endif |
} |
} |
|
|
int generic_gauss_elim_hensel_dalg(MAT mat,DP *mb,MAT *nmmat,Q *dn,int **rindp,int **cindp) |
int generic_gauss_elim_hensel_dalg(MAT mat,DP *mb,MAT *nmmat,Q *dn,int **rindp,int **cindp) |
Line 1947 int generic_gauss_elim_hensel_dalg(MAT mat,DP *mb,MAT |
|
Line 1966 int generic_gauss_elim_hensel_dalg(MAT mat,DP *mb,MAT |
|
} |
} |
} |
} |
} |
} |
|
#if defined(__MINGW32__) || defined(__MINGW64__) |
|
fflush(stderr); |
|
#endif |
} |
} |
|
|
int f4_nocheck; |
int f4_nocheck; |