From: Oz N. T. <na...@gm...> - 2011-06-14 07:11:38
|
On Tue, Jun 14, 2011 at 5:12 AM, Dominique Orban <dom...@gm...> wrote: > > Hi Oz, > I just pushed improved wrappers for the iterative solvers to the 'suitesparse' branch of Pysparse (which is meant to become trunk in the near future). The wrappers are now more streamlined. For some reason, Bicgstab and Gmres were left out when we wrote those wrappers. I just added them. If you look in itsolvers_utils.py, you will see an example towards the bottom of the file. The wrappers accept both ll_mat and PysparseMatrix input. This is the better way to call the iterative solvers. There shouldn't be much of a difference in speed. > What do you mean by "why is your code not there"? > > -- > Dominique Hi Dominique, I looked at the code in itsolvers_utils.py, which uses the classes and objects concept. It is not unfamiliar for me to write code like that, but imho it adds a layer of complexity I prefer to avoid. What I meant in my question is from the two ways I used the matrices, is there some which gives better performance ? >What do you mean by "why is your code not there"? About code in scipy: it is good to have both projects, I just thought if you merge your code into scipy, it will reach bigger audiences. Thanks a lot for you answers, Oz |