% GENFIG_7C Reproduces figure 7c.
%
% USAGE     
%
%       GENFIG_7C(color, display, filename)
%
% INPUTS    
%
%       color     Color (true) or grayscale (false) (optional, default = true)
%       display   Displays (true) or does not display (false) the figure (optional, default = true)
%       filename  Saves the figure in eps format under the name 'filename' (optional, default = not saved)
%
% COPYRIGHT
%
%       This file is part of the Matlab code provided for the following
%       reproducible paper:
%
%       Olivier Roy and Martin Vetterli,
%       "Dimensionality Reduction for Distributed Estimation in the Infinite Dimensional Regime",
%       vol. 54, no. 4, pp. 1655-1669, April 2008.
%
%       This program is free software; you can redistribute it and/or modify it  
%       under the terms of the GNU General Public License as published by the     
%       Free Software Foundation; either version 2 of the License, or (at your    
%       option) any later version. This software is distributed in the hope that  
%       it will be useful, but without any warranty; without even the implied     
%       warranty of merchantability or fitness for a particular purpose.          
%       See the GNU General Public License for more details                       
%       (enclosed in the file GPL).   
% 
%       GNU General Public License,
%       Copyright (C) 2008,
%       Audiovisual Communications Laboratory (LCAV),
%       Ecole Polytechnique F閐閞ale de Lausanne (EPFL),
%       CH-1015 Lausanne.
%
% COMMENTS
%
%       Author: Olivier Roy
%       Latest modifications: April 2, 2008.
function genfig_7c(color, display, filename)

% We set the default values for the input arguments
save = true;
if nargin < 3
    save = false;
end
if nargin < 2
    display = true;
end
if nargin < 1
    color = true;
end

% We compute the infinite block-length distortion (IBLD) for different values of rho and alpha
rho_val = [0:0.001:0.999];
alpha = 0.1;
IBLD_central_val = (1/2)*(1 - (2/pi)*atan(((1+rho_val.^2)./(1-rho_val.^2))*tan(pi*alpha/2)));
IBLD_partial_val = 1 + (alpha/2)*(1-rho_val.^2)./(1+rho_val.^2) - (2/pi)*atan(((1+rho_val.^2)./(1-rho_val.^2))*tan(pi*alpha/2));

% We plot the figure
figure;
if (color)
    plot(rho_val, IBLD_partial_val-IBLD_central_val, 'b-');
else
    plot(rho_val, IBLD_partial_val-IBLD_central_val, 'k-');
end
box on
grid on

if save
    xlabel('xLabel');
    ylabel('yLabel');
    title('title');
    saveas(gcf, filename,'psc2');
end
if display
    xlabel('\rho');
    ylabel('Loss');
else
    close gcf
end