Re: [petsc-users] SuperLU_dist issue in 3.7.4 failure of repeated calls to MatLoad() or MatMPIAIJSetPreallocation() with the same matrix

[Anton]

On 10/24/16 8:21 PM, Hong wrote:

 Anton :
If replacing superlu_dist with mumps, does your code work?

yes

Hong

    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 <ba...@mcs.anl.gov
    <mailto:ba...@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