summaryrefslogtreecommitdiffstats
path: root/src/ffts_real_nd.c
diff options
context:
space:
mode:
author= <henry.gouk@gmail.com>2013-06-02 23:44:28 +1200
committer= <henry.gouk@gmail.com>2013-06-02 23:44:28 +1200
commit32b31f840aea0a79b028c3f77a254a08adaa1298 (patch)
tree0d1b831271c8b3b38c77a7081b458a73a5b1ab67 /src/ffts_real_nd.c
parent31dfafd8c4a4141ad3afa44d70a861bed07ca5eb (diff)
downloadffts-32b31f840aea0a79b028c3f77a254a08adaa1298.zip
ffts-32b31f840aea0a79b028c3f77a254a08adaa1298.tar.gz
Fixed #5
Diffstat (limited to 'src/ffts_real_nd.c')
-rw-r--r--src/ffts_real_nd.c10
1 files changed, 5 insertions, 5 deletions
diff --git a/src/ffts_real_nd.c b/src/ffts_real_nd.c
index fe0b402..bf46254 100644
--- a/src/ffts_real_nd.c
+++ b/src/ffts_real_nd.c
@@ -79,9 +79,9 @@ void ffts_execute_nd_real(ffts_plan_t *p, const void * in, void * out) {
size_t i,j;
for(i=0;i<p->Ns[0];i++) {
- p->plans[0]->transform(p->plans[0], din + (i * p->Ms[0]), buf + (i * p->Ms[0]));
+ p->plans[0]->transform(p->plans[0], din + (i * p->Ms[0]), buf + (i * (p->Ms[0] / 2 + 1)));
}
- ffts_scalar_transpose(buf, dout, p->Ms[0], p->Ns[0], p->transpose_buf);
+ ffts_scalar_transpose(buf, dout, p->Ms[0] / 2 + 1, p->Ns[0], p->transpose_buf);
for(i=1;i<p->rank;i++) {
for(j=0;j<p->Ns[i];j++) {
@@ -159,9 +159,9 @@ ffts_plan_t *ffts_init_nd_real(int rank, size_t *Ns, int sign) {
p->Ns[i] = p->Ns[i] / 2 + 1;
}
}else{
- for(i=0;i<rank-1;i++) {
- p->Ms[i] = p->Ms[i] / 2 + 1;
- }
+ for(i=0;i<rank-1;i++) {
+ p->Ms[i] = p->Ms[i] / 2 + 1;
+ }
}
p->transpose_buf = valloc(sizeof(float) * 2 * 8 * 8);
OpenPOWER on IntegriCloud