I'll do it in a few days if noone else does, but right now I'm out of town for Thanksgiving and don't have access to my AMD computer.
On Thu, Nov 25, 2010 at 10:06 AM, Don Clugston <[email protected]>wrote: > With regard to the AMD vs Intel bug: > > The previous results narrowed it down significantly, but I'm going to > need more info to fix this. > Please, somebody with AMD add this to line to > std.math.internal.gammafunction line 624, add an import > core.stdc.stdio; in the file, > > if( x == 0.0L ) { > di = 0.5; > x = x0 + di * (x1 - x0); > if( x == 0.0 ) > goto under; > } > y = betaIncomplete( a, b, x ); > + printf("%d %La %La %La %La %La %La %La\n", i, x, y, di, > yl, yh, a, b); > yp = (x1 - x0)/(x1 + x0); > if( fabs(yp) < dithresh ) > goto newt; > yp = (y-y0)/y0; > if( fabs(yp) < dithresh ) > goto newt; > > compile Phobos, and then run this program: > > import std.mathspecial; > import std.stdio; > > void main() > { > writefln("%a", betaIncompleteInverse(0x1.ff1275ae5b939bcap-41, 4.6713e18, > 0.0813601)); > } > _______________________________________________ > phobos mailing list > [email protected] > http://lists.puremagic.com/mailman/listinfo/phobos >
_______________________________________________ phobos mailing list [email protected] http://lists.puremagic.com/mailman/listinfo/phobos
