【发布时间】:2015-06-16 00:07:37
【问题描述】:
大家好,提前致谢。 我试图在 C 中为 Gauss-Legendre 求积写一个函数 - 积分的数值方法 - 我必须将它应用于 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 代表原始间隔的开始和结束,x 代表 x 的表格 gauss-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;
}
我不会在这里发布 main 函数,因为它只是打印和说明不重要的返回内容。 我的问题是它一直在发送消息:
"EP3.c:140:4: 警告:从不兼容的指针类型传递‘adapt’的参数 3 [默认启用] 适应(a, b, &ptrt, &ptrm, x[i]);
EP3.c:16:6:注意:预期为“double *”,但参数类型为“double **” 无效适应(双a,双b,双*ptrt,双*ptrm,双x)“
和另一个指针相同。 我不明白为什么编译器会理解这个指针是“指向指针的指针”
我已经通过代码查看了数十次,但找不到原因。 如有错误,第一次来这里抱歉=)
【问题讨论】:
-
因为您使用指向 ptrt 的指针调用 adapt()。从 &ptrt 中删除 & 符号,它将是你想要的。
-
@PaulR 哎呀,你是对的。我的错,我错过了。 hayesti,你说得对。 :)
-
非常感谢,现在我明白了.......
标签: c pointers numerical-methods