I agree with Barry here. When I first started
playing with petsc, I found the "python" directory confusing and
frustrating for quite some time. I think it was just
because of the name!
Simon.
On Mon, 10 Dec 2007 10:16:13 -0600
Barry Smith wrote:
>
> I sent email about a week ago and no on
On Thu, 03 Nov 2005 17:49:46 -0600
Matthew Knepley wrote:
>
> Simon Burton writes:
>
> This was patched in 2.3.0, but maybe you have an unpatched version. It
> is
> easy to check. The symbol comes from
We have been tracking your bitkeeper repository. I don't
gramming (with all these
different matrix types). That really excites me, but i guess it's beside the
point right now.
I'm cc'ing the petsc-dev list, hope thats OK.
ciao!
Simon.
--
Simon Burton, B.Sc.
Licensed PO Box 8066
ANU Canberra 2601
Australia
Ph. 61 02 6249 6940
http://arrowtheory.com
_SeqDense_SeqDense()?
>
> Did you use inplace mat product in your python test?
We ended up wrapping BLASgemm_ and using it directly.
This call provides a kind of MatMatMultAdd (except the result
gets stored in the add matrix) which we need. But obviously
this will not work for the other ma
so, do you have testing code for it? I can test it
> once I see the MatMatMult_seqdense.
We tested it using python. I see also there is a test in
src/mat/examples/tests/ex94.c .
Simon.
--
Simon Burton, B.Sc.
Licensed PO Box 8066
ANU Canberra 2601
Australia
Ph. 61 02 6249 6940
http://arrowtheory.com
Hi,
I sent through a patch to implement MatMatMult* for seqdense matrices.
Here it is again; is there anything wrong with it ? It just calls through to
BLASgemm_ .
thanks,
Simon.
--
Simon Burton, B.Sc.
Licensed PO Box 8066
ANU Canberra 2601
Australia
Ph. 61 02 6249 6940
http
x27;s invarient under linear isometries).
Simon.
On Thu, 18 Aug 2005 23:12:21 -0500 (CDT)
Barry Smith wrote:
>
> In terms of exp( -1/2\sigma^{2} ||x_{i} - x_{j}||_{2}^{2})
> what are they?
>
>Thanks
>
>Barry
>
>
> On Fri, 19 Aug 2005, Simon Burto
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
se) != 2-d array , the are completely different
> beasts mathematically
>
>
>
> On Fri, 19 Aug 2005, Simon Burton wrote:
>
> >
> > 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
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
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 ?
Simo
numarray do the work, but this seems doomed to failure because
numarray is inherently dense.
Simon.
[1]: http://stsdas.stsci.edu/numarray/numarray-1.3.html/node35.html
--
Simon Burton, B.Sc.
Licensed PO Box 8066
ANU Canberra 2601
Australia
Ph. 61 02 6249 6940
http://arrowtheory.com
12 matches
Mail list logo