function varargout = comparesessions(varargin)
% Implements Compare Sessions GUI
% COMPARESESSIONS M-file for comparesessions.fig
%      COMPARESESSIONS, by itself, creates a new COMPARESESSIONS or raises the existing
%      singleton*.
%
%      H = COMPARESESSIONS returns the handle to a new COMPARESESSIONS or the handle to
%      the existing singleton*.
%
%      COMPARESESSIONS('CALLBACK',hObject,eventData,handles,...) calls the local
%      function named CALLBACK in COMPARESESSIONS.M with the given input arguments.
%
%      COMPARESESSIONS('Property','Value',...) creates a new COMPARESESSIONS or raises the
%      existing singleton*.  Starting from the left, property value pairs are
%      applied to the GUI before comparesessions_OpeningFunction gets called.  An
%      unrecognized property name or invalid value makes property application
%      stop.  All inputs are passed to comparesessions_OpeningFcn via varargin.
%
%      *See GUI Options on GUIDE's Tools menu.  Choose "GUI allows only one
%      instance to run (singleton)".
%
% See also: GUIDE, GUIDATA, GUIHANDLES

% Edit the above text to modify the response to help comparesessions

% Last Modified by GUIDE v2.5 22-May-2003 12:54:40
% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name',       mfilename, ...
                   'gui_Singleton',  gui_Singleton, ...
                   'gui_OpeningFcn', @comparesessions_OpeningFcn, ...
                   'gui_OutputFcn',  @comparesessions_OutputFcn, ...
                   'gui_LayoutFcn',  [] , ...
                   'gui_Callback',   []);
if nargin & isstr(varargin{1})
    gui_State.gui_Callback = str2func(varargin{1});
end

if nargout
    [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
    gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT


% --- Executes just before comparesessions is made visible.
function comparesessions_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject    handle to figure
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
% varargin   command line arguments to comparesessions (see VARARGIN)

% Choose default command line output for comparesessions
handles.output = hObject;
% Set the background of the figure to White
set(hObject,'Color','white');
% Set the toolbar for the figure
pronytoolbar(hObject,'on');
% Flag to check if the GUI is picking up data from file or from workspace
handles.flag_menuopen=0;
handles.idennz_mat=0;
handles.Index_Selected=0;
% Update handles structure
guidata(hObject, handles);

% UIWAIT makes comparesessions wait for user response (see UIRESUME)
% uiwait(handles.figure1);



% --- Outputs from this function are returned to the command line.
function varargout = comparesessions_OutputFcn(hObject, eventdata, handles)
% varargout  cell array for returning output args (see VARARGOUT);
% hObject    handle to figure
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Get default command line output from handles structure
varargout{1} = handles.output;
% --- Executes during object creation, after setting all properties.function lstbx_session_CreateFcn(hObject, eventdata, handles)% hObject    handle to lstbx_session (see GCBO)% eventdata  reserved - to be defined in a future version of MATLAB% handles    empty - handles not created until after all CreateFcns called% Hint: listbox controls usually have a white background on Windows.%       See ISPC and COMPUTER.if ispc    set(hObject,'BackgroundColor','white');else    set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor'));end
% --- Executes on selection change in lstbx_session.function lstbx_session_Callback(hObject, eventdata, handles)% hObject    handle to lstbx_session (see GCBO)% eventdata  reserved - to be defined in a future version of MATLAB% handles    structure with handles and user data (see GUIDATA)% Hints: contents = get(hObject,'String') returns lstbx_session contents as cell array%        contents{get(hObject,'Value')} returns selected item from lstbx_session% --- Executes on button press in push_close.function push_close_Callback(hObject, eventdata, handles)% hObject    handle to push_close (see GCBO)% eventdata  reserved - to be defined in a future version of MATLAB% handles    structure with handles and user data (see GUIDATA)close(gcf);% --------------------------------------------------------------------function menu_compareset_Callback(hObject, eventdata, handles)% hObject    handle to menu_compareset (see GCBO)% eventdata  reserved - to be defined in a future version of MATLAB% handles    structure with handles and user data (see GUIDATA)% --------------------------------------------------------------------function menu_open_Callback(hObject, eventdata, handles)% hObject    handle to menu_open (see GCBO)% eventdata  reserved - to be defined in a future version of MATLAB% handles    structure with handles and user data (see GUIDATA)% Flag to check if the GUI is picking up data from file or from workspace

% Update handles structure

[filename, pathname] = uigetfile( ...
    {'*.cmp', 'All CMP-Files (*.cmp)'},...
    'Open Prony Comparator File');
% If "Cancel" is selected then return
if isequal([filename,pathname],[0,0])
    handles.flag_menuopen=0;
    guidata(hObject, handles);
    return
    % Otherwise load the file.
else
    cd (pathname);
    CmpData=load(eval('filename'),'-mat');
    assignin('base','CmpData',CmpData);
    handles.flag_menuopen=1;
    guidata(hObject, handles);
end
setstatus(gcbf,'Select Sessions from Listbox and Choose an Option from "Plots"');
set(handles.ppmenu_plots,'enable','on');
maxNum=evalin('base','CmpData.maxNum');
SaveData=evalin('base','CmpData.SaveData');
[promptstr]=  prepare_comparelist(maxNum,SaveData);
set(handles.lstbx_session,'String',promptstr)% --------------------------------------------------------------------function menu_close_Callback(hObject, eventdata, handles)% hObject    handle to menu_close (see GCBO)% eventdata  reserved - to be defined in a future version of MATLAB% handles    structure with handles and user data (see GUIDATA)close(gcf);% --------------------------------------------------------------------function menu_view_Callback(hObject, eventdata, handles)% hObject    handle to menu_view (see GCBO)% eventdata  reserved - to be defined in a future version of MATLAB% handles    structure with handles and user data (see GUIDATA)% --------------------------------------------------------------------function menu_load_Callback(hObject, eventdata, handles)% hObject    handle to menu_load (see GCBO)% eventdata  reserved - to be defined in a future version of MATLAB% handles    structure with handles and user data (see GUIDATA)
setstatus(gcbf,'Select Sessions from Listbox and Choose an Option from "Plots"');
set(handles.ppmenu_plots,'enable','on');
maxNum=evalin('base','maxNum');
SaveData=evalin('base','SaveData');
[promptstr]=  prepare_comparelist(maxNum,SaveData);
set(handles.lstbx_session,'ForegroundColor','black');
set(handles.lstbx_session,'String',promptstr);

% --- Executes during object creation, after setting all properties.function ppmenu_plots_CreateFcn(hObject, eventdata, handles)% hObject    handle to ppmenu_plots (see GCBO)% eventdata  reserved - to be defined in a future version of MATLAB% handles    empty - handles not created until after all CreateFcns called% Hint: popupmenu controls usually have a white background on Windows.%       See ISPC and COMPUTER.if ispc    set(hObject,'BackgroundColor','white');else    set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor'));end% --- Executes on selection change in ppmenu_plots.function ppmenu_plots_Callback(hObject, eventdata, handles)% hObject    handle to ppmenu_plots (see GCBO)% eventdata  reserved - to be defined in a future version of MATLAB% handles    structure with handles and user data (see GUIDATA)% Hints: contents = get(hObject,'String') returns ppmenu_plots contents as cell array%        contents{get(hObject,'Value')} returns selected item from ppmenu_plots
% --- Executes on button press in push_compareplots.
function push_compareplots_Callback(hObject, eventdata, handles)
% hObject    handle to push_compareplots (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

    %evalin('base','clc');
    if(handles.flag_menuopen==1)
        SaveData=evalin('base','CmpData.SaveData');
    else
        SaveData=evalin('base','SaveData');
    end
    Index_Selected=get(handles.lstbx_session,'Value');
    handles.Index_Selected=Index_Selected;
    guidata(hObject, handles);
    % Get the column size of Index_Selected
    Run_limit=size(Index_Selected,2);
    % if nothing is selected
    if (isempty(Index_Selected))
        errordlg('Please Select the Sessions to be Compared','No Session Selected !!');
        whitebg('white'); 
    end
    
    for index=1:Run_limit
        Run=Index_Selected(index);
        mode_sizes(index)=SaveData(Run).modes;
    end
        max_modes=max(mode_sizes);
        
    
val=get(handles.ppmenu_plots,'Value');
switch val
    case 1 % poles plot
        
        for index=1:Run_limit
            Run=Index_Selected(index);
            [dp_x,dp_y,ai,a_list,SUB_IND,p,energy]= preparecomparedata(SaveData,Run);
            [inummz,idennz]=residuez(a_list(SUB_IND),p(SUB_IND),1);
            [z,p,k] = tf2zp(inummz,idennz);
            p_size=size(p,1); % Get the column size of p
            diff_size=max_modes-p_size;
            if (p_size<max_modes)
                dummy(1:diff_size,1)=-9;   % Create dummy poles of -9
                dummy_rowsize=size(dummy,1);
                if(dummy_rowsize>diff_size)
                    dummy_new=dummy(1:diff_size);
                else
                    dummy_new=dummy(1:diff_size);
                end
                p1=vertcat(p,dummy_new);
            else
                p1=p;
            end
            idennz_mat(:,index)=p1;
        end
        zplane(-9,idennz_mat,{'w^','x'});  % dummy zero at -9
        axis([-1.1 1.1 -1.1 1.1]);xlabel('Real Part'), ylabel('Imaginary Part');
        title('Poles Plot');
        legend_poles(handles,Index_Selected,Run_limit);
        handles.idennz_mat=idennz_mat;
        % Update handles structure
        guidata(hObject, handles);
        
        
    case 2 % Squared Error
        S=['m:';'b-';'y:';'r-';'g:'];
        row_S=size(S,1);
        
        for index=1:Run_limit
            colorindex=mod(index,row_S);% limit colorindex to the row size of S
            % avoid zero value of colorindex, it will happen when row_
            if (colorindex==0)
                colorindex=row_S;
            end
            Run=Index_Selected(index);
            [dp_x,dp_y,ai,a_list,SUB_IND,p,energy]= preparecomparedata(SaveData,Run);
            seerror=(dp_y(:)-ai(:)).^2;
            axes(handles.axes_pole);
            plot(dp_x,seerror,S(colorindex));
            hold on;
        end
         xlabel('Time'), ylabel('Squared Error')
         title('Squared Error Plot');
         legend_seerror(handles,Index_Selected,Run_limit);
         hold off;
            
     case 3 % Residue Plot
         S=['m:';'b-';'y:';'r-';'g:'];
         row_S=size(S,1);
         for index=1:Run_limit
            colorindex=mod(index,row_S);% limit colorindex to the row size of S
            % avoid zero value of colorindex, it will happen when row_
            if (colorindex==0)
                colorindex=row_S;
            end
            Run=Index_Selected(index);
            [dp_x,dp_y,ai,a_list,SUB_IND,p,energy]= preparecomparedata(SaveData,Run);
            SUB_N=size(SUB_IND,2);
            aa_list=abs(2*a_list(SUB_IND'));
            axes(handles.axes_pole);