On Tue, 7 Mar 2006, Patrick Wallace wrote:
> With Myles Standish of JPL I'm trying to find out why the TDB-TT routine
> they're using gives distinctly different answers on different platforms.
> I suspect it isn't just rounding but something more pathological. If any
> of you can compile this program on any Fortrans you have and let me know
> the number, it would be immensely useful.
I get -65.497099 for three different setups.
On my PowerBook G4 running OS X 10.4.5 with the G95 compiler:
[pippin:~] bradc% uname -a
Darwin pippin.local 8.5.0 Darwin Kernel Version 8.5.0: Sun Jan 22 10:38:46
PST 2006; root:xnu-792.6.61.obj~1/RELEASE_PPC Power Macintosh powerpc
[pippin:~] bradc% g95 -v -o bum bum.f
Driving: g95 -v -o bum bum.f -lf95 -lm
Using built-in specs.
Target:
Configured with: /Users/andy/g95/osx/gcc.osx/configure
--enable-languages=c
Thread model: posix
gcc version 4.0.1 (g95!) Feb 27 2006
/local-star/bin/../lib/gcc-lib/powerpc-apple-darwin6.8/4.0.1//f951 bum.f
-fPIC -quiet -dumpbase bum.f -auxbase bum -version -o /var/tmp//ccGmJlfH.s
G95 Fortran 95 version 4.0.1 (g95!) Feb 27 2006 (powerpc-apple-darwin6.8)
compiled by GNU C version 3.1 20020420 (prerelease).
GGC heuristics: --param ggc-min-expand=73 --param ggc-min-heapsize=81920
as -arch ppc -o /var/tmp//cc6hvj1a.o /var/tmp//ccGmJlfH.s
ld -dynamic -arch ppc -weak_reference_mismatches non-weak -o bum -lcrt1.o
/local-star/bin/../lib/gcc-lib/powerpc-apple-darwin6.8/4.0.1//crt2.o
-L/local-star/bin/../lib/gcc-lib/powerpc-apple-darwin6.8/4.0.1/
-L/local-star/bin/../lib/gcc-lib/powerpc-apple-darwin6.8/4.0.1
-L/usr/lib/gcc//
-L/local-star/bin/../lib/gcc-lib/powerpc-apple-darwin6.8/4.0.1///
-L/local-star/bin/../lib/gcc-lib/powerpc-apple-darwin6.8/4.0.1/
-L/usr/lib// /var/tmp//cc6hvj1a.o -lf95 -lm -lgcc -lgcc_eh -lSystem
[pippin:~] bradc% ./bum
FB2DT -65.497099
On my Linux PC (Intel Pentium 4) running WBEL4 with the f77 compiler:
akule % uname -a
Linux akule 2.6.9-5.0.5.ELsmp #1 SMP Thu May 19 19:24:44 CDT 2005 i686
i686 i386 GNU/Linux
akule % f77 -v -o bum bum.f
Driving: f77 -v -o bum bum.f -lfrtbegin -lg2c -lm -shared-libgcc
Reading specs from /usr/lib/gcc/i386-redhat-linux/3.4.3/specs
Configured with: ../configure --prefix=/usr --mandir=/usr/share/man
--infodir=/usr/share/info --enable-shared --enable-threads=posix
--disable-checking --with-system-zlib --enable-__cxa_atexit
--disable-libunwind-exceptions --enable-java-awt=gtk
--host=i386-redhat-linux
Thread model: posix
gcc version 3.4.3 20041212 (Red Hat 3.4.3-9.EL4)
/usr/libexec/gcc/i386-redhat-linux/3.4.3/f771 bum.f -quiet -dumpbase
bum.f -auxbase bum -version -o /tmp/ccBT3EjF.s
GNU F77 version 3.4.3 20041212 (Red Hat 3.4.3-9.EL4) (i386-redhat-linux)
compiled by GNU C version 3.4.3 20041212 (Red Hat 3.4.3-9.EL4).
GGC heuristics: --param ggc-min-expand=100 --param ggc-min-heapsize=131072
as -V -Qy -o /tmp/cccIREKI.o /tmp/ccBT3EjF.s
GNU assembler version 2.15.92.0.2 (i386-redhat-linux) using BFD version
2.15.92.0.2 20040927
/usr/libexec/gcc/i386-redhat-linux/3.4.3/collect2 --eh-frame-hdr -m
elf_i386 -dynamic-linker /lib/ld-linux.so.2 -o bum
/usr/lib/gcc/i386-redhat-linux/3.4.3/../../../crt1.o
/usr/lib/gcc/i386-redhat-linux/3.4.3/../../../crti.o
/usr/lib/gcc/i386-redhat-linux/3.4.3/crtbegin.o
-L/usr/lib/gcc/i386-redhat-linux/3.4.3
-L/usr/lib/gcc/i386-redhat-linux/3.4.3
-L/usr/lib/gcc/i386-redhat-linux/3.4.3/../../.. /tmp/cccIREKI.o -lfrtbegin
-lg2c -lm -lgcc_s -lgcc -lc -lgcc_s -lgcc
/usr/lib/gcc/i386-redhat-linux/3.4.3/crtend.o
/usr/lib/gcc/i386-redhat-linux/3.4.3/../../../crtn.o
akule % ./bum
FB2DT -65.497099
Same computer using g95:
akule % g95 -v -o bum95 bum.f
Driving: g95 -v -o bum95 bum.f -lf95 -lm
Using built-in specs.
Target:
Configured with: /home/andy/g95.build/linux/gcc.linux/configure
--enable-languages=c
Thread model: posix
gcc version 4.0.1 (g95!) Sep 5 2005
/local/progs/g95-install/bin/../lib/gcc-lib/i686-pc-linux-gnu/4.0.1//f951
bum.f -quiet -dumpbase bum.f -mtune=pentiumpro -auxbase bum -version -o
/tmp/ccGRk5pd.s
G95 Fortran 95 version 4.0.1 (g95!) Sep 5 2005 (i686-pc-linux-gnu)
compiled by GNU C version 2.96 20000731 (Red Hat Linux 7.1
2.96-98).
GGC heuristics: --param ggc-min-expand=100 --param ggc-min-heapsize=131072
as --traditional-format -V -Qy -o /tmp/ccbQinph.o /tmp/ccGRk5pd.s
GNU assembler version 2.15.92.0.2 (i386-redhat-linux) using BFD version
2.15.92.0.2 20040927
ld -m elf_i386 -dynamic-linker /lib/ld-linux.so.2 -o bum95
/usr/lib///crt1.o /usr/lib///crti.o
/local/progs/g95-install/bin/../lib/gcc-lib/i686-pc-linux-gnu/4.0.1//crtbegin.o
-L/local/progs/g95-install/bin/../lib/gcc-lib/i686-pc-linux-gnu/4.0.1/
-L/local/progs/g95-install/bin/../lib/gcc-lib/i686-pc-linux-gnu/4.0.1
-L/usr/lib/gcc// -L/lib//
-L/local/progs/g95-install/bin/../lib/gcc-lib/i686-pc-linux-gnu/4.0.1///
-L/local/progs/g95-install/bin/../lib/gcc-lib/i686-pc-linux-gnu/4.0.1/
-L/lib// -L/usr/lib// /tmp/ccbQinph.o -lf95 -lm -lgcc -lc -lgcc
/local/progs/g95-install/bin/../lib/gcc-lib/i686-pc-linux-gnu/4.0.1//crtend.o
/usr/lib///crtn.o
akule % ./bum95
FB2DT -65.497099
|