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>
