如何把这段代码改成OMP多线程代码？谢谢

[situjunc]

额，有一个三维数组h[2][100][100]，其实就是两个二维数组h[0][99][99]和h[1][99][99]。
其中h[0][49][49]=100,h[1][49][49]=100。剩下的所有数据初始化为20。执行下面的迭代代码后，就会数据就会改变。其实就是热量散布问题啦.....
现在要怎么把这段顺序改成#pragma omp parallel 这种并行代码呢？谢谢

       for (iteration = 0; iteration < 999; iteration++) {
		for (i = 1; i < 99; i++){
			 for (j = 1; j < 99; j++){
			 	h[next][i][j] = 0.25 * (h[current][i-1][j] + h[current][i+1][j] + h[current][i][j-1] + h[current][i][j+1]);
			 }
		}
		current = next; // swap values of current and next
		next = 1 - current;
	}