[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
About unpreconditioned residuals in Left Preconditioned GMRES
- To: petsc-users@xxxxxxxxxxx
- Subject: About unpreconditioned residuals in Left Preconditioned GMRES
- From: jiaxun hou <jiaxun_hou@xxxxxxxxxxxx>
- Date: Sat, 2 Feb 2008 20:49:03 +0800 (CST)
- Domainkey-signature: a=rsa-sha1; q=dns; c=nofws; s=s1024; d=yahoo.com.cn; h=X-YMail-OSG:Received:Date:From:Subject:To:MIME-Version:Content-Type:Content-Transfer-Encoding:Message-ID; b=odiIUu09PClEl/EIoC3quVIFMGKpRbtR1tDFqOQ7pT/HZC9r0B+DgJAFms1tn1GMbIezGoNzhTmjTgkuyPYoIj5qs52iD0SgJYy0KTz49m4l/BA/Dl+3wtuIx0l5bSHBGSAiFbDOjQIHQj/eTZIPgNcnsmOXQejKYHthfmyrB60=;
- Reply-to: petsc-users@xxxxxxxxxxx
- Sender: owner-petsc-users@xxxxxxxxxxx
Hi everyone,
I want to use the Left Preconditioned GMRES to solve a linear system, and the stopping criterion must be based on the actual residuals (b-Ax). But the GMRES codes of PETSc seems to use the preconditioned residuals (B^-1(b-Ax)) only. In addition, when I set KSPSetNormType(ksp,KSP_NORM_UNPRECONDITIONED), I receive the error message: "Currently can use GMRES with only preconditioned residual (right preconditioning not coded)". So, is there any way to set stopping criterion based on the actual residuals?
Best regards,
Jiaxun
雅虎邮箱传递新年祝福,个性贺卡送亲朋!