|
[Sponsors] |
May 7, 2004, 18:48 |
ARCHITECTURE/OS and round-offs
|
#1 |
Guest
Posts: n/a
|
I was using the cyclic reduction procedure for inverting directly( without iterations) the block-triagonal matrix arising out of a Poisson Equation.
I found the results to the same code are slightly different on my desktop PC running Windows NT and and on a somewhat DATED cluster of 3 UNIX-machines running Sun's Solaris... The machines are probably of CRAY make.. My desktop PC gives errors of 10 power -5 whereas the unix machine(s) give a higher error 10 power -4.. I would appreciate any idea or comments on this.. Incidentally, when I use the SUN's porprietory fortran compiler f77 and g77 on the UNIX system, I get identical results for this direct solution method. |
|
May 8, 2004, 11:39 |
Re: ARCHITECTURE/OS and round-offs
|
#2 |
Guest
Posts: n/a
|
What word length (single or double precision?) you using? Are the single and double precision words the same length on both systems?
I'd not be surprised at the difference you're seeing if you're using single (32-bit words)precision on Windows and single on the UNIX boxes. Whether using single or double, I'd expect the chips to have different rounding in their internal operations. Try double precision (real*64 ...) on both systems. |
|
May 8, 2004, 23:42 |
Re: ARCHITECTURE/OS and round-offs
|
#3 |
Guest
Posts: n/a
|
In Fortran double precision is real*8. But real*64?
|
|
May 9, 2004, 01:25 |
Re: ARCHITECTURE/OS and round-offs
|
#4 |
Guest
Posts: n/a
|
I used quadruple precision in the Sun solaris system using 'f77 -r8' ..the errors got better...not identical but of comparable magnitude with the windows output...i would check later if there are some accidental real*4 in my code
|
|
May 9, 2004, 14:06 |
Re: ARCHITECTURE/OS and round-offs
|
#5 |
Guest
Posts: n/a
|
Would you believe "poetic license?"
Or, perhaps, careless proof reading? |
|
May 10, 2004, 04:49 |
Re: ARCHITECTURE/OS and round-offs
|
#6 |
Guest
Posts: n/a
|
My comments:
Differences will always exist between different architectures, operating systems and compilers. Only worry about them if they cause a substantial change in your answers (not your error terms or residuals). If they do, your solution method is inappropriate. - Steve |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
turbulent round jet | jet | OpenFOAM | 2 | March 5, 2012 22:42 |
Two round jet modeling | shg74 | FLUENT | 0 | February 2, 2010 10:47 |
turbulent round jet | jet | OpenFOAM Post-Processing | 0 | January 3, 2010 19:28 |
Meshing round orifice on physical wall | Dora | Main CFD Forum | 0 | August 28, 2006 17:10 |
Transient behaviour of a round jet | Clifford Arnold | Main CFD Forum | 15 | November 10, 1998 16:47 |