为什么NumPy.exp在C

时间:2017-01-05 07:12:32

标签: python c numpy exponent

我正在将一些Python代码转换为C代码。

下面的Python NumPy exp (6.12323399574e-17-1j)k=1的复数输出l=4

numpy.exp(-2.0*1j*np.pi*k/l)

我将它转换为C代码,如下所示。但输出为1.000000

#include <complex.h>
#define PI 3.1415926535897932384626434
exp(-2.0*I*PI*k/l)

我错过了什么?

2 个答案:

答案 0 :(得分:2)

您必须使用cimagcreal来打印数据。

C版:

#include <stdio.h>
#include <complex.h>
#include <tgmath.h>


int main(){
    int k=1;
    int l=4;
    double PI = acos(-1);
    double complex z = exp(-2.0*I*PI*k/l);
    printf(" %.1f%+.1fj\n", creal(z), cimag(z));
    return 0;
}

输出:

0.0-1.0j

Python版本:

import numpy as np

k=1
l=4
z = np.exp(-2.0*1j*np.pi*k/l)
print(z)

输出:

6.12323399574e-17-1j

答案 1 :(得分:0)

这里有正确的C代码来打印你的答案:

    x = cexp(-2.0*I*PI*0.25);
    printf("%f + i%f\n", creal(x), cimag(x));