推 yoco315:Leon 大大駕臨 <(_ _)> 08/15 19:34
※ 引述《Lepton (輕子)》之銘言：
: 我這邊做的事情基本上就是迭代計算
: vector<vector<double> > V, VNew ,rho;
: void Jacobi()
: {
: #pragma omp parallel for
: for (int i = 1; i <= L; i++)
: #pragma omp parallel for
: for (int j = 1; j <= L; j++)
: VNew[i][j] = 0.25 * (V[i - 1][j] + V[i + 1][j] +
: V[i][j - 1] + V[i][j + 1] +
: h * h * rho[i][j]);
: }
....
在最佳化之前, 我習慣是去 check algorithm 上面的解法.
你的問題是 parallel Jacobi,
這個問題非常古老, 有很多 multi-threading 的方式去做.
www.cfm.brown.edu/people/gk/AM258/Handouts/Web_paraJacobi.ppt
這個問題應該是可以找到 Cuda 的範例.
--
趙客縵胡纓，吾鉤霜雪明。銀鞍照白馬，颯沓如流星。
十步殺一人，千里不留行。是了拂衣去，深藏身與名。
閑過信陵飲，脫劍膝前橫。將炙啖朱亥，持觴勸侯贏。
三杯吐然諾，五嶽倒為輕。眼花耳熱後，意氣素霓生。
就趙揮金錘，邯鄲先震驚。千秋二壯士，烜赫大梁城。
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 96.41.13.29