[petsc-users] MatRARt

[Hong Zhang]

John:

>
> > Can you test MatRARt on
> > petsc-3.3/src/mat/examples/tests/ex94 with attached matrix data file by
> >
> > ./ex94 -f0 medium -f1 medium
> >
> > Do it crash?
>
> Yes, it does crash in the same exact fashion as the code I sent.  This
> is 3.3-p2.  petsc-dev also crashes, but I haven't pulled changes in
> about a month (and the server is down right now).
>

Hmm, ex94 runs well on my tests.
May I have your ~petsc/configure.log?

>
> > From flops analysis and experiments, MatRARt() is not as efficient as
> > MatPtAP,
> > likely costs twice as MatPtAP.
>
> I can just as easily construct P, so this isn't a big deal.  How would
> MatPtAP compare to a hypothetical MatRAP?  Is there a reason this
> routine doesn't exist?
>

Using petsc aij matrix format, sparse inner-product used in
A*Rt (MatMatTransposeMult)
costs twice as A*P (MatMatMult).
Out-product used in Pt*A (MatTransposeMatMult) requires extensive
data movement. I'm planning to implement
MatRAP (MatMatMatMult), but have not been able to work on it.
Suggestions are welcome.

Hong
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120719/36db808c/attachment.html>