回复: about combining the GMRES and smooth methods
Barry Smith
bsmith at mcs.anl.gov
Sat Jan 27 22:43:28 CST 2007
If you want help with why a compile fails you need to at least
send the entire output when the compile fails.
Barry
On Sun, 28 Jan 2007, jiaxun hou wrote:
> But I already had another preconditioner. I want to study the behave of the smooth method in the GMRES with my preconditioner. How can I do?
Thanks!
-------------------------------------------------------------
----- ÔʼÓʼþ ----
·¢¼þÈË£º Barry Smith <bsmith at mcs.anl.gov>
ÊÕ¼þÈË£º petsc-users at mcs.anl.gov
ÒÑ·¢ËÍ£º 2007/1/28(ÖÜÈÕ), ÉÏÎç11:28:35
Ö÷Ì⣺ Re: about combining the GMRES and smooth methods
It seems extremely likely that you do not want to do this.
If you want to use SOR or SSOR as a preconditioner for GMRES simply
use -pc_type sor; see the manual page for PCSOR for other options.
Barry
On Sun, 28 Jan 2007, jiaxun hou wrote:
> Hi,
I want to add a smooth method in the GMRES routine. I do it in the function KSPSolve_GMRES of /src/ksp/ksp/impls/gmres/gmres.c. The code likes that:
PetscErrorCode KSPSolve_GMRES(KSP ksp) {
.....
ierr = GMREScycle(&its,ksp);CHKERRQ(ierr);
ierr = MatRelax(ksp->pc->mat, ksp->vec_rhs,1,SOR_FORWARD_SWEEP,0,1,1,ksp->vec_sol);CHKERRQ(ierr);
.....
}
But compile fail. Is it the right way?
Thanks.
___________________________________________________________
ÑÅ»¢Ãâ·ÑÓÊÏä-3.5GÈÝÁ¿£¬20M¸½¼þ
http://cn.mail.yahoo.com/
___________________________________________________________
ÇÀ×¢ÑÅ»¢Ãâ·ÑÓÊÏä-3.5GÈÝÁ¿£¬20M¸½¼þ£¡
http://cn.mail.yahoo.com
More information about the petsc-users
mailing list