Return to b7.f CVS log

Up to [Power 6/32 Unix Tahoe 4.2BSD] / cci / usr / src / usr.bin / f77 / f77pass1 / tst

Power 6/32 Unix version 1.21

c	benchmark program number 7

c	statistics

	integer array(10)
	data array/500, 1000, 2000, 4000, 100, 200, 400, 800, 1600, 3200/

	do 1 i = 1, 30000
1	call stat(array,10)
	end

	subroutine stat(arr,nel)
	integer    arr(nel)
	integer sum,mean
	integer gtmean

	sum = 0
	do 2 i=1,nel
2	sum = sum + arr(i)
	mean = gtmean(arr,nel,sum)
	call getmed(arr,nel,sum)
	end

	integer function gtmean(arr,nel,sum)
	integer    arr(nel),sum

	mean = 0
	do 3 i=1,nel
3	mean = mean + arr(i) * i
	gtmean = mean / sum
	return
	end

	subroutine getmed(arr,nel,sum)
	integer    arr(nel),sum
	integer	partsm

	partsm = 0
	do 4 i=1,nel
	partsm = partsm + arr(i)
	if(partsm .gt. rshift(sum,2))goto 5
4	continue
5	return
	end

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.