好久没用C了,居然报个奇怪的问题....急用,求各位高手帮忙

[iamltd]

数据是这样的
float ans[j][mdu+2];当然实际不会这样定义的哈.我拿malloc分配的
float te;

这样处理了一下以后
*(ans + i*(mdu+2) + j) -= *(ans + i*(mdu+2) + aft) * *(ans + be*(mdu+2) + j) / te;
想把太小的数处理成0
if(fabs(*(ans + i*(mdu+2) + j))<1e-6)*(ans + i*(mdu+2) + j)=0;
结果运行到这句话老报错.Floating point error:Domain.

到底是什么错啊??


急用,求各位高手帮忙

[aa1298]

加一句试试：float *p = ans;

[coyprightbao]

mark!

[dudu妈]

up

[ocpbase]

if(fabs(*(ans + i*(mdu+2) + j))<1e-6)*(ans + i*(mdu+2) + j)=0;

少了括号?
赋值号本意是等价号?没说清楚.

[积木]

我看是最后的那个=号闹得

[qingyuan18]

加一句试试：
float *p = ans;

[pacman2000]

看得莫名奇妙。。。把if的条件改清楚一点吧。

[iamltd]

if(fabs(*(ans + i*(mdu+2) + j))<1e-6) *(ans + i*(mdu+2) + j)=0;
后面那个*(ans + i*(mdu+2) + j)确实是想访问ans[i][j]的值.因为我没有声明成2维数组,用ans[i][j]会报错.

IF的意思是,如果*(ans + i*(mdu+2) + j)接近0,就把它改成0.

各位大大帮帮忙啊

[iamltd]

to xunknown(爱你是我一生中理想)

访问一个2维数组,a[i][j]和*(a+i*维数+j)该是一样的吧?

[pacman2000]

在if之前把这个*(ans + i*(mdu+2) + j)显示出来看看，究竟是多少。

[eduhf_123]

将if(fabs(*(ans + i*(mdu+2) + j))<1e-6) *(ans + i*(mdu+2) + j)=0;改为如下：
(int)temp=fabs( *(ans+i*(mdu+2)+j) )<1e-6 ? ( *(ans + i*(mdu+2) + j)=0 ):0;
试试看……

[Flood1984]

Floating point error:Domain.
除数为0

[xunknown]

二维数组指针访问的方式：
int two[N][M];
访问元素two[i][j]可以使用指针在方式访问，其中 0<= i <N,0<=j<M.
two[i][j] == *( two + i*j );

[xunknown]

上面搞错了。
应该是：
two[i][j] = *( two + i*N + j );

[llmingjie]

if(fabs(*(ans + i*(mdu+2) + j))<1e-6) * (ans + i*(mdu+2) + j) = 0;

这句话怎么理解？会不会是问题

[Leaveye]

这很可能与运行时的值有关。或者发赋值代码出来，或者在这方面仔细查一下。

[xunknown]

*(ans + i*(mdu+2) + j) ==???????????
你这个是什么意思啊？？？
你是想取ans[i][mdu+2]这个元素的值吗？
为什么要+j ?

[xunknown]

if(fabs(*(ans + i*(mdu+2) + j))<1e-6)*(ans + i*(mdu+2) + j)=0;
而且你这句话本来就是错误的！
简单点就是：
if ( fabs(K) < zero ) * (A) = 0

其中
K = *(ans + i*(mdu+2) + j);
zero = 1e-6;
A = (ans + i*(mdu+2) + j);

然后再简单点：
if M*A = 0
其中：
M = fabs( K ) < zero ;

////

是不是你的表达式有错误啊？

[maddevil]

float ans[j][mdu+2];当然实际不会这样定义的哈.我拿malloc分配的
*(ans + i*(mdu+2) + j)
>> ans 是什么类型??

[iamltd]

多谢各位大大哈.终于知道哪儿错了
我没有include <math.h>
居然编译不报错!!!

汗.....