Okay, I don't actually need it. I ran it until it finally completed (45 minutes for 1 million nonzeros) and it didn't seem to help my preconditioner.
Thanks for looking into it. Seems like you could leave it in as a Developer function and note that it is not efficient. > On Mon, Nov 29, 2010 at 4:34 PM, Jed Brown <jed at 59a2.org> wrote: > >> On Mon, Nov 29, 2010 at 23:32, Jed Brown <jed at 59a2.org> wrote: >> >>> On Mon, Nov 29, 2010 at 23:30, Barry Smith <bsmith at mcs.anl.gov> wrote: >>> >>>> who wrote this routine >>> >>> >>> changeset: 19:eb3da2ea4a05 >>> parent: 16:4babff491f33 >>> user: knepley at spinetta.mcs.anl.gov >>> date: Thu Oct 11 06:00:09 2001 -0600 >>> summary: bk-changeset-1.21.1.1 >>> >> >> MatMultConstrained also still exists from this changeset, is marked as >> "beginner", but has zero implementations in PETSc. >> > > I thought Barry ripped that out long ago. > > Matt > > >> Jed >> > > > > -- > What most experimenters take for granted before they begin their > experiments > is infinitely more interesting than any results to which their experiments > lead. > -- Norbert Wiener >
