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

Rename restrictions

parent e5c7ec7d
No related branches found
No related tags found
No related merge requests found
...@@ -113,23 +113,23 @@ typename Functional::SmallVector minimise( ...@@ -113,23 +113,23 @@ typename Functional::SmallVector minimise(
typename Functional::SmallVector tmp; typename Functional::SmallVector tmp;
J.A.mv(descDir, tmp); // Av J.A.mv(descDir, tmp); // Av
double const rest_A = tmp * descDir; // <Av,v> double const JRestA = tmp * descDir; // <Av,v>
J.A.mv(x, tmp); // Au J.A.mv(x, tmp); // Au
double const rest_b = (J.b - tmp) * descDir; // <b-Au,v> double const JRestb = (J.b - tmp) * descDir; // <b-Au,v>
typedef MyNonlinearity<Functional::SmallVector::dimension, typedef MyNonlinearity<Functional::SmallVector::dimension,
typename Functional::FunctionType> MyNonlinearityType; typename Functional::FunctionType> MyNonlinearityType;
MyNonlinearityType phi; MyNonlinearityType phi;
typedef DirectionalConvexFunction<MyNonlinearityType> typedef DirectionalConvexFunction<MyNonlinearityType>
MyDirectionalConvexFunctionType; MyDirectionalConvexFunctionType;
MyDirectionalConvexFunctionType rest(rest_A, rest_b, phi, x, descDir); MyDirectionalConvexFunctionType JRest(JRestA, JRestb, phi, x, descDir);
// FIXME: default values are used // FIXME: default values are used
Bisection bisection; Bisection bisection;
int count; int count;
// FIXME: does x_old = 1 make any sense?! // FIXME: does x_old = 1 make any sense?!
double const m = bisection.minimize(rest, 0.0, 1.0, count); double const m = bisection.minimize(JRest, 0.0, 1.0, count);
Dune::dverb << "Number of iterations in the bisection method: " << count Dune::dverb << "Number of iterations in the bisection method: " << count
<< std::endl; << std::endl;
; ;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment