diff --git a/dune/tectonic/circularconvexfunction.hh b/dune/tectonic/circularconvexfunction.hh
index c7fdc54736498f71e40e0882fecf29e3fc33c202..ba4f4071ea70fafee2ed07cd9e93977d871ba786 100644
--- a/dune/tectonic/circularconvexfunction.hh
+++ b/dune/tectonic/circularconvexfunction.hh
@@ -42,8 +42,8 @@ template <class NonlinearityType> class CircularConvexFunction {
   }
 
   void domain(Interval<double> &domain) const {
-    domain[0] = -M_PI;
-    domain[1] = M_PI;
+    domain[0] = -2 * M_PI;
+    domain[1] = 2 * M_PI;
   }
 
   void cartesian(double m, VectorType &y) const {
diff --git a/src/test-circle.cc b/src/test-circle.cc
index 341b702a019a0cdd5bc00be8a046bfdc16e67e03..19bbb3f6151b3a93ddbe10220bbc329ae47fa5f5 100644
--- a/src/test-circle.cc
+++ b/src/test-circle.cc
@@ -57,19 +57,16 @@ int main() {
   }
 
   double const intervals = 10000;
-  double const iterations = 2;
   for (size_t i = 0; i < intervals; ++i) {
     double alpha = i / (double)intervals * 2 * M_PI;
     SmallVector x;
     x[0] = scale * std::sin(alpha);
     x[1] = scale * std::cos(alpha);
-    SmallVector descDir;
 
-    for (int i = 0; i < iterations; ++i) {
-      descDir[0] = x[1];
-      descDir[1] = -x[0];
-      tangentialMinimisation(J, x, descDir, bisection);
-    }
+    SmallVector descDir;
+    descDir[0] = x[1];
+    descDir[1] = -x[0];
+    tangentialMinimisation(J, x, descDir, bisection);
 
     bool minimum_hit(false);
     for (auto const &minimum : minima) {