diff options
author | = <henry.gouk@gmail.com> | 2013-06-02 23:44:28 +1200 |
---|---|---|
committer | = <henry.gouk@gmail.com> | 2013-06-02 23:44:28 +1200 |
commit | 32b31f840aea0a79b028c3f77a254a08adaa1298 (patch) | |
tree | 0d1b831271c8b3b38c77a7081b458a73a5b1ab67 /src/ffts_real_nd.c | |
parent | 31dfafd8c4a4141ad3afa44d70a861bed07ca5eb (diff) | |
download | ffts-32b31f840aea0a79b028c3f77a254a08adaa1298.zip ffts-32b31f840aea0a79b028c3f77a254a08adaa1298.tar.gz |
Fixed #5
Diffstat (limited to 'src/ffts_real_nd.c')
-rw-r--r-- | src/ffts_real_nd.c | 10 |
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); |