MueLu or Amesos2 behavior changed between 12.10.1 and develop, such that preconditioner setup now fails with a singular matrix (Was: "Amesos2 KLU2 Crash")
Created by: dlkarnit
I'm experiencing a crash while attempting to solve a 3D poisson problem with Belos/MueLu/Tpetra built from the develop branch that I retrieved yesterday afternoon. Previously this code would crash at exit related to issue #1889 (closed) when It was fully integrated into our product. This crash though is now occurring within our unit tests. I can provide some sample code describing our problem setup that, however, does not reproduce this behavior. I've trace the crash to a NULL field in the KLU numeric factorization step, specifically:
Amesos2_KLU2_def.hpp:327 at this point the data_.numeric_ pointer is 0x0 The state of data_ is: {symbolic_ = 0x561e4d0, numeric_ = 0x0, common_ = {tol = 0.001, memgrow = 1.2, initmem_amd = 1.2, initmem = 10, maxwork = 0, btf = 1, ordering = 0, scale = 2, malloc_memory = 0x40e580 malloc@plt, realloc_memory = 0x40d8e0 <realloc@p lt>, free_memory = 0x40f1b0 free@plt, calloc_memory = 0x40e540 calloc@plt, user_order = 0x0, user_data = 0x0, halt_if_singular = 1, status = 1, nrealloc = 0, structural_rank = 16, numerical_rank = 14, singular_col = 6, noffdiag = 0, fl ops = -1, rcond = -1, condest = -1, rgrowth = -1, work = 0, memusage = 420, mempeak = 3404}}
I've also attached the "verbose: high" output from the preconditioner parameters and the stack trace.