C++ calculation DFT