Floating-point consistent cross-verification methodology for reproducible and interoperable DDA solvers with fair benchmarking

Abstract

The discrete dipole approximation (DDA) is a widely used and versatile numerical method for solving electromagnetic scattering by arbitrarily shaped objects. Despite its popularity, quantitative comparisons between independent implementations remain challenging due to differences in linear-system conventions, solver settings, and default numerical parameters. In this work, we introduce a unified software-assisted methodology for cross-verification and benchmarking of three major open-source DDA solvers: DDSCAT, ADDA, and IFDDA. We demonstrate how machine-precision agreement can be achieved across implementations by aligning all free parameters and provide practical equivalence tables enabling reproducible and interoperable simulations. Using this methodology, we perform systematic CPU and GPU performance comparisons covering OpenMP, MPI, and CUDA/OpenCL parallelization. Beyond benchmarking, our approach serves as a practical guide for configuring consistent DDA simulations and for understanding how precision, solver choice, and hardware architecture affect runtime, scalability, and accuracy in computational light-scattering studies. The software package also supports regression testing and bitwise reproducibility verification for future code releases.

References in corpus (7)

• The discrete dipole approximation: an overview and recent developments
• Radiative Torques on Interstellar Grains: I. Superthermal Spinup
• Experimental and theoretical study of light scattering by individual mature red blood cells by use of scanning flow cytometry and discrete dipole approximation
• Convergence of the discrete dipole approximation. I. Theoretical analysis
• Reproducibility, accuracy and performance of the Feltor code and library on parallel computer architectures
• Vector Bessel beams: general classification and scattering simulations
• Simulating electron energy-loss spectroscopy and cathodoluminescence for particles in arbitrary host medium using the discrete dipole approximation