So, I'm trying to write some macros that use pow(X, V).
I'm trying to make it past the "domain error" when pow(0,0) occurs.

So I tried:
#local PowX = select (X * V * (1-X), 0, 0, pow(X, V) );
But that gave me the same error, presumably because select evaluates the full
expression even though the threshold switch won't implement the >0 result.

So I finally did the following:

#declare E = 0.0000001;
#local PowX = select ((X * V * (1-X))-E, 0, 0, pow(X+E, V) );

#local Pow1minX = select ((1-X)*(N-V), 0, 0, pow((1-X),(N-V)) );

At some point X=1, and so 1-1=0, and you'd expect to get the same domain error,
but that doesn't happen.

Is there a better, more elegant way to do this?

Is this a floating point thing again - out at the 15th decimal place?

Isn't there a pow function that tolerates / passes through pow(0,0) as 0, like
the magical atan2 ?

I must say that it get get REALLY frustrating to get bogged down in fixing
things that you wouldn't have thought needed fixing, just so you can get back to
the point of working on the ACTUAL math that you started trying to fix.

Post a reply to this message