| 95 | |
| 96 | |
| 97 | void dMakeRandomMatrix (dReal *A, int n, int m, dReal range) |
| 98 | { |
| 99 | int i,j; |
| 100 | int skip = dPAD(m); |
| 101 | dSetZero (A,n*skip); |
| 102 | for (i=0; i<n; i++) { |
| 103 | for (j=0; j<m; j++) A[i*skip+j] = (dRandReal()*REAL(2.0)-REAL(1.0))*range; |
| 104 | } |
| 105 | } |
| 106 | |
| 107 | |
| 108 | void dClearUpperTriangle (dReal *A, int n) |