From: Dominique Orban <dominique.orban@gm...>  20080128 21:34:10

Hello, I wrote an addition to the UMFPack module of PySparse that lets users retrieve the details of the factorization. Here is an example: from pysparse import spmatrix, umfpack A = spmatrix.ll_mat_from_mtx('bcsstk11.mtx') LU = umfpack.factorize(A) (L, U, P, Q, R, do_recip) = LU.lu() The help of lu() is as follows:  lu(...) self.lu() Returns L and U factors, permutation and scaling information. Use: (L, U, P, Q, R, do_recip) = self.lu(). The original matrix A is factorized into L U = P R A Q where L is unit lower triangular, U is upper triangular, P and Q are permutation matrices, R is a rowscaling diagonal matrix such that the ith row of A has been divided by R[i] if do_recip = True, the ith row of A has been multiplied by R[i] if do_recip = False. L and U are returned as ll_mat sparse matrices. P, Q and R are returned as NumPy arrays.  I would very much like to see this function included in PySparse. I will email the code directly to Daniel Wheeler. Comments welcome. Cheers, Dominique 
From: Daniel Wheeler <daniel.wheeler2@gm...>  20080129 15:52:48

I committed the changes and bumped the version number. I'll hold off on a release unless there is a pressing need. Cheers On Jan 28, 2008 4:34 PM, Dominique Orban <dominique.orban@...> wrote: > Hello, > > I wrote an addition to the UMFPack module of PySparse that lets users > retrieve the details of the factorization. Here is an example: > > from pysparse import spmatrix, umfpack > A = spmatrix.ll_mat_from_mtx('bcsstk11.mtx') > LU = umfpack.factorize(A) > (L, U, P, Q, R, do_recip) = LU.lu() > > The help of lu() is as follows: > >  > lu(...) > self.lu() > Returns L and U factors, permutation and scaling information. > > Use: (L, U, P, Q, R, do_recip) = self.lu(). > The original matrix A is factorized into > L U = P R A Q > where L is unit lower triangular, > U is upper triangular, > P and Q are permutation matrices, > R is a rowscaling diagonal matrix such that > the ith row of A has been divided by R[i] if do_recip = True, > the ith row of A has been multiplied by R[i] if do_recip = False. > > L and U are returned as ll_mat sparse matrices. > P, Q and R are returned as NumPy arrays. >  > > I would very much like to see this function included in PySparse. I > will email the code directly to Daniel Wheeler. > > Comments welcome. > > Cheers, > Dominique > >  > This SF.net email is sponsored by: Microsoft > Defy all challenges. Microsoft(R) Visual Studio 2008. > http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/ > _______________________________________________ > Pysparseusers mailing list > Pysparseusers@... > https://lists.sourceforge.net/lists/listinfo/pysparseusers >  Daniel Wheeler 
From: Daniel Wheeler <daniel.wheeler2@gm...>  20080129 18:25:09

I think that was initiated before my time using pysparse. It is not a branch but a new project. We have wraped a lot of the pysparse functionality. Take a look at the following. <http://matforge.org/fipy/browser/trunk/fipy/tools/pysparseMatrix.py>; and <http://matforge.org/fipy/browser/trunk/fipy/solvers/pysparse/linearCGSSolver.py>; Maybe some of these classes should be added back in. Cheers On Jan 29, 2008 11:25 AM, Dominique Orban <dominique.orban@...> wrote: > Ok, thanks. > > I had a question about the development of PySparse. I saw a branch > called 'pysparse2' in the CVS tree. At first glance, it seemed like an > effort to have nice Python wrappers around the C modules. Is there > still any effort going into that? > > Cheers, > Dominique > > > > On 1/29/08, Daniel Wheeler <daniel.wheeler2@...> wrote: > > I committed the changes and bumped the version number. I'll hold off > > on a release > > unless there is a pressing need. > > > > Cheers > > > > On Jan 28, 2008 4:34 PM, Dominique Orban <dominique.orban@...> wrote: > > > Hello, > > > > > > I wrote an addition to the UMFPack module of PySparse that lets users > > > retrieve the details of the factorization. Here is an example: > > > > > > from pysparse import spmatrix, umfpack > > > A = spmatrix.ll_mat_from_mtx('bcsstk11.mtx') > > > LU = umfpack.factorize(A) > > > (L, U, P, Q, R, do_recip) = LU.lu() > > > > > > The help of lu() is as follows: > > > > > >  > > > lu(...) > > > self.lu() > > > Returns L and U factors, permutation and scaling information. > > > > > > Use: (L, U, P, Q, R, do_recip) = self.lu(). > > > The original matrix A is factorized into > > > L U = P R A Q > > > where L is unit lower triangular, > > > U is upper triangular, > > > P and Q are permutation matrices, > > > R is a rowscaling diagonal matrix such that > > > the ith row of A has been divided by R[i] if do_recip = True, > > > the ith row of A has been multiplied by R[i] if do_recip = False. > > > > > > L and U are returned as ll_mat sparse matrices. > > > P, Q and R are returned as NumPy arrays. > > >  > > > > > > I would very much like to see this function included in PySparse. I > > > will email the code directly to Daniel Wheeler. > > > > > > Comments welcome. > > > > > > Cheers, > > > Dominique > > > > > >  > > > This SF.net email is sponsored by: Microsoft > > > Defy all challenges. Microsoft(R) Visual Studio 2008. > > > http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/ > > > _______________________________________________ > > > Pysparseusers mailing list > > > Pysparseusers@... > > > https://lists.sourceforge.net/lists/listinfo/pysparseusers > > > > > > > > > > >  > > Daniel Wheeler > > > >  > > This SF.net email is sponsored by: Microsoft > > Defy all challenges. Microsoft(R) Visual Studio 2008. > > http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/ > > _______________________________________________ > > Pysparsedevelopers mailing list > > Pysparsedevelopers@... > > https://lists.sourceforge.net/lists/listinfo/pysparsedevelopers > > >  Daniel Wheeler 