diff --git a/src/mynonlinearity.hh b/src/mynonlinearity.hh index 7e2bda65122340323b739f6893e534290a8ef209..b0fb20d57813633f77e2d634dbf29300bc2b3562 100644 --- a/src/mynonlinearity.hh +++ b/src/mynonlinearity.hh @@ -15,8 +15,8 @@ namespace Dune { template <int dimension> class MyNonlinearity { public: - typedef Dune::FieldVector<double, dimension> VectorType; - typedef Dune::FieldMatrix<double, dimension, dimension> MatrixType; + typedef FieldVector<double, dimension> VectorType; + typedef FieldMatrix<double, dimension, dimension> MatrixType; MyNonlinearity(NiceFunction const &func) : func_(func) {} diff --git a/src/samplefunctional.hh b/src/samplefunctional.hh index 9403f9ea639f9f4a906cb31a5549c2b5c3b03713..8e9db33b07dde806443075b07dec776c91c69763 100644 --- a/src/samplefunctional.hh +++ b/src/samplefunctional.hh @@ -16,8 +16,8 @@ namespace Dune { template <int dimension> class SampleFunctional { public: - typedef Dune::FieldVector<double, dimension> SmallVector; - typedef Dune::FieldMatrix<double, dimension, dimension> SmallMatrix; + typedef FieldVector<double, dimension> SmallVector; + typedef FieldMatrix<double, dimension, dimension> SmallMatrix; typedef MyNonlinearity<dimension> NonlinearityType; @@ -64,19 +64,19 @@ template <int dimension> class SampleFunctional { return; } else if (pgx >= 0 && mgx >= 0) { ret = pg; - Dune::dverb << "## Directional derivative (as per scalar product w/ " - "semigradient): " << -(ret * mg) - << " (coordinates of the restriction)" << std::endl; + dverb << "## Directional derivative (as per scalar product w/ " + "semigradient): " << -(ret * mg) + << " (coordinates of the restriction)" << std::endl; } else if (pgx <= 0 && mgx <= 0) { ret = mg; - Dune::dverb << "## Directional derivative (as per scalar product w/ " - "semigradient): " << -(ret * pg) - << " (coordinates of the restriction)" << std::endl; + dverb << "## Directional derivative (as per scalar product w/ " + "semigradient): " << -(ret * pg) + << " (coordinates of the restriction)" << std::endl; } else { ret = project(smoothGradient(x), x); - Dune::dverb << "## Directional derivative (as per scalar product w/ " - "semigradient): " << -(ret * ret) - << " (coordinates of the restriction)" << std::endl; + dverb << "## Directional derivative (as per scalar product w/ " + "semigradient): " << -(ret * ret) + << " (coordinates of the restriction)" << std::endl; } ret *= -1; } @@ -156,7 +156,7 @@ void minimise(const Functional J, const typename Functional::SmallVector x, { // Debug Interval<double> D; JRest.subDiff(0, D); - Dune::dverb + dverb << "## Directional derivative (as per subdifferential of restriction): " << D[1] << " (coordinates of the restriction)" << std::endl; assert(D[1] <= @@ -168,8 +168,8 @@ void minimise(const Functional J, const typename Functional::SmallVector x, int count; // FIXME: The value of x_old should not matter if the factor is 1.0, correct? double const stepsize = bisection.minimize(JRest, 0.0, 1.0, count); - Dune::dverb << "Number of iterations in the bisection method: " << count - << std::endl; + dverb << "Number of iterations in the bisection method: " << count + << std::endl; ; corr = descDir;