version 1.7, 2000/08/22 05:03:56 |
version 1.8, 2000/09/21 09:19:25 |
|
|
* 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.6 2000/08/21 08:31:18 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/builtin/array.c,v 1.7 2000/08/22 05:03:56 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "base.h" |
#include "base.h" |
|
|
t = mat[i]; |
t = mat[i]; |
if ( i != j && (a = t[j]) ) |
if ( i != j && (a = t[j]) ) |
for ( k = j, a = md - a; k <= n; k++ ) { |
for ( k = j, a = md - a; k <= n; k++ ) { |
|
unsigned int tk; |
/* t[k] = dmar(pivot[k],a,t[k],md); */ |
/* t[k] = dmar(pivot[k],a,t[k],md); */ |
DMAR(pivot[k],a,t[k],md,t[k]) |
DMAR(pivot[k],a,t[k],md,tk) |
|
t[k] = tk; |
} |
} |
} |
} |
} |
} |
|
|
DMAR(inv,m,0,md,t[k]) |
DMAR(inv,m,0,md,t[k]) |
for ( j = k+1, m = md - t[k]; j < col; j++ ) |
for ( j = k+1, m = md - t[k]; j < col; j++ ) |
if ( pivot[j] ) { |
if ( pivot[j] ) { |
DMAR(m,pivot[j],t[j],md,t[j]) |
unsigned int tj; |
|
|
|
DMAR(m,pivot[j],t[j],md,tj) |
|
t[j] = tj; |
} |
} |
} |
} |
} |
} |
Line 1517 int **rinfo,**cinfo; |
|
Line 1522 int **rinfo,**cinfo; |
|
DMAR(inv,m,0,md,t[k]) |
DMAR(inv,m,0,md,t[k]) |
for ( j = k+1, m = md - t[k]; j < col; j++ ) |
for ( j = k+1, m = md - t[k]; j < col; j++ ) |
if ( pivot[j] ) { |
if ( pivot[j] ) { |
DMAR(m,pivot[j],t[j],md,t[j]) |
unsigned int tj; |
|
DMAR(m,pivot[j],t[j],md,tj) |
|
t[j] = tj; |
} |
} |
} |
} |
} |
} |
|
|
|
|
for ( i = 0; i < row; i++ ) { |
for ( i = 0; i < row; i++ ) { |
for ( j = 0; j < col; j++ ) { |
for ( j = 0; j < col; j++ ) { |
printnum(mat[i][j]); printf(" "); |
printnum((Num)mat[i][j]); printf(" "); |
} |
} |
printf("\n"); |
printf("\n"); |
} |
} |