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

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

[Bruno IGLESIAS]

Hi all,                                                               
                                                                      
	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.    
                                                                      
	The calculation time for the test is 56 hours. Under HPUX, it is     
about 30 hours.                                                       
                                                                      
	A link is available in our page                                      
(http://mkhppa1.esiee.fr/en/result.html) showing the comparison       
between the client and the official results.                          
                                                                      
Bye                                                                   
	Bruno IGLESIAS                                                       

Re: [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/

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

[Matthew Wilcox]

On Fri, Mar 09, 2001 at 05:43:48PM +0000, Bruno IGLESIAS wrote:
> 	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.    

i'd really like to try this client myself.  Do you have the source available
for ftp somewhere?  i can't see it on the seti site.

-- 
Revolutions do not require corporate support.