<div dir="ltr">LBFGSB optimizer number of iterations printout error<div><br></div><div style>Hi all,</div><div style><br></div><div style>Im printing out values of metric following Bspline registrations, using an observer connected to the optimizer.</div>
<div style><br></div><div style>However when  using  boptimizer-&gt;SetMaximumNumberOfIterations( 2 ); </div><div style><br></div><div style>I actually get 8 values printed out: see below.  The command iteration update code is also attched.  At first I thought it was because I use multiple observers during each multiresolution step but also removeallobservers() between each multi layer step.  Any suggestions appreciated?</div>
<div style><br></div><div style><br></div><div style><div>Starting Registration with Medium resolution transform</div><div>1.07399e+06   </div><div>bspline</div><div>1.06241e+06   </div><div>bspline</div><div>1.01962e+06   </div>
<div>bspline</div><div>1.01962e+06   </div><div>bspline</div><div>895448   </div><div>bspline</div><div>895448   </div><div>bspline</div><div>878673   </div><div>bspline</div><div>878673   </div><div>bspline</div><div><br>
</div></div><div style><br></div><div style><div> </div><div>//This is for B spline</div><div>  </div><div>  class BCommandIterationUpdate : public itk::Command</div><div>{</div><div>public:</div><div>  typedef  BCommandIterationUpdate   Selfb;</div>
<div>  typedef  itk::Command            Superclassb;</div><div>  typedef itk::SmartPointer&lt;Selfb&gt;   Pointerb;</div><div>  itkNewMacro( Selfb );</div><div>protected:</div><div>  BCommandIterationUpdate() {};</div><div>
public:</div><div>   typedef itk::LBFGSBOptimizer       bOptimizerType;</div><div>  typedef   const bOptimizerType *                  bOptimizerPointer;</div><div><br></div><div>  void Execute(itk::Object *caller, const itk::EventObject &amp; event)</div>
<div>    {</div><div>    Execute( (const itk::Object *)caller, event);</div><div>    }</div><div><br></div><div>  void Execute(const itk::Object * object, const itk::EventObject &amp; event)</div><div>    {</div><div>    bOptimizerPointer boptimizer =</div>
<div>                      dynamic_cast&lt; bOptimizerPointer &gt;( object );</div><div>    if( ! itk::IterationEvent().CheckEvent( &amp;event ) )</div><div>      {</div><div>      return;</div><div>      }</div><div>     //std::cout &lt;&lt; optimizer-&gt;GetCurrentIteration() &lt;&lt; &quot;   &quot;;</div>
<div>     std::cout &lt;&lt; boptimizer-&gt;GetValue() &lt;&lt; &quot;   &quot; &lt;&lt;std::endl;</div><div>     std::cout &lt;&lt; &quot;bspline&quot; &lt;&lt;std::endl;</div><div>     //std::cout &lt;&lt; optimizer-&gt;GetCurrentPosition() &lt;&lt; std::endl;</div>
<div>    }</div><div>};</div><div>  </div></div><div style><br></div><div style><br></div><div style><br></div><div style><br></div><div style><br></div></div>