| 418 |
|
rmx_identity(const int dim, const int n) |
| 419 |
|
{ |
| 420 |
|
RMATRIX *rid = rmx_alloc(dim, dim, n); |
| 421 |
< |
int i; |
| 421 |
> |
int i, k; |
| 422 |
|
|
| 423 |
|
if (rid == NULL) |
| 424 |
|
return(NULL); |
| 425 |
< |
memset(rid->mtx, 0, sizeof(rid->mtx[0])*dim*dim); |
| 425 |
> |
memset(rid->mtx, 0, sizeof(rid->mtx[0])*n*dim*dim); |
| 426 |
|
for (i = dim; i--; ) |
| 427 |
< |
rmx_lval(rid,i,i,0) = 1; |
| 428 |
< |
for (i = n; --i; ) |
| 429 |
< |
memcpy(rid->mtx+i*(dim*dim), rid->mtx, |
| 430 |
< |
sizeof(rid->mtx[0])*dim*dim); |
| 427 |
> |
for (k = n; k--; ) |
| 428 |
> |
rmx_lval(rid,i,i,k) = 1; |
| 429 |
|
return(rid); |
| 430 |
|
} |
| 431 |
|
|