* G77 DOUBLE/SINGLE PRECISION PUZZLE
@ 2002-07-18 17:58 stephen.buckley
2002-07-19 5:04 ` Toon Moene
0 siblings, 1 reply; 3+ messages in thread
From: stephen.buckley @ 2002-07-18 17:58 UTC (permalink / raw)
To: gcc-help
Dear Helpers
I am using Force2.0.7 which is an IDE for g77 running under MS windows.
What I was doing...
I tried an algorithm in a text for solving a partial
differential equation numerically.
What happened..
running the program at work [we have Windows NT 4.0 installed]
the answer came out right on the nose--the result column gave
identical numbers to those in the text.
Thought I'd run same program at home [I have to do a lot of
work at home too] and the result column came out really crook--
nothing like the correct solution. Went over the code and
couldn't see anything wrong. Tried same program on one of my
children's computers and got the same result.
Sitting there wondering what to do next when the thought
occurred to me to try to see what happens for a small t
increment [pde is in x and t--t from 0.0 to 0.5]. Tried a small
step and solution looked OK. A little more and it started to
drift--by time I went to 60 steps [t=0.3] it was getting way off.
So I thought hang it let's see what happens if I try double
precision--and guess what--answer came in right on the nose like
at work.
My home computer [and my children's] are running windows 98.
So does NT 4.0 give more precision for g77 than windows 98?
It's not a serious problem for me as I can work around it by
using double precision at home but it is a little puzzling.
Anyone got any thoughts on this??
kind regards
Stephen
------------------------------------------------
Join Excite! - http://www.excite.com
The most personalized portal on the Web!
^ permalink raw reply [flat|nested] 3+ messages in thread
* Re: G77 DOUBLE/SINGLE PRECISION PUZZLE
2002-07-18 17:58 G77 DOUBLE/SINGLE PRECISION PUZZLE stephen.buckley
@ 2002-07-19 5:04 ` Toon Moene
0 siblings, 0 replies; 3+ messages in thread
From: Toon Moene @ 2002-07-19 5:04 UTC (permalink / raw)
To: stephen.buckley; +Cc: gcc-help
stephen.buckley@excite.com wrote:
>
> Dear Helpers
>
> I am using Force2.0.7 which is an IDE for g77 running under MS windows.
> What I was doing...
> I tried an algorithm in a text for solving a partial
> differential equation numerically.
> What happened..
> running the program at work [we have Windows NT 4.0 installed]
> the answer came out right on the nose--the result column gave
> identical numbers to those in the text.
> Thought I'd run same program at home [I have to do a lot of
> work at home too] and the result column came out really crook--
> nothing like the correct solution. Went over the code and
> couldn't see anything wrong. Tried same program on one of my
> children's computers and got the same result.
> Sitting there wondering what to do next when the thought
> occurred to me to try to see what happens for a small t
> increment [pde is in x and t--t from 0.0 to 0.5]. Tried a small
> step and solution looked OK. A little more and it started to
> drift--by time I went to 60 steps [t=0.3] it was getting way off.
> So I thought hang it let's see what happens if I try double
> precision--and guess what--answer came in right on the nose like
> at work.
> My home computer [and my children's] are running windows 98.
> So does NT 4.0 give more precision for g77 than windows 98?
>
> It's not a serious problem for me as I can work around it by
> using double precision at home but it is a little puzzling.
>
> Anyone got any thoughts on this??
No, not really ...
Is it an option to send this program to me so that I could try it on Linux
?
Thanks in advance,
--
Toon Moene - mailto:toon@moene.indiv.nluug.nl - phoneto: +31 346 214290
Saturnushof 14, 3738 XG Maartensdijk, The Netherlands
Maintainer, GNU Fortran 77: http://gcc.gnu.org/onlinedocs/g77_news.html
Join GNU Fortran 95: http://g95.sourceforge.net/ (under construction)
^ permalink raw reply [flat|nested] 3+ messages in thread
* Re: G77 DOUBLE/SINGLE PRECISION PUZZLE
@ 2002-07-19 6:17 Bjorn R. Jensen
0 siblings, 0 replies; 3+ messages in thread
From: Bjorn R. Jensen @ 2002-07-19 6:17 UTC (permalink / raw)
To: gcc-help
Hi there,
I think, the problem is merely a symptom of using too low a precision
solving your differential equation. I doubt, there is significant
difference in the precision of floats on various windows platforms.
Your
differential equation most likely has degenerate or near degenerate
regions in its
Hilbert space, so where you end up would be fairly random, if the
precision
is too low.
You could try to map out the surfaces of the relevant eigenvalues in
the regions,
where the calculations start to diverge.
Yours sincerely,
Bjorn
^ permalink raw reply [flat|nested] 3+ messages in thread
end of thread, other threads:[~2002-07-19 13:17 UTC | newest]
Thread overview: 3+ messages (download: mbox.gz / follow: Atom feed)
-- links below jump to the message on this page --
2002-07-18 17:58 G77 DOUBLE/SINGLE PRECISION PUZZLE stephen.buckley
2002-07-19 5:04 ` Toon Moene
2002-07-19 6:17 Bjorn R. Jensen
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).