Problems VPATH, pow

[Davies, Daniel]

Jeff,

Thanks for your very quick response!

Thanks also for pushing me into making a simple test case.  I discovered
along the way that there was a tangle of dependencies between my RTOS
and newlib that resulted in using the RTOS's buggy version of modf.

When that was straightened out, I discovered that pow(0,2) works fine,
but pow(0,0.5) fails.  That's what had actually triggered my first bug,
but I generalized the problem incorrectly based on my mis-reading of the
pow source code.  Sorry!

#include <stdio.h>
#include <math.h>

int
main(int argc, char **argv)
{
    printf("pow (1,2) is %f\n", pow(1,2));

    printf("pow (0,2) is %f\n", pow(0,2));

    printf("pow (0,1) is %f\n", pow(0,1));

    printf("pow (0,0.5) is %f\n", pow(0,0.5));

    printf("pow (0,0.1) is %f\n", pow(0,0.1));

    printf("pow (4,0.5) is %f\n", pow(4,0.5));

    printf("pow (2,2) is %f\n", pow(2,2));

    return 0;
}

Produces the following output on my system (powerpc-unknown-eabialtivec,
an embedded board.  The boards has a 7447A.  The code is compiled with
gcc-4.0.1 using hard floating point).

pow (1,2) is 1.000000
pow (0,2) is 0.000000
pow (0,1) is 0.000000
pow (0,0.5) is 1.000000		<wrong>
pow (0,0.1) is 1.000000		<wrong>
pow (4,0.5) is 2.000000
pow (2,2) is 4.000000

Dan

 

-----Original Message-----
From: Jeff Johnston [mailto:jjohnstn@redhat.com] 
Sent: Tuesday, August 30, 2005 10:11 AM
To: Davies, Daniel
Cc: newlib@sources.redhat.com
Subject: Re: Problems VPATH, pow

Davies, Daniel wrote:
> Sorry for the subject-less message - I hit the send button too soon. 
> 
> -----Original Message-----
> From: Davies, Daniel
> Sent: Monday, August 29, 2005 1:57 PM
> To: 'newlib@sources.redhat.com'
> Cc: Davies, Daniel
> Subject: 
> 
> First, let me thank everyone involved with newlib, it has certainly 
> made my life easier!
>  
> I'm cross compiling newlib 1.13.0 for powerpc-unknown-eabialtivec on 
> i386-pc-solaris2.10.  I'm using version 3.80 of gnu make and gcc
4.0.1.
>  
> Two nits:
>  
> 1. This is probably a solairs-ism.  The VPATH declaration in 
> newlib-1.13.0/libgloss/rs6000/Makefile.in is
>  
> VPATH = @srcdir@ @srcdir@/..
>  
> This does not work in Solairs.
>  
> VPATH = @srcdir@:@srcdir@/..
>  
> (replaced space with colon) works.
>  

I have changed this over to use colon as gnu make accepts either space
or colon.

> 2. pow(0, y) = 1 for all y > 0, e.g. 0^2 = 1.  It looks like pow 
> doesn't check for this case, takes the log of 0 (which is undefined), 
> gets a number very close to 0 instead of an approximation to minus 
> infinity, later takes exp of this number and gets 1.
>

Can you provide a test that demonstrates your problem?  I have a simple
test which when run on mn10300 and linux, works as expected (result of
pow (0,2) is 0).

#include <stdio.h>
#include <math.h>

int main (void)
{
   printf ("pow (1,2) is %g\n", pow (1, 2));
   printf ("pow (0,2) is %g\n", pow (0, 2));
   printf ("pow (2,2) is %g\n", pow (2, 2));
   return 0;
}
~
bash-3.00$ ./a.out
pow (1,2) is 1
pow (0,2) is 0
pow (2,2) is 4


> Dan