c - 警告 : assignment from incompatible pointer type [-Wincompatible-pointer-types] qt1. 函数 = &H1

Question

我使用 C 编程语言中的 gsl 库集成了一个复杂的函数。在这段代码中，我必须使用我已成功完成的指针声明两个变量。但是当我在主函数中传递这些变量时遇到问题。

请帮帮我。

我的代码写在这里:

struct har{

double t;
double k;
double x;
};

#include"str.h"
struct har H1( void  * params, float q , int p )

    {

    struct har  *p_params = (void *)params;


    float  mu=2;
    double x =p_params ->x ;
    double t = p_params -> t;
    double  k = p_params -> k;
    //printf("%d",k);
    struct har  H1 =  {x*cos(t*x)/(pow((mu*t*k),2)+pow(x,2))};
    return H1;
    }


 double H (double x,void * params )
{
double e;
double t = *(double *) params;
//printf("%f\n",t);
//struct har t1,z;
//double t = z.params ->t1;
double H  = (pow(e,-x)/x) ;//I(x,&t)*(sin(x*t));  
return H;
}


#include<stdio.h>
#include <math.h>
#include <gsl/gsl_integration.h>
#include"H.h"
#include"H1.h"

int
main (void)
{
  gsl_integration_workspace * w
    = gsl_integration_workspace_alloc (10000);

  struct har t,k;

  double x, qtr, qbartr, qdottr, qddotr,q7r, qatr, qbtr, qt1r, error;
  double expected = -4.0;
  double a1 = 1e-14;
  double a= 150;//150;
  double pi = 3.141;
  double b = 1.05263;
  double mu = (2*pi)/b;
  double q0=0;
  double p0=6.5;
  double om=7.07;
  double tau=0.141;
  double gamma = 0.180;

for(int m=2; m<4; m++)
{


 float  t = 4.14937*(m-1);


 gsl_function qt;
 qt.function = &H;
 qt.params =&t;
 gsl_integration_qags (&qt, a1, a, 0,1e-7, 10000, w, &qtr, &error);
// printf ("qtresult = % .18f\n\n", qtr);

 for(int k=1; k<=2; k++)
    {
      {
struct har item = {x,t,k};
struct har *p_params = &item;


       double  gama=gamma/((mu*k*tau)+1);
//      printf("gama = % .18f\n", gama);

       gsl_function qt1;

       qt1.function = &H1;
       qt1.params = &p_params;
    // qt1.params = &k;
       gsl_integration_qags (&qt1, a1, a, 0,1e-7, 10000, w, &qt1r, &error);
      printf("qt1result= % .18f\n", qt1r);


     }

}


  gsl_integration_workspace_free (w);

  return 0;
}

在上面的代码中，我有四个文件，即 code.c、str.h、H1.h 和 H.h。其中 H1 和 H 包含两个用户定义的函数，str 是结构文件，code.c 是我要编译的主要文件。

感谢您的时间和观察。

问候

拉梅什

Answer 1

字段qt1.function类型为double (*)(double, void *) 。函数H1您尝试分配给此字段的类型为 struct har (*)( void *, float, int) 。

这些类型不兼容，因为参数数量、参数类型和返回类型都不同。现在，您有两个名为 p 的参数。和q从未被使用过。另外，struct har您构建为返回值甚至没有设置所有字段，只是第一个字段，那么为什么您还要返回一个结构体？

为了H1与 gsl_function 一起工作，必须接受 double和一个 void *并返回 double 。 x值对应于double这是传入的，因此您的结构中不需要它。由于您也只是构建一个要返回的值，因此这就是您返回的值。

struct har{
    double t;
    double k;
};

double H1( double x, void  *params)
{
    struct har *p_params = params;   // no need to cast to/from void *

    float mu=2;
    double t = p_params->t;
    double k = p_params->k;

    return  x*cos(t*x)/(pow((mu*t*k),2)+pow(x,2));
}

此外，您没有设置 qt1 的值正确地:

qt1.params = &p_params;

自 p_params是 struct har * ，它的地址是 struct har ** 。然而，您的函数期望 struct har *因此，将其更改为结构体的地址，而不是指向它的指针的地址:

qt1.params = &item;

当您设置 qt.params 时，您也会遇到不匹配的情况。 。你给它一个float * ，但是你的函数H期待 double * 。只需更改t的类型即可匹配:

double t = 4.14937*(m-1);