<div dir="ltr"><div>Carol,</div><div><br></div><div>I tested qcg with petsc/src/ksp/ksp/examples/tutorials/ex2.c and found it crashes with np=2, pc=bjacobi.<br></div><div><br></div><div>Then I found in <a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/KSP/KSPQCG.html">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/KSP/KSPQCG.html</a></div><div><br></div><div><span style="color:rgb(0,0,0);font-family:Times;font-size:medium">Notes: This is rarely used directly</span><br></div><div><h3 style="color:rgb(0,0,0);font-family:Times"><font color="#CC3333">Currently we allow symmetric preconditioning with the following scaling matrices</font></h3><a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/PC/PCNONE.html#PCNONE" style="font-family:Times;font-size:medium">PCNONE</a><span style="color:rgb(0,0,0);font-family:Times;font-size:medium">: D = Identity matrix </span><a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/PC/PCJACOBI.html#PCJACOBI" style="font-family:Times;font-size:medium">PCJACOBI</a><span style="color:rgb(0,0,0);font-family:Times;font-size:medium">: D = diag [d_1, d_2, ...., d_n], where d_i = sqrt(H[i,i]) </span><a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/PC/PCICC.html#PCICC" style="font-family:Times;font-size:medium">PCICC</a><span style="color:rgb(0,0,0);font-family:Times;font-size:medium">: D = L^T, implemented with forward and backward solves. Here L is an incomplete Cholesky factor of H.</span></div><div><br></div><div>Hong</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jun 29, 2015 at 9:22 AM,  <span dir="ltr"><<a href="mailto:Carol.Brickley@awe.co.uk" target="_blank">Carol.Brickley@awe.co.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">All,<br>
<br>
I am trying to use KSPSolve for a QCG solver using petsc 3.4.3 and it gives me the error of "PC does not have left symmetric apply!". Further reading on the web pages suggest I need to use KSPSetPCSide to set up a symmetric side for the solver. However, this does not seem to work. I am doing the following in an F90 code:<br>
<br>
call KSPSetPCType(ksp_ptr%ksp, 'qcg', petsc_err)<br>
call KSPSetPCSide(ksp_ptr%ksp, PC_SYMMETRIC, petsc_err)<br>
<br>
but I still get the same failure message when I call KSPSolve(ksp_ptr%ksp..)  later on.<br>
<br>
Any ideas as to why this is not working?<br>
<br>
Carol<br>
<br>
<br>
<br>
<br>
<br>
Dr Carol Brickley<br>
BSc,PhD,ARCS,DIC,MBCS<br>
<br>
Senior Software Engineer<br>
Applied Computer Science<br>
DS+T,<br>
AWE<br>
Aldermaston<br>
Reading<br>
Berkshire<br>
RG7 4PR<br>
<br>
Direct: 0118 9855035<br>
<br>
<br>
<br>
___________________________________________________<br>
____________________________<br>
<br>
The information in this email and in any attachment(s) is<br>
commercial in confidence. If you are not the named addressee(s)<br>
or<br>
if you receive this email in error then any distribution, copying or<br>
use of this communication or the information in it is strictly<br>
prohibited. Please notify us immediately by email at<br>
admin.internet(at)<a href="http://awe.co.uk" rel="noreferrer" target="_blank">awe.co.uk</a>, and then delete this message from<br>
your computer. While attachments are virus checked, AWE plc<br>
does not accept any liability in respect of any virus which is not<br>
detected.<br>
<br>
AWE Plc<br>
Registered in England and Wales<br>
Registration No 02763902<br>
AWE, Aldermaston, Reading, RG7 4PR<br>
</blockquote></div><br></div>