xdsopl/robot36
1
0
Fork 0
mirror of https://github.com/xdsopl/robot36.git synced 2025-12-06 07:12:07 +01:00
Code Issues Actions Packages Projects Releases Wiki Activity

fixed bug introduced with "replaced radix2 input with float"

Browse source 
it really sucks not having a complex data type .. damn you c99
This commit is contained in:
Ahmet Inan 2015-01-15 16:23:03 +01:00
parent 5fc7dfc458
commit c0f911f1e9
1 changed files with 7 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

14
app/src/main/rs/radix2.rsh
View file

@ -23,18 +23,18 @@ limitations under the License.
static void radix2(complex_t *out, float *in, int N, int S, int L)
{
if (1 == N) {
out[0] = in[0];
out[0] = complex(in[0], 0.0f);
return;
} else if (2 == N) {
out[0] = in[0] + in[S];
out[1] = in[0] - in[S];
out[0] = complex(in[0] + in[S], 0.0f);
out[1] = complex(in[0] - in[S], 0.0f);
return;
} else if (4 == N) {
complex_t w = radix2_z[1 << L];
out[0] = in[0] + in[S] + in[2 * S] + in[3 * S];
out[1] = in[0] + w * in[S] - in[2 * S] - w * in[3 * S];
out[2] = in[0] - in[S] + in[2 * S] - in[3 * S];
out[3] = in[0] - w * in[S] - in[2 * S] + w * in[3 * S];
out[0] = complex(in[0] + in[S] + in[2 * S] + in[3 * S], 0.0f);
out[1] = complex(in[0] - in[2 * S], 0.0f) + w * (in[S] - in[3 * S]);
out[2] = complex(in[0] - in[S] + in[2 * S] - in[3 * S], 0.0f);
out[3] = complex(in[0] - in[2 * S], 0.0f) + w * (in[3 * S] - in[S]);
return;
}
radix2(out, in, N / 2, 2 * S, L + 1);