diff --git a/src/bisection-example-new.cc b/src/bisection-example-new.cc index 49e7321bc502b688961900650778b21b0d8faa37..b9f09477c23b7a998cceb6d4c4f4981b0f7f82a9 100644 --- a/src/bisection-example-new.cc +++ b/src/bisection-example-new.cc @@ -104,6 +104,8 @@ typename Functional::SmallVector minimise( if (descDir == typename Functional::SmallVector(0.0)) return typename Functional::SmallVector(0.0); + // {{{ Construct a restriction of J to the line x + t * descDir + /* We have 1/2 <A(u+xv),u+xv>-<b,u+xv> = 1/2 <Av,v> x^2 - <b-Au,v> x + <1/2 Au-b,u> @@ -124,6 +126,7 @@ typename Functional::SmallVector minimise( typedef DirectionalConvexFunction<MyNonlinearityType> MyDirectionalConvexFunctionType; MyDirectionalConvexFunctionType JRest(JRestA, JRestb, phi, x, descDir); + // }}} // FIXME: default values are used Bisection bisection;