← Back to team overview

dolfin team mailing list archive

petsc backend and periodic boundary conditions by PeriodicBC

 


hello!

i observed some strange behavior when i run the periodic poisson demo
(/pde/periodic). i use PETSC as backend and when it comes to applying
the bc's i get strange messages all over std output, here is a cut-out
from that:


[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: MatRestoreRow() line 242 in src/mat/interface/matrix.c
[0]PETSC ERROR: --------------------- Error Message
------------------------------------
[0]PETSC ERROR: Object is in wrong state!
[0]PETSC ERROR: Not for unassembled matrix!
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 2.3.3, Patch 15, Tue Sep 23
10:02:49 CDT 2008 HG revision: 31306062cd1a6f6a2496fccb4878f485c9b91760
[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: Unknown Name on a x64-ompi named jokull by priesen Wed
Jan 21 14:37:50 2009
[0]PETSC ERROR: Libraries linked from
/home/priesen/softw/petsc-2.3.3-p15/lib/x64-ompi
[0]PETSC ERROR: Configure run at Fri Nov 28 17:20:38 2008
[0]PETSC ERROR: Configure options --with-shared
--with-mpi-dir=/home/priesen/num/openmpi
--with-blas-lapack-dir=/usr/pack/atlas-3.7.24-ds --download-ml
--download-hypre --download-umfpack --download-mumps
--download-scalapack --download-blacs
[0]PETSC ERROR:

and and and...

this stuff occurs exactly between the "Applying periodic boundary
conditions" and the "Applying periodic boundary conditions (finished)"
and then the program finishes successfully, though i do not trust the
solution.
Do you ever had something similar to this?
Is there a problem that i can not use periodic boundary conditions
together with PETSC or did i possibly miss-configured either PETSC or
dolfin?


Thanks for your help,

Patrick





Follow ups