在 gsl 中分配复数值

我正在尝试使用 GSLcomplex numbers, complex vectors and complex matrices在我的项目中。我使用的是VS2010，我在中添加了库的地址配置属性>C/C++>常规>其他包含目录。但我有一个愚蠢的问题。据我了解，我不能使用=分配两个gsl_complex, gsl_vector_complex or gsl_matrix_complex对彼此。

对于向量我必须使用gsl_vector_complex_set对于矩阵gsl_matrix_complex_set。但是对于gsl_complex，我只发现GSL_SET_COMPLEX其中我应该将实部和虚部分别作为两个参数给出：

GSL_SET_COMPLEX (zp, real, imaginary)

在我的代码中我有这样的功能：

gsl_complex cmx, cmx2;
void vector_complex_exp(gsl_vector_complex *v)
{
    for (i = 0; i < v->size; i++)
    {
        gsl_vector_complex_set(v, i, gsl_complex_exp(gsl_vector_complex_get(v, i)));
    }
}

使用这个，我得到以下错误：

错误 LNK1120：2 未解析的外部引用。

错误 LNK2001：无法解析的外部符号“_hypot”。

错误 LNK2001：无法解析的外部符号“_log1p”。

错误 LNK2001：无法解析的外部符号“_log1p”。

我不明白这些错误背后的原因。但我这样重写我的代码：

void vector_complex_exp(gsl_vector_complex *v)
{
    for (i = 0; i < v->size; i++)
    {
        cmx = gsl_vector_complex_get(v, i);
        //cmx2 = gsl_complex_exp(cmx);
        gsl_vector_complex_set(v, i, cmx2);
    }
}

这里当 for 的第二行被注释时，没有错误。但是当我取消注释时，我得到以下信息：

错误 LNK1120：2 个未解析的外部引用。

错误 LNK2001：无法解析的外部符号“_log1p”。

错误 LNK2019：在函数“_gsl_complex_div”中引用未解析的外部符号“_hypot”。

错误 LNK2019：在函数“_gsl_complex_logabs”中引用未解析的外部符号“_log1p”。

我没有_gsl_complex_div or _gsl_complex_logabs我的代码中的函数。所以我很确定问题出在这里的分配上。但我也不能在这里使用 GSL_SET_COMPLEX 。

有人可以帮我弄这个吗？难道真的没有办法直接给gsl_complex赋值吗？

如果您在这里发布所有代码会更好，因此，我立即使用了以下代码最低的例子GSL 的。我做了一些小改变：

#include <stdio.h>
#include <gsl/gsl_math.h>
#include <gsl/gsl_eigen.h>

int main(void)
{
double data[] = { -1.0, 1.0, -1.0, 1.0,
    -8.0, 4.0, -2.0, 1.0,
    27.0, 9.0, 3.0, 1.0,
    64.0, 16.0, 4.0, 1.0 };

gsl_matrix_view m
    = gsl_matrix_view_array(data, 4, 4);

gsl_vector_complex *eval = gsl_vector_complex_alloc(4);
gsl_matrix_complex *evec = gsl_matrix_complex_alloc(4, 4);

gsl_eigen_nonsymmv_workspace * w = gsl_eigen_nonsymmv_alloc(4);

gsl_eigen_nonsymmv(&m.matrix, eval, evec, w);

gsl_eigen_nonsymmv_free(w);

gsl_eigen_nonsymmv_sort(eval, evec, GSL_EIGEN_SORT_ABS_DESC);

{
    int i, j;

    for (i = 0; i < 4; i++)
    {
        gsl_complex eval_i
            = gsl_vector_complex_get(eval, i);
        gsl_vector_complex_view evec_i
            = gsl_matrix_complex_column(evec, i);

        printf("\n eigenvalue = %g + %gi\n",
            GSL_REAL(eval_i), GSL_IMAG(eval_i));
        printf(" eigenvector = \n");
        for (j = 0; j < 4; ++j)
        {
            gsl_complex z =
                gsl_vector_complex_get(&evec_i.vector, j);
            printf("        %g + %gi\n", GSL_REAL(z), GSL_IMAG(z));
        }
    }
}

gsl_vector_complex_free(eval);
gsl_matrix_complex_free(evec);

system("pause");

return 0;
}

该代码的输出：（从红色箭头和下面来看，与 GSL 示例中的预期输出不匹配） 要获得我的输出，您需要 IDE（我使用 Visual Studio 2015）：

1. 进入“your_application”属性页 -> 配置属性 -> VC++ 目录 -> （右窗格）行可执行目录 type: C：\用户\...（您的 GSL 构建目录路径）...\gsl\Release;$(可执行路径)
2. 同上。排队包括目录 type: C：\用户\...（您的 GSL 构建目录路径）...\gsl;$(包含路径)
3. 同上。排队图书馆目录 type: C：\用户\...（您的 GSL 构建目录路径）...\gsl\Release;$(库路径)
4. 下面，在左窗格中，选择 C/C++ -> Peprocessor -> (右窗格)预处理器定义 type: WIN32;_DEBUG;_CONSOLE;GSL_DLL;%(预处理器定义)（我使用调试模式，创建空控制台应用程序）。保存设置（按“应用”和“确定”按钮）
5. 复制并放入应用程序项目文件夹的 Debug 目录中gsl.dll and dll from C：\用户\...（您的 GSL 构建目录路径）...\gsl\发布目录
6. 构建您的应用程序并运行它
7. NOTE!:一开始，最好使用目标应用程序的编译器重建 GSL - 这样工作就能得到保证。

祝你好运！