[petsc-users] SuperLU_dist issue in 3.7.4 failure of repeated calls to MatLoad() or MatMPIAIJSetPreallocation() with the same matrix
Anton Popov
popov at uni-mainz.de
Mon Oct 24 12:09:37 CDT 2016
On 10/24/2016 05:47 PM, Hong wrote:
> Barry,
> Your change indeed fixed the error of his testing code.
> As Satish tested, on your branch, ex16 runs smooth.
>
> I do not understand why on maint or master branch, ex16 creases inside
> superlu_dist, but not with mumps.
>
I also confirm that ex16 runs fine with latest fix, but unfortunately
not my code.
This is something to be expected, since my code preallocates once in the
beginning. So there is no way it can be affected by multiple
preallocations. Subsequently I only do matrix assembly, that makes sure
structure doesn't change (set to get error otherwise).
Summary: we don't have a simple test code to debug superlu issue anymore.
Anton
> Hong
>
> On Mon, Oct 24, 2016 at 9:34 AM, Satish Balay <balay at mcs.anl.gov
> <mailto:balay at mcs.anl.gov>> wrote:
>
> On Mon, 24 Oct 2016, Barry Smith wrote:
>
> >
> > > [Or perhaps Hong is using a different test code and is observing bugs
> > > with superlu_dist interface..]
> >
> > She states that her test does a NEW MatCreate() for each
> matrix load (I cut and pasted it in the email I just sent). The
> bug I fixed was only related to using the SAME matrix from one
> MatLoad() in another MatLoad().
>
> Ah - ok.. Sorry - wasn't thinking clearly :(
>
> Satish
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20161024/1023efc1/attachment.html>
More information about the petsc-users
mailing list