ufuncs, iterators

2005-08-18 Thread Matthew Knepley
Simon Burton simon at arrowtheory.com writes: Hi, Along the lines of python's numarray [1], we need some way of operating pointwise (and inner/outer operations) on Mat/Vec objects. In particular, we need things like the following: (a) v=add.reduce(m) (sum along rows/cols of a Mat to

ufuncs, iterators

2005-08-18 Thread Barry Smith
Simon, There is not because what does MatSet() mean for a sparse matrix? Make it dense, just set all the current nonzeros, error? You tell us, I could live with either of the later two but don't really like the first. Barry BTW: possibly related note, are you using dense matrices

ufuncs, iterators

2005-08-18 Thread Barry Smith
On Fri, 19 Aug 2005, S V N Vishwanathan wrote: Hi! BTW: possibly related note, are you using dense matrices sometimes to represent just 2-arrays; that is, not as representations of linear operators. If so, I do not think this is the correct approach! Conceptually PETSc Mat's are

ufuncs, iterators

2005-08-18 Thread Barry Smith
What is x1, x2 and ip? Barry On Fri, 19 Aug 2005, S V N Vishwanathan wrote: Hi! I am not sure I understand the fine difference. As far as we are concerned, all the operations which we are doing (point wise addition, addition, multiplication etc.) are on the linear operator.