[petsc-dev] The errors when running petsc/src/snes/examples/tutorials/ex19.c in PETSc

Peter Brune prbrune at gmail.com
Thu May 9 14:33:57 CDT 2013


The fix is now in next.

- Peter


On Thu, May 9, 2013 at 3:09 AM, Lulu Liu <lulu.liu at kaust.edu.sa> wrote:

> I replace  the file "nasm.c" in branch *next* with "nasm.c" in branch *
> prbrune/snes-nasmview*, and it works now!
>
> Thank you!
>
>
> On Wed, May 8, 2013 at 11:12 PM, Peter Brune <prbrune at gmail.com> wrote:
>
>> I've made -snes_view for NASM a bit more full featured.  Once again, I
>> cannot reproduce your problem, but it needed to be fixed up anyways.
>>
>> It's in branch prbrune/snes-nasmview.  Trying it (to see if your problem
>> still exists) would be appreciated.
>>
>> - Peter
>>
>>
>> On Tue, May 7, 2013 at 2:05 PM, Lulu Liu <lulu.liu at kaust.edu.sa> wrote:
>>
>>> I used next branch and got these errors. However, I pulled master and
>>> test ex19.c, and got the following errors:
>>>  0 SNES Function norm 5381.22
>>>   1 SNES Function norm 3104.04
>>>   2 SNES Function norm 777.712
>>>   3 SNES Function norm 45.9649
>>>   4 SNES Function norm 0.195115
>>>   5 SNES Function norm 6.28208e-06
>>> SNES Object: 2 MPI processes
>>>   type: newtonls
>>>   maximum iterations=50, maximum function evaluations=10000
>>>   tolerances: relative=1e-08, absolute=1e-50, solution=1e-08
>>>   total number of linear solver iterations=47
>>>   total number of function evaluations=0
>>>   SNESLineSearch Object:   2 MPI processes
>>>     type: bt
>>>       interpolation: cubic
>>>       alpha=1.000000e-04
>>>     maxstep=1.000000e+08, minlambda=1.000000e-12
>>>     tolerances: relative=1.000000e-08, absolute=1.000000e-15,
>>> lambda=1.000000e-08
>>>     maximum iterations=40
>>>   SNES Object:  (npc_)   2 MPI processes
>>>     type: nasm
>>>       Nonlinear Additive Schwarz: total subdomain blocks = 2
>>>       Nonlinear Additive Schwarz: restriction/interpolation type - BASIC
>>>       Nonlinear Additive Schwarz: subSNES iterations: 275635578 subKSP
>>> iterations: -2128747920
>>>       [0] number of local blocks = 1
>>>       [1] number of local blocks = 1
>>>       Local SNES objects:
>>>       SNES Object:      (npc_sub_)       1 MPI processes
>>>         type: newtonls
>>>         maximum iterations=50, maximum function evaluations=10000
>>>         tolerances: relative=1e-08, absolute=1e-50, solution=1e-08
>>>         total number of linear solver iterations=1
>>>         total number of function evaluations=2
>>>         SNESLineSearch Object:        (npc_sub_)         1 MPI processes
>>>           type: bt
>>>             interpolation: cubic
>>>             alpha=1.000000e-04
>>>           maxstep=1.000000e+08, minlambda=1.000000e-12
>>>           tolerances: relative=1.000000e-08, absolute=1.000000e-15,
>>> lambda=1.000000e-08
>>>           maximum iterations=40
>>>         KSP Object:        (npc_sub_)         1 MPI processes
>>>           type: preonly
>>>           maximum iterations=10000, initial guess is zero
>>>           tolerances:  relative=1e-05, absolute=1e-50, divergence=10000
>>>           left preconditioning
>>>           using NONE norm type for convergence test
>>>         PC Object:        (npc_sub_)         1 MPI processes
>>>           type: lu
>>>             LU: out-of-place factorization
>>>             tolerance for zero pivot 2.22045e-14
>>>             matrix ordering: nd
>>>             factor fill ratio given 5, needed 2.95194
>>>               Factored matrix follows:
>>>                 Matrix Object:                 1 MPI processes
>>>                   type: seqaij
>>>                   rows=468, cols=468, bs=4
>>>                   package used to perform factorization: petsc
>>>                   total: nonzeros=25552, allocated nonzeros=25552
>>>                   total number of mallocs used during MatSetValues calls
>>> =0
>>>                     using I-node routines: found 117 nodes, limit used
>>> is 5
>>>           linear system matrix = precond matrix:
>>>           Matrix Object:           1 MPI processes
>>>             type: seqaij
>>>             rows=468, cols=468, bs=4
>>>             total: nonzeros=8656, allocated nonzeros=8656
>>>             total number of mallocs used during MatSetValues calls =0
>>>               using I-node routines: found 117 nodes, limit used is 5
>>> [1]PETSC ERROR: PetscCommDuplicate() line 183 in
>>> /Users/liul/soft/petsc/src/sys/objects/tagm.c
>>> [1]PETSC ERROR: PetscSynchronizedFlush() line 450 in
>>> /Users/liul/soft/petsc/src/sys/fileio/mprint.c
>>> [1]PETSC ERROR: PetscViewerFlush_ASCII() line 135 in
>>> /Users/liul/soft/petsc/src/sys/classes/viewer/impls/ascii/filev.c
>>> [1]PETSC ERROR: PetscViewerFlush() line 30 in
>>> /Users/liul/soft/petsc/src/sys/classes/viewer/interface/flush.c
>>> [1]PETSC ERROR: SNESView_NASM() line 222 in
>>> /Users/liul/soft/petsc/src/snes/impls/nasm/nasm.c
>>>  [1]PETSC ERROR: SNESView() line 251 in
>>> /Users/liul/soft/petsc/src/snes/interface/snes.c
>>> [1]PETSC ERROR: SNESView() line 336 in
>>> /Users/liul/soft/petsc/src/snes/interface/snes.c
>>> [1]PETSC ERROR: SNESSolve() line 3677 in
>>> /Users/liul/soft/petsc/src/snes/interface/snes.c
>>>  [1]PETSC ERROR: main() line 157 in src/snes/examples/tutorials/ex19.c
>>> application called MPI_Abort(MPI_COMM_WORLD, 872626446) - process 1
>>> [cli_1]: aborting job:
>>> application called MPI_Abort(MPI_COMM_WORLD, 872626446) - process 1
>>>
>>>
>>>
>>> On Tue, May 7, 2013 at 9:14 PM, Peter Brune <prbrune at gmail.com> wrote:
>>>
>>>> I can't reproduce.  Have you pulled master recently?  Those line
>>>> numbers don't correspond to anything, and that file hasn't changed very
>>>> much in a while.
>>>>
>>>> - Peter
>>>>
>>>>
>>>> On Tue, May 7, 2013 at 11:31 AM, Lulu Liu <lulu.liu at kaust.edu.sa>wrote:
>>>>
>>>>> mpirun -n 2  ./ex19 -da_refine 2 -da_overlap 2 -snes_monitor_short
>>>>> -snes_type aspin  -grashof 4e4 -lidvelocity 100   -snes_view
>>>>>
>>>>>  0 SNES Function norm 5381.22
>>>>>   1 SNES Function norm 3104.04
>>>>>   2 SNES Function norm 777.712
>>>>>   3 SNES Function norm 45.9649
>>>>>   4 SNES Function norm 0.195115
>>>>>   5 SNES Function norm 6.28208e-06
>>>>> SNES Object: 2 MPI processes
>>>>>   type: newtonls
>>>>>   maximum iterations=50, maximum function evaluations=10000
>>>>>   tolerances: relative=1e-08, absolute=1e-50, solution=1e-08
>>>>>   total number of linear solver iterations=47
>>>>>   total number of function evaluations=0
>>>>>   SNESLineSearch Object:   2 MPI processes
>>>>>     type: bt
>>>>>       interpolation: cubic
>>>>>       alpha=1.000000e-04
>>>>>     maxstep=1.000000e+08, minlambda=1.000000e-12
>>>>>     tolerances: relative=1.000000e-08, absolute=1.000000e-15,
>>>>> lambda=1.000000e-08
>>>>>     maximum iterations=40
>>>>>   SNES Object:  (npc_)   2 MPI processes
>>>>>     type: nasm
>>>>>       Nonlinear Additive Schwarz: total subdomain blocks = 2
>>>>>       Nonlinear Additive Schwarz: restriction/interpolation type -
>>>>> BASIC
>>>>>       Nonlinear Additive Schwarz: subSNES iterations: 114022506 subKSP
>>>>> iterations: 564006000
>>>>>       [0] number of local blocks = 1
>>>>>       [1] number of local blocks = 1
>>>>>       Local SNES objects:
>>>>>       SNES Object:      (npc_sub_)       1 MPI processes
>>>>>         type: newtonls
>>>>>         maximum iterations=50, maximum function evaluations=10000
>>>>>         tolerances: relative=1e-08, absolute=1e-50, solution=1e-08
>>>>>         total number of linear solver iterations=1
>>>>>         total number of function evaluations=3
>>>>>         SNESLineSearch Object:        (npc_sub_)         1 MPI
>>>>> processes
>>>>>           type: bt
>>>>>             interpolation: cubic
>>>>>             alpha=1.000000e-04
>>>>>           maxstep=1.000000e+08, minlambda=1.000000e-12
>>>>>           tolerances: relative=1.000000e-08, absolute=1.000000e-15,
>>>>> lambda=1.000000e-08
>>>>>           maximum iterations=40
>>>>>         KSP Object:        (npc_sub_)         1 MPI processes
>>>>>           type: preonly
>>>>>           maximum iterations=10000, initial guess is zero
>>>>>           tolerances:  relative=1e-05, absolute=1e-50, divergence=10000
>>>>>           left preconditioning
>>>>>           using NONE norm type for convergence test
>>>>>         PC Object:        (npc_sub_)         1 MPI processes
>>>>>           type: lu
>>>>>             LU: out-of-place factorization
>>>>>             tolerance for zero pivot 2.22045e-14
>>>>>             matrix ordering: nd
>>>>>             factor fill ratio given 5, needed 2.95194
>>>>>               Factored matrix follows:
>>>>>                 Matrix Object:                 1 MPI processes
>>>>>                   type: seqaij
>>>>>                   rows=468, cols=468, bs=4
>>>>>                   package used to perform factorization: petsc
>>>>>                   total: nonzeros=25552, allocated nonzeros=25552
>>>>>                   total number of mallocs used during MatSetValues
>>>>> calls =0
>>>>>                     using I-node routines: found 117 nodes, limit used
>>>>> is 5
>>>>>           linear system matrix = precond matrix:
>>>>>           Matrix Object:           1 MPI processes
>>>>>             type: seqaij
>>>>>             rows=468, cols=468, bs=4
>>>>>             total: nonzeros=8656, allocated nonzeros=8656
>>>>>             total number of mallocs used during MatSetValues calls =0
>>>>>               using I-node routines: found 117 nodes, limit used is 5
>>>>> [1]PETSC ERROR: PetscCommDuplicate() line 183 in
>>>>> /Users/liul/soft/petsc/src/sys/objects/tagm.c
>>>>> [1]PETSC ERROR: PetscSynchronizedFlush() line 450 in
>>>>> /Users/liul/soft/petsc/src/sys/fileio/mprint.c
>>>>> [1]PETSC ERROR: PetscViewerFlush_ASCII() line 135 in
>>>>> /Users/liul/soft/petsc/src/sys/classes/viewer/impls/ascii/filev.c
>>>>> [1]PETSC ERROR: PetscViewerFlush() line 30 in
>>>>> /Users/liul/soft/petsc/src/sys/classes/viewer/interface/flush.c
>>>>> [1]PETSC ERROR: SNESView_NASM() line 240 in
>>>>> /Users/liul/soft/petsc/src/snes/impls/nasm/nasm.c
>>>>> [1]PETSC ERROR: SNESView() line 251 in
>>>>> /Users/liul/soft/petsc/src/snes/interface/snes.c
>>>>> [1]PETSC ERROR: SNESView() line 336 in
>>>>> /Users/liul/soft/petsc/src/snes/interface/snes.c
>>>>> [1]PETSC ERROR: SNESSolve() line 3804 in
>>>>> /Users/liul/soft/petsc/src/snes/interface/snes.c
>>>>> [1]PETSC ERROR: main() line 157 in src/snes/examples/tutorials/ex19.c
>>>>> application called MPI_Abort(MPI_COMM_WORLD, 872626446) - process 1
>>>>> [cli_1]: aborting job:
>>>>> application called MPI_Abort(MPI_COMM_WORLD, 872626446) - process 1
>>>>>  --
>>>>> Best wishes,
>>>>> Lulu Liu
>>>>> Applied Mathematics and Computational Science
>>>>> King Abdullah University of Science and Technology
>>>>> Tel:+966-0544701599
>>>>>
>>>>> ------------------------------
>>>>> This message and its contents, including attachments are intended
>>>>> solely for the original recipient. If you are not the intended recipient or
>>>>> have received this message in error, please notify me immediately and
>>>>> delete this message from your computer system. Any unauthorized use or
>>>>> distribution is prohibited. Please consider the environment before printing
>>>>> this email.
>>>>
>>>>
>>>>
>>>
>>>
>>> --
>>> Best wishes,
>>> Lulu Liu
>>> Applied Mathematics and Computational Science
>>> King Abdullah University of Science and Technology
>>> Tel:+966-0544701599
>>>
>>> ------------------------------
>>> This message and its contents, including attachments are intended solely
>>> for the original recipient. If you are not the intended recipient or have
>>> received this message in error, please notify me immediately and delete
>>> this message from your computer system. Any unauthorized use or
>>> distribution is prohibited. Please consider the environment before printing
>>> this email.
>>>
>>
>>
>
>
> --
> Best wishes,
> Lulu Liu
> Applied Mathematics and Computational Science
> King Abdullah University of Science and Technology
> Tel:+966-0544701599
>
> ------------------------------
> This message and its contents, including attachments are intended solely
> for the original recipient. If you are not the intended recipient or have
> received this message in error, please notify me immediately and delete
> this message from your computer system. Any unauthorized use or
> distribution is prohibited. Please consider the environment before printing
> this email.
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130509/53ee2ad9/attachment.html>


More information about the petsc-dev mailing list