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 ([email protected]) 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 --------- _______________________________________________ Getfem-users mailing list [email protected] https://mail.gna.org/listinfo/getfem-users
