dolfin team mailing list archive

[Andy Ray Terrel <aterrel@xxxxxxxxxxxx>]

Hi Patrick,

Are you using the development version of FEniCS or still 0.8.1 with
dolfin-grade2? I'm trying to update my library at the moment and the
periodic demo works without petsc on the dev version, so you might be
able to work on it with that.  It looks like somewhere in the code is
using MatRestoreRow with a modified array instead of MatSetValues, but
not having the dev environment setup with Petsc I haven't found the
problem yet.

-- Andy

On Wed, Jan 21, 2009 at 7:56 AM, Patrick Riesen <priesen@xxxxxxxxxxxxxxx> wrote:
>
> 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
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev
>