Print

Print


Dear Lars,

I have seen that particular crash before, and then it has been caused by a “nonsensical” topup output file. Can you please send me the d1_1001_b0_topup_fieldcoef.nii.gz file (off mailbase) and I can take a look at it?

Jesper

> On 27 Jul 2017, at 16:59, Lars M. Rimol <[log in to unmask]> wrote:
> 
> Hi,
> 
> I'm trying to run eddy using these parameters:
> 
> $ cat d1_1001_acqParams.txt
> 0 1 0 0.1
> 0 -1 0 0.1
> 
> $ cat ../_eddy_index.txt
> 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
> 
> Here is the command and the error message:
> 
> $ eddy --imain=d1_1001_DTI.nii.gz --mask=d1_1001_topup_mask.nii.gz --acqp=d1_1001_acqParams.txt --index=../_eddy_index.txt --bvecs=d1_1001.bvec --bvals=d1_1001.bval --topup=d1_1001_b0_topup --out=d1_DTI_eddycorrected
> 
> Running Register
> Segmentation violation, Address not mapped, Offending address = 0x41bca6ed0
> eddy SPLINTERPOLATOR::Splinterpolator<float>::value_at(double const*) const [0x523c93]
> eddy NEWIMAGE::volume<float>::splineinterpolate(float, float, float) const [0x56838e]
> eddy NEWIMAGE::volume<float>::interpolate(float, float, float) const [0x568ae0]
> eddy void NEWIMAGE::raw_general_transform<float>(NEWIMAGE::volume<float> const&, NEWMAT::Matrix const&, NEWIMAGE::volume4D<float> const&, std::vector<int, std::allocator<int> > const&, std::vector<int, std::allocator<int> > const&, NEWMAT::Matrix const*, NEWMAT::Matrix const*, NEWIMAGE::volume<float>&, NEWIMAGE::volume4D<float>&, NEWIMAGE::volume<char>*) [0x48f843]
> eddy void NEWIMAGE::apply_warp<float>(NEWIMAGE::volume<float> const&, NEWMAT::Matrix const&, NEWIMAGE::volume4D<float>, NEWMAT::Matrix const&, NEWMAT::Matrix const&, NEWIMAGE::volume<float>&, NEWIMAGE::volume<char>&) [0x49e58e]
> eddy EDDY::EddyUtils::transform_model_to_scan_space(NEWIMAGE::volume<float> const&, EDDY::ECScan const&, boost::shared_ptr<NEWIMAGE::volume<float> const>, bool, NEWIMAGE::volume<float>&, NEWIMAGE::volume<float>*, NEWIMAGE::volume4D<float>*) [0x498625]
> eddy EDDY::EddyUtils::param_update(NEWIMAGE::volume<float> const&, boost::shared_ptr<NEWIMAGE::volume<float> const>, NEWIMAGE::volume<float> const&, EDDY::Parameters, bool, EDDY::ECScan&, NEWMAT::ColumnVector*) [0x49a35b]
> eddy ) [0x46f909] [n
> eddy Register(EDDY::EddyCommandLineOptions const&, EDDY::ScanType, EDDY::ECScanManager&, unsigned int, NEWMAT::Matrix&, NEWMAT::Matrix&) [0x47066b]
> eddy main [0x4713dd]
> /lib64/libc.so.6 __libc_start_main [0x7f92d461cd1d]
> eddy ) [0x46e3d9] [�F
> 
> 
> I've tried running it also with the eddy index as one column. Any idea why this doesn't work? 
> 
> Thank you!