dolfin team mailing list archive

["Garth N. Wells" <gnw20@xxxxxxxxx>]

Patrick Riesen wrote:
> hello Garth
>
> thanks, that works fine with the demo and but switches to a krylov
> iterative solver because i only have umfpack via petsc. the krylov does
> not converge on the dolfin-grade2, a LU is favored but petsc obviously
> inhibits me to use periodic bc, what about another way to access another
> LU solver? 

Install UMFPACK if you can. It's available as a package for most systems.

Garth

point dolfin to the umfpack (in the externalpackages folder
> of my petsc) directly or use trilinois?
>
> patrick
>
> Garth N. Wells wrote:
> > Patrick Riesen wrote:
> > > Hello Andy
> > >
> > > Thanks for your reply! I am using dolfin 0.8.1 together with the 
> > > dolfin-grade2 library. I did not yet update to keep consistent with 
> > > the grade2lib. otherwise there might be a possibility to switch of the 
> > > petsc backend in the periodic demo or do i have to recompile dolfin to 
> > > not using petsc?
> > Try
> >
> >     dolfin_set("linear algebra backend", "uBLAS");
> >
> > Garth
> >
> > > patrick
> > >
> > >
> > > Andy Ray Terrel wrote:
> > > > 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
> > > _______________________________________________
> > > DOLFIN-dev mailing list
> > > DOLFIN-dev@xxxxxxxxxx
> > > http://www.fenics.org/mailman/listinfo/dolfin-dev
>
>
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev