[petsc-dev] ex43_3 stokes problem FAIL
Mark F. Adams
mark.adams at columbia.edu
Tue May 15 14:48:42 CDT 2012
I do not seem to be able to reproduce these errors. This is what I see:
~/Codes/petsc-dev/src/ksp/ksp/examples/tutorials>make runex49
30c30
< 28 KSP Residual norm 0.000350713
---
> 28 KSP Residual norm 0.000350714
41c41
< 39 KSP Residual norm 9.38531e-05
---
> 39 KSP Residual norm 9.38532e-05
43c43
< 41 KSP Residual norm 5.55156e-05
---
> 41 KSP Residual norm 5.55157e-05
/Users/markadams/Codes/petsc-dev/src/ksp/ksp/examples/tutorials
Possible problem with with ex49_1, diffs above
=========================================
~/Codes/petsc-dev/src/ksp/ksp/examples/tutorials>make runex49_3
~/Codes/petsc-dev/src/ksp/ksp/examples/tutorials>make runex49_2
~/Codes/petsc-dev/src/ksp/ksp/examples/tutorials>
Mark
On May 11, 2012, at 10:34 PM, Barry Smith wrote:
>
> Same problem with ex49, what the fuck are you guys thinking?
>
>
> [2]PETSC ERROR: --------------------- Error Message ------------------------------------
>> [2]PETSC ERROR: Arguments are incompatible!
>> [2]PETSC ERROR: Local size 315 not compatible with block size 2!
>> [2]PETSC ERROR: ------------------------------------------------------------------------
>> [2]PETSC ERROR: Petsc Development HG revision: 5df5da2fbee0018d6214063ab12ea49c0847b4f8 HG Date: Fri May 11 13:26:59 2012 -0500
>> [2]PETSC ERROR: See docs/changes/index.html for recent updates.
>> [2]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
>> [2]PETSC ERROR: See docs/index.html for manual pages.
>> [2]PETSC ERROR: ------------------------------------------------------------------------
>> [2]PETSC ERROR: ./ex49 on a arch-gnu named barry-smiths-macbook-pro.local by barrysmith Fri May 11 18:47:27 2012
>> [2]PETSC ERROR: Libraries linked from /Users/barrysmith/Src/petsc-dev/arch-gnu/lib
>> [2]PETSC ERROR: Configure run at Fri May 11 18:30:02 2012
>> [2]PETSC ERROR: Configure options --download-blacs --download-fftw --download-hypre --download-libyaml --download-metis --download-ml --download-mpich --download-mumps --download-parmetis --download-ptscotch --download-scalapack --download-superlu --download-superlu_dist --download-yaml --with-ams-dir=/Users/barrysmith/Src/ams-dev --with-java --with-openmp --with-server --with-shared-libraries PETSC_ARCH=arch-gnu
>> [2]PETSC ERROR: ------------------------------------------------------------------------
>> [2]PETSC ERROR: PetscLayoutSetBlockSize() line 459 in /Users/barrysmith/Src/petsc-dev/src/vec/vec/impls/mpi/pmap.c
>> [2]PETSC ERROR: MatSetBlockSizes() line 6843 in /Users/barrysmith/Src/petsc-dev/src/mat/interface/matrix.c
>> [2]PETSC ERROR: MatGetSubMatrices_MPIAIJ_Local() line 1312 in /Users/barrysmith/Src/petsc-dev/src/mat/impls/aij/mpi/mpiov.c
>> [2]PETSC ERROR: MatGetSubMatrices_MPIAIJ() line 786 in /Users/barrysmith/Src/petsc-dev/src/mat/impls/aij/mpi/mpiov.c
>> [2]PETSC ERROR: MatGetSubMatrices() line 6498 in /Users/barrysmith/Src/petsc-dev/src/mat/interface/matrix.c
>> [2]PETSC ERROR: MatGetSubMatrix_MPIAIJ_Private() line 3554 in /Users/barrysmith/Src/petsc-dev/src/mat/impls/aij/mpi/mpiaij.c
>> [2]PETSC ERROR: MatGetSubMatrix_MPIAIJ() line 3515 in /Users/barrysmith/Src/petsc-dev/src/mat/impls/aij/mpi/mpiaij.c
>> [2]PETSC ERROR: MatGetSubMatrix() line 7344 in /Users/barrysmith/Src/petsc-dev/src/mat/interface/matrix.c
>> [2]PETSC ERROR: DMDABCApplySymmetricCompression() line 1348 in src/ksp/ksp/examples/tutorials/ex49.c
>> [2]PETSC ERROR: solve_elasticity_2d() line 1080 in src/ksp/ksp/examples/tutorials/ex49.c
>> [2]PETSC ERROR: main() line 1138 in src/ksp/ksp/examples/tutorials/ex49.c
>> application called MPI_Abort(MPI_COMM_WORLD, 75) - process 2
>> [cli_2]: aborting job:
>
> On May 11, 2012, at 9:31 PM, Barry Smith wrote:
>
>>
>> Please run make alltests before pushing changes!
>>
>> Who ever pushed this broken code please fix.
>> runex43_3:
>> -@${MPIEXEC} -n 4 ./ex43 -stokes_ksp_type gcr -stokes_ksp_gcr_restart 60 -stokes_ksp_norm_type unpreconditioned -stokes_ksp_rtol 1e-8 -c_str 3 -sinker_eta0 1.0 -sinker_eta1 100 -sinker_dx 0.4 -sinker_dy 0.3 -mx 128 -my 128 -stokes_ksp_monitor_short -stokes_pc_type mg -stokes_mg_levels_pc_type fieldsplit -stokes_pc_mg_galerkin -stokes_mg_levels_pc_fieldsplit_block_size 3 -stokes_mg_levels_pc_fieldsplit_0_fields 0,1 -stokes_mg_levels_pc_fieldsplit_1_fields 2 -stokes_mg_levels_fieldsplit_0_pc_type sor -stokes_mg_levels_fieldsplit_1_pc_type sor -stokes_mg_levels_ksp_type chebyshev -stokes_mg_levels_ksp_max_it 1 -stokes_mg_levels_ksp_chebyshev_estimate_eigenvalues 0,0.2,0,1.1 -stokes_pc_mg_levels 4 -stokes_ksp_view > ex43_3.tmp 2>&1; \
>> ${DIFF} output/ex43_3.out ex43_3.tmp || echo ${PWD} "\nPossible problem with with ex43_3, diffs above \n========================================="; \
>> ${RM} -f ex43_3.tmp
>>
>> [2]PETSC ERROR: --------------------- Error Message ------------------------------------
>>> [2]PETSC ERROR: Arguments are incompatible!
>>> [2]PETSC ERROR: Local size 8320 not compatible with block size 3!
>>> [2]PETSC ERROR: ------------------------------------------------------------------------
>>> [2]PETSC ERROR: Petsc Development HG revision: 5df5da2fbee0018d6214063ab12ea49c0847b4f8 HG Date: Fri May 11 13:26:59 2012 -0500
>>> [2]PETSC ERROR: See docs/changes/index.html for recent updates.
>>> [2]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
>>> [2]PETSC ERROR: See docs/index.html for manual pages.
>>> [2]PETSC ERROR: ------------------------------------------------------------------------
>>> [2]PETSC ERROR: ./ex43 on a arch-gnu named barry-smiths-macbook-pro.local by barrysmith Fri May 11 21:28:34 2012
>>> [2]PETSC ERROR: Libraries linked from /Users/barrysmith/Src/petsc-dev/arch-gnu/lib
>>> [2]PETSC ERROR: Configure run at Fri May 11 18:30:02 2012
>>> [2]PETSC ERROR: Configure options --download-blacs --download-fftw --download-hypre --download-libyaml --download-metis --download-ml --download-mpich --download-mumps --download-parmetis --download-ptscotch --download-scalapack --download-superlu --download-superlu_dist --download-yaml --with-ams-dir=/Users/barrysmith/Src/ams-dev --with-java --with-openmp --with-server --with-shared-libraries PETSC_ARCH=arch-gnu
>>> [2]PETSC ERROR: ------------------------------------------------------------------------
>>> [2]PETSC ERROR: PetscLayoutSetBlockSize() line 459 in /Users/barrysmith/Src/petsc-dev/src/vec/vec/impls/mpi/pmap.c
>>> [2]PETSC ERROR: MatSetBlockSizes() line 6843 in /Users/barrysmith/Src/petsc-dev/src/mat/interface/matrix.c
>>> [2]PETSC ERROR: MatGetSubMatrices_MPIAIJ_Local() line 1312 in /Users/barrysmith/Src/petsc-dev/src/mat/impls/aij/mpi/mpiov.c
>>> [2]PETSC ERROR: MatGetSubMatrices_MPIAIJ() line 786 in /Users/barrysmith/Src/petsc-dev/src/mat/impls/aij/mpi/mpiov.c
>>> [2]PETSC ERROR: MatGetSubMatrices() line 6498 in /Users/barrysmith/Src/petsc-dev/src/mat/interface/matrix.c
>>> [2]PETSC ERROR: MatGetSubMatrix_MPIAIJ_Private() line 3554 in /Users/barrysmith/Src/petsc-dev/src/mat/impls/aij/mpi/mpiaij.c
>>> [2]PETSC ERROR: MatGetSubMatrix_MPIAIJ() line 3515 in /Users/barrysmith/Src/petsc-dev/src/mat/impls/aij/mpi/mpiaij.c
>>> [2]PETSC ERROR: MatGetSubMatrix() line 7344 in /Users/barrysmith/Src/petsc-dev/src/mat/interface/matrix.c
>>> [2]PETSC ERROR: PCSetUp_FieldSplit() line 445 in /Users/barrysmith/Src/petsc-dev/src/ksp/pc/impls/fieldsplit/fieldsplit.c
>>> [2]PETSC ERROR: PCSetUp() line 832 in /Users/barrysmith/Src/petsc-dev/src/ksp/pc/interface/precon.c
>>> [2]PETSC ERROR: KSPSetUp() line 278 in /Users/barrysmith/Src/petsc-dev/src/ksp/ksp/interface/itfunc.c
>>> [2]PETSC ERROR: PCSetUp_MG() line 678 in /Users/barrysmith/Src/petsc-dev/src/ksp/pc/impls/mg/mg.c
>>> [2]PETSC ERROR: PCSetUp() line 832 in /Users/barrysmith/Src/petsc-dev/src/ksp/pc/interface/precon.c
>>> [2]PETSC ERROR: KSPSetUp() line 278 in /Users/barrysmith/Src/petsc-dev/src/ksp/ksp/interface/itfunc.c
>>> [2]PETSC ERROR: KSPSolve() line 402 in /Users/barrysmith/Src/petsc-dev/src/ksp/ksp/interface/itfunc.c
>>> [2]PETSC ERROR: solve_stokes_2d_coupled() line 1502 in src/ksp/ksp/examples/tutorials/ex43.c
>>> [2]PETSC ERROR: main() line 1621 in src/ksp/ksp/examples/tutorials/ex43.c
>>> application called MPI_Abort(MPI_COMM_WORLD, 75) - process 2
>>> [cli_2]: aborting job:
>>> application called MPI_Abort(MPI_COMM_WORLD, 75) - process 2
>> /Users/barrysmith/Src/petsc-dev/src/ksp/ksp/examples/tutorials
>> Possible problem with with ex43_3, diffs above
>> =========================================
>>
>>
>
>
More information about the petsc-dev
mailing list