%% analyze net
%-------------

% gain field / receptive field analysis
for lay = 1:2,
    xx = -45:45; % extend of receptive field
    ee = -45:15:45; % different eye positions
    Nxx = length(xx);
    Nee = length(ee);
    TAin = [repmat(xx',Nee,1) reshape(repmat(ee',1,Nxx)',Nxx*Nee,1)]';
    Nana = Nxx*Nee;
    TAout = exp(-(repmat(x',1,Nana)-repmat(TAin(1,:) - TAin(2,:),Ni,1)).^2./10.^2./2);
    Pdana{1,1} = exp(-(repmat(x',1,Nana)-repmat(TAin(1,:),Ni,1)).^2./10.^2./2);
    Pdana{2,1} = [(TAin(2,:)+50)/100; (TAin(2,:)-50)/100];
    ZA = sim(net,Pdana); % simulate all conditions

    figure % plot results (1-D)
    for i = 1:Nee,
        for j = 1:9,
            subplot(3,3,j); hold on;
            plot(xx,ZA{lay}(j,(i-1)*Nxx+(1:Nxx)),'color', i/Nee*[1 1 1]);
            if (i == 1 & j == 8), xlabel('Position (deg)'); end
            if (i == 1 & j == 4), ylabel('Activation (a.u.)'); end
            if (i == 1 & j == 2), title(['\bf{Receptive field changes with eye position, layer ' num2str(lay) '}']); end
            if (i == Nee & j == 3), legend([repmat('eye = ',Nee,1) num2str(ee') repmat(' deg',Nee,1)]); end
        end
    end
    
    figure % plot results (2-D)
    [XX,YY] = meshgrid(ee,xx);
    for j = 1:9,
        subplot(3,3,j);
        surf(XX,YY,reshape(ZA{lay}(j,:),size(XX)));
        view([0 0 1]); shading interp
        if (j == 8), xlabel('Eye position (deg)'); end
        if (j == 4), ylabel('Target position (deg)'); end
        if (j == 2), title(['\bf{Receptive field changes with eye position, layer ' num2str(lay) '}']); end
    end
end