X-Git-Url: https://bilbo.iut-bm.univ-fcomte.fr/and/gitweb/GMRES2stage.git/blobdiff_plain/704a29e1a3accbeb491a424bc8e02cb8a54dd3b2..5019b4b212ec23d6c70ce59e5389525c741d7629:/code/ex15.c diff --git a/code/ex15.c b/code/ex15.c index f62763f..1fffcbe 100644 --- a/code/ex15.c +++ b/code/ex15.c @@ -56,7 +56,7 @@ int KrylovMinimize(Mat A, Vec b, Vec x) { //Variables PetscScalar gamma, alpha, oldgamma, beta; - PetscReal norm=20, Eprecision=1e-8, cgprec=1e-40; + PetscReal norm=20, Eprecision=1e-6, cgprec=1e-40; PetscInt giter=0, ColS=12, col=0, Emaxiter=50000000, iter=0, iterations=15, Iiter=0; PetscErrorCode ierr; PetscScalar T1, T2; @@ -233,7 +233,7 @@ int KrylovMinimizeLSQR(Mat A, Vec b, Vec x) { //Variables PetscScalar alpha, beta; - PetscReal norm=20, Eprecision=1e-8, tol=1e-40; + PetscReal norm=20, Eprecision=1e-6, tol=1e-40; PetscInt giter=0, ColS=12, col=0, Emaxiter=50000000, iter=0, iterations=20, Iiter=0; PetscErrorCode ierr; PetscScalar T1, T2; @@ -310,7 +310,7 @@ int KrylovMinimizeLSQR(Mat A, Vec b, Vec x) { //Initializations // ierr = KSPGMRESSetRestart(ksp, 16); CHKERRQ(ierr); - ierr = KSPSetTolerances(ksp, 1e-10, 1e-10, PETSC_DEFAULT, 16); CHKERRQ(ierr); + ierr = KSPSetTolerances(ksp, 1e-13, 1e-13, PETSC_DEFAULT, 16); CHKERRQ(ierr); ierr = KSPSetInitialGuessNonzero(ksp, PETSC_TRUE); CHKERRQ(ierr); @@ -572,7 +572,7 @@ int main(int argc,char **args) to set various options. */ ierr = KSPGetPC(ksp,&pc);CHKERRQ(ierr); - ierr = KSPSetTolerances(ksp,1e-9,1e-9,PETSC_DEFAULT,5000000);CHKERRQ(ierr); + ierr = KSPSetTolerances(ksp,1e-7,1e-7,PETSC_DEFAULT,5000000);CHKERRQ(ierr); /* Set runtime options, e.g.,