fftw_execute_dft_r2c(planlos, windowedData, output);
}
+// idea: if match, then sxx * syy = sxy * sxy
+
double vectorDot(fftw_complex *v1, fftw_complex *v2, size_t length)
{
size_t i;
double sum = 0;
for(i = 0; i < length; ++i)
sum += cabs(v1[i]) * cabs(v2[i]);
+ for(i = 0; i < length; ++i)
+ sum += 0.001 * creal(v1[i]) * creal(v2[i]);
+ for(i = 0; i < length; ++i)
+ sum += 0.001 * cimag(v1[i]) * cimag(v2[i]);
return sum;
}