ufuncs, iterators

On Thu, 18 Aug 2005 22:49:29 -0500 (CDT) Barry Smith wrote: > > What is x1, x2 and ip? > >Barry x1 and x2 are 2-arrays; their rows are the 'sample' vectors. ip is the matrix of all inner products from x1 and x2. Simon. -- Simon Burton, B.Sc. Licensed PO Box 8066 ANU Canberra 2601 Aus

ufuncs, iterators

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. Barry> Certainly pointwise addition is the same as adding two operators Barry> together

ufuncs, iterators

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 linear operators I think it would be a big mistake

ufuncs, iterators

On Thu, 18 Aug 2005 20:50:12 -0500 (CDT) Barry Smith wrote: > > 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

ufuncs, iterators

Is there a way to simulate a MatSet (like VecSet) ? Simon. -- Simon Burton, B.Sc. Licensed PO Box 8066 ANU Canberra 2601 Australia Ph. 61 02 6249 6940 http://arrowtheory.com

ufuncs, iterators

On Thu, 18 Aug 2005 19:57:02 -0500 Matthew Knepley wrote: > > > (d) m=exp(m)(pointwise exp) > > We could add pointwise operations just like the VecPointwise*(). > > Matt Yes, we need MatPointwiseMult aswell. Should I make a start on this ? Simon. -- Simon Burton, B

Fortran 90 interfaces with bfort?

Bill, Do you think it would be possible to enhance bfort to generate Fortran 90 interface definitions? Something like PetscErrorCode MatFactorInfoInitialize(MatFactorInfo info) becomes interface subroutine MatFactorInfoInitialize(info,ierr) MatFactorInfo :: info(

ufuncs, iterators

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 produce a Vec) (b) m=add.outer(v1,v2) (sum of all el

ufuncs, iterators

After sleeping on it, it may be ok to include these methods in Mat as methods for "filling up matrices" so long as in the end you end up with a Mat that you then use as an operator. But I'd still like to see/understand a little more of the "construction" process. Classically one would do th

ufuncs, iterators

Simon Burton writes: > On Thu, 18 Aug 2005 19:57:02 -0500 > Matthew Knepley wrote: > >> >> > (d) m=exp(m) (pointwise exp) >> >> We could add pointwise operations just like the VecPointwise*(). >> >> Matt > > Yes, we need MatPointwiseMult aswell. > > Should I make a start

ufuncs, iterators

S V N Vishwanathan writes: > 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 linear operator