Re: [Pysparse-developers] [Pysparse-users] Factors in Umfpack: contribution to PySparse

[Daniel W. <dan...@gm...>]

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 <dom...@gm...> 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 <dan...@gm...> 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 <dom...@gm...> 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 row-scaling diagonal matrix such that
> > >            the i-th row of A has been divided    by R[i] if do_recip = True,
> > >            the i-th 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 e-mail 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/
> > > _______________________________________________
> > > Pysparse-users mailing list
> > > Pys...@li...
> > > https://lists.sourceforge.net/lists/listinfo/pysparse-users
> > >
> >
> >
> >
> > --
> > 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/
> > _______________________________________________
> > Pysparse-developers mailing list
> > Pys...@li...
> > https://lists.sourceforge.net/lists/listinfo/pysparse-developers
> >
>



-- 
Daniel Wheeler