Re: What's the deal with C99?

Dann Corbit wrote:
3. If 128 bit platforms are also included, it would mean yet another version. Of course, the Cephes collection handles all of this nicely. I think it would be good to make the Cepehs collection a "Reference" implementation which can be selected or ignored because it is so fully complete.


Cephes is very good. Lcc-win uses Cephes internally. I have rewritten
the 352 bit part of it in assembler, extended some functions, and
integrated them into the package.

[snip]


In C++ there is a nice work around, since I can use QD:
http://crd.lbl.gov/~dhbailey/mpdist/
to create my own double-double and quad-double libraries to get lots more precision if I want it but using C it would be a real pain. Without a hardware type, intermediate precisions like that become an ugly mess of library calls instead of the simple, fundamental algorithms that they ought to be.

Not in lcc-win since lcc-win accepts operator overloading,
precisely because of the points you make above.

--
jacob navia
jacob at jacob point remcomp point fr
logiciels/informatique
http://www.cs.virginia.edu/~lcc-win32
.

Relevant Pages

  • Re: Gauss Hypergeometric Function 2F1(a,b;c;x)
    ... I have programmed the Gauss Hypergeometric Function 2F1for all possible real arguments and would like to check it with ... some benchmarks or a particularly challenging set of parameters. ... Results should be in extended or double precision format. ... jacob at jacob point remcomp point fr ...
    (sci.math.num-analysis)
  • Re: Multi precision floating point
    ... to go beyond single precision, ... A whole three minutes to advertise your compiler, jacob? ... If, on the other hand, you want to write multiple-precision floating-point as a learning exercise, or you forsee ever needing to use your code on an implementation other than jacob's, or you simply want to keep the choice of switching vendors, then it is perfectly possible to write multiple-precision floating-point in Standard C. ... That is because those are libraries, whereas qfloat is a language extension, which adds to the C syntax - something libraries cannot do. ...
    (comp.lang.c)
  • Re: Trouble with integer floating point conversion
    ... For my purposes, however, it will be better to choose the constants such ... Actually for the most part of it calculation with higher precision are very ... jacob at jacob point remcomp point fr ...
    (comp.lang.c)
  • Re: Calculating log gamma
    ... of Bernoulli numbers, and I can't get any better precision ... -- John Harper, School of Mathematics, Statistics and Computer Science, Victoria University, PO Box 600, Wellington 6140, New Zealand ... I obtain 132 digits precision now. ... jacob navia ...
    (sci.math.num-analysis)
  • Re: Trouble with Cody&Waite exponentition x**y
    ... Have you looked at the powimplementations in the Cephes collection ... for the precision that you are after? ... Do they suffer from the same defect? ...
    (sci.math.num-analysis)