Small improvements to the amg_hierarchy. Use splu direct solver and
use max_coarse = 10000.