fortran90 - Fortran behavior of a tiny program -


i new in fortran90 (30 minutes ago...) , have program:

program example1     implicit none     real (kind=8) :: x,y,z     x = 3.d0     y = 2.d-1     z = x + y      print *, "y = ", y     print *, "x = ", x     print *, "z = ", z  end program example1

but when run with:

gfortran example1.f90 ./a.out

the output is:

y =   0.20000000000000001      x =    3.0000000000000000      z =    3.2000000000000002

why not 3.2000000000000000 ??? doing wrong? why y has 1 in last digit?? , why z has 2 in last digit?? sorry if dumb question, dont understand doing wrong...

thanks!!

there's absolutely nothing wrong program. issue has real's inability represent decimals precisely, without error. problem numbers not composed of negative powers of 2 must represented approximately. that's why there small error in 16-th decimal place. more information representation of reals take @ article on wikipedia. another great article on same subject.

if replace 0.2 0.25, problem go away, because 0.25 2 ^ -2.


Comments

Post a Comment

Popular posts from this blog

SPSS keyboard combination alters encoding -

my question not have programming practical issue on spss. accident keyword combination cannot remember because did not intend press changed encoding in spss , text in data in different language english shown small boxes. nevertheless, when kind of analysis performed different language text shown in output. does know reason? tried change encoding did not produce result. first, @ edit > options , see whether setting unicode or locale. might fix it. if have text in multiple languages, want unicode. otherwise, may have specific locale set wrong. can change set locale, need know locale want use.
Read more

Socket.connect doesn't throw exception in Android -

if put nonsense url no exception thrown , none of rest of code executed not rest of asynctask called method connects. try { socketcliente.connect(new inetsocketaddress("f", port), 2000); } catch (unknownhostexception e) { geterror("host doesn't exist"); return -1; } catch (ioexception e) { geterror("could not connect: host down"); return -1; } add catch statement. catch ( exception e ) { log.d(tag, e.getmessage(),e); toast.maketext(getapplicationcontext(), "unexpected error:" + e.getmessage(), toast.length_long).show(); } this write log , pop toast tell what's going on.
Read more

iphone - How do I keep MDScrollView from truncating my row headers and making my cells look bad? -

Image
here's screenshot show describing: and relevant code: - (nsinteger)spreadview:(mdspreadview *)aspreadview numberofcolumnsinsection:(nsinteger)section{ return self.viewmodel.matrix.width - 1; } - (nsinteger)spreadview:(mdspreadview *)aspreadview numberofrowsinsection:(nsinteger)section{ return self.viewmodel.matrix.height - 1; } -(mdspreadviewcell *)spreadview:(mdspreadview *)aspreadview cellforrowatindexpath:(mdindexpath *)rowpath forcolumnatindexpath:(mdindexpath *)columnpath{ mdspreadviewcell *cell = [aspreadview dequeuereusablecellwithidentifier:@"infocell"]; if (cell == nil){ cell = [[[mdspreadviewcell alloc] initwithstyle:mdspreadviewcellstyledefault reuseidentifier:@"infocell"] autorelease]; } cell.textlabel.text = ((dtvendor *)self.vendors[rowpath.row]).name; return cell; } -(cgfloat)spreadview:(mdspreadview *)aspreadview heightforrowheaderinsection:(nsinteger)rowsection{ return 200; } -(cgfloat)spreadview:(mdspreadview...
Read more