Parallel 3D finite-difference time-domain method on multi-GPU systems. (English) Zbl 1218.78146

Summary: Finite-difference time-domain (FDTD) is a popular but computational intensive method to solve Maxwell’s equations for electrical and optical devices simulation. This paper presents implementations of three-dimensional FDTD with convolutional perfect match layer (CPML) absorbing boundary conditions on graphics processing unit (GPU). Electromagnetic fields in Yee cells are calculated in parallel millions of threads arranged as a grid of blocks with compute unified device architecture (CUDA) programming model and considerable speedup factors are obtained versus sequential CPU code. We extend the parallel algorithm to multiple GPUs in order to solve electrically large structures. Asynchronous memory copy scheme is used in data exchange procedure to improve the computation efficiency. We successfully use this technique to simulate pointwise source radiation and validate the result by comparison to high precision computation, which shows favorable agreements. With four commodity GTX295 graphics cards on a single personal computer, more than 4000 million Yee cells can be updated in one second, which is hundreds of times faster than traditional CPU computation.


78M20 Finite difference methods applied to problems in optics and electromagnetic theory
65Y05 Parallel numerical computation
65Y10 Numerical algorithms for specific classes of architectures
Full Text: DOI


[1] Yee K., IEEE Trans. Antennas Propag. 14 pp 302– · Zbl 1155.78304
[2] DOI: 10.1109/22.55775
[3] DOI: 10.1364/OE.14.011796
[4] Wenhua Y., IEEE Trans. Adv. Packag. 30 pp 335–
[5] DOI: 10.2528/PIER07120703
[6] Sano K., Interdisci. Inf. Sci. 15 pp 67–
[7] DOI: 10.1109/MAP.2005.1608730
[8] DOI: 10.1016/j.jcp.2009.03.018 · Zbl 1167.82347
[9] DOI: 10.1109/TMAG.2009.2012614
[10] DOI: 10.1155/2009/308606 · Zbl 05761165
[11] DOI: 10.1016/j.cpc.2010.05.005 · Zbl 1219.82020
[12] DOI: 10.1016/j.jcp.2010.02.009 · Zbl 1214.68477
[13] DOI: 10.1109/TMAG.2010.2043655
[14] Mur G., IEEE Trans. Electromagn. Compat. 23 pp 377–
[15] Liao Z. P., Scientia SinicQ (Series A) pp 1063–
[16] DOI: 10.1006/jcph.1994.1159 · Zbl 0814.65129
[17] DOI: 10.1109/8.546249
[18] DOI: 10.1002/1098-2760(20001205)27:5<334::AID-MOP14>3.0.CO;2-A
[19] A. Taflove and S. C. Hagness, Computational Electrodynamics: The Finite-Difference Time-domain method (Artech House Norwood, Boston, 2005) p. 305. · Zbl 0840.65126
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.