[Getfem-users] Gmm++ problem with superlu interface.

[Dawid Midura]

Hi all,

I have the following code that gives me trouble:

  typedef gmm::wsvector<std::complex<double> > complex_sparse_wvector;

  typedef gmm::col_matrix<complex_sparse_wvector> complex_sparse_wmatrix;

  

  complex_sparse_wvector F(N);

  complex_sparse_wmatrix K(N,N);

  ...

  // Assemble stiffness matrix K, and load vector F.

  ...

  double condition;

  vector<complex<double> > U(N);

  gmm::SuperLU_solve(K, U, F, condition);

I compile with appropriate flags for superlu, when I run the code I get the following:

============================================

|      An error has been detected !!!      |

============================================

Error in /usr/include/gmm/gmm_superlu_interface.h, line 244 : 

SuperLU solve failed: info=-2

From what I have read about SuperLU, info=-2 means that there is a problem with the second argument the the SuperLU_gssvx which is matrix K in the SuperLU format SuperMatrix, I am not using getfem to assemble K or F, I have written my own procedures to do that. I know K is not singular, when I make my assembly procedures work with dense matrices and vectors I don't have any problems. Any help would be greatly appreciated. 

Thanks,

Dawid M.