MShaw

January 17th, 2010, 04:00 PM

Could some please compile this and tell me the last value of n and i. I want to see if it works.

program sieve_t

implicit none

integer*1 s(1000000), offset (10), sequence

integer i, j, n

n=0

sequence=4

do i=1, 10

offset(i)=0

enddo

offset(1)=2

offset(2)=4

offset(3)=2

do i=1, 1000000

s(i) = 1

enddo

do i=2, 1000000

if (s(i).eq.1) then

do j=2, (1000000/i)

s(i*j)=0

enddo

endif

enddo

do i=2, 1000000

if (s(i).eq.1) then

do j=1,sequence-1

if (s(i+offset(j)).ne.1) goto 10

enddo

n=n+1

write(*,*) n, i,i+2,i+6,i+8

endif

10 continue

enddo

end

Thanks in advance.

program sieve_t

implicit none

integer*1 s(1000000), offset (10), sequence

integer i, j, n

n=0

sequence=4

do i=1, 10

offset(i)=0

enddo

offset(1)=2

offset(2)=4

offset(3)=2

do i=1, 1000000

s(i) = 1

enddo

do i=2, 1000000

if (s(i).eq.1) then

do j=2, (1000000/i)

s(i*j)=0

enddo

endif

enddo

do i=2, 1000000

if (s(i).eq.1) then

do j=1,sequence-1

if (s(i+offset(j)).ne.1) goto 10

enddo

n=n+1

write(*,*) n, i,i+2,i+6,i+8

endif

10 continue

enddo

end

Thanks in advance.