diff --git a/dune/upscaling/UpscalerBase.hpp b/dune/upscaling/UpscalerBase.hpp
index 64f633f..def0416 100644
--- a/dune/upscaling/UpscalerBase.hpp
+++ b/dune/upscaling/UpscalerBase.hpp
@@ -80,9 +80,12 @@ namespace Dune
                   double perm_threshold,
                   double z_tolerance = 0.0,
                   double residual_tolerance = 1e-8,
+                  int linsolver_maxit = 0,
+                  double linsolver_prolongate_factor =1.6,
                   int linsolver_verbosity = 0,
                   int linsolver_type = 1,
-                  bool twodim_hack = false);
+                  bool twodim_hack = false,
+                  int linsolver_smooth_steps=2);
 
 	/// Access the grid.
 	const GridType& grid() const;
@@ -130,8 +133,11 @@ namespace Dune
 	BoundaryConditionType bctype_;
 	bool twodim_hack_;
 	double residual_tolerance_;
+        int linsolver_maxit_;
+        double linsolver_prolongate_factor_;
 	int linsolver_verbosity_;
         int linsolver_type_;
+        int linsolver_smooth_steps_;
 
 	GridType grid_;
 	GridInterface ginterf_;
diff --git a/dune/upscaling/UpscalerBase_impl.hpp b/dune/upscaling/UpscalerBase_impl.hpp
index f277d28..ca4aef0 100644
--- a/dune/upscaling/UpscalerBase_impl.hpp
+++ b/dune/upscaling/UpscalerBase_impl.hpp
@@ -47,8 +47,11 @@ namespace Dune
 	: bctype_(Fixed),
 	  twodim_hack_(false),
 	  residual_tolerance_(1e-8),
+	  linsolver_maxit_(0),
+	  linsolver_prolongate_factor_(1.6),
 	  linsolver_verbosity_(0),
-          linsolver_type_(1)
+          linsolver_type_(1),
+          linsolver_smooth_steps_(2)
     {
     }
 
@@ -84,6 +87,9 @@ namespace Dune
 	residual_tolerance_ = param.getDefault("residual_tolerance", residual_tolerance_);
 	linsolver_verbosity_ = param.getDefault("linsolver_verbosity", linsolver_verbosity_);
         linsolver_type_ = param.getDefault("linsolver_type", linsolver_type_);
+        linsolver_maxit_ = param.getDefault("linsolver_max_iterations", linsolver_maxit_);
+        linsolver_prolongate_factor_ = param.getDefault("linsolver_prolongate_factor", linsolver_prolongate_factor_);
+        linsolver_smooth_steps_ = param.getDefault("linsolver_smooth_steps", linsolver_smooth_steps_);
 
         // Ensure sufficient grid support for requested boundary
         // condition type.
@@ -125,14 +131,20 @@ namespace Dune
                                                   double perm_threshold,
                                                   double z_tolerance,
                                                   double residual_tolerance,
+                                                  int linsolver_maxit,
+                                                  double linsolver_prolongate_factor,
                                                   int linsolver_verbosity,
                                                   int linsolver_type,
-                                                  bool twodim_hack)
+                                                  bool twodim_hack,
+                                                  int linsolver_smooth_steps)
     {
 	bctype_ = bctype;
 	residual_tolerance_ = residual_tolerance;
 	linsolver_verbosity_ = linsolver_verbosity;
         linsolver_type_ = linsolver_type;
+	    linsolver_maxit_ = linsolver_maxit;
+         linsolver_prolongate_factor_ = linsolver_prolongate_factor;
+        linsolver_smooth_steps_ = linsolver_smooth_steps;
 	twodim_hack_ = twodim_hack;
 
 	// Faking some parameters depending on bc type.
@@ -231,7 +243,9 @@ namespace Dune
 	    // Run pressure solver.
             bool same_matrix = (bctype_ != Fixed) && (pdd != 0);
 	    flow_solver_.solve(fluid, sat, bcond_, src, residual_tolerance_,
-                               linsolver_verbosity_, linsolver_type_, same_matrix);
+                               linsolver_maxit_, linsolver_prolongate_factor_,
+                               linsolver_verbosity_, 
+                               linsolver_type_, same_matrix, linsolver_smooth_steps_);
             double max_mod = flow_solver_.postProcessFluxes();
             std::cout << "Max mod = " << max_mod << std::endl;
 
@@ -291,8 +305,8 @@ namespace Dune
 			if (canon_bid - 1 == 2*flow_dir) {
 			    ++num_side1;
 			    if (flow_dir == pdrop_dir && flux > 0.0) {
-				std::cerr << "Flow may be in wrong direction at bid: " << f->boundaryId()
-					  << " Magnitude: " << std::fabs(flux) << std::endl;
+			      std::cerr << "Flow may be in wrong direction at bid: " << f->boundaryId()<<" (canonical: "<<canon_bid
+					  << ") Magnitude: " << std::fabs(flux) << std::endl;
 				// THROW("Detected outflow at entry face: " << face);
 			    }
 			    side1_flux += flux*norm_comp;
diff --git a/examples/Makefile.am b/examples/Makefile.am
index b9f7024..b8205d9 100644
--- a/examples/Makefile.am
+++ b/examples/Makefile.am
@@ -13,6 +13,10 @@ upscale_cap					\
 upscale_cond					\
 upscale_elasticity				\
 upscale_perm					\
+upscale_relperm_ssor                            \
+upscale_relperm_bgs                             \
+upscale_relperm_aniso                           \
+upscale_relperm_aniso_ssor                      \
 upscale_relperm					\
 upscale_relpermvisc                             \
 upscale_steadystate_implicit
@@ -44,6 +48,18 @@ upscale_perm_SOURCES = upscale_perm.cpp
 
 upscale_relperm_SOURCES = upscale_relperm.cpp
 
+upscale_relperm_ssor_SOURCES = upscale_relperm.cpp
+upscale_relperm_ssor_CXXFLAGS = -DSMOOTHER_ILU=0 $(AM_CPPFLAGS)
+
+upscale_relperm_bgs_SOURCES = upscale_relperm.cpp
+upscale_relperm_bgs_CXXFLAGS = -DSMOOTHER_ILU=0 -DSMOOTHER_BGS=1 $(AM_CPPFLAGS)
+
+upscale_relperm_aniso_SOURCES = upscale_relperm.cpp
+upscale_relperm_aniso_CXXFLAGS = -DANISOTROPIC_3D=1
+
+upscale_relperm_aniso_ssor_SOURCES = upscale_relperm.cpp
+upscale_relperm_aniso_ssor_CXXFLAGS = -DANISOTROPIC_3D=1 -DSMOOTHER_ILU=0 $(AM_CPPFLAGS)
+
 upscale_relpermvisc_SOURCES = upscale_relpermvisc.cpp
 
 upscale_steadystate_implicit_SOURCES = upscale_steadystate_implicit.cpp
diff --git a/examples/upscale_relperm.cpp b/examples/upscale_relperm.cpp
index 69803ff..2a88fd4 100644
--- a/examples/upscale_relperm.cpp
+++ b/examples/upscale_relperm.cpp
@@ -246,6 +246,8 @@ int main(int varnum, char** vararg)
                                                               // give so small contributions near endpoints.
    options.insert(make_pair("linsolver_tolerance", "1e-12"));  // residual tolerance for linear solver
    options.insert(make_pair("linsolver_verbosity", "0"));     // verbosity level for linear solver
+   options.insert(make_pair("linsolver_max_iterations", "0"));         // Maximum number of iterations allow, specify 0 for default
+   options.insert(make_pair("linsolver_prolongate_factor", "1.6")); // Factor to scale the prolongate coarse grid correction,
    options.insert(make_pair("linsolver_type",      "1"));     // type of linear solver: 0 = ILU/BiCGStab, 1 = AMG/CG
    options.insert(make_pair("fluids",              "ow")); // wheater upscaling for oil/water (ow) or gas/oil (go)
    options.insert(make_pair("krowxswirr",          "-1")); // relative permeability in x direction of oil in corresponding oil/water system
@@ -253,7 +255,7 @@ int main(int varnum, char** vararg)
    options.insert(make_pair("krowzswirr",          "-1")); // relative permeability in z direction of oil in corresponding oil/water system
    options.insert(make_pair("doEclipseCheck",      "true")); // Check if minimum relpermvalues in input are zero (specify critical saturations)
    options.insert(make_pair("critRelpermThresh",   "1e-6")); // Threshold for setting minimum relperm to 0 (thus specify critical saturations)
-   
+   options.insert(make_pair("linsolver_smooth_steps", "2")); // Number of pre and postsmoothing steps for AMG
 
    // Conversion factor, multiply mD numbers with this to get m² numbers
    const double milliDarcyToSqMetre = 9.869233e-16;
@@ -941,11 +943,15 @@ int main(int varnum, char** vararg)
    double linsolver_tolerance = atof(options["linsolver_tolerance"].c_str());
    int linsolver_verbosity = atoi(options["linsolver_verbosity"].c_str());
    int linsolver_type = atoi(options["linsolver_type"].c_str());
+   int linsolver_maxit = atoi(options["linsolver_max_iterations"].c_str()); 
+   int smooth_steps = atoi(options["linsolver_smooth_steps"].c_str());
+   double linsolver_prolongate_factor = atof(options["linsolver_prolongate_factor"].c_str());
    bool twodim_hack = false;
    eclParser.convertToSI();
    upscaler.init(eclParser, boundaryCondition,
                  Opm::unit::convert::from(minPerm, Opm::prefix::milli*Opm::unit::darcy),
-                 ztol, linsolver_tolerance, linsolver_verbosity, linsolver_type, twodim_hack);
+                 ztol, linsolver_tolerance, linsolver_maxit, linsolver_prolongate_factor,
+                 linsolver_verbosity, linsolver_type, twodim_hack, smooth_steps);
 
    finish = clock();   timeused_tesselation = (double(finish)-double(start))/CLOCKS_PER_SEC;
    if (isMaster) cout << " (" << timeused_tesselation <<" secs)" << endl;