[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 <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/ee02b1cf/attachment.html>