<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<div class="moz-cite-prefix">Hi,<br>
<br>
Part of my code is :<br>
<br>
<b><i>do ...<br>
<br>
call
MatSetValues(A_mat,1,II,1,JJ,big_A(ijk,kk),ADD_VALUES,ierr)</i></b><b><i><br>
</i></b><b><i> </i></b><b><i><br>
</i></b><b><i>end do</i></b><b><i> </i></b><b><i><br>
</i></b><b><i><br>
</i></b><b><i> call
MatAssemblyBegin(A_mat,MAT_FINAL_ASSEMBLY,ierr)</i></b><b><i><br>
</i></b><b><i><br>
</i></b><b><i> call
MatAssemblyEnd(A_mat,MAT_FINAL_ASSEMBLY,ierr)</i></b><b><i><br>
</i></b><b><i> </i></b><b><i><br>
</i></b><b><i> call
MatSetOption(A_mat,MAT_STRUCTURALLY_SYMMETRIC,PETSC_TRUE,ierr)</i></b><b><i><br>
</i></b><b><i><br>
</i></b><b><i> call
MatSetOption(A_mat,MAT_NEW_NONZERO_LOCATIONS,PETSC_FALSE,ierr)</i></b><b><i><br>
</i></b><b><i> </i></b><b><i><br>
</i></b><b><i> call
KSPSetOperators(ksp,A_mat,A_mat,SAME_NONZERO_PATTERN,ierr)</i></b><b><i><br>
</i></b><b><i><br>
</i></b><b><i> call KSPGetPC(ksp,pc,ierr) </i></b><b><i><br>
</i></b><b><i> </i></b><b><i><br>
</i></b><b><i> call KSPSetOptionsPrefix(ksp,"poisson_",ierr)</i></b><b><i><br>
</i></b><b><i><br>
</i></b><b><i> ksptype=KSPBCGS</i></b><b><i><br>
</i></b><b><i><br>
</i></b><b><i> call KSPSetType(ksp,ksptype,ierr)</i></b><b><i><br>
</i></b><b><i><br>
</i></b><b><i> !call PCSetType(pc,PCGAMG,ierr)</i></b><b><i><br>
</i></b><b><i><br>
</i></b><b><i><br>
</i></b><b><i> call KSPSetFromOptions(ksp,ierr)</i></b><b><i><br>
</i></b><b><i><br>
</i></b><b><i> tol=1.e-5</i></b><b><i><br>
</i></b><b><i><br>
</i></b><b><i> call
KSPSetTolerances(ksp,tol,PETSC_DEFAULT_DOUBLE_PRECISION,PETSC_DEFAULT_DOUBLE_PRECISION,PETSC_DEFAULT_INTEGER,ierr)</i></b><b><i><br>
</i></b><b><i><br>
do ...<br>
<br>
</i></b><b><i>call
VecSetValue(b_rhs,II,q_p(ijk),INSERT_VALUES,ierr)</i></b><b><i><br>
<br>
end do<br>
</i></b><b><i><br>
</i></b><i><b>c</b></i><b><i>all VecAssemblyBegin(b_rhs,ierr)</i></b><b><i><br>
</i></b><b><i><br>
</i></b><b><i>call VecAssemblyEnd(b_rhs,ierr)</i></b><b><i><br>
</i></b><b><i><br>
</i></b><b><i>call VecAssemblyBegin(xx,ierr)</i></b><b><i><br>
</i></b><b><i><br>
</i></b><b><i>call VecAssemblyEnd(xx,ierr)</i></b><b><i><br>
</i></b><b><i><br>
</i></b><b><i>call KSPSolve(ksp,b_rhs,xx,ierr)</i></b><b><i> -
hang<br>
</i></b><b><i><br>
</i></b><b><i>call KSPGetConvergedReason(ksp,reason,ierr)</i></b><br>
<br>
It hangs at <b><i>KSPSolve</i></b>, when used with the option :<br>
<br>
-poisson_pc_gamg_agg_nsmooths 1 -poisson_pc_type gamg<br>
<pre class="moz-signature" cols="72">Yours sincerely,
TAY wee-beng</pre>
On 20/7/2013 8:36 PM, Jed Brown wrote:<br>
</div>
<blockquote cite="mid:87fvv984f8.fsf@mcs.anl.gov" type="cite">
<pre wrap="">-poisson_pc_gamg_agg_nsmooths 1 -poisson_pc_type gamg</pre>
</blockquote>
<br>
</body>
</html>