Re: [Libmesh-users] Dense Petsc Matrix

[Paul T. B. <pb...@ic...>]

If you have MUMPS installed, you can do:

-ksp_type preonly -pc_type lu -mat_type aijmumps

or if you have SuperLu_Dist

-mat_type superlu_dist

both *should* work in parallel with PETSc (unless you compiled MUMPS 
with serial only). MUMPS is not as error catching friendly as SuperLU 
for zero pivots, etc, but both are reasonably fast sparse-direct solvers.

The key is having a call to the PETSc routine SetMatFromOptions which 
will set the MatType from the command line. If you don't have anything, 
it just uses what you've set already, or the defaults. This is very 
useful for debugging and quickly switching between direct and iterative 
solvers. Not sure if libMesh has this call or not, but you can probably 
work it in yourself without too much trouble.

FWIW.

Paul

Andrea Hawkins wrote:
> Ah. Ok. Well, I guess I"ll just run it in serial for now. And if I need it
> in parallel, I'll add the option in PetscMatrix and see what happens! =)
>
> Thanks!
> Andrea
>
> On Fri, Mar 13, 2009 at 1:16 PM, John Peterson <jwp...@gm...> wrote:
>
>   
>> On Fri, Mar 13, 2009 at 1:15 PM, Andrea Hawkins <and...@gm...>
>> wrote:
>>     
>>> Well, the suggestion PETSc gives for performing a direct solve is to use
>>>       
>> the
>>     
>>> options
>>>
>>>  -ksp_type preonly -pc_type lu
>>>
>>> But, when I use it I get the error:
>>>
>>> [0]PETSC ERROR: [1]PETSC ERROR: --------------------- Error Message
>>> ------------------------------------
>>> [2]PETSC ERROR: --------------------- Error Message
>>> ------------------------------------
>>> [1]PETSC ERROR: No support for this operation for this object type!
>>> [2]PETSC ERROR: No support for this operation for this object type!
>>> [1]PETSC ERROR: [2]PETSC ERROR: Matrix type mpiaij  symbolic LU!
>>> [1]PETSC ERROR: Matrix type mpiaij  symbolic LU!
>>> [2]PETSC ERROR:
>>>       
>> Sorry, I forgot to mention: it only works in serial.
>>
>> --
>> John
>>
>>     
> ------------------------------------------------------------------------------
> Apps built with the Adobe(R) Flex(R) framework and Flex Builder(TM) are
> powering Web 2.0 with engaging, cross-platform capabilities. Quickly and
> easily build your RIAs with Flex Builder, the Eclipse(TM)based development
> software that enables intelligent coding and step-through debugging.
> Download the free 60 day trial. http://p.sf.net/sfu/www-adobe-com
> _______________________________________________
> Libmesh-users mailing list
> Lib...@li...
> https://lists.sourceforge.net/lists/listinfo/libmesh-users
>
>