[petsc-dev] [pflotran-dev: 976] Re: [pflotran-users: 998] Distorted structured grid

[Paolo Orsini]

Hi Jed,

I tried the MatConvert updated version, (from MPIAdj type to MPIAij type)
but it seems still slow as before.

Anyway, don't worry to much. I abandoned the idea of the conversion.
I am forming the adjacency matrix directly in MPI AIJ format, so I don't
have to do the conversion.
I also realised that the Dual Graph is formed in parallel (in the current
PFLOTRAN implementation), so I can't use the MatMatTransposeMult (which
works only for sequential matrix).

I will use instead MatCreatTranspose and MAtMatMult, then filtering the
resulting matrix row by row and from the DualGraph.
It will require more memory, but it should work. Do you see any problem
with this?

Paolo


On Fri, Nov 15, 2013 at 10:07 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:

> Jed Brown <jedbrown at mcs.anl.gov> writes:
>
> > Jed Brown <jedbrown at mcs.anl.gov> writes:
> >> Thanks for the reminder, I forgot that we were missing a specialized
> >> function.  What I recommended should work and be fast with 'next'.  Can
> >> you try that?
> >
> > Hang on, I messed something up.  Will fix and reply here.
>
> Okay, done.  Please check it out in 'next'.
>
>
> https://bitbucket.org/petsc/petsc/wiki/Home#markdown-header-following-petsc-development
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20131118/f6715fe7/attachment.html>