In my project, we had to migrate a few old Fortran codes from one server to another server. in that process, i am able to compile the Fortran code and create the executable file. but, i get a small error in the in the floating point numbers. I have tried to compile the codes in different ways. But am not able to succeed.
below is a small piece of the comparison between the outputs from both the servers.(just look at the last fp number)
Compare: (<) new_server_swtst.3o (115221 bytes)
with : (>) orig_server_swtst.3o (115221 bytes)
< ag 0.40934E-04 0.3708229293E-09
> ag 0.40934E-04 0.3708229117E-09
The new server to which we are moving the Fortran code is a "Sun Fire[tm] V880 Server" system and it uses the IEEE floating point system & and the old system is a "Sun Blade[tm] 150 Workstation". This is all that i know abt the systems.
& about the software package, it is a program called EQ3/6 that does a few chemical calculations. according to the documentation, it says "Portability to machines with odd word lengths (other than 16, 32, or 64 bits), or machines which lack compilers fully consistent with the FORTRAN 77 standard, may be quite difficult. " but, both of the servers have 128 bit word lengths. somehow they have installed it correctly in the older server.
both the servers have a GNU fortran 77 compiler. & the new server has a Sun's f77,f90 & f95 compiler.
I would be glad if anyone can help me out with this problem, as i am new to both Fortran and Unix server stuff.
(Software Developer, ASU SF)
Fortran Floating point problem
Thu, 01/06/2005 - 00:20