[Insight-developers] problem with itk::LinearInterpolateImageFunction

[pj gouttenoire]

Le 29/06/2012 11:15, pj gouttenoire a écrit :
> Hello,
>
> i use itk::RegularStepGradientDescentBaseOptimizer for computing path 
> between structures, the first version is doing in a for loop and one 
> by one and it works, but it's very long so i use openMP for parallel 
> computation. I made the directive fot the for loop and i launch it in 
> multprocessor and i have this seg fault
>
> Program received signal SIGSEGV, Segmentation fault.
>
> [Switching to Thread 1126189376 (LWP 17584)]
> 0x00000000004be7d4 in 
> itk::LinearInterpolateImageFunction<itk::Image<double, 3u>, 
> double>::EvaluateAtContinuousIndex (
> this=0x2aaabc028150, index=@0x43202c70)
> at 
> /data/id19/bones/Creatis/All/thirdparty_install/include/InsightToolkit/Common/itkLinearInterpolateImageFunction.txx:150
> 150 value += static_cast<RealType>( this->GetInputImage()->GetPixel( 
> neighIndex ) ) * overlap;
> (gdb) bt
> #0 0x00000000004be7d4 in 
> itk::LinearInterpolateImageFunction<itk::Image<double, 3u>, 
> double>::EvaluateAtContinuousIndex (
> this=0x2aaabc028150, index=@0x43202c70)
> at 
> /data/id19/bones/Creatis/All/thirdparty_install/include/InsightToolkit/Common/itkLinearInterpolateImageFunction.txx:150
> #1 0x00000000004bcc9e in 
> itk::InterpolateImageFunction<itk::Image<double, 3u>, 
> double>::Evaluate (this=0x2aaabc028150,
> point=@0x43202d00)
> at 
> /data/id19/bones/Creatis/All/thirdparty_install/include/InsightToolkit/Common/itkInterpolateImageFunction.h:93
> #2 0x000000000046c406 in 
> itk::PhysicalCentralDifferenceImageFunction<itk::Image<double, 3u>, 
> double>::Evaluate (
> this=0x2aaabc027360, point=@0x43202db0)
> at 
> /data/id19/bones2/maria/linux/source/lcNetwork/Src/FastMarching/itkPhysicalCentralDifferenceImageFunction.txx:65
> #3 0x00000000004b8dbe in 
> itk::SingleImageCostFunction<itk::Image<double, 3u> >::GetDerivative 
> (this=0x2aaab8033a20,
> parameters=@0x2aaab9740790, derivative=@0x2aaab97407e0)
> at 
> /data/id19/bones2/maria/linux/source/lcNetwork/Src/FastMarching/itkSingleImageCostFunction.txx:136
> #4 0x00000000004627f2 in 
> itk::SingleValuedCostFunction::GetValueAndDerivative 
> (this=0x2aaab8033a20,
> parameters=@0x2aaab9740790, value=@0x2aaab9740818, 
> derivative=@0x2aaab97407e0)
> at 
> /data/id19/bones/Creatis/All/thirdparty_install/include/InsightToolkit/Numerics/itkSingleValuedCostFunction.h:76
> #5 0x00002aaaabd4fc57 in 
> itk::RegularStepGradientDescentBaseOptimizer::ResumeOptimization ()
> from 
> /data/id19/bones/Creatis/All/thirdparty_install/lib/InsightToolkit/libITKNumerics.so.3.20
> #6 0x00002aaaabd503b9 in 
> itk::RegularStepGradientDescentBaseOptimizer::StartOptimization ()
> from 
> /data/id19/bones/Creatis/All/thirdparty_install/lib/InsightToolkit/libITKNumerics.so.3.20
> #7 0x00000000004ba2ce in 
> itk::ArrivalFunctionToPathFilter<itk::Image<double, 3u>, 
> itk::PolyLineParametricPath<3u> >::GenerateData (this=0x2aaabc026f60)
> at 
> /data/id19/bones2/maria/linux/source/lcNetwork/Src/FastMarching/itkArrivalFunctionToPathFilter.txx:199
> #8 0x00000000004ba752 in 
> itk::SpeedFunctionToPathFilter<itk::Image<double, 3u>, 
> itk::PolyLineParametricPath<3u> >::GenerateData (this=0x2aaabc026f60)
> at 
> /data/id19/bones2/maria/linux/source/lcNetwork/Src/FastMarching/itkSpeedFunctionToPathFilter.txx:185
> #9 0x00002aaab313e7af in itk::ProcessObject::UpdateOutputData ()
> from 
> /data/id19/bones/Creatis/All/thirdparty_install/lib/InsightToolkit/libITKCommon.so.3.20
> #10 0x00002aaab30fb192 in itk::DataObject::UpdateOutputData ()
> from 
> /data/id19/bones/Creatis/All/thirdparty_install/lib/InsightToolkit/libITKCommon.so.3.20
> #11 0x00002aaab30fb0b6 in itk::DataObject::Update ()
> from 
> /data/id19/bones/Creatis/All/thirdparty_install/lib/InsightToolkit/libITKCommon.so.3.20
> #12 0x00002aaab313dd51 in itk::ProcessObject::Update ()
> from 
> /data/id19/bones/Creatis/All/thirdparty_install/lib/InsightToolkit/libITKCommon.so.3.20
> #13 0x00000000004605da in calculFastMarching (startfile=
> {static npos = 18446744073709551615, _M_dataplus = 
> {<std::allocator<char>> = {<__gnu_cxx::new_allocator<char>> = {<No 
> data fields>}, <No data fields>}, _M_p = 0x43204050 "�q'\n"}}, endfile=
> {static npos = 18446744073709551615, _M_dataplus = 
> {<std::allocator<char>> = {<__gnu_cxx::new_allocator<char>> = {<No 
> data fields>}, <No data fields>}, _M_p = 0x43204040 "\210r'\n"}}, 
> startindex=1, debug=1, advanced=false,
> vote=0.050000000000000003, inputfile=
> {static npos = 18446744073709551615, _M_dataplus = 
> {<std::allocator<char>> = {<__gnu_cxx::new_allocator<char>> = {<No 
> data fields>}, <No data fields>}, _M_p = 0x43204030 "�p'\n"}}, 
> addImage={m_Pointer = 0x43204020})
> at 
> /data/id19/bones2/maria/linux/source/lcNetwork/Apps/FastMarching/fastmarchingTestMP.cxx:306
> #14 0x000000000046210c in main.omp_fn.0 (.omp_data_i=0x7fff042f66b0)
> at 
> /data/id19/bones2/maria/linux/source/lcNetwork/Apps/FastMarching/fastmarchingTestMP.cxx:191
> #15 0x00002aaab59ce4ac in GOMP_single_copy_start () from 
> /usr/lib64/libgomp.so.1
> #16 0x00002aaab479d2e7 in start_thread () from /lib64/libpthread.so.0
> #17 0x0000003caf4ce3bd in clone () from /lib64/libc.so.6
>
>
> _______________________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
>
> Kitware offers ITK Training Courses, for more information visit:
> http://kitware.com/products/protraining.php
>
> Please keep messages on-topic and check the ITK FAQ at:
> http://www.itk.org/Wiki/ITK_FAQ
>
> Follow this link to subscribe/unsubscribe:
> http://www.itk.org/mailman/listinfo/insight-developers
>
i'm using itk 3.20