vb_load_current_ext

 load estimated current for extra dipoles

 --- Syntax
 function [Jext,Pointlist] = vb_load_current_ext(currfile,ave_mode,ix_trial)

 --- Input
  currfile  - current file name

 --- Output
 Jext     : estimated current for extra dipoles
 Pointlist: Information of extra dipoles

 --- History
 2008-08-21 Taku Yoshioka

 Copyright (C) 2011, ATR All Rights Reserved.
 License : New BSD License(see VBMEG_LICENSE.txt)

0001 function [Jext,Plist] = vb_load_current_ext(currfile,ave_mode,ix_trial,verbose)
0002 % load estimated current for extra dipoles
0003 %
0004 % --- Syntax
0005 % function [Jext,Pointlist] = vb_load_current_ext(currfile,ave_mode,ix_trial)
0006 %
0007 % --- Input
0008 %  currfile  - current file name
0009 %
0010 % --- Output
0011 % Jext     : estimated current for extra dipoles
0012 % Pointlist: Information of extra dipoles
0013 %
0014 % --- History
0015 % 2008-08-21 Taku Yoshioka
0016 %
0017 % Copyright (C) 2011, ATR All Rights Reserved.
0018 % License : New BSD License(see VBMEG_LICENSE.txt)
0019 
0020 const = vb_define_verbose;
0021 VERBOSE_LEVEL_INFO = const.VERBOSE_LEVEL_INFO;
0022 
0023 if nargin<4, verbose = 2; end
0024 if nargin<3, ix_trial = []; end
0025 if nargin<2 | isempty(ave_mode), ave_mode = true; end
0026 
0027 load(currfile,'Pointlist');
0028 Plist = Pointlist;
0029 
0030 if ave_mode==OFF
0031   % trial mode
0032   load(currfile,'Jinfo');
0033   
0034   if isfield(Jinfo,'jactdir'),
0035     jactdir = [fileparts(currfile) '/' Jinfo.jactdir];
0036     if exist(jactdir,'dir')==0,
0037       error('No trial data directory')
0038     end
0039     
0040     % Session & Trial number list
0041     session = [];
0042     ntrial = [];
0043     
0044     for n=1:Jinfo.Nsession
0045       Ntry = Jinfo.Ntrial(n);
0046       session = [session, repmat(n,1,Ntry)];
0047       ntrial  = [ntrial, 1:Ntry];
0048     end
0049     
0050     % Trial loop
0051     if isempty(ix_trial), ix_trial = 1:length(ntrial); end
0052     Ntrial = length(ix_trial);
0053 
0054     for n=1:Ntrial
0055       nt = ix_trial(n);
0056       vb_disp(sprintf('Session: %04d , Trial: %04d', ...
0057                       session(nt), ntrial(nt)), VERBOSE_LEVEL_INFO);
0058       %if verbose==1,
0059       %  fprintf('.')
0060       %elseif verbose==2,
0061       %  fprintf('Session: %04d , Trial: %04d\n',session(nt), ntrial(nt))
0062       %end
0063       
0064       fname = sprintf('data_s%04dt%04d', session(nt), ntrial(nt));
0065       load([jactdir '/' fname],'Jext');    
0066       Jall(:,:,n) = Jext;
0067     end
0068     
0069     Jext = Jall;
0070   else
0071     load(currfile,'Jext');
0072     if ~isempty(ix_trial), Jext = Jext(:,:,ix_trial); end
0073   end
0074 else
0075   load(currfile,'Jext');
0076 end