| 761 |
|
return(1); |
| 762 |
|
} |
| 763 |
|
if (rm->nrows == rm->ncols) { /* square matrix case */ |
| 764 |
< |
for (i = rm->nrows; i--; ) |
| 765 |
< |
for (j = rm->ncols; j--; ) { |
| 766 |
< |
if (i == j) continue; |
| 764 |
> |
for (i = rm->nrows; --i > 0; ) |
| 765 |
> |
for (j = i; j-- > 0; ) { |
| 766 |
|
memcpy(val, rmx_val(rm,i,j), |
| 767 |
|
sizeof(rmx_dtype)*rm->ncomp); |
| 768 |
|
memcpy(rmx_lval(rm,i,j), rmx_val(rm,j,i), |
| 769 |
|
sizeof(rmx_dtype)*rm->ncomp); |
| 770 |
< |
memcpy(rmx_val(rm,j,i), val, |
| 770 |
> |
memcpy(rmx_lval(rm,j,i), val, |
| 771 |
|
sizeof(rmx_dtype)*rm->ncomp); |
| 772 |
|
} |
| 773 |
< |
return(1); |
| 773 |
> |
return(1); |
| 774 |
|
} |
| 775 |
|
#define bmbyte(r,c) bmap[((r)*rm->ncols+(c))>>3] |
| 776 |
|
#define bmbit(r,c) (1 << ((r)*rm->ncols+(c) & 7)) |
