summaryrefslogtreecommitdiff
path: root/test-chill/test-cases/examples/cuda-chill/mriq-fh.c
diff options
context:
space:
mode:
authordhuth <derickhuth@gmail.com>2014-09-17 18:09:29 -0600
committerdhuth <derickhuth@gmail.com>2014-09-17 18:09:29 -0600
commit600fa18324c21a162c50c40ae5f00c899a41dd24 (patch)
treed399a8ea49c71a85abf5c07cb96b24676df32a0a /test-chill/test-cases/examples/cuda-chill/mriq-fh.c
parenta2bd0557344bbd8d06e94814abd409f552b0efec (diff)
downloadchill-600fa18324c21a162c50c40ae5f00c899a41dd24.tar.gz
chill-600fa18324c21a162c50c40ae5f00c899a41dd24.tar.bz2
chill-600fa18324c21a162c50c40ae5f00c899a41dd24.zip
removed submodule, added test-chill
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;
+ }
+
+ }
+}
+