[parisc-linux] Bug 72 ( Floating point problems )

[David Huggins-Daines]

Bruno IGLESIAS <iglesiab@esiee.fr> writes:

> 	The test of the Seti@client reveals lots of problems. When the       
> calculation doesn't result in an illegal instruction (cf A500), the   
> results are false. The kernel running was compiled on a B132L and the 
> xc seems to be an old one. We now try with an newer kernel and xc.    

That's likely to be either:

 1) a compiler bug.
 2) a problem with saving/restoring floating point state in the kernel.
 3) a consequence of the mostly non-existant floating point completion
    handler in the kernel. (see FIXME in arch/parisc/kernel/traps.c)

...in increasing order of likeliness.

A whole bunch of instructions are supposed to be trapped and emulated
in software, most notably the rounding instructions.  We catch
co-processor traps and handle exceptions, but we don't (yet) emulate
anything, so the results of these instructions are undefined.

Part of the problem is that the details of which instructions are to
be emulated are not in the architecture manuals, but rather in the
chip manuals, and we only have these for the 7100 and 7300LC, as far
as I know.

What probably needs to be done is someone needs to take small pieces
of numeric code and track down all the operations that are throwing
traps or giving incorrect results, and fix them all individually :(

> 	The calculation time for the test is 56 hours. Under HPUX, it is     
> about 30 hours.                                                       

That's probably to be expected for floating-point code.  HP/UX has
lots of extra, super-optimized millicode to handle various low-level
FP operations.

-- 
David Huggins-Daines              -                     dhd@pobox.com
                      http://www.pobox.com/~dhd/