[petsc-users] MatRARt

Jed Brown jedbrown at mcs.anl.gov
Thu Jul 19 08:31:38 CDT 2012


On Thu, Jul 19, 2012 at 8:10 AM, John Fettig <john.fettig at gmail.com> wrote:

> 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).
>

You can use this mirror

https://bitbucket.org/petsc/petsc-dev


> > 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?
>

Time and calling convenience. It's also not clear that it can be
implemented significantly more efficiently than R(AP).
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120719/9aacb881/attachment.html>


More information about the petsc-users mailing list