// calculate the real and imaginary results
// scaling appropriately
- real = (q1 - q2 * cosine) / scalingFactor;
- imag = (q2 * sine) / scalingFactor;
+ real = (q1 * cosine - q2) / scalingFactor;
+ imag = (q1 * sine) / scalingFactor;
magnitude = sqrtf(real*real + imag*imag);
return magnitude;