<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:SimSun;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:SimSun;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:"\@SimSun";
panose-1:2 1 6 0 3 1 1 1 1 1;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal">Hi, I’m trying to use PETSc ksp to solve a linear equation with the following option: “-ksp_type preonly -pc_type lu -pc_factor_mat_ordering_type amd”<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I got an error message below. Does anyone have an idea how it comes from?<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR:<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: Running KSP of preonly doesn't make sense with nonzero initial guess<o:p></o:p></p>
<p class="MsoNormal"> you probably want a KSP type of Richardson<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: Petsc Release Version 3.6.0, Jun, 09, 2015<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: Unknown Name on a olympus-openmpi_f-gnu-cxx-complex-opt named olympus.local by d3m956 Wed Jul 22 11:37:26 2015<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: Configure options PETSC_ARCH=olympus-openmpi_f-gnu-cxx-complex-opt --with-prefix=./ --with-mpi=1 --with-cc=mpicc --with-fc=mpif90 --with-cxx=mpicxx --with-c++-support=1 --with-c-support=0 --with-c-language=C++ --with-fortran=0
--with-scalar-type=complex --with-fortran-kernels=generic --download-superlu_dist --download-superlu --download-parmetis --download-metis --download-f2cblaslapack=1 --download-suitesparse --with-clanguage=c++ --with-shared-libraries=0 --with-x=0 --with-mpirun=mpirun
--with-mpiexec=mpiexec --with-debugging=0<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: #1 KSPSolve_PREONLY() line 23 in /pic/projects/gridpack/software_f/petsc-3.6.0/src/ksp/ksp/impls/preonly/preonly.c<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: #2 KSPSolve() line 604 in /pic/projects/gridpack/software_f/petsc-3.6.0/src/ksp/ksp/interface/itfunc.c<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: #3 p_solve() line 120 in /pic/projects/gridpack/gridpack-Github/src/math/petsc/petsc_linear_solver_implementation.hpp<o:p></o:p></p>
<p class="MsoNormal">terminate called after throwing an instance of 'gridpack::math::PETScException'<o:p></o:p></p>
<p class="MsoNormal"> what(): PETSc error (0): Error detected in C PETSc<o:p></o:p></p>
<p class="MsoNormal">[olympus:27408] *** Process received signal ***<o:p></o:p></p>
<p class="MsoNormal">[olympus:27408] Signal: Aborted (6)<o:p></o:p></p>
<p class="MsoNormal">[olympus:27408] Signal code: (-6)<o:p></o:p></p>
<p class="MsoNormal">[olympus:27408] [ 0] /lib64/libpthread.so.0 [0x3d2920eca0]<o:p></o:p></p>
<p class="MsoNormal">[olympus:27408] [ 1] /lib64/libc.so.6(gsignal+0x35) [0x3d2862ffc5]<o:p></o:p></p>
<p class="MsoNormal">[olympus:27408] [ 2] /lib64/libc.so.6(abort+0x110) [0x3d28631a70]<o:p></o:p></p>
<p class="MsoNormal">[olympus:27408] [ 3] /share/apps/gcc/4.6.2/lib64/libstdc++.so.6(_ZN9__gnu_cxx27__verbose_terminate_handlerEv+0x11d) [0x7f4aa53705bd]<o:p></o:p></p>
<p class="MsoNormal">[olympus:27408] [ 4] /share/apps/gcc/4.6.2/lib64/libstdc++.so.6 [0x7f4aa536e766]<o:p></o:p></p>
<p class="MsoNormal">[olympus:27408] [ 5] /share/apps/gcc/4.6.2/lib64/libstdc++.so.6 [0x7f4aa536e793]<o:p></o:p></p>
<p class="MsoNormal">[olympus:27408] [ 6] /share/apps/gcc/4.6.2/lib64/libstdc++.so.6 [0x7f4aa536e88e]<o:p></o:p></p>
<p class="MsoNormal">[olympus:27408] [ 7] ./ds2.x(_ZNK8gridpack4math31PETScLinearSolverImplementationISt7complexIdEiE7p_solveERKNS0_7VectorTIS3_iEERS6_+0x4cf) [0x74634f]<o:p></o:p></p>
<p class="MsoNormal">[olympus:27408] [ 8] ./ds2.x(_ZNK8gridpack4math25BaseLinearSolverInterfaceISt7complexIdEiE5solveERKNS0_7VectorTIS3_iEERS6_+0x37) [0x65e84f]<o:p></o:p></p>
<p class="MsoNormal">[olympus:27408] [ 9] ./ds2.x(_ZNK8gridpack4math13LinearSolverTISt7complexIdEiE7p_solveERKNS0_7VectorTIS3_iEERS6_+0x3b) [0x68fb89]<o:p></o:p></p>
<p class="MsoNormal">[olympus:27408] [10] ./ds2.x(_ZNK8gridpack4math25BaseLinearSolverInterfaceISt7complexIdEiE5solveERKNS0_7VectorTIS3_iEERS6_+0x37) [0x65e84f]<o:p></o:p></p>
<p class="MsoNormal">[olympus:27408] [11] ./ds2.x(_ZN8gridpack18dynamic_simulation6DS2App5solveENS0_9DS2Branch5EventE+0x8f0) [0x656386]<o:p></o:p></p>
<p class="MsoNormal">[olympus:27408] [12] ./ds2.x(main+0x41d) [0x64d9b1]<o:p></o:p></p>
<p class="MsoNormal">[olympus:27408] [13] /lib64/libc.so.6(__libc_start_main+0xf4) [0x3d2861d9f4]<o:p></o:p></p>
<p class="MsoNormal">[olympus:27408] [14] ./ds2.x [0x64d4dd]<o:p></o:p></p>
<p class="MsoNormal">[olympus:27408] *** End of error message ***<o:p></o:p></p>
<p class="MsoNormal">Abort<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thanks,<o:p></o:p></p>
<p class="MsoNormal">Shuangshuang<o:p></o:p></p>
</div>
</body>
</html>