diff options
author | Derick Huth <derickhuth@gmail.com> | 2015-09-24 11:34:04 -0600 |
---|---|---|
committer | Derick Huth <derickhuth@gmail.com> | 2015-09-24 11:34:04 -0600 |
commit | 99c062c028c7f4e94fb38cde50772cfd3ea5ad3b (patch) | |
tree | 2d86a6d5a9b7dadd88ee83f9fc0f576a4ce1d451 /test-chill/test-cases/examples/chill/qr.c | |
parent | c285135eb903c31cd221f90f03e288a6b67770cd (diff) | |
download | chill-99c062c028c7f4e94fb38cde50772cfd3ea5ad3b.tar.gz chill-99c062c028c7f4e94fb38cde50772cfd3ea5ad3b.tar.bz2 chill-99c062c028c7f4e94fb38cde50772cfd3ea5ad3b.zip |
v0.2.1
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; -} |