On 2008-08-05 09:59+0100 Andrew Ross wrote: > On Sat, Aug 02, 2008 at 08:54:07AM -0700, Alan Irwin wrote: >> On 2008-08-02 10:09+0100 Andrew Ross wrote: >> >>> C and f95 versions of example 21 give identical results here for me. >>> Either this is a 64-bit issue, or it is a difference between different >>> versions of gfortran (I'm using 4.2.3). Can you test further or narrow >>> down where the differences are? >> >> Hi Andrew: >> >> The numerical differences are limited, but where they occur the differences >> are substantial [...]
> Alan, > > A first initial look at the postscript shows differences even on the > first page, suggesting perhaps a slightly different random distribution > of numbers to start with. I've made a few of the constants in the > example explicitly of kind plflt, including the max / min values for the > range. This might make a difference. Can you check? If I can elimiate > these differences it might make a difference to the subsequent pages. If > not, at least I have narrowed down the options. Hi Andrew: Your instincts were absolutely right. As the result of your f95 changes the f95 and C results for example 21 now agree exactly! You should also make the corresponding changes in the f77 example (which currently agrees with the old f95 version done without your recent change). Alan __________________________ Alan W. Irwin Astronomical research affiliation with Department of Physics and Astronomy, University of Victoria (astrowww.phys.uvic.ca). Programming affiliations with the FreeEOS equation-of-state implementation for stellar interiors (freeeos.sf.net); PLplot scientific plotting software package (plplot.org); the libLASi project (unifont.org/lasi); the Loads of Linux Links project (loll.sf.net); and the Linux Brochure Project (lbproject.sf.net). __________________________ Linux-powered Science __________________________ ------------------------------------------------------------------------- This SF.Net email is sponsored by the Moblin Your Move Developer's challenge Build the coolest Linux based applications with Moblin SDK & win great prizes Grand prize is a trip for two to an Open Source event anywhere in the world http://moblin-contest.org/redirect.php?banner_id=100&url=/ _______________________________________________ Plplot-devel mailing list Plplot-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/plplot-devel
