Skip to content
Snippets Groups Projects
Commit 4f03644d authored by Elias Pipping's avatar Elias Pipping Committed by Elias Pipping
Browse files

Clean up formatting

parent 63428372
No related branches found
No related tags found
No related merge requests found
...@@ -4,6 +4,8 @@ ...@@ -4,6 +4,8 @@
#include <cassert> #include <cassert>
#include <boost/format.hpp>
#include <dune/common/exceptions.hh> #include <dune/common/exceptions.hh>
#include <dune/common/shared_ptr.hh> #include <dune/common/shared_ptr.hh>
#include <dune/common/stdstreams.hh> #include <dune/common/stdstreams.hh>
...@@ -24,10 +26,10 @@ double functionTester(Dune::SampleFunctional<dim> J, ...@@ -24,10 +26,10 @@ double functionTester(Dune::SampleFunctional<dim> J,
1e-12, // requiredResidual: ? 1e-12, // requiredResidual: ?
true, // fastQuadratic true, // fastQuadratic
0); // safety: acceptance factor for inexact minimization 0); // safety: acceptance factor for inexact minimization
std::cout << "Old value: J(...) = " << J(start) << std::endl; double const original = J(start);
Dune::minimise(J, start, runs, bisection); Dune::minimise(J, start, runs, bisection);
double const final = J(start); double const final = J(start);
std::cout << "Final value J(...) = " << final << std::endl; std::cout << boost::format("%8g -> %e") % original % final << std::endl;
return final; return final;
} }
...@@ -160,7 +162,6 @@ void testSampleFunctionNonsmooth() { ...@@ -160,7 +162,6 @@ void testSampleFunctionNonsmooth() {
functionTester(J, start, 6); functionTester(J, start, 6);
} }
std::cout << std::endl;
{ {
start = b; start = b;
start /= (b.two_norm() - 1e-12); // Make sure the norm is above 1; start /= (b.two_norm() - 1e-12); // Make sure the norm is above 1;
...@@ -211,7 +212,6 @@ void testTrivialFunction() { ...@@ -211,7 +212,6 @@ void testTrivialFunction() {
assert(error.two_norm() < 1e-10); assert(error.two_norm() < 1e-10);
double const ret1 = functionTester(J, start, 6); double const ret1 = functionTester(J, start, 6);
std::cout << std::endl;
// Something random // Something random
start[0] = 279; start[0] = 279;
...@@ -479,28 +479,38 @@ void testSteepFunction() { ...@@ -479,28 +479,38 @@ void testSteepFunction() {
int main() { int main() {
try { try {
std::cout << "testIdentity:" << std::endl;
testIdentity(); testIdentity();
std::cout << std::endl;
std::cout << std::endl << "testSampleFunction:" << std::endl;
testSampleFunction(); testSampleFunction();
std::cout << std::endl;
std::cout << std::endl << "testSampleFunctionNonsmooth:" << std::endl;
testSampleFunctionNonsmooth(); testSampleFunctionNonsmooth();
std::cout << std::endl;
std::cout << std::endl << "testTrivialFunction:" << std::endl;
testTrivialFunction(); testTrivialFunction();
std::cout << std::endl;
std::cout << std::endl << "testHorribleFunction:" << std::endl;
testHorribleFunction(); testHorribleFunction();
std::cout << std::endl;
std::cout << std::endl << "testHorribleFunctionLogarithmic:" << std::endl;
testHorribleFunctionLogarithmic(); testHorribleFunctionLogarithmic();
std::cout << std::endl;
std::cout << std::endl << "testSampleFunction2:" << std::endl;
testSampleFunction2(); testSampleFunction2();
std::cout << std::endl;
std::cout << std::endl << "testSampleFunction3D:" << std::endl;
testSampleFunction3D(); testSampleFunction3D();
std::cout << std::endl;
std::cout << std::endl << "testSampleFunctionSteep1:" << std::endl;
testSampleFunctionSteep1(); testSampleFunctionSteep1();
std::cout << std::endl;
std::cout << std::endl << "testSampleFunctionSteep2:" << std::endl;
testSampleFunctionSteep2(); testSampleFunctionSteep2();
std::cout << std::endl;
std::cout << std::endl << "testSteepFunction:" << std::endl;
testSteepFunction(); testSteepFunction();
return 0;
} }
catch (Dune::Exception &e) { catch (Dune::Exception &e) {
Dune::derr << "Dune reported error: " << e << std::endl; Dune::derr << "Dune reported error: " << e << std::endl;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment