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.