Relay-Version: version B 2.10 5/3/83; site utzoo.UUCP Posting-Version: version B 2.10.1 6/24/83; site astrovax.UUCP Path: utzoo!watmath!clyde!burl!ulysses!allegra!princeton!astrovax!gam From: gam@astrovax.UUCP (Gary Mamon) Newsgroups: net.bugs.4bsd,net.lang.f77 Subject: Another F77 bug in 4.2BSD UNIX Message-ID: <238@astrovax.UUCP> Date: Thu, 15-Mar-84 17:04:18 EST Article-I.D.: astrovax.238 Posted: Thu Mar 15 17:04:18 1984 Date-Received: Fri, 16-Mar-84 03:01:21 EST Organization: Princeton Univ. Astrophysics Lines: 31 Subject: F77 exponentiation. Index: /mnt/gam/BUGS/sumpow.f Description: 4.2 F77 mishandles the following in single precision: z = x1**alpha + x2**alpha + ... + xlast**alpha where alpha is a floating point constant, and z and x sub i are real variables. F77 returns for z the value N*xlast**alpha, where N is the number of terms in the sum. Furthermore F77 mishandles the following in both single and double precision. z = (x1-y1)**alpha + (x2-y2)**alpha where again alpha is a floating point constant, and z, the x sub i and the y sub i are real (real*8 in double precision) variables. Fix: Unknown. Comment: This is one of a few 4.2 BSD F77 bugs already discovered at our site. The majority of Fortran users here are now using the 4.1BSD F77 compiler.