Constness

src/test-circle.cc

@@ -26,12 +26,12 @@ int main() {
   typedef Functional::SmallMatrix SmallMatrix;
   typedef Functional::SmallVector SmallVector;
 
-  SmallMatrix A = { { 4, 1.5 }, { 1.5, 3 } };
-  SmallVector b = { 1, 2 };
+  SmallMatrix const A = { { 4, 1.5 }, { 1.5, 3 } };
+  SmallVector const b = { 1, 2 };
 
   auto const f = Dune::make_shared<Dune::SampleFunction<2> const>();
   auto const phi = Dune::make_shared<Functional::NonlinearityType const>(f);
-  Functional J(A, b, phi);
+  Functional const J(A, b, phi);
 
   Bisection const bisection(0.0, 1.0, 1e-12, false, 0);
 

src/test-gradient-horrible-logarithmic.cc

@@ -19,12 +19,12 @@ int main() {
   typedef Functional::SmallMatrix SmallMatrix;
   typedef Functional::SmallVector SmallVector;
 
-  SmallMatrix A = { { 3, 1.5 }, { 1.5, 4 } };
-  SmallVector b = { 1, 2 };
+  SmallMatrix const A = { { 3, 1.5 }, { 1.5, 4 } };
+  SmallVector const b = { 1, 2 };
 
   auto const f = Dune::make_shared<Dune::HorribleFunctionLogarithmic const>();
   auto const phi = Dune::make_shared<Functional::NonlinearityType const>(f);
-  Functional J(A, b, phi);
+  Functional const J(A, b, phi);
 
   double ret1;
   {

src/test-gradient-horrible.cc

@@ -19,12 +19,12 @@ int main() {
   typedef Functional::SmallMatrix SmallMatrix;
   typedef Functional::SmallVector SmallVector;
 
-  SmallMatrix A = { { 3, 1.5 }, { 1.5, 4 } };
-  SmallVector b = { 1, 2 };
+  SmallMatrix const A = { { 3, 1.5 }, { 1.5, 4 } };
+  SmallVector const b = { 1, 2 };
 
   auto const f = Dune::make_shared<Dune::HorribleFunction const>();
   auto const phi = Dune::make_shared<Functional::NonlinearityType const>(f);
-  Functional J(A, b, phi);
+  Functional const J(A, b, phi);
 
   double ret1;
   {

src/test-gradient-identity.cc

@@ -19,12 +19,12 @@ int main() {
   typedef Functional::SmallMatrix SmallMatrix;
   typedef Functional::SmallVector SmallVector;
 
-  SmallMatrix A = { { 3, 1.5 }, { 1.5, 4 } };
-  SmallVector b = { 1, 2 };
+  SmallMatrix const A = { { 3, 1.5 }, { 1.5, 4 } };
+  SmallVector const b = { 1, 2 };
 
   auto const f = Dune::make_shared<Dune::LinearFunction const>(1);
   auto const phi = Dune::make_shared<Functional::NonlinearityType const>(f);
-  Functional J(A, b, phi);
+  Functional const J(A, b, phi);
 
   double ret1;
   {

src/test-gradient-parabola.cc

@@ -21,15 +21,15 @@ int main() {
   typedef Functional::SmallMatrix SmallMatrix;
   typedef Functional::SmallVector SmallVector;
 
-  SmallMatrix A = { { 3, 1.5 }, { 1.5, 4 } };
-  SmallVector b = { 1, 2 };
+  SmallMatrix const A = { { 3, 1.5 }, { 1.5, 4 } };
+  SmallVector const b = { 1, 2 };
 
   // |x|^2 as the nonlinearity is the same as having no nonlinearity
   // but twice the identity matrix added to A. In other words, we're
   // solving A + 2*id = b
   auto const f = Dune::make_shared<Dune::Parabola const>();
   auto const phi = Dune::make_shared<Functional::NonlinearityType const>(f);
-  Functional J(A, b, phi);
+  Functional const J(A, b, phi);
 
   SmallVector solution = { 4.0 / 37.0, 34.0 / 111.0 }; // Analytic solution
   SmallMatrix M = A;

src/test-gradient-sample-3d.cc

@@ -19,12 +19,12 @@ int main() {
   typedef Functional::SmallMatrix SmallMatrix;
   typedef Functional::SmallVector SmallVector;
 
-  SmallMatrix A = { { 3, 1.5, 0 }, { 1.5, 4, 1.5 }, { 0, 1.5, 5 } };
-  SmallVector b = { 1, 2, 3 };
+  SmallMatrix const A = { { 3, 1.5, 0 }, { 1.5, 4, 1.5 }, { 0, 1.5, 5 } };
+  SmallVector const b = { 1, 2, 3 };
 
   auto const f = Dune::make_shared<Dune::SampleFunction<2> const>();
   auto const phi = Dune::make_shared<Functional::NonlinearityType const>(f);
-  Functional J(A, b, phi);
+  Functional const J(A, b, phi);
 
   double ret1;
   {

src/test-gradient-sample-nonsmooth.cc

@@ -20,12 +20,12 @@ int main() {
   typedef Functional::SmallMatrix SmallMatrix;
   typedef Functional::SmallVector SmallVector;
 
-  SmallMatrix A = { { 3, 1.5 }, { 1.5, 4 } };
-  SmallVector b = { 1, 2 };
+  SmallMatrix const A = { { 3, 1.5 }, { 1.5, 4 } };
+  SmallVector const b = { 1, 2 };
 
   auto const f = Dune::make_shared<Dune::SampleFunction<2> const>();
   auto const phi = Dune::make_shared<Functional::NonlinearityType const>(f);
-  Functional J(A, b, phi);
+  Functional const J(A, b, phi);
 
   SmallVector start;
   /*

src/test-gradient-sample-steep.cc

@@ -19,12 +19,12 @@ int main() {
   typedef Functional::SmallMatrix SmallMatrix;
   typedef Functional::SmallVector SmallVector;
 
-  SmallMatrix A = { { 1, 0 }, { 0, 1 } };
-  SmallVector b = { 1, 2 };
+  SmallMatrix const A = { { 1, 0 }, { 0, 1 } };
+  SmallVector const b = { 1, 2 };
 
   auto const f = Dune::make_shared<Dune::SampleFunction<2> const>();
   auto const phi = Dune::make_shared<Functional::NonlinearityType const>(f);
-  Functional J(A, b, phi);
+  Functional const J(A, b, phi);
 
   double ret1;
   {

src/test-gradient-sample-steep2.cc

@@ -19,12 +19,12 @@ int main() {
   typedef Functional::SmallMatrix SmallMatrix;
   typedef Functional::SmallVector SmallVector;
 
-  SmallMatrix A = { { 1, 0 }, { 0, 1 } };
-  SmallVector b = { 1, 2.5 };
+  SmallMatrix const A = { { 1, 0 }, { 0, 1 } };
+  SmallVector const b = { 1, 2.5 };
 
   auto const f = Dune::make_shared<Dune::SampleFunction<2> const>();
   auto const phi = Dune::make_shared<Functional::NonlinearityType const>(f);
-  Functional J(A, b, phi);
+  Functional const J(A, b, phi);
 
   double ret1;
   {

src/test-gradient-sample-verysteep.cc

@@ -19,12 +19,12 @@ int main() {
   typedef Functional::SmallMatrix SmallMatrix;
   typedef Functional::SmallVector SmallVector;
 
-  SmallMatrix A = { { 1, 0 }, { 0, 1 } };
-  SmallVector b = { 1, 2.5 };
+  SmallMatrix const A = { { 1, 0 }, { 0, 1 } };
+  SmallVector const b = { 1, 2.5 };
 
   auto const f = Dune::make_shared<Dune::SampleFunction<100> const>();
   auto const phi = Dune::make_shared<Functional::NonlinearityType const>(f);
-  Functional J(A, b, phi);
+  Functional const J(A, b, phi);
 
   double ret1;
   {

src/test-gradient-sample.cc

@@ -21,12 +21,12 @@ int main() {
   typedef Functional::SmallMatrix SmallMatrix;
   typedef Functional::SmallVector SmallVector;
 
-  SmallMatrix A = { { 3, 1.5 }, { 1.5, 4 } };
-  SmallVector b = { 1, 2 };
+  SmallMatrix const A = { { 3, 1.5 }, { 1.5, 4 } };
+  SmallVector const b = { 1, 2 };
 
   auto const f = Dune::make_shared<Dune::SampleFunction<2> const>();
   auto const phi = Dune::make_shared<Functional::NonlinearityType const>(f);
-  Functional J(A, b, phi);
+  Functional const J(A, b, phi);
 
   /*
     j(x)

src/test-gradient-sample2.cc

@@ -21,12 +21,12 @@ int main() {
   typedef Functional::SmallMatrix SmallMatrix;
   typedef Functional::SmallVector SmallVector;
 
-  SmallMatrix A = { { 1, 0 }, { 0, 1 } };
-  SmallVector b = { 1, 1 };
+  SmallMatrix const A = { { 1, 0 }, { 0, 1 } };
+  SmallVector const b = { 1, 1 };
 
   auto const f = Dune::make_shared<Dune::SampleFunction<2> const>();
   auto const phi = Dune::make_shared<Functional::NonlinearityType const>(f);
-  Functional J(A, b, phi);
+  Functional const J(A, b, phi);
 
   double ret1;
   {

src/test-gradient-trivial.cc

@@ -20,12 +20,12 @@ int main() {
   typedef Functional::SmallMatrix SmallMatrix;
   typedef Functional::SmallVector SmallVector;
 
-  SmallMatrix A = { { 3, 1.5 }, { 1.5, 4 } };
-  SmallVector b = { 1, 2 };
+  SmallMatrix const A = { { 3, 1.5 }, { 1.5, 4 } };
+  SmallVector const b = { 1, 2 };
 
   auto const f = Dune::make_shared<Dune::TrivialFunction const>();
   auto const phi = Dune::make_shared<Functional::NonlinearityType const>(f);
-  Functional J(A, b, phi);
+  Functional const J(A, b, phi);
 
   SmallVector solution = { 4.0 / 39.0, 6.0 / 13.0 }; // Analytic solution