  Topological defects form at cosmological phase transitions by the Kibble
mechanism. Cosmic strings and superstrings can lead to particularly interesting
astrophysical and cosmological consequences, but this study is is currently
limited by the availability of accurate numerical simulations, which in turn is
bottlenecked by hardware resources and computation time. Aiming to eliminate
this bottleneck, in recent work we introduced and validated a GPU-accelerated
evolution code for local Abelian-Higgs strings networks. While this leads to
significant gains in speed, it is still limited by the physical memory
available on a graphical accelerator. Here we report on a further step towards
our main goal, by implementing and validating a multiple GPU extension of the
earlier code, and further demonstrate its good scalability, both in terms of
strong and weak scaling. A $8192^3$ production run, using $4096$ GPUs, runs in
$33.2$ minutes of wall clock time on the Piz Daint supercomputer.
