From 8a6af80b29277f6ff22e64ff6615d3ea0edf6434 Mon Sep 17 00:00:00 2001
From: Pietro Incardona <incardon@mpi-cbg.de>
Date: Wed, 12 Jul 2017 19:02:10 +0200
Subject: [PATCH] Fixing update site for the dev version
---
example/Numerics/Vortex_in_cell/main_vic_petsc.cpp | 14 ++++++++++++--
update | 2 +-
2 files changed, 13 insertions(+), 3 deletions(-)
diff --git a/example/Numerics/Vortex_in_cell/main_vic_petsc.cpp b/example/Numerics/Vortex_in_cell/main_vic_petsc.cpp
index f0360eeb6..7905066ae 100644
--- a/example/Numerics/Vortex_in_cell/main_vic_petsc.cpp
+++ b/example/Numerics/Vortex_in_cell/main_vic_petsc.cpp
@@ -477,18 +477,26 @@ void helmotz_hodge_projection(grid_type & gr, const Box<3,float> & domain, petsc
//! \cond [solve_petsc] \endcond
+ timer tm_solve;
if (init == true)
{
// Set the conjugate-gradient as method to solve the linear solver
solver.setSolver(KSPBCGS);
// Set the absolute tolerance to determine that the method is converged
- solver.setAbsTol(0.1);
+ solver.setAbsTol(0.001);
// Set the maximum number of iterations
solver.setMaxIter(500);
- timer tm_solve;
+ solver.setPreconditioner(PCHYPRE_BOOMERAMG);
+ solver.setPreconditionerAMG_nl(6);
+ solver.setPreconditionerAMG_maxit(1);
+ solver.setPreconditionerAMG_relax("SOR/Jacobi");
+ solver.setPreconditionerAMG_cycleType("V",0,4);
+ solver.setPreconditionerAMG_coarsenNodalType(0);
+ solver.setPreconditionerAMG_coarsen("HMIS");
+
tm_solve.start();
// Give to the solver A and b, return x, the solution
@@ -498,7 +506,9 @@ void helmotz_hodge_projection(grid_type & gr, const Box<3,float> & domain, petsc
}
else
{
+ tm_solve.start();
solver.solve(x_,fd.getB());
+ tm_solve.stop();
}
// copy the solution x to the grid psi
diff --git a/update b/update
index c2f6a8c22..71880b9be 100644
--- a/update
+++ b/update
@@ -1 +1 @@
-http://openfpm.mpi-cbg.de/upload/update_0.8
+http://openfpm.mpi-cbg.de/upload/update_dev
--
GitLab