Проблемы параллелизма с gccgo
Я только что приземлился в GO-сфере, поэтому, как обычно, я пытаюсь реализовать некоторые классики на новом языке, чтобы провести некоторое сравнение.
Я испытываю странную ситуацию с симуляцией Монте-Карло для PI, используя параллелизм:
- когда я компилирую код, используя
gc
для любого количества событий (nevents
) Я обеспечил одновременное выполнение всех подпрограмм и все работает как положено (даже если "параллелизм не параллелизм", я знаю) - когда я компилирую код, используя
gccgo
работает как указано выше, когдаnevents
меньше или равно ~10.000.000, но для больших значений каждая процедура запускается, когда предыдущая почти завершена, и это поведение тем более очевидно, чем большеnvalues
получает (~100.000.000)
Что мне не хватает? Будь хорошим, это моя первая неделя с Go!
Вот код:
package main
import ("os"
"fmt"
"log"
"strconv"
"math"
"math/rand"
"time"
)
var (ncpus int = 2
threadspercpu int = 2
ngophers int = ncpus*threadspercpu
rayu float64 = 1.0
)
func montecarlo(returnvalues chan int, gevents int, id int) {
// debug
fmt.Println("goroutine started")
randgen := rand.New(rand.NewSource(time.Now().UnixNano() +
int64(id*id + 1)))
insidecircle := 0
for i := 0; i < gevents; i++ {
if (math.Pow(randgen.Float64(), 2.0) +
math.Pow(randgen.Float64(), 2.0)) < rayu {
insidecircle++
}
}
// debug {
time.Sleep(time.Duration(id)*time.Second)
fmt.Println("goroutine finished")
// }
returnvalues <- insidecircle
}
func main() {
arg := os.Args[1]
nevents64, err := strconv.ParseInt(arg, 10, 32)
if (err != nil || nevents64 < 1) {
log.Fatalf("Invalid number of cases provided: '%s'", arg)
}
nevents := int(nevents64)
returnvalues := make(chan int)
gopherevents := nevents / ngophers
for i := 0; i < ngophers-1; i++ {
go montecarlo(returnvalues, gopherevents, i)
}
remaining := gopherevents + (nevents % ngophers)
go montecarlo(returnvalues, remaining, ngophers-1)
insidecircle_events := 0
for i := 0; i < ngophers; i++ {
insidecircle_events += <- returnvalues
}
fmt.Println("pi: ", (float64(insidecircle_events) /
float64(nevents) * 4))
}
Я добавил задержку для облегчения отладки.