X. Sherry Li

If only b (and x) are complex, it's more efficient to solve 2 real linear systems: A x_r = b_r (real part of b) A x_i = b_i (imaginary part...

"inconsistent" means the two solutions are not identical? Have you checked the error : ||X_matlab - X_superlu|| / ||X_matlab|| ?

I don't know what goes wrong. If you see file SRC/ssp_blas2.c, it includes slu_sdefs.h, which contains the prototype definition: extern int strsv_(char*, char*, char*, int*, float*, int*, float*, int*); But...

CMakeLists.txt does use find_package(BLAS). What more needs to be done? If you have suggestions, please set up a PR.

It dies in iterative refinement, after 2 iterations. Can you turn off refinement? In fact, your first solution already looks good -- backward error 'berr' is small: ( 0) .....

Sorry, it seems it passes iterative refinement (you can still safely turn it off.) It probably dies in the final redistribution of the solution vector X. Which version (or commit)...

Hard to tell what my be wrong. In this line: 2 0x00000000000c7dd8 dScatter_B3d() /home/q8j/GIT/superlu_dist/SRC/dnrformat_loc3d.c:565 matCopy(A3d->m_loc, nrhs, B, ldb, Btmp, A3d->m_loc); Can you do some detective work, printing nrhs, ldb, A3d->m_loc,...

The numbers all look correct to me. It dies at line 34, in the loop. Do you know what is "i" at that point?

The copy destination 'Dst' is just the right-hand side of the original linear system. Somehow it got interference from your program. Can you run the standalone pddrive3d(), set a breakpoint...

Do you know how the vector Tpetra::MultiVector is stored? Perhaps it contains some structure. i.e., not a flat vector. @jwillenbring @egboman