get_current_space

 function J = get_current_space(F,B)

 Return the spacial distribution of the power of the estimated 
 current for the dipoles associated with the inverse filter. 
 2004-05-31 Taku Yoshioka

 --- Input parameters
 F: Inverse filter, I-by-N (sparse) matrix. 
 B: MEG signal, N-by-T matrix. 
 
 (I: number of dipoles)
 (N: number of sensors)
 (T: number of time steps)

 --- Output values
 J : Spacial distribution of the estimated current. 


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

0001 function J = get_current_space(F,B)
0002 %
0003 % function J = get_current_space(F,B)
0004 %
0005 % Return the spacial distribution of the power of the estimated
0006 % current for the dipoles associated with the inverse filter.
0007 % 2004-05-31 Taku Yoshioka
0008 %
0009 % --- Input parameters
0010 % F: Inverse filter, I-by-N (sparse) matrix.
0011 % B: MEG signal, N-by-T matrix.
0012 %
0013 % (I: number of dipoles)
0014 % (N: number of sensors)
0015 % (T: number of time steps)
0016 %
0017 % --- Output values
0018 % J : Spacial distribution of the estimated current.
0019 %
0020 %
0021 % Copyright (C) 2011, ATR All Rights Reserved.
0022 % License : New BSD License(see VBMEG_LICENSE.txt)
0023 
0024 J = sum((F*(B*B')).*F,2);
0025