Hi,
I am solving simple poisson equation using IPDG. I have two questions when I
try to model my equation on step 40. I hope someone can help me with this.
1. I get this error when I say make_flux_sparsity_pattern.
An error occurred in line <3174> of file
</gpfs/runtime/opt/dealii/7.0.0/include/deal.II/dofs/dof_accessor.templates.h>
in function
void dealii::DoFCellAccessor<DH>::get_dof_indices(std::vector<unsigned
int>&) const [with DH = dealii::DoFHandler<2, 2>]
The violated condition was:
this->is_artificial() == false
The name and call sequence of the exception was:
ExcMessage ("Can't ask for DoF indices on artificial cells.")
Additional Information:
Can't ask for DoF indices on artificial cells.
Stacktrace:
-----------
#0 ./Poisson: dealii::DoFCellAccessor<dealii::DoFHandler<2, 2>
>::get_dof_indices(std::vector<unsigned int, std::allocator<unsigned int>
>&) const
#1 /gpfs/runtime/opt/dealii/7.0.0/lib/libdeal_II.g.so.7.0.0: void
dealii::DoFTools::make_flux_sparsity_pattern<dealii::DoFHandler<2, 2>,
dealii::CompressedSimpleSparsityPattern>(dealii::DoFHandler<2, 2> const&,
dealii::CompressedSimpleSparsityPattern&, dealii::ConstraintMatrix const&,
bool, unsigned int)
#2 /gpfs/runtime/opt/dealii/7.0.0/lib/libdeal_II.g.so.7.0.0: void
dealii::DoFTools::make_flux_sparsity_pattern<dealii::DoFHandler<2, 2>,
dealii::CompressedSimpleSparsityPattern>(dealii::DoFHandler<2, 2> const&,
dealii::CompressedSimpleSparsityPattern&)
#3 ./Poisson: Poisson<2>::setup_system()
#4 ./Poisson: Poisson<2>::run()
#5 ./Poisson: main
2. On a related note, I don't use the sparsity pattern and just reinit
system_matrix using dof_handler.n_dofs(). But when I run this it assembles
fine and gives me the right results when I run with 1 processor but when I
use more than 1 processor, I get the following error:
[1]PETSC ERROR:
------------------------------------------------------------------------
[1]PETSC ERROR: MatILUFactorSymbolic_SeqAIJ() line 1248 in
src/mat/impls/aij/seq/aijfact.c
[1]PETSC ERROR: MatILUFactorSymbolic() line 5314 in
src/mat/interface/matrix.c
[1]PETSC ERROR: PCSetUp_ILU() line 293 in src/ksp/pc/impls/factor/ilu/ilu.c
[1]PETSC ERROR: PCSetUp() line 794 in src/ksp/pc/interface/precon.c
[1]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c
[1]PETSC ERROR: PCSetUpOnBlocks_BJacobi_Singleblock() line 753 in
src/ksp/pc/impls/bjacobi/bjacobi.c
[1]PETSC ERROR: PCSetUpOnBlocks() line 827 in src/ksp/pc/interface/precon.c
[1]PETSC ERROR: KSPSetUpOnBlocks() line 159 in
src/ksp/ksp/interface/itfunc.c
[1]PETSC ERROR: KSPSolve() line 354 in src/ksp/ksp/interface/itfunc.c
[0]PETSC ERROR: --------------------- Error Message
------------------------------------
[0]PETSC ERROR: Nonconforming object sizes!
[0]PETSC ERROR: Incompatible partition of A (4500) and xx (2250)!
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 3.0.0, Patch 12, Tue Mar 16 23:20:08
CDT 2010
[0]PETSC ERROR: See docs/changes/index.html for recent updates.
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[0]PETSC ERROR: See docs/index.html for manual pages.
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: ./Poisson on a linux-gnu named smp009 by gk Thu May 19
14:59:07 2011
[0]PETSC ERROR: Libraries linked from
/gpfs/runtime/opt/petsc/3.0.0-p12/linux-gnu-cxx-opt/lib
[0]PETSC ERROR: Configure run at Fri May 13 12:16:12 2011
[0]PETSC ERROR: Configure options --download-superlu_dist=1
--download-hypre=1 --with-debugging=0 --with-precision=double --with-dynamic
--with-shared=1 --with-cc=mpicc --with-fc=mpif77 --with-cxx=mpiCC
--with-clanguage=c++
--with-blas-lib=/gpfs/runtime/opt/gotoblas2/1.13/lib/libgoto2.so
--with-lapack-lib=/gpfs/runtime/opt/gotoblas2/1.13/lib/libgoto2.so
--with-mumps-dir=/gpfs/runtime/opt/mumps/4.9.2
--with-parmetis-dir=/gpfs/runtime/opt/parmetis/3.1.1
--with-scalapack-dir=/gpfs/runtime/opt/gotoblas2/1.13
--with-blacs-lib="[/gpfs/runtime/opt/gotoblas2/1.13/lib/libblacsCinit.a,/gpfs/runtime/opt/gotoblas2/1.13/lib/libblacs.a]"
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: MatMult_MPIAIJ() line 720 in src/mat/impls/aij/mpi/mpiaij.c
[0]PETSC ERROR: MatMult() line 1784 in src/mat/interface/matrix.c
[0]PETSC ERROR: KSPInitialResidual() line 54 in
src/ksp/ksp/interface/itres.c
[0]PETSC ERROR: KSPSolve_GMRES() line 241 in src/ksp/ksp/impls/gmres/gmres.c
[0]PETSC ERROR: KSPSolve() line 385 in src/ksp/ksp/interface/itfunc.c
terminate called after throwing an instance of
'dealii::PETScWrappers::SolverBase::ExcPETScError'
terminate called after throwing an instance of
'dealii::PETScWrappers::SolverBase::ExcPETScError'
what(): --------------------------------------------------------
An error occurred in line <123> of file
</gpfs/runtime/opt/dealii/7.0.0/source/lac/petsc_solver.cc> in function
void dealii::PETScWrappers::SolverBase::solve(const
dealii::PETScWrappers::MatrixBase&, dealii::PETScWrappers::VectorBase&,
const dealii::PETScWrappers::VectorBase&, const
dealii::PETScWrappers::PreconditionerBase&)
The violated condition was:
ierr == 0
The name and call sequence of the exception was:
ExcPETScError(ierr)
Additional Information:
An error with error number 73 occured while calling a PETSc function
--------------------------------------------------------
what(): --------------------------------------------------------
An error occurred in line <123> of file
</gpfs/runtime/opt/dealii/7.0.0/source/lac/petsc_solver.cc> in function
void dealii::PETScWrappers::SolverBase::solve(const
dealii::PETScWrappers::MatrixBase&, dealii::PETScWrappers::VectorBase&,
const dealii::PETScWrappers::VectorBase&, const
dealii::PETScWrappers::PreconditionerBase&)
The violated condition was:
ierr == 0
The name and call sequence of the exception was:
ExcPETScError(ierr)
Additional Information:
An error with error number 60 occured while calling a PETSc function
--------------------------------------------------------
[smp009:25583] *** Process received signal ***
[smp009:25583] Signal: Aborted (6)
[smp009:25583] Signal code: (-6)
[smp009:25582] *** Process received signal ***
[smp009:25582] Signal: Aborted (6)
[smp009:25582] Signal code: (-6)
[smp009:25583] [ 0] /lib64/libpthread.so.0 [0x37f680eb10]
[smp009:25583] [ 1] /lib64/libc.so.6(gsignal+0x35) [0x37f6030265]
[smp009:25583] [ 2] /lib64/libc.so.6(abort+0x110) [0x37f6031d10]
[smp009:25583] [ 3]
/gpfs/runtime/opt/gcc/4.5.1/lib64/libstdc++.so.6(_ZN9__gnu_cxx27__verbose_terminate_handlerEv+0x11d)
[0x2ad2d1b4e6fd]
[smp009:25583] [ 4] /gpfs/runtime/opt/gcc/4.5.1/lib64/libstdc++.so.6
[0x2ad2d1b4c916]
[smp009:25583] [ 5] /gpfs/runtime/opt/gcc/4.5.1/lib64/libstdc++.so.6
[0x2ad2d1b4c943]
[smp009:25583] [ 6] /gpfs/runtime/opt/gcc/4.5.1/lib64/libstdc++.so.6
[0x2ad2d1b4ca4e]
[smp009:25583] [ 7] /gpfs/runtime/opt/dealii/7.0.0/lib/libdeal_II.g.so.7.0.0
[0x2ad2caa8c6e5]
[smp009:25583] [ 8]
/gpfs/runtime/opt/dealii/7.0.0/lib/libdeal_II.g.so.7.0.0(_ZN6dealii13PETScWrappers10SolverBase5solveERKNS0_10MatrixBaseERNS0_10VectorBaseERKS5_RKNS0_18PreconditionerBaseE+0x42d)
[0x2ad2caa8b231]
[smp009:25583] [ 9] ./Poisson(_ZN14PoissonLi2EE5solveEi+0x12d) [0x4430dd]
[smp009:25583] [10] ./Poisson(_ZN14PoissonILi2EE3runEv+0x39b) [0x43f31b]
[smp009:25583] [11] ./Poisson(main+0x73) [0x437947]
[smp009:25583] [12] /lib64/libc.so.6(__libc_start_main+0xf4) [0x37f601d994]
[smp009:25583] [13] ./Poisson [0x437819]
[smp009:25583] *** End of error message ***
[smp009:25582] [ 0] /lib64/libpthread.so.0 [0x37f680eb10]
[smp009:25582] [ 1] /lib64/libc.so.6(gsignal+0x35) [0x37f6030265]
[smp009:25582] [ 2] /lib64/libc.so.6(abort+0x110) [0x37f6031d10]
[smp009:25582] [ 3]
/gpfs/runtime/opt/gcc/4.5.1/lib64/libstdc++.so.6(_ZN9__gnu_cxx27__verbose_terminate_handlerEv+0x11d)
[0x2b03b3ad66fd]
[smp009:25582] [ 4] /gpfs/runtime/opt/gcc/4.5.1/lib64/libstdc++.so.6
[0x2b03b3ad4916]
[smp009:25582] [ 5] /gpfs/runtime/opt/gcc/4.5.1/lib64/libstdc++.so.6
[0x2b03b3ad4943]
[smp009:25582] [ 6] /gpfs/runtime/opt/gcc/4.5.1/lib64/libstdc++.so.6
[0x2b03b3ad4a4e]
[smp009:25582] [ 7] /gpfs/runtime/opt/dealii/7.0.0/lib/libdeal_II.g.so.7.0.0
[0x2b03aca146e5]
[smp009:25582] [ 8]
/gpfs/runtime/opt/dealii/7.0.0/lib/libdeal_II.g.so.7.0.0(_ZN6dealii13PETScWrappers10SolverBase5solveERKNS0_10MatrixBaseERNS0_10VectorBaseERKS5_RKNS0_18PreconditionerBaseE+0x42d)
[0x2b03aca13231]
[smp009:25582] [ 9] ./Poisson(_ZN14PoissonILi2EE5solveEi+0x12d) [0x4430dd]
[smp009:25582] [10] ./Poisson(_ZN14PoissonILi2EE3runEv+0x39b) [0x43f31b]
[smp009:25582] [11] ./Poisson(main+0x73) [0x437947]
[smp009:25582] [12] /lib64/libc.so.6(__libc_start_main+0xf4) [0x37f601d994]
[smp009:25582] [13] ./Poisson [0x437819]
--
Thanks,
gk
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
