| 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"); |
| } |
} |