[petsc-users] PETSC with SuperLU

Fri Apr 1 01:44:42 CDT 2011

On Fri, Apr 1, 2011 at 05:21, Ormiston, Scott J.
<SJ_Ormiston at umanitoba.ca>wrote:

> I rebuilt PETSc with SuperLU and ParMETIS, and then I tried to run ex15f.F
> (it was working with a previous build of PETSc).
>

SuperLU is not the same package as SuperLU_Dist (don't ask me why they
organize software that way). You also misspelled the command line option. If
you run the example below in serial with

  -pc_factor_mat_solver_package superlu

it should work (using SuperLU instead of PETSc's native direct solver). To
use SuperLU_Dist (which works in parallel), configure with
--download-superlu_dist and run with -pc_factor_mat_solver_package
superlu_dist. You do not have to touch the source code.

>
> I ran
>
> % mpiexec -n 4 ex15f -pc_type lu -pc_factor_mat_solver_type superlu_dist
>
> but this generated the following error messages:
>
>
> ====================================================================================
> [0]PETSC ERROR: --------------------- Error Message
> ------------------------------------
> [0]PETSC ERROR: No support for this operation for this object type!
> [0]PETSC ERROR: Matrix format mpiaij does not have a built-in PETSc direct
> solver!
> [0]PETSC ERROR:
> ------------------------------------------------------------------------
> [0]PETSC ERROR: Petsc Release Version 3.1.0, Patch 6, Tue Nov 16 17:02:32
> CST 2010
> [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: ex15f on a linux-gnu named mecfd02 by sormist Thu Mar 31
> 21:15:34 2011
> [0]PETSC ERROR: Libraries linked from
> /home/mecfd/common/sw/petsc-3.1-p6/linux-gnu-c-debug/lib
> [0]PETSC ERROR: Configure run at Thu Mar 31 19:59:09 2011
> [0]PETSC ERROR: Configure options
> --with-mpi-dir=/home/mecfd/common/sw/openmpi-1.4.3 --download-superlu_dist
> --download-parmetis
> [0]PETSC ERROR:
> ------------------------------------------------------------------------
> [0]PETSC ERROR: MatGetFactor() line 3644 in src/mat/interface/matrix.c
> [0]PETSC ERROR: PCSetUp_LU() line 133 in src/ksp/pc/impls/factor/lu/lu.c
> [0]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c
> [0]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c
> [0]PETSC ERROR: KSPSolve() line 353 in src/ksp/ksp/interface/itfunc.c
> [1]PETSC ERROR: --------------------- Error Message
> ------------------------------------
> [1]PETSC ERROR: No support for this operation for this object type!
> [1]PETSC ERROR: Matrix format mpiaij does not have a built-in PETSc direct
> solver!
> [1]PETSC ERROR:
> ------------------------------------------------------------------------
> [1]PETSC ERROR: Petsc Release Version 3.1.0, Patch 6, Tue Nov 16 17:02:32
> CST 2010
> [1]PETSC ERROR: See docs/changes/index.html for recent updates.
> [1]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
> [1]PETSC ERROR: See docs/index.html for manual pages.
> [1]PETSC ERROR:
> ------------------------------------------------------------------------
> [1]PETSC ERROR: ex15f on a linux-gnu named mecfd02 by sormist Thu Mar 31
> 21:15:34 2011
> [1]PETSC ERROR: Libraries linked from
> /home/mecfd/common/sw/petsc-3.1-p6/linux-gnu-c-debug/lib
> [1]PETSC ERROR: Configure run at Thu Mar 31 19:59:09 2011
> [1]PETSC ERROR: Configure options
> --with-mpi-dir=/home/mecfd/common/sw/openmpi-1.4.3 --download-superlu_dist
> --download-parmetis
> [1]PETSC ERROR:
> ------------------------------------------------------------------------
> [1]PETSC ERROR: MatGetFactor() line 3644 in src/mat/interface/matrix.c
> [1]PETSC ERROR: PCSetUp_LU() line 133 in src/ksp/pc/impls/factor/lu/lu.c
> [1]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c
> [1]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c
> [1]PETSC ERROR: KSPSolve() line 353 in src/ksp/ksp/interface/itfunc.c
> [2]PETSC ERROR: --------------------- Error Message
> ------------------------------------
> [2]PETSC ERROR: No support for this operation for this object type!
> [2]PETSC ERROR: Matrix format mpiaij does not have a built-in PETSc direct
> solver!
> [2]PETSC ERROR:
> ------------------------------------------------------------------------
> [2]PETSC ERROR: Petsc Release Version 3.1.0, Patch 6, Tue Nov 16 17:02:32
> CST 2010
> [2]PETSC ERROR: See docs/changes/index.html for recent updates.
> [2]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
> [2]PETSC ERROR: See docs/index.html for manual pages.
> [2]PETSC ERROR:
> ------------------------------------------------------------------------
> [2]PETSC ERROR: ex15f on a linux-gnu named mecfd02 by sormist Thu Mar 31
> 21:15:34 2011
> [2]PETSC ERROR: Libraries linked from
> /home/mecfd/common/sw/petsc-3.1-p6/linux-gnu-c-debug/lib
> [2]PETSC ERROR: Configure run at Thu Mar 31 19:59:09 2011
> [2]PETSC ERROR: Configure options
> --with-mpi-dir=/home/mecfd/common/sw/openmpi-1.4.3 --download-superlu_dist
> --download-parmetis
> [2]PETSC ERROR:
> ------------------------------------------------------------------------
> [2]PETSC ERROR: MatGetFactor() line 3644 in src/mat/interface/matrix.c
> [2]PETSC ERROR: PCSetUp_LU() line 133 in src/ksp/pc/impls/factor/lu/lu.c
> [2]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c
> [2]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c
> [2]PETSC ERROR: KSPSolve() line 353 in src/ksp/ksp/interface/itfunc.c
> [3]PETSC ERROR: --------------------- Error Message
> ------------------------------------
> [3]PETSC ERROR: No support for this operation for this object type!
> [3]PETSC ERROR: Matrix format mpiaij does not have a built-in PETSc direct
> solver!
> [3]PETSC ERROR:
> ------------------------------------------------------------------------
> [3]PETSC ERROR: Petsc Release Version 3.1.0, Patch 6, Tue Nov 16 17:02:32
> CST 2010
> [3]PETSC ERROR: See docs/changes/index.html for recent updates.
> [3]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
> [3]PETSC ERROR: See docs/index.html for manual pages.
> [3]PETSC ERROR:
> ------------------------------------------------------------------------
> [3]PETSC ERROR: ex15f on a linux-gnu named mecfd02 by sormist Thu Mar 31
> 21:15:34 2011
> [3]PETSC ERROR: Libraries linked from
> /home/mecfd/common/sw/petsc-3.1-p6/linux-gnu-c-debug/lib
> [3]PETSC ERROR: Configure run at Thu Mar 31 19:59:09 2011
> [3]PETSC ERROR: Configure options
> --with-mpi-dir=/home/mecfd/common/sw/openmpi-1.4.3 --download-superlu_dist
> --download-parmetis
> [3]PETSC ERROR:
> ------------------------------------------------------------------------
> [3]PETSC ERROR: MatGetFactor() line 3644 in src/mat/interface/matrix.c
> [3]PETSC ERROR: PCSetUp_LU() line 133 in src/ksp/pc/impls/factor/lu/lu.c
> [3]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c
> [3]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c
> [3]PETSC ERROR: KSPSolve() line 353 in src/ksp/ksp/interface/itfunc.c
> Norm of error 2.0050E+02 iterations     0
> WARNING! There are options you set that were not used!
> WARNING! could be spelling mistake, etc!
> Option left: name:-pc_factor_mat_solver_type value: superlu_dist
>
> ====================================================================================
>
> Does this require changing the code in ex15f.F or is there something else
> that I am doing wrong?
>
> Scott Ormiston
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110401/334b2aee/attachment-0001.htm>