Search papers, labs, and topics across Lattice.
This paper introduces a GPU-initiated discrete simulated bifurcation (dSB) architecture leveraging NVIDIA DOCA GPUNetIO to overcome communication bottlenecks in GPU-based optimization for both low-latency requests and large, dense models. For resident models, a persistent service enables single-thread-block solves and batched transmissions, achieving median round-trip latencies as low as 189 $\mu$s under high concurrent load. A streaming solver handles models exceeding device memory by retaining dynamical state and reusing coupling tiles, enabling evaluation of ten-million-variable dense binary matrices at 307 Gb/s.
GPU-native optimization achieves sub-200 $\mu$s latency for millions of requests per second and streams multi-terabyte models through gigabyte buffers.
GPU-based optimization faces two communication bottlenecks: coordinating frequent requests and delivering dense models that exceed device memory. We present a discrete simulated bifurcation (dSB) architecture that addresses both through NVIDIA DOCA GPUNetIO. For resident models, a persistent service receives field updates, executes each solve within one GPU thread block, and returns the result. Exact integer coupling sums, GPU work queues, and batched transmission keep the receive--solve--reply path on the device without a dedicated CPU data-path core. In comparisons with socket-based servers using the same solver, the largest latency gains occur under concurrent load. As the offered load increases from 400 to 800 thousand requests per second, median round-trip latency rises by only 6\%. At the highest tested load, median and 99th-percentile latencies are 189 and 218~$\mu$s, compared with 288 and 609~$\mu$s for the tuned persistent CPU proxy across repeated runs. For models larger than device memory, a streaming solver retains dynamical state on the GPU and reuses incoming coupling tiles across replicas. It evaluates ten-million-variable dense binary matrices at approximately 307~Gb/s, consuming a 12.5-TB logical matrix through a 64-MiB packet buffer. Ground-state recovery on planted instances and agreement with reference executions verify the computation. Together, the two modes scale dSB to concurrent requests and dense models beyond GPU memory.