2008-07-14 16:13:59 -06:00
|
|
|
// $G $D/$F.go && $L $F.$A && ./$A.out
|
|
|
|
|
|
|
|
// Copyright 2009 The Go Authors. All rights reserved.
|
|
|
|
// Use of this source code is governed by a BSD-style
|
|
|
|
// license that can be found in the LICENSE file.
|
|
|
|
|
2010-02-23 22:21:16 -07:00
|
|
|
// Generate primes up to 100 using channels, checking the results.
|
|
|
|
// This sieve consists of a linear chain of divisibility filters,
|
|
|
|
// equivalent to trial-dividing each n by all primes p ≤ n.
|
2008-07-14 16:13:59 -06:00
|
|
|
|
|
|
|
package main
|
|
|
|
|
|
|
|
// Send the sequence 2, 3, 4, ... to channel 'ch'.
|
2009-01-20 15:40:40 -07:00
|
|
|
func Generate(ch chan<- int) {
|
2008-07-15 21:52:07 -06:00
|
|
|
for i := 2; ; i++ {
|
2008-09-16 20:33:40 -06:00
|
|
|
ch <- i // Send 'i' to channel 'ch'.
|
2008-07-14 16:13:59 -06:00
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
// Copy the values from channel 'in' to channel 'out',
|
|
|
|
// removing those divisible by 'prime'.
|
2009-01-20 15:40:40 -07:00
|
|
|
func Filter(in <-chan int, out chan<- int, prime int) {
|
2009-03-20 12:32:58 -06:00
|
|
|
for i := range in { // Loop over values received from 'in'.
|
2008-07-14 16:13:59 -06:00
|
|
|
if i % prime != 0 {
|
2008-09-16 20:33:40 -06:00
|
|
|
out <- i // Send 'i' to channel 'out'.
|
2008-07-14 16:13:59 -06:00
|
|
|
}
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
// The prime sieve: Daisy-chain Filter processes together.
|
2009-01-20 15:40:40 -07:00
|
|
|
func Sieve(primes chan<- int) {
|
2009-01-06 16:19:02 -07:00
|
|
|
ch := make(chan int); // Create a new channel.
|
2008-07-15 21:52:07 -06:00
|
|
|
go Generate(ch); // Start Generate() as a subprocess.
|
2008-07-14 16:13:59 -06:00
|
|
|
for {
|
2009-03-20 12:32:58 -06:00
|
|
|
// Note that ch is different on each iteration.
|
2008-07-15 21:52:07 -06:00
|
|
|
prime := <-ch;
|
2008-09-16 20:33:40 -06:00
|
|
|
primes <- prime;
|
2009-01-06 16:19:02 -07:00
|
|
|
ch1 := make(chan int);
|
2008-07-14 16:13:59 -06:00
|
|
|
go Filter(ch, ch1, prime);
|
|
|
|
ch = ch1
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
func main() {
|
2009-01-06 16:19:02 -07:00
|
|
|
primes := make(chan int);
|
2008-07-14 16:13:59 -06:00
|
|
|
go Sieve(primes);
|
2009-03-03 09:39:12 -07:00
|
|
|
a := []int{2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97};
|
2008-09-16 20:33:40 -06:00
|
|
|
for i := 0; i < len(a); i++ {
|
2009-03-20 12:32:58 -06:00
|
|
|
if x := <-primes; x != a[i] { panic(x, " != ", a[i]) }
|
2008-09-16 20:33:40 -06:00
|
|
|
}
|
2008-07-14 16:13:59 -06:00
|
|
|
}
|