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

[Yves Renard]

Dear Dawid,

The test program helmholtz.cc in the directory 'tests' use complex sparse 
matrices and superlu and it works.

You can also export your matrix and compute the smallest eigen values in 
Matlab to see if it is singular.

Yves.







On jeudi 3 novembre 2011, Dawid Midura wrote:
> 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.


-- 

  Yves Renard (address@hidden)       tel : (33) 04.72.43.87.08
  Pole de Mathematiques, INSA-Lyon             fax : (33) 04.72.43.85.29
  20, rue Albert Einstein
  69621 Villeurbanne Cedex, FRANCE
  http://math.univ-lyon1.fr/~renard

---------