Back to the main page.

Bug 409 - ft_compute_leadfield crashes without vol.skin_surface field for bemcp BEM

Status CLOSED INVALID
Reported 2011-01-18 11:28:00 +0100
Modified 2011-01-26 15:42:15 +0100
Product: FieldTrip
Component: forward
Version: unspecified
Hardware: PC
Operating System: Windows
Importance: P1 critical
Assigned to: Robert Oostenveld
URL:
Tags:
Depends on:
Blocks:
See also:

Vladimir Litvak - 2011-01-18 11:28:17 +0100

A new field was added to the bemcp volume (skin_surface) and ft_compute_leadfield crashes without it. This is not good as SPM saves pre-computed vols on disk. Backward compatibility should be ensured.


Robert Oostenveld - 2011-01-20 11:26:41 +0100

Hi Vladimir ft_compute_leadfield does not use it, but ft_prepare_vol_sens does. why is this part (line 73) not executed properly? % determine the skin compartment if ~isfield(vol, 'skin_surface') if isfield(vol, 'bnd') vol.skin_surface = find_outermost_boundary(vol.bnd); elseif isfield(vol, 'r') && length(vol.r)<=4 [dum, vol.skin_surface] = max(vol.r); end end and then in line 345 it once more does if ~isfield(vol, 'skin_surface') vol.skin_surface = find_outermost_boundary(vol.bnd); fprintf('determining skin compartment (%d)\n', vol.skin_surface); end


Vladimir Litvak - 2011-01-25 16:58:55 +0100

Hi Robert, The bug seems to be gone now after updating again. Vladimir (In reply to comment #1) > Hi Vladimir > > ft_compute_leadfield does not use it, but ft_prepare_vol_sens does. > > why is this part (line 73) not executed properly? > > % determine the skin compartment > if ~isfield(vol, 'skin_surface') > if isfield(vol, 'bnd') > vol.skin_surface = find_outermost_boundary(vol.bnd); > elseif isfield(vol, 'r') && length(vol.r)<=4 > [dum, vol.skin_surface] = max(vol.r); > end > end > > and then in line 345 it once more does > > if ~isfield(vol, 'skin_surface') > vol.skin_surface = find_outermost_boundary(vol.bnd); > fprintf('determining skin compartment (%d)\n', vol.skin_surface); > end


Robert Oostenveld - 2011-01-25 23:10:59 +0100

if it cannot be reproduced, then lets close it. best, Robert


Robert Oostenveld - 2011-01-26 15:41:05 +0100

I closed all the bugs that I resolved earlier today. These have been discussed at the Fieldtrip meeting and should not show up any more at the next meeting. If you feel that one of these bugs should not have been closed, you have to reopen it.


Robert Oostenveld - 2011-01-26 15:42:15 +0100

<<once more, I forgot to actually change the status to CLOSED>> I closed all the bugs that I resolved earlier today. These have been discussed at the Fieldtrip meeting and should not show up any more at the next meeting. If you feel that one of these bugs should not have been closed, you have to reopen it.