|
|
Power 6/32 Unix version 1.21
c benchmark program number 7 c statistics common /a/ array(10) integer array c data array/500, 1000, 2000, 4000, 100, 200, 400, 800, 1600, 3200/ do 1 i = 1, 30000 1 call stat(10) end subroutine stat(nel) common /a/ arr(10) integer sum,mean,arr integer gtmean sum = 0 do 2 i=1,nel 2 sum = sum + arr(i) mean = gtmean(nel,sum) call getmed(nel,sum) end integer function gtmean(nel,sum) common /a/ arr(10) integer sum,arr mean = 0 do 3 i=1,nel 3 mean = mean + arr(i) * i gtmean = mean / sum return end subroutine getmed(nel,sum) common /a/ arr(10) integer sum,arr 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.