  Gauge covariant smearing based on the 3D lattice Laplacian can be used to
create extended operators that have better overlap with hadronic ground states.
For staggered quarks, we make use of two-link parallel transport to preserve
taste properties. We have implemented the procedure in QUDA. We present the
performance of this code on the NVIDIA A100 GPUs in Indiana University's Big
Red 200 supercomputer and on the AMD MI250X GPUs in Oak Ridge Leadership
Computer Facility's (OLCF's) Crusher and discuss its scalability. We also study
the performance improvement from using NVSHMEM on OLCF's Summit. Reusing
precomputed two-link products for all sources and sinks, it reduces the total
smearing time for a baryon correlator measurement by a factor of 100-120 as
compared with the original MILC code and reduces the overall time by 60-70%.
