Question

大家好，并提前致谢。我试图用高斯 - 勒让德正交的C语言编写一个函数 - 积分的数值方法 - 我必须将它应用于2个案例，其中一个区间被转换为[-1,1]，另一个案例的区间又被分割分四部分。由于两者都使用相同的表格数字作为权重和横坐标，我已经为两者制作了数组，它们之间的唯一区别是分割（或重新划分）间隔的功能。问题是，对于每个我将在我想要集成的函数（在erf的情况下）和不同的乘数上有不同的应用程序（我不会详细说明，这是该方法的内在数学部分）我有一个功能：

void adapt(double a, double b, double *ptrt, double *ptrm, double x)
{
   *ptrt=1/2*(x*(b-a)+a+b);
   *ptrm=(b-a)/2;
}

其中a和b代表beggingin和原始区间的结尾，x代表x的tabled gauss-legendre值，指针将返回我需要在正交函数的公式中使用的数字：< / p>

double gauss(double a, double b)
{
int i;
double I, t, m,  p[16],  x[16];
double *ptrt, *ptrm;
ptrt=&t;
ptrm=&m;
I=t=m=0.0;
x[1]=0.0000000000000000; x[2]=-0.2011940939974345; x[3]=0.2011940939974345;
x[4]=-0.3941513470775634; x[5]=0.3941513470775634; x[6]=-0.5709721726085388;
x[7]=0.5709721726085388; x[8]=-0.7244177313601701; x[9]=0.7244177313601701;
x[10]=-0.8482065834104272; x[11]=0.8482065834104272; x[12]=-0.9372733924007060 ;
x[13]=0.9372733924007060; x[14]=-0.9879925180204854; x[15]=0.9879925180204854 ;

p[1]=0.2025782419255613; p[2]=0.1984314853271116; p[3]=0.1984314853271116;
p[4]=0.1861610000155622; p[5]=0.1861610000155622; p[6]=0.1662692058169939;
p[7]=0.1662692058169939; p[8]=0.1395706779261543; p[9]=0.1395706779261543;
p[10]=0.1071592204671719; p[11]=0.1071592204671719; p[12]=0.0703660474881081;
p[13]=0.0703660474881081; p[14]=0.0307532419961173; p[15]=0.0307532419961173;

for(i=1;i<=15;i++)
    {
            adapt(a, b, &ptrt, &ptrm, x[i]);
            if(a=0.0 && b=XF) I=I+C*m*(p[i]*dfunc(t));
            else I=I+C*m*XF/2*(p[i]*dfunc(XF/2*(t+1));

    }
return I;
}

我不会在这里发布主要功能，因为它只是打印并说出不重要的回报。 我的问题它不断发出消息：

“EP3.c：140：4：警告：从不兼容的指针类型传递'adapt'的参数3 [默认启用] 适应（a，b，＆amp; ptrt，＆amp; ptrm，x [i]）;

EP3.c：16：6：注意：预期'double *'但参数类型为'double **' void adapt（double a，double b，double * ptrt，double * ptrm，double x）“

和另一个指针相同。 我不明白为什么编译器理解这个指针是“指针指针”

我已经通过代码看了几十次，无法找到原因。对不起，第一次来这里=）

Answer 1

ptrm指向m(double)的地址。因此ptrm属于double *类型。在函数调用中，您传递ptrm的地址，使其变为double **类型。因此，在函数调用中只传递pointers。

 adapt(a, b, ptrt, ptrm, x[i]);

预期'double *'但参数类型为'double **'

1 个答案: