summaryrefslogtreecommitdiff
path: root/test-chill/test-cases/examples/chill/qr.c
diff options
context:
space:
mode:
authorDerick Huth <derickhuth@gmail.com>2015-09-24 12:22:41 -0600
committerDerick Huth <derickhuth@gmail.com>2015-09-24 12:22:41 -0600
commit4631ad76927d433da5d55c3c373a1dfd0f74c9d4 (patch)
treef8dcba88576ec95e403f0c14efd80e970f30a260 /test-chill/test-cases/examples/chill/qr.c
parent6eb2b89896da66a77d0dcdf2d72b98c122826949 (diff)
parent0cff3f9a3c4ccd434900162ebef4bd814850f481 (diff)
downloadchill-4631ad76927d433da5d55c3c373a1dfd0f74c9d4.tar.gz
chill-4631ad76927d433da5d55c3c373a1dfd0f74c9d4.tar.bz2
chill-4631ad76927d433da5d55c3c373a1dfd0f74c9d4.zip
Merge pull request #7 from dhuth/master
V0.2.1
Diffstat (limited to 'test-chill/test-cases/examples/chill/qr.c')
-rw-r--r--test-chill/test-cases/examples/chill/qr.c44
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;
-}