By the way, I'm calling the gmm++ version of SuperLU_solve, that is I include gmm/gmm_superlu_interface.h header file, not getfem_superlu.h. I just looked at getfem_superlu.h file, and in it is says that the gmm version of superlu interface is not being used, I guess getfem++ provides its own? Maybe if I try to use getfem superlu interface the problem will go away.
Dawid M. On Thu, Nov 3, 2011 at 1:01 PM, Yves Renard <[email protected]>wrote: > > > 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 >
_______________________________________________ Getfem-users mailing list [email protected] https://mail.gna.org/listinfo/getfem-users
