Question

无数的GSL函数将其结果作为指针返回到第一个参数中。例如

int gsl_matrix_get_col（gsl_vector * v，const gsl_matrix * m，size_t j）

我的编程水平非常低，但我被告知这样的事情是不可能的局部变量（在函数结束时删除），但可以使用指针，只要它们被调用函数声明和正确分配。我觉得很奇怪，指针和普通变量之间应该存在这样的根本差异，但是我试图将这个结果存储在一个简单的GSL程序的变量中，我想要一个函数（ fetch_eigenvalue（） ）输出两件事。我失败了。我的计划如下：

#include <math.h>
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_blas.h>


/* Parameters */
#define N 3


int CREATE_MATRIX_AND_VECTOR(gsl_matrix *m, gsl_vector *v);
double fetch_eigenvalue(gsl_matrix *M, gsl_vector *v, double *rescos);

int main()
{
gsl_matrix *unit_matrix = gsl_matrix_calloc(N, N); //soon to be unity 
gsl_vector *v = gsl_vector_calloc(N);          //soon to be unit x vector
double *outcos = (double*)malloc(sizeof(double) );

printf("**********************************************\n");
CREATE_MATRIX_AND_VECTOR(unit_matrix, v);
fetch_eigenvalue(unit_matrix, v, outcos);   
printf("==IN MAIN: outcos = %e\n", *outcos);

free((void *)outcos);
gsl_vector_free(v);
gsl_matrix_free(unit_matrix);
printf("**********************************************\n");
return(0);
}

int CREATE_MATRIX_AND_VECTOR(gsl_matrix * m, gsl_vector *v)
{
int i;

    for (i = 0; i < N; i++)
    {
        gsl_matrix_set(m, i, i, 1.0);
    }

    gsl_vector_set(v, 0, 1.0);

return(0);
}

double fetch_eigenvalue(gsl_matrix *M, gsl_vector *v, double *rescos)   //fetches eigenvalue, if Mv is parallel to v within accuracy gvaccu
                                    //rescos is the cosine of the angle between Mv and v
{

int i,lv;
double v0, v1, cos;
double result;
double vnorm, pnorm;
double rdot;
lv = v->size;

double gvaccu = 1e-10;


            gsl_vector *prod = gsl_vector_calloc(lv);
            gsl_matrix_get_row(prod, M, 0);

            if(gsl_blas_dnrm2(prod)==0.0)
            {
                result = 0.0;
            }
            else
            {
                gsl_blas_dgemv( CblasNoTrans,1.0, M, v, 0.0, prod);
                gsl_blas_ddot(prod, v, &rdot);
                pnorm = gsl_blas_dnrm2(prod);
                vnorm = gsl_blas_dnrm2(v);
                cos = rdot/pnorm/vnorm;
                cos = fabs(cos);
                rescos = &cos;

                if(fabs(cos -1.0) > gvaccu)
                {
                    result = -1.0;
                }
                else
                {
                v0 = gsl_vector_get(v,0);
                v1 = gsl_vector_get(prod,0);
                result = v1/v0;
                }
            }
        printf("==IN FETCH_EV: COS = %e\n", cos);//print cheat!!
        printf("==IN FETCH_EV: RESCOS = %e\n", *rescos);//print cheat!!
gsl_vector_free(prod);
return(result);     
}

我运行它并得到以下输出：

ludi@ludi-M17xR4:~/Desktop/Healpix_3.20$ g++ -o wrong_output wrong_output.c -L. -L/sw/lib -I/sw/include -lgsl -lblas && ./wrong_output 
**********************************************
==IN FETCH_EV: COS = 1.000000e+00
==IN FETCH_EV: RESCOS = 1.000000e+00
==IN MAIN: outcos = 0.000000e+00
**********************************************
ludi@ludi-M17xR4:~/Desktop/Healpix_3.20$

所以，调用者main（）对fetch_eigenvalue（）中发生的事情一无所知，尽管我使用了指针。我究竟做错了什么？我有一种感觉，我误解了一些非常重要的东西。

Answer 1

我总结了您对rescos函数中参数fetch_eigenvalue的处理方式：

double fetch_eigenvalue(gsl_matrix *M, gsl_vector *v, double *rescos)
{
    double cos;
    // some code
    rescos = &cos;
    // some code
    return(result);
}

此处您未修改rescos指向的双精度值，您正在修改变量rescos本身，这是使用的变量outcos的副本在你的主要。

您要在fetch_eigenvalue中执行的操作是将cos的值复制到rescos指向的变量中：

double fetch_eigenvalue(gsl_matrix *M, gsl_vector *v, double *rescos)
{
    double cos;
    // some code
    *rescos = cos;
    // some code
    return(result);
}

编辑：正如其他答案所述，尽可能避免使用malloc，在这里你可以：

double outcos;
fetch_eigenvalue(unit_matrix, v, &outcos);

Answer 2

我怀疑这是因为语句rescos = &cos;将rescos保存到本地变量cos的地址中。但是，此变量的范围仅为本地，因此您无法在main()中使用它。我想你想要做的是：

在rescos = &cos;函数中将*rescos = cos;更改为fetch_eigenvalue，以便将cos的值存储在rescos指向的地址< / p>
仅在double outcos;函数中使用main()，即不使用指针
将fetch_eigenvalue称为fetch_eigenvalue(unit_matrix, v, &outcos);

Answer 3

我对GSL一无所知，但它似乎是一个使用C风格界面的库。要从函数设置值，它们使用指针。您似乎还不知道如何使用这样的API，所以这里有一些提示。

陈述

double *outcos = (double*)malloc(sizeof(double) );
...
fetch_eigenvalue(unit_matrix, v, outcos);

不是您想要使用此类API的方式。相反，您只需定义一个double变量，并在调用中使用operator的地址：

double outcos;
...
fetch_eigenvalue(unit_matrix, v, &outcos);

此外，在您的方法中，要分配值，请使用

        cos = fabs(cos);
        rescos = &cos;

但

        cos = fabs(cos);
        *rescos = cos;

将值赋给指向的变量，而不是指针。希望这会有所帮助。

在函数参数

3 个答案: