{"id":"arxiv/2604.26441","url":"https://arxiv.org/abs/2604.26441","tags":["cs.CE","math.NA"],"score":0,"title":"A Matrix-Free Galerkin Multigrid Solver and Failure-Mode Screen for Single-GPU 3D SIMP Linear Systems","domain":"cs.CE","source":"arxiv","authors":["Shaoliang Yang","Jun Wang","Yunsheng Wang"],"reviews":0,"subject":"cs.CE","summary":"Large 3D SIMP studies require repeated elasticity solves for density-dependent operators whose finest matrices are expensive to assemble and whose conditioning degrades under high contrast. We study this linear-solver layer rather than claiming end-to-end optimization acceleration. The solver builds a matrix-free Galerkin geometric multigrid (GMG) hierarchy around a fused fine operator: the finest level remains matrix-free, the first coarse level is assembled by local Galerkin aggregation, and deeper levels use sparse Galerkin products. The practical default is FP32-GMG; BF16 is evaluated as…","category":"research","comments":0,"disputed":false,"division":"cs","mentions":0,"discussed":0,"published":"2026-04-29","reactions":{}}