TR1 Math

[Paolo Carlini]

Paolo Carlini wrote:

> Hi Ed,
>
>> I have to accept that the TR1 math won't be done for 4.2.  However, I 
>> have enough done that I think it would be worthwhile to report on the 
>> status and get some feedback and suggestions and maybe some help. 
>
> I'm looking into the patch draft. A few immediate, rather obvious 
> comments for the next iteration.

A couple of additional sparse comments:
10- M_PI is a *double* constant, as such not suited for long double 
computations (the glibc math.h header is pretty clear about that)
11- At this point, certainly I'm not going to discuss the merit of each 
specific algorithm (maybe a little bit after some bibliographic 
citations...), but it seems to me that all those replicated __max_iter 
cannot be the final solution: at minimum the specific values should be 
discussed a bit and changed to commonized macros, I think.

Paolo.