Ravi, I need more info for debugging. Can you provide a simple stand alone code and matrices in petsc binary format that reproduce the error?
MatTransposeMatMult() for mpiaij is a newly developed subroutine - less than one month old and not well tested yet :-( I used petsc-dev/src/mat/examples/tests/ex94.c for testing. Thanks, Hong On Thu, Jan 12, 2012 at 9:17 PM, Mark F. Adams <mark.adams at columbia.edu>wrote: > It looks like the problem is in MatTransposeMatMult and Hong (cc'ed) is > working on it. > > I'm hoping that your output will be enough for Hong to figure this out but > I could not reproduce this problem with any of my tests. > > If Hong can not figure this out then we will need to get the matrix from > you to reproduce this. > > Mark > > > On Jan 12, 2012, at 6:25 PM, Ravi Kannan wrote: > > Hi Mark,**** > ** ** > Any luck with the gamg bug fix?**** > ** ** > Thanks,**** > Ravi.**** > ** ** > *From:* petsc-dev-bounces at mcs.anl.gov [mailto: > petsc-dev-bounces at mcs.anl.gov] *On Behalf Of *Mark F. Adams > *Sent:* Wednesday, January 11, 2012 1:54 PM > *To:* For users of the development version of PETSc > *Subject:* Re: [petsc-dev] boomerAmg scalability**** > ** ** > This seems to be dying earlier than it was last week, so it looks like a > new bug in MatTransposeMatMult.**** > ** ** > Mark**** > ** ** > On Jan 11, 2012, at 1:59 PM, Matthew Knepley wrote:**** > > > **** > On Wed, Jan 11, 2012 at 12:23 PM, Ravi Kannan <rxk at cfdrc.com> wrote:**** > > Hi Mark,**** > **** > I downloaded the dev version again. This time, the program crashes even > earlier. Attached is the serial and parallel info outputs.**** > **** > Could you kindly take a look.**** > > ** ** > It looks like this is a problem with MatMatMult(). Can you try to > reproduce this using KSP ex10? You put**** > your matrix in binary format and use -pc_type gamg. Then you can send us > the matrix and we can track**** > it down. Or are you running an example there?**** > ** ** > Thanks,**** > ** ** > Matt**** > ** ** > **** > > **** > Thanks,**** > Ravi.**** > **** > *From:* petsc-dev-bounces at mcs.anl.gov [mailto: > petsc-dev-bounces at mcs.anl.gov] *On Behalf Of *Mark F. Adams > *Sent:* Monday, January 09, 2012 3:08 PM**** > > *To:* For users of the development version of PETSc > *Subject:* Re: [petsc-dev] boomerAmg scalability**** > ** ** > **** > Yes its all checked it, just pull from dev.**** > Mark**** > **** > On Jan 9, 2012, at 2:54 PM, Ravi Kannan wrote:**** > > ** ** > Hi Mark,**** > **** > Thanks for your efforts.**** > **** > Do I need to do the install from scratch once again? Or some particular > files (check out gamg.c for instance)?**** > **** > Thanks,**** > Ravi.**** > **** > *From:* petsc-dev-bounces at mcs.anl.gov [mailto: > petsc-dev-bounces at mcs.anl.gov] *On Behalf Of *Mark F. Adams > *Sent:* Friday, January 06, 2012 10:30 AM > *To:* For users of the development version of PETSc > *Subject:* Re: [petsc-dev] boomerAmg scalability**** > **** > I think I found the problem. You will need to use petsc-dev to get the > fix.**** > **** > Mark**** > **** > On Jan 6, 2012, at 8:55 AM, Mark F. Adams wrote:**** > > > > **** > Ravi, I forgot but you can just use -ksp_view_binary to output the matrix > data (two files). You could run it with two procs and a Jacobi solver to > get it past the solve, where it writes the matrix (I believe).**** > Mark**** > **** > On Jan 5, 2012, at 6:19 PM, Ravi Kannan wrote:**** > > > > **** > Just send in another email with the attachment.**** > **** > *From:* petsc-dev-bounces at mcs.anl.gov [mailto: > petsc-dev-bounces at mcs.anl.gov] *On Behalf Of *Jed Brown > *Sent:* Thursday, January 05, 2012 5:15 PM > *To:* For users of the development version of PETSc > *Subject:* Re: [petsc-dev] boomerAmg scalability**** > **** > On Thu, Jan 5, 2012 at 17:12, Ravi Kannan <rxk at cfdrc.com> wrote:**** > I have attached the verbose+info outputs for both the serial and the > parallel (2 partitions). NOTE: the serial output at some location says > PC=Jacobi! Is it implicitly converting the PC to a Jacobi?**** > **** > Looks like you forgot the attachment.**** > **** > **** > **** > > > > **** > ** ** > -- > 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**** > ** ** > > > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120112/fe84d368/attachment.html>
