[petsc-users] run direct linear solver in parallel

[Hong Zhang]

Note: spooles has been out of support from its developers for more
than 10 years.
Recommend to use superlu_dist or mumps.

Hong

On Fri, Dec 3, 2010 at 4:22 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:
>
>  You are using an ANCIENT version of PETSc but using documentation from a much newer version. You should upgrade to petsc-3.1 immediately then life will be easy :-)
>
>
>   Barry
>
> On Dec 3, 2010, at 4:19 PM, Xiangdong Liang wrote:
>
>> Hi everyone,
>>
>> I am wondering how I can run the direct solver in parallel. I can run
>> my program in a single processor with direct linear solver by
>>
>> ./foo.out  -ksp_type preonly -pc_type lu -pc_factor_mat_solver_package spooles
>>
>> However, when I try to run it with mpi:
>>
>> mpirun.openmpi -np 2 ./foo.out -ksp_type preonly -pc_type lu
>> -pc_factor_mat_solver_package spooles
>>
>> I got error like this:
>>
>> [0]PETSC ERROR: --------------------- Error Message
>> ------------------------------------
>> [0]PETSC ERROR: No support for this operation for this object type!
>> [0]PETSC ERROR: Matrix type mpiaij  symbolic LU!
>>
>> [0]PETSC ERROR: Libraries linked from
>> /home/hazelsct/petsc-2.3.3/lib/linux-gnu-c-opt
>> [0]PETSC ERROR: Configure run at Mon Jun 30 14:37:52 2008
>> [0]PETSC ERROR: Configure options --with-shared --with-dynamic
>> --with-debugging=0 --useThreads 0 --with-mpi-dir=/usr/lib/openmpi
>> --with-mpi-shared=1 --with-blas-lib=-lblas --with-lapack-lib=-llapack
>> --with-umfpack=1 --with-umfpack-include=/usr/include/suitesparse
>> --with-umfpack-lib="[/usr/lib/libumfpack.so,/usr/lib/libamd.so]"
>> --with-superlu=1 --with-superlu-include=/usr/include/superlu
>> --with-superlu-lib=/usr/lib/libsuperlu.so --with-spooles=1
>> --with-spooles-include=/usr/include/spooles
>> --with-spooles-lib=/usr/lib/libspooles.so --with-hypre=1
>> --with-hypre-dir=/usr --with-babel=1 --with-babel-dir=/usr
>> [0]PETSC ERROR:
>> ------------------------------------------------------------------------
>> [0]PETSC ERROR: MatLUFactorSymbolic() line 2174 in src/mat/interface/matrix.c
>> [0]PETSC ERROR: PCSetUp_LU() line 257 in src/ksp/pc/impls/factor/lu/lu.c
>> -------------------------------------------------------
>>
>> Would you like to tell me where I am doing wrong? I appreciate your help.
>>
>> Xiangdong
>
>