GPU加速模拟混合自旋的Ising模型.pdf

V01．13No．26 第13卷第26期2013年9月科学技术与工程 Sep．2013 Science and Technology @2013 1671—1815(2013)26—7720—06 Engineenng Sci．Tech．Engrg． 崔沛东蔡静魏亮 张伟 (暨南大学理工学院物理学系，广州510632) 摘要基于传统的Monte unifieddevice Carlo方法，利用CUDA(compute 024X1 Y。、Y。。与同时期的CPU算法相比，基于GPU的并行模拟算法大大提高了算法效率，对于尺寸为1024的混合自旋的 Ising模型，在正方晶格上实现了40倍的加速比，在三角晶格上实现了4l倍的加速比。 关键词GPU加速 混合自旋Ising模型 中图法分类号0414．2； 文献标志码A Ising模型…是统计物理中的经典模型，单一自了基于GPU的Metropolis算法¨J。 旋Ising模型的临界行为得到了深入地研究。近年 现分别使用CPU和GPU算法，模拟了二维正 来，越来越多的学者开始研究混合自旋的Ising模方晶格的混合Ising模型，验证了结果的正确性，而 型。国内的张国民、杨传章等人对混合自旋的Ising且对三角晶格上的混合Ising模型进行了探索，观察 到磁化强度，磁化率，比热随温度的变化关系，并使 模型进行了探索，推测了Metropolis算法同样可以 用4阶宾得累积量拟合出不同浓度下的临界温度， 应用于混合自旋的Ising模型心』。2009年M．Godoy 计算出了临界指数Y。、Y。。相对于同时期的CPU算 等人研究了正方晶格上混合自旋(S=1，S=1／2)的 法，基于GPU的算法大大提高了计算效率，实现了 Ising模型，其中自旋S=1和自旋S=1／2的粒子数 40多倍的加速比。 相等，并通过Binder累积量得到了相应模型的临界 温度L3 J。2012年D．S．Cambui等人用Metropolis算 1 GPU架构以及研究模型 法模拟了二维正方晶格上的混合Ising模型，得出了 不同自旋比例的临界温度Ho，但其精度较低，同时 1．1 GPU和CUDA架构 GPU是一种专门用来处理在个人电脑、工作站 Metropolis算法需要大量的重复性计算，在CPU运 算耗时较长，效率较低。 或游戏机上那些图像运算工作的微处理器。随着 近年来GPU迅速发展，已经不再局限于图形渲 也越来越多地应用于计算领域，掀起了一场通用计 染的工作，CUDA的开发为GPU提供了更为方便的 算的革命。相比于同时期的CPU，GPU的浮点运算 编程环境，使得本身具有强大运算能力的GPU应用 于越来越多的新领域，如物理学、金融学、医学等， 能力更强，高显存带宽和强大的处理能力提供了很