समस्या कास्टिंग एसटीएल जटिल <double> fftw_complex को

Question

FFTW मैनुअल says कि अपनी fftw_complex प्रकार बिट एसटीएल में std::complex<double> वर्ग के लिए संगत है। लेकिन यह मेरे लिए काम नहीं करता:

error: invalid cast from type ‘std::complex<double>’ to type ‘double [2]’ 

क्या मैं गलत कर रहा हूँ:

#include <complex> 
#include <fftw3.h> 
int main() 
{ 
    std::complex<double> x(1,0); 
    fftw_complex fx; 
    fx = reinterpret_cast<fftw_complex>(x); 
} 

यह मैं एक त्रुटि देता है?

Answer 1

अपने कोड को फिर से लिखना:

#include <complex> 
#include <fftw3.h> 
int main() 
{ 
    std::complex<double> x(1,0); 
    fftw_complex fx; 
    memcpy(&fx, &x, sizeof(fftw_complex)); 
} 

हर संकलक मैं memcpy बाहर अनुकूलित करेंगे का उपयोग किया है, क्योंकि यह एक निश्चित अर्थात संकलन समय पर, डेटा की मात्रा को कॉपी कर रहा है।

यह pointer aliasing issues से बचाता है।

संपादित करें: तुम भी इस प्रकार एक संघ का उपयोग कर सख्त अलियासिंग मुद्दों से बचने कर सकते हैं:

#include <complex> 
#include <fftw3.h> 
int main() 
{ 
    union stdfftw 
    { 
     std::complex<double> stdc; 
     fftw_complex   fftw; 
    }; 
    std::complex<double> x(1,0); 
    stdfftw u; 
    u.stdc = x; 
    fftw_complex fx = u.fftw; 
} 

हालांकि सख्ती से इस C99 नियम (सी के बारे में ++ सुनिश्चित नहीं हैं) एक संघ का एक अलग सदस्य से पढ़ने के रूप में टूट रहे हैं लिखे गए एक को भी अपरिभाषित है। यह हालांकि अधिकांश कंपाइलरों पर काम करता है। व्यक्तिगत रूप से मैं अपनी मूल विधि पसंद करता हूं।

Answer 2

reinterpret_cast केवल पॉइंटर्स और संदर्भों के लिए काम करता है। तो तुम यह करने के लिए होगा:

#include <complex> 
#include <fftw3.h> 
int main() 
{ 
    std::complex<double> x(1,0); 
    fftw_complex fx(*reinterpret_cast<fftw_complex*>(&x)); 
} 

इसमें यह माना जाता fftw_complex एक प्रति निर्माता है। सख्त एलियासिंग के साथ समस्याओं से बचने के लिए, Goz's solution को प्राथमिकता दी जानी चाहिए। इस प्रकार

Answer 3

fftw_complex और सी 99 और सी ++ जटिल प्रकारों की बिट-संगतता के पीछे विचार यह नहीं है कि उन्हें आसानी से एक दूसरे से बनाया जा सकता है, लेकिन एफएफटीडब्ल्यू में सभी फ़ंक्शन जो fftw_complex पर पॉइंटर्स लेते हैं, वे भी C++ std :: जटिल।

std::vector<std::complex<double> > a1, a2; 
.... 
.... 
fftw_plan_dft(N, reinterpret_cast<fftw_complex*>(&a1[0]), 
       reinterpret_cast<fftw_complex*>(&a2[0]), 
       FFTW_FORWARD, FFTW_ESTIMATE); 
....