ok, i'll do that, thank you
patrick
Garth N. Wells wrote:
>
>
> 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
>>>>> <[email protected]> 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
>>>>>> [email protected]
>>>>>> http://www.fenics.org/mailman/listinfo/dolfin-dev
>>>>>>
>>>> _______________________________________________
>>>> DOLFIN-dev mailing list
>>>> [email protected]
>>>> http://www.fenics.org/mailman/listinfo/dolfin-dev
>>>
>>
>>
>> _______________________________________________
>> DOLFIN-dev mailing list
>> [email protected]
>> http://www.fenics.org/mailman/listinfo/dolfin-dev
_______________________________________________
DOLFIN-dev mailing list
[email protected]
http://www.fenics.org/mailman/listinfo/dolfin-dev
