summaryrefslogtreecommitdiff
path: root/test-chill/test-cases/examples/cuda-chill/mriq-fh.c
diff options
context:
space:
mode:
authorDerick Huth <derickhuth@gmail.com>2016-01-18 15:43:52 -0700
committerDerick Huth <derickhuth@gmail.com>2016-01-18 15:43:52 -0700
commit983749787ee0dc1beb1107873e8a13ebdaeba576 (patch)
treee9bff337b4d5582b87ad2edc25baa4d3b0c163fa /test-chill/test-cases/examples/cuda-chill/mriq-fh.c
parent0cff3f9a3c4ccd434900162ebef4bd814850f481 (diff)
downloadchill-983749787ee0dc1beb1107873e8a13ebdaeba576.tar.gz
chill-983749787ee0dc1beb1107873e8a13ebdaeba576.tar.bz2
chill-983749787ee0dc1beb1107873e8a13ebdaeba576.zip
restore test suite
Diffstat (limited to 'test-chill/test-cases/examples/cuda-chill/mriq-fh.c')
-rwxr-xr-xtest-chill/test-cases/examples/cuda-chill/mriq-fh.c38
1 files changed, 38 insertions, 0 deletions
diff --git a/test-chill/test-cases/examples/cuda-chill/mriq-fh.c b/test-chill/test-cases/examples/cuda-chill/mriq-fh.c
new file mode 100755
index 0000000..1e924b7
--- /dev/null
+++ b/test-chill/test-cases/examples/cuda-chill/mriq-fh.c
@@ -0,0 +1,38 @@
+#define X 32768
+#define K 256
+struct kValues {
+ float Kx;
+ float Ky;
+ float Kz;
+ float PhiMag;
+};
+extern float sin(float);
+extern float cos(float);
+
+void mriFH_cpu(float *rPhi,float *rRho,float *iRho, float *iPhi, float *rD, float *iD, float *kx, float *ky, float *kz, float *dx, float *dy, float *dz, float *rFHref, float *iFHref)
+{
+
+ float rfh;
+ float ifh;
+ float exp;
+ float cArg;
+ float sArg;
+ //float rRho[K];
+ //float iRho[K];
+ unsigned int k;
+ unsigned int x;
+
+
+ for (x = 0; x < X; ++x) {
+ for (k = 0; k < K; ++k) {
+
+ exp = 2 * 3.14159 * (kx[k]* dx[x] + ky[k]* dy[x] + kz[k]* dz[x]);
+ cArg = cos(exp);
+ sArg = sin(exp);
+ rFHref[x] += rRho[k]* cArg - iRho[k]* sArg;
+ iFHref[x] += iRho[k]*cArg + rRho[k]*sArg;
+ }
+
+ }
+}
+