0001 function cross_trial = fft2csdm(fft_trialdata) 0002 0003 cross_trial = zeros(fix(size(fft_trialdata,1)),(size(fft_trialdata,2).^2+size(fft_trialdata,2))/2); 0004 icount = 1; 0005 0006 if strcmp(computer, 'MAC2') 0007 for i = 1:size(fft_trialdata,2) 0008 for j = i:size(fft_trialdata,2) 0009 cross_trial(:,icount) = ppc_cmplx_mult(fft_trial(1:fix(size(fft_trialdata,1)),i),conj(fft_trial(1:size(fft_trialdata,1),j))); 0010 icount = icount+1; 0011 end; 0012 end; 0013 else 0014 for i = 1:size(fft_trialdata,2) 0015 for j = i:size(fft_trialdata,2) 0016 cross_trial(:,icount) = fft_trialdata(1:fix(size(fft_trialdata,1)),i).*conj(fft_trialdata(1:fix(size(fft_trialdata,1)),j)); 0017 icount = icount+1; 0018 end; 0019 end; 0020 end 0021