Displaying 1 result from an estimated 1 matches for "ixmin".
Did you mean:
i_min
1999 Oct 08
1
floor(NaN) problem fixed in massdist.c (PR#291)
...Then ix = floor(xpos) returns an unexpected result.
On linux/i386, ix become -2147483648 when xpos = NaN(0x80000000007a2).
However ix become 0 on Linux/Alpha when xpos=NaN. When xpos = NaN, it
shouldn't go into the following "if" statement. But it goes inside on
Alpha/Linux (because ixmin is 0), and set y[0] and y[1] to NaN.
Following is extracted from massdist.c
for(i=0 ; i<*nx ; i++) {
xpos = (x[i] - *xlow) / xdelta;
ix = floor(xpos);
fx = xpos - ix;
if(ixmin <= ix && ix <= ixmax) {
y[ix] += (1.0 - fx);
y[ix + 1] += fx;...