[petsc-users] MUMPS solver crash with petsc-3.2
Gong Ding
gdiso at ustc.edu
Wed Dec 7 00:48:09 CST 2011
With a small matrix, mumps of AIX reports that
Entering DMUMPS driver with JOB, N, NZ = 1 25 105
DMUMPS 4.10.0
L U Solver for unsymmetric matrices
Type of parallelism: Working host
****** ANALYSIS STEP ********
Resetting candidate strategy to 0 because NSLAVES=1
Entering analysis phase with ...
N NZ LIW INFO(1)
25 105 287 0
Matrix entries: IRN() ICN()
1 1 1 2 1 3
2 1 2 2 2 4
2 5 3 1 3 3
3 4 3 6 4 2
4 3 4 4 4 7
4 8 5 2 5 5
5 7 5 9 6 3
6 6 6 8 6 10
7 4 7 5 7 7
7 11 7 12 8 4
8 6 8 8 8 12
8 13 9 5 9 9
9 11 9 14 10 6
10 10 10 13 10 15
11 7 11 9 11 11
11 16 11 17 12 7
12 8 12 12 12 17
12 18 13 8 13 10
13 13 13 18 13 19
14 9 14 14 14 16
15 10 15 15 15 19
16 11 16 14 16 16
16 20 17 11 17 12
17 17 17 20 17 21
18 12 18 13 18 18
18 21 18 22 19 13
19 15 19 19 19 22
20 16 20 17 20 20
20 23 21 17 21 18
21 21 21 23 21 24
22 18 22 19 22 22
22 24 23 20 23 21
23 23 23 25 24 21
24 22 24 24 24 25
25 23 25 24 25 25
... Structural symmetry (in percent)= 100
Density: NBdense, Average, Median = 0 3 3
... No column permutation
Ordering based on AMF
IKEEP(.,1)= 14 16 15 17 18 24 13 9 11 5
25 19 23 10 4 12 8 7 6 22
21 20 3 2 1
IKEEP(.,2)= 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 2 7 0 0
0 0 0 0 0
IKEEP(.,3)= 0 1 1 0 1 1 0 0 0 0
1 1 0 0 1 1 4 5 0 1
1 3 4 2 5
PTRAR(.,1)= 14 16 15 17 17 18 13 9 11 5
18 18 18 10 4 12 8 7 6 18
18 18 3 2 1
PTRAR(.,2)= 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0
PTRAR(.,3)= 3 4 4 3 4 4 5 5 5 3
4 4 5 3 4 4 7 7 0 1
1 1 0 0 0
NFSIZ(.) = 3 4 4 7 0 0 5 5 4 4
0 7 0 3 3 4 5 5 4 0
0 0 4 4 3
FILS (.) = 0 -3 -1 5 -8 11 0 0 -14 -15
-23 22 6 0 0 -9 0 0 -10 13
20 21 -24 -25 0
FRERE(.) = -3 -4 -2 -12 26 26 2 16 -16 -19
26 0 26 -9 -10 7 4 17 18 26
26 26 19 -23 -24
ERROR: 0031-250 task 0: Segmentation fault
The output is exactly the same as on linux, except the error.
On Linux, the naxt line is
Leaving analysis phase with ...
INFOG(1) = 0
INFOG(2) = 0
-- (20) Number of entries in factors (estim.) = 185
-- (3) Storage of factors (REAL, estimated) = 185
-- (4) Storage of factors (INT , estimated) = 372
-- (5) Maximum frontal size (estimated) = 7
-- (6) Number of nodes in the tree = 18
-- (32) Type of analysis effectively used = 1
-- (7) Ordering option effectively used = 2
ICNTL(6) Maximum transversal option = 0
ICNTL(7) Pivot order option = 7
Percentage of memory relaxation (effective) = 20
Number of level 2 nodes = 0
Number of split nodes = 0
I compiled petsc-dev with debug , however, my code failed with
[0] PetscCommDuplicate(): Using internal PETSc communicator 0 5
[0] PetscCommDuplicate(): Using internal PETSc communicator 0 5
[0] PetscCommDuplicate(): Using internal PETSc communicator 0 5
[0] PetscCommDuplicate(): Duplicating a communicator 1 6 max tags = 1073741823
The assert subroutine failed: context >= _min_context, file
/project/sprelhya/build/rhyas002a/src/ppe/poe/include/mpi_cc_inlines.h, line 314
Any suggestion?
> '-mat_mumps_icntl_4 4' which will display mumps info
>
>
>
> Use a debugger to see where the program terminates.
>
>
>
> Hong
>
>
>
> On Tue, Dec 6, 2011 at 10:54 PM, Gong Ding <gdiso at ustc.edu> wrote:
>
> > I can only analysis the core file at present.
>
> > It reported
>
> > Program terminated with signal 11, Segmentation fault.
>
> > #0 0x000000010217319c in dmumps_462 ()
>
> >
>
> > Run program with -info, the output is
>
> >
>
> > Entering DMUMPS driver with JOB, N, NZ = 1 49 217
>
> >
>
> > DMUMPS 4.10.0
>
> > L U Solver for unsymmetric matrices
>
> > Type of parallelism: Working host
>
> >
>
> > ****** ANALYSIS STEP ********
>
> >
>
> > Resetting candidate strategy to 0 because NSLAVES=1
>
> >
>
> > ... Structural symmetry (in percent)= 100
>
> > Density: NBdense, Average, Median = 0 3 3
>
> > ... No column permutation
>
> > Ordering based on AMF
>
> >
>
> > Program terminated with signal SIGSEGV, Segmentation fault.
>
> >
>
> > The example in src/ksp/ksp/examples/tests/ex3 run ok -- it solves a small
>
> > linear problem.
>
> > Also, if I reduce the matrix size to ~20x20, no crash happen.
>
> >
>
> > I'll try to provide more information.
>
> >
>
> >
>
> >
>
> >
>
> >
>
> > On Tue, Dec 6, 2011 at 12:52 AM, Gong Ding <gdiso at ustc.edu> wrote:
>
> >>
>
> >> Bad news.
>
> >> MUMPS still crash on AIX/PPC.
>
> >> I tried petsc-3.1, 3.2 and dev. All of them crash with ERROR: 0031-250
>
> >> task 0: Segmentation fault.
>
> >> I guess the interface of petsc-mumps has some problem, i.e. alignment.
>
> >
>
> >
>
> > Please get a stack trace with the debugger.
>
> >
>
> > Matt
>
> >
>
> >>>
>
> >>> Yes, petsc-dev works. Please merge the update to petsc-3.2.
>
> >>
>
> >>
>
> >> Parmetis 4.0 was released after PETSc 3.2 so it will not be updated to
>
> >> work with it. In general, software projects don't update a previous release
>
> >> on a dependent package update. Please use petsc-dev or wait until PETSc 3.3.
>
> >
>
> >
>
> >
>
> >
>
> > --
>
> > 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
>
>
More information about the petsc-users
mailing list