ryuon
並列化(OpenMP)。
ev-dhのOpenMP版:
さてさて、こいつは reduction を使えばいいのだな (と昨日の勉強から)、 とやってみるが、文句を言われる。
調べてみたら、 reduction に配列は使えない(Cの場合)だと……ガックシ。
API version 3.0 (PDF) でもそう書いてあった。 p.99 lines 32-33:Aggregate types (including arrays), pointer types and reference types may not appear in a reduction clause.
ここに OpenMP を使うには、したがって、 (shared で済むように)多少冗長に書き直す必要があるが、 最終的な効率はどんなもんなのかなあ。
10/30/2008: 素朴な高速化に戻る。