For what it’s worth, I had the same issue, I switched to an iterative solver (i.e. opts[f"{option_prefix}pc_type"] = "ksp") and I no longer have the problem.
I think LU is just too greedy memory-wise for large systems.
1 Like