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 _______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
