Dear Tsai, The determinant of a large sparse matrix is useless 99.99% of the time because except if the eigenvalues are very close from 1, either the determinant is too small (represented by zero, even if all the eigenvalues are nonzero) or too large (> 1E308). So, for instance, if you want to test the invertibility of a system, the computation of the determinant is not appropriate. A better approach is for instance compute the smallest singluar value with a Lanczos or Von Mises algorithm. In fact, depending on the use of the determinant you need, you better have to try to find an indirect mean not to compute the determinant ... Yves. Le 19/11/2012 15:26, Finite Tsai a écrit : > > Dear Yves, > How can I use the getfem_superlu.h for obtaining the determinant of a > sparse matrix? Any comment of yours will be highly appreciated. > > best regard > CC Tsai > > > _______________________________________________ > Getfem-users mailing list > [email protected] > https://mail.gna.org/listinfo/getfem-users -- 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
