大家好,谢谢大家。我试图用C写一个函数,用于Gauss求积--积分的数值方法--我必须将它应用于2种情况,一种是将区间转换为-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代表乞讨和原始间隔的结尾,x代表提交的高斯-legendre值,指针将返回我在求积函数的公式中需要使用的数字:
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:警告:从默认启用的不兼容指针类型(a,b,&ptrt,&ptrm,xi)传递‘adapt’的参数3;
EP3.c:16:6:注:预期的‘双重*’但参数是‘双**’类型的无效适应(双a,双b,双*ptrt,双*ptrm,双x)“
另一个指针也是一样。我不明白为什么编译器会理解这个指针是“指向指针的指针”。
我看了几十遍这些代码,找不到原因。对不起,我第一次犯了错误
发布于 2015-06-15 15:57:58
ptrm指向m(double)的地址。因此,ptrm属于double *型。在函数调用中,您将传递ptrm的地址,因此它将成为double **类型。因此,在您的函数调用中,只传递pointers。
adapt(a, b, ptrt, ptrm, x[i]);https://stackoverflow.com/questions/30849632
复制相似问题