[Insight-users] [bug] segfault when enabling verbose in Examples/Patented/IterativeClosestPoint3.cxx

[Ramón Casero Cañas]

Maybe this is something that I'm doing wrong, but I believe this is a 
bug, possibly somewhere between 
Utilities/vxl/core/vnl/algo/vnl_levenberg_marquardt.* and 
Utilities/vxl/core/vnl/vnl_nonlinear_minimizer.h.

ITK v3.14.0
Platform: Ubuntu GNU/Linux 64 bit

Building Examples/Patented/IterativeClosestPoint3.cxx, and running it 
with the provided example data works

<SCREEN>
$ ./IterativeClosestPoint3 
~/Desktop/InsightToolkit-3.14.0/Examples/Data/IterativeClosestPointFixedPoints.txt 
~/Desktop/InsightToolkit-3.14.0/Examples/Data/IterativeClosestPointMovingPoints.txt
Number of fixed Points = 8
Number of moving Points = 8
Solution = [-2, -2]
</SCREEN>


* To reproduce the error:

Enable the verbose flag in vnl_nonlinear_minimizer.h, adding line

<CODE>
    optimizer->GetOptimizer()->set_verbose( true );
</CODE>

right after

<CODE>
   // Optimizer Type
   typedef itk::LevenbergMarquardtOptimizer OptimizerType;

   OptimizerType::Pointer      optimizer     = OptimizerType::New();
   optimizer->SetUseCostFunctionGradient(false);
</CODE>

Recompiling and running now throws a segfault

<SCREEN>
$ ./IterativeClosestPoint3 
~/Desktop/InsightToolkit-3.14.0/Examples/Data/IterativeClosestPointFixedPoints.txt 
~/Desktop/InsightToolkit-3.14.0/Examples/Data/IterativeClosestPointMovingPoints.txt 

Number of fixed Points = 8
Number of moving Points = 8
Segmentation fault
</SCREEN>

Cheers,

Ramon.


-- 
Ramón Casero Cañas, DPhil

Computational Biology, Computing Laboratory
University of Oxford
Wolfson Building, Parks Rd
Oxford OX1 3QD

tlf     +44 (0) 1865 610807
web     http://web.comlab.ox.ac.uk/people/Ramon.CaseroCanas
photos  http://www.flickr.com/photos/rcasero/