Multi-gpu问题
域在z方向划分为num_gpus段,其中num_gpus表示可用GPU的数量,然后每个GPU相应地负责一个大小为nx×ny×(nz/num_gpus)的子域。
虽然整个域在主机端表示,但gpu只存储它们的子域。由于更新GPU k(k = 0,…,num_gpus -1)上的磁场和电场分量需要分别来自设备k+1和k-1的场值,因此我们需要在多GPU域分解边界上增加额外的nx×ny细胞来存储这些数据。
1.这些鬼细胞需要先复制到CPU上再进行GPU之间的传递,会增加时间
解决:使用Open-MP API将一个CPU线程分配给一个GPU,这样每个设备在主机上都有自己的上下文
//如果没有在第一个设备上操作
//将电场边界单元从设备复制到主机(顶部)
//在读取更新的幽灵单元之前与其他线程同步
#pragama omp barrier
//如果没有在最后一个设备上运行
//将电场幽灵单元从主机复制到设备
//调用内核,以计算磁场分量
kH<<<dimGrid, dimBlock>>>