diff options
| author | Tuowen Zhao <ztuowen@gmail.com> | 2016-09-18 21:39:45 +0000 |
|---|---|---|
| committer | Tuowen Zhao <ztuowen@gmail.com> | 2016-09-18 21:39:45 +0000 |
| commit | f255f2498da1fd985ad1ed79362580bbf4675723 (patch) | |
| tree | 3cd02e9d3147054820c58aacd1e5d3d336eea0d7 /test-chill/test-cases/examples/chill/qr.c | |
| parent | 7233faacb7990a6c3a40d2435ede88d7725dfc6e (diff) | |
| download | chill-f255f2498da1fd985ad1ed79362580bbf4675723.tar.gz chill-f255f2498da1fd985ad1ed79362580bbf4675723.tar.bz2 chill-f255f2498da1fd985ad1ed79362580bbf4675723.zip | |
rm test-chill, create doc subdir
Diffstat (limited to 'test-chill/test-cases/examples/chill/qr.c')
| -rw-r--r-- | test-chill/test-cases/examples/chill/qr.c | 44 |
1 files changed, 0 insertions, 44 deletions
diff --git a/test-chill/test-cases/examples/chill/qr.c b/test-chill/test-cases/examples/chill/qr.c deleted file mode 100644 index 8d18b72..0000000 --- a/test-chill/test-cases/examples/chill/qr.c +++ /dev/null @@ -1,44 +0,0 @@ -#include <math.h> - -int main() { - - int M, N; - float** A; - float *s; - float *Rdiag; - float *nrm; - int i, j, k; - float t; - for (k = 0; k < N; k++) { - nrm[k] = 0; - - for (i = k; i < M; i++) - nrm[k] = sqrt(nrm[k] * nrm[k] + A[i][k] * A[i][k]); - //t = A[k][k]; - - //if (t < 0) - // nrm[k] = -nrm[k]; - for (i = k; i < M; i++) - A[i][k] = A[i][k] / nrm[k]; - - A[k][k] = A[k][k] + 1; - - for (j = k + 1; j < N; j++) { - s[j] = 0; //S6 - - for (i = k; i < M; i++) - s[j] = s[j] + A[i][k] * A[i][j]; //S7 - - s[j] = -s[j] / A[k][k]; //S8 - - for (i = k; i < M; i++) - A[i][j] = A[i][j] + s[j] * A[i][k]; //S9 - - } - - Rdiag[k] = -nrm[k]; - - } - - return 0; -} |