<div dir="ltr"><div class="gmail_quote">On Fri, Apr 1, 2011 at 05:21, Ormiston, Scott J. <span dir="ltr"><<a href="mailto:SJ_Ormiston@umanitoba.ca">SJ_Ormiston@umanitoba.ca</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
I rebuilt PETSc with SuperLU and ParMETIS, and then I tried to run ex15f.F (it was working with a previous build of PETSc).<br></blockquote><div><br></div><div>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</div>
<div><br></div><div> -pc_factor_mat_solver_package superlu</div><div><br></div><div>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.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<br>
I ran<br>
<br>
% mpiexec -n 4 ex15f -pc_type lu -pc_factor_mat_solver_type superlu_dist<br>
<br>
but this generated the following error messages:<br>
<br>
====================================================================================<br>
[0]PETSC ERROR: --------------------- Error Message ------------------------------------<br>
[0]PETSC ERROR: No support for this operation for this object type!<br>
[0]PETSC ERROR: Matrix format mpiaij does not have a built-in PETSc direct solver!<br>
[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[0]PETSC ERROR: Petsc Release Version 3.1.0, Patch 6, Tue Nov 16 17:02:32 CST 2010<br>
[0]PETSC ERROR: See docs/changes/index.html for recent updates.<br>
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<br>
[0]PETSC ERROR: See docs/index.html for manual pages.<br>
[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[0]PETSC ERROR: ex15f on a linux-gnu named mecfd02 by sormist Thu Mar 31 21:15:34 2011<br>
[0]PETSC ERROR: Libraries linked from /home/mecfd/common/sw/petsc-3.1-p6/linux-gnu-c-debug/lib<br>
[0]PETSC ERROR: Configure run at Thu Mar 31 19:59:09 2011<br>
[0]PETSC ERROR: Configure options --with-mpi-dir=/home/mecfd/common/sw/openmpi-1.4.3 --download-superlu_dist --download-parmetis<br>
[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[0]PETSC ERROR: MatGetFactor() line 3644 in src/mat/interface/matrix.c<br>
[0]PETSC ERROR: PCSetUp_LU() line 133 in src/ksp/pc/impls/factor/lu/lu.c<br>
[0]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c<br>
[0]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c<br>
[0]PETSC ERROR: KSPSolve() line 353 in src/ksp/ksp/interface/itfunc.c<br>
[1]PETSC ERROR: --------------------- Error Message ------------------------------------<br>
[1]PETSC ERROR: No support for this operation for this object type!<br>
[1]PETSC ERROR: Matrix format mpiaij does not have a built-in PETSc direct solver!<br>
[1]PETSC ERROR: ------------------------------------------------------------------------<br>
[1]PETSC ERROR: Petsc Release Version 3.1.0, Patch 6, Tue Nov 16 17:02:32 CST 2010<br>
[1]PETSC ERROR: See docs/changes/index.html for recent updates.<br>
[1]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<br>
[1]PETSC ERROR: See docs/index.html for manual pages.<br>
[1]PETSC ERROR: ------------------------------------------------------------------------<br>
[1]PETSC ERROR: ex15f on a linux-gnu named mecfd02 by sormist Thu Mar 31 21:15:34 2011<br>
[1]PETSC ERROR: Libraries linked from /home/mecfd/common/sw/petsc-3.1-p6/linux-gnu-c-debug/lib<br>
[1]PETSC ERROR: Configure run at Thu Mar 31 19:59:09 2011<br>
[1]PETSC ERROR: Configure options --with-mpi-dir=/home/mecfd/common/sw/openmpi-1.4.3 --download-superlu_dist --download-parmetis<br>
[1]PETSC ERROR: ------------------------------------------------------------------------<br>
[1]PETSC ERROR: MatGetFactor() line 3644 in src/mat/interface/matrix.c<br>
[1]PETSC ERROR: PCSetUp_LU() line 133 in src/ksp/pc/impls/factor/lu/lu.c<br>
[1]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c<br>
[1]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c<br>
[1]PETSC ERROR: KSPSolve() line 353 in src/ksp/ksp/interface/itfunc.c<br>
[2]PETSC ERROR: --------------------- Error Message ------------------------------------<br>
[2]PETSC ERROR: No support for this operation for this object type!<br>
[2]PETSC ERROR: Matrix format mpiaij does not have a built-in PETSc direct solver!<br>
[2]PETSC ERROR: ------------------------------------------------------------------------<br>
[2]PETSC ERROR: Petsc Release Version 3.1.0, Patch 6, Tue Nov 16 17:02:32 CST 2010<br>
[2]PETSC ERROR: See docs/changes/index.html for recent updates.<br>
[2]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<br>
[2]PETSC ERROR: See docs/index.html for manual pages.<br>
[2]PETSC ERROR: ------------------------------------------------------------------------<br>
[2]PETSC ERROR: ex15f on a linux-gnu named mecfd02 by sormist Thu Mar 31 21:15:34 2011<br>
[2]PETSC ERROR: Libraries linked from /home/mecfd/common/sw/petsc-3.1-p6/linux-gnu-c-debug/lib<br>
[2]PETSC ERROR: Configure run at Thu Mar 31 19:59:09 2011<br>
[2]PETSC ERROR: Configure options --with-mpi-dir=/home/mecfd/common/sw/openmpi-1.4.3 --download-superlu_dist --download-parmetis<br>
[2]PETSC ERROR: ------------------------------------------------------------------------<br>
[2]PETSC ERROR: MatGetFactor() line 3644 in src/mat/interface/matrix.c<br>
[2]PETSC ERROR: PCSetUp_LU() line 133 in src/ksp/pc/impls/factor/lu/lu.c<br>
[2]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c<br>
[2]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c<br>
[2]PETSC ERROR: KSPSolve() line 353 in src/ksp/ksp/interface/itfunc.c<br>
[3]PETSC ERROR: --------------------- Error Message ------------------------------------<br>
[3]PETSC ERROR: No support for this operation for this object type!<br>
[3]PETSC ERROR: Matrix format mpiaij does not have a built-in PETSc direct solver!<br>
[3]PETSC ERROR: ------------------------------------------------------------------------<br>
[3]PETSC ERROR: Petsc Release Version 3.1.0, Patch 6, Tue Nov 16 17:02:32 CST 2010<br>
[3]PETSC ERROR: See docs/changes/index.html for recent updates.<br>
[3]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<br>
[3]PETSC ERROR: See docs/index.html for manual pages.<br>
[3]PETSC ERROR: ------------------------------------------------------------------------<br>
[3]PETSC ERROR: ex15f on a linux-gnu named mecfd02 by sormist Thu Mar 31 21:15:34 2011<br>
[3]PETSC ERROR: Libraries linked from /home/mecfd/common/sw/petsc-3.1-p6/linux-gnu-c-debug/lib<br>
[3]PETSC ERROR: Configure run at Thu Mar 31 19:59:09 2011<br>
[3]PETSC ERROR: Configure options --with-mpi-dir=/home/mecfd/common/sw/openmpi-1.4.3 --download-superlu_dist --download-parmetis<br>
[3]PETSC ERROR: ------------------------------------------------------------------------<br>
[3]PETSC ERROR: MatGetFactor() line 3644 in src/mat/interface/matrix.c<br>
[3]PETSC ERROR: PCSetUp_LU() line 133 in src/ksp/pc/impls/factor/lu/lu.c<br>
[3]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c<br>
[3]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c<br>
[3]PETSC ERROR: KSPSolve() line 353 in src/ksp/ksp/interface/itfunc.c<br>
Norm of error 2.0050E+02 iterations 0<br>
WARNING! There are options you set that were not used!<br>
WARNING! could be spelling mistake, etc!<br>
Option left: name:-pc_factor_mat_solver_type value: superlu_dist<br>
====================================================================================<br>
<br>
Does this require changing the code in ex15f.F or is there something else that I am doing wrong?<br><font color="#888888">
<br>
Scott Ormiston<br>
</font></blockquote></div><br></div>