From 5019b4b212ec23d6c70ce59e5389525c741d7629 Mon Sep 17 00:00:00 2001 From: raphael couturier Date: Sun, 17 Aug 2014 17:05:22 +0200 Subject: [PATCH] ajout ex29 --- code/ex29.c | 9 +++++---- 1 file changed, 5 insertions(+), 4 deletions(-) diff --git a/code/ex29.c b/code/ex29.c index 910237b..f647a9b 100644 --- a/code/ex29.c +++ b/code/ex29.c @@ -1,6 +1,7 @@ -// /home/couturie/work/petsc-3.5.1/arch-linux2-c-debug/bin/mpirun -np 3 ex29 -da_grid_x 600 -da_grid_y 600 -ksp_type fgmres +// /home/couturie/work/petsc-3.5.1/arch-linux2-c-debug/bin/mpirun -np 3 ex29 -da_grid_x 600 -da_grid_y 600 -ksp_type fgmres -bc_type neumann + @@ -64,7 +65,7 @@ int KrylovMinimize(Mat A, Vec b, Vec x) { //Variables PetscScalar gamma, alpha, oldgamma, beta; - PetscReal norm=20, Eprecision=1e-7, 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; @@ -241,7 +242,7 @@ int KrylovMinimizeLSQR(Mat A, Vec b, Vec x) { //Variables PetscScalar alpha, beta; - PetscReal norm=20, Eprecision=1e-7, 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; @@ -501,7 +502,7 @@ int main(int argc,char **argv) ierr = KSPSetFromOptions(ksp);CHKERRQ(ierr); ierr = KSPSetUp(ksp);CHKERRQ(ierr); - ierr = KSPSetTolerances(ksp, 1e-7, 1e-7, PETSC_DEFAULT, 50000000); CHKERRQ(ierr); + ierr = KSPSetTolerances(ksp, 1e-6, 1e-6, PETSC_DEFAULT, 50000000); CHKERRQ(ierr); T1 = MPI_Wtime(); ierr = KSPSolve(ksp,NULL,NULL);CHKERRQ(ierr); T2 = MPI_Wtime(); -- 2.39.5