将std :: complex <double>转换为__complex128

时间:2017-01-21 05:30:54

标签: c++ floating-point precision quad

我正在尝试在GCC中使用quadmath库。我有一个复杂的double值,我想将其转换为相应的四精度复数__complex128。以下是最小(非)工作示例:

#include <quadmath.h>
#include <complex>
#include <stdio.h>
using namespace std::complex_literals;

int main(){
   std::complex<double> x = 1 + 2i;
   std::printf("x = %5.5g + %5.5g\n", x.real(), x.imag()); 
   __complex128 y = 2+2i;
   y = x;
  return 0;
}

当我尝试使用

编译此代码时
 g++ test.cpp -lquadmath -o test

我收到以下错误:

test.cpp:10:6: error: cannot convert 'std::complex<double>' to '__complex128 {aka __complex__ __float128}' in assignment
y = x;

如果我尝试使用显式类型转换来替换赋值行,

y = (__complex128) x;

我收到类似的错误

test.cpp:10:21: error: invalid cast from type 'std::complex<double>' to type '__complex128 {aka __complex__ __float128}'
y = (__complex128) x;

如何在这两种类型之间进行转换?

2 个答案:

答案 0 :(得分:2)

我猜您使用的是GCC,在这种情况下,您可以使用__real____imag__扩展来设置__complex128的各个组成部分:

__complex128 y;
__real__ y = x.real();
__imag__ y = x.imag();

这也适用于Clang for __complex64(Clang还不支持__complex128)。

答案 1 :(得分:0)

我必须假设这里存在某种类型的兼容性问题,因为据我所知,__complex__非常古老(见https://gcc.gnu.org/onlinedocs/gcc/Complex.html)。作为解决此问题的方法,您可以尝试:

y = 1.0i;
y *= x.imag();
y += x.real();