diff --git a/dune/tectonic/minimisation.hh b/dune/tectonic/minimisation.hh index ad079b027c187c60e4567603527e6bea8c748e08..7d26a015f5f6b3f0631f8171150ba2926a6cfce9 100644 --- a/dune/tectonic/minimisation.hh +++ b/dune/tectonic/minimisation.hh @@ -19,11 +19,11 @@ void descentMinimisation(Functional const &J, using SmallVector = typename Functional::SmallVector; using LocalNonlinearityType = typename Functional::NonlinearityType; - // {{{ Construct a restriction of J to the line x + t * descDir + // {{{ Construct a restriction of J to the line x + t * v /* 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> + 1/2 <A(x+tv),x+tv>-<b,x+tv> = 1/2 <Av,v> t^2 - <b-Ax,v> t + <1/2 Ax-b,x> since A is symmetric. */