function varargout = Urun(varargin)
% URUN MATLAB code for Urun.fig
%      URUN, by itself, creates a new URUN or raises the existing
%      singleton*.
%
%      H = URUN returns the handle to a new URUN or the handle to
%      the existing singleton*.
%
%      URUN('CALLBACK',hObject,eventData,handles,...) calls the local
%      function named CALLBACK in URUN.M with the given input arguments.
%
%      URUN('Property','Value',...) creates a new URUN or raises the
%      existing singleton*.  Starting from the left, property value pairs are
%      applied to the GUI before Urun_OpeningFcn gets called.  An
%      unrecognized property name or invalid value makes property application
%      stop.  All inputs are passed to Urun_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 Urun

% Last Modified by GUIDE v2.5 10-Jan-2014 14:44:15

% Begin initialization code - DO NOT EDIT
warning('off','all');
gui_Singleton = 1;
gui_State = struct('gui_Name',       mfilename, ...
    'gui_Singleton',  gui_Singleton, ...
    'gui_OpeningFcn', @Urun_OpeningFcn, ...
    'gui_OutputFcn',  @Urun_OutputFcn, ...
    'gui_LayoutFcn',  [] , ...
    'gui_Callback',   []);
if nargin && ischar(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 Urun is made visible.
function Urun_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 Urun (see VARARGIN)

% Choose default command line output for Urun
handles.output = hObject;
handles.jobs={};
handles.run=0;
% Update handles structure
guidata(hObject, handles);

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


% --- Outputs from this function are returned to the command line.
function varargout = Urun_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;


% --------------------------------------------------------------------
function open_job_ClickedCallback(hObject, eventdata, handles)
% hObject    handle to open_job (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
[fils,paths]=uigetfile({'*.dat','Ufreckles input files';'*.ufr','Ufreckles input ascii files'},'Open data sets');
%if ~iscell(fils)
%   fils={fils};
%end
handles.jobs{1}=fils;
handles.jobs{2}=paths;
%  nbjobs=size(handles.jobs,1);
%  for id=1:length(fils)
%      handles.jobs{nbjobs+id,1}=fils{id};
%      handles.jobs{nbjobs+id,2}=paths;
%  end
guidata(hObject, handles);
set(handles.job_list,'String',handles.jobs{1});
set(handles.runing_text,'String',sprintf('No job running...'))


% --- Executes on selection change in job_list.

% --- Executes during object creation, after setting all properties.
function job_list_CreateFcn(hObject, eventdata, handles)
% hObject    handle to job_list (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 && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end
function run_live(hrun,param,model)
global phix phiy
set(hrun,'String','Stop');
run=1;
nmod=0;
clear functions
param.onflight=1;
nscale=1;
model.nscale=1;
LoadParameters(param);
LoadParameters(model,nmod);
ReferenceImage(nmod);
LoadMask(nmod);
filres=param.result_file;
iscale=nscale;
switch model.basis
    case 'fem'
LoadMeshes(nmod);
LoadMat(nmod);
load(fullfile('TMP','0_mesh_0'),'Nnodes','Nelems','xo','yo','conn','elt','ng','rflag','rint');
tips=zeros(size(model.zone,2),2);
cracks=[];
Ks=[];
U=zeros(2*prod(Nnodes),1);
param.ulive=1;
save(filres,'U','Ks','tips','cracks','Nnodes','Nelems','xo','yo','param','model','nmod','conn','elt','rint','ng','rflag','-v7.3');
CreateBasisFunction(iscale,nmod);
ComputeGradFPhi(iscale,nmod);
CreateGradBasisFunction(iscale,nmod);
AssembleCorrelationOperator(iscale,nmod);
    case 'uni'
        gages=model.zone(2,:);
        assert(length(gages)==1)
param.ulive=1;ig=1;
    gage=gages{ig};
    roi=param.roi;
    sizeim=[roi(2)-roi(1),roi(4)-roi(3)]+1;
    xp=gage(1:4,1)+roi(1)-1;
    yp=gage(1:4,2)+roi(3)-1;
    xo=round(xp);
    yo=round(yp);
elt=4;
zo=1;
conn=1:4;
Nnodes=[length(xo),1,1];
Nelems=[length(elt),1,1];
        Smesh=[max(xo)-min(xo),max(yo)-min(yo)];
        Vmesh=[max(xo)-min(xo),max(yo)-min(yo),max(zo)-min(zo)];
ns=ones(1,3);
    rflag=1;
    rint=false;
    ng=0;
    selected=[];
      U1=zeros(6,1);
      U=zeros(2*prod(Nnodes),1);
      Up=cell(length(gages),1);
      for ig=1:length(gages)
          Up{ig}=zeros(9,1);
      end
    [Yi, Xi]=meshgrid(1:(sizeim(2)),1:(sizeim(1)));
    mask=inpolygon(Xi,Yi,xp([1:length(xp),1])-roi(1)+1,yp([1:length(xp),1])-roi(3)+1);
              save(filres,'U','Up','Nnodes','Nelems','xo','yo','param','model','nmod','conn','elt','rint','ng','rflag','-v7.3');
        mask=diag(sparse(mask(:)));
        save(fullfile('TMP',sprintf('%d_mask_%d',nmod,iscale-1)),'mask','-append');
CreateBasisFunction(iscale,nmod);
ComputeGradFPhi(iscale,nmod);
AssembleCorrelationOperator(iscale,nmod);
                    save(fullfile('TMP',sprintf('%d_mesh_%d',nmod,iscale-1)),'rflag','rint','xo','yo','zo','elt','conn','Nnodes','Nelems','Smesh','Vmesh','selected','ng','ns');

        try
            delete([strrep(filres,'.res',''),'-gage.csv']);
        end
        dlmwrite([strrep(filres,'.res',''),'-gage.csv'],zeros(1,4));
       
end
        switch model.basis
            case 'uni'
               indp=xo-roi(1)+1+(yo-roi(3)+1-1)*sizeim(1);
    phixo=phix(indp,:);
    phiyo=phiy(indp,:);
     
                
        end
[~,fil,ext]=fileparts(param.reference_image);
hmaster=findobj('Tag',param.umaster);
for ih=1:numel(hmaster)
    if ~isempty(strfind(hmaster(ih).Name,'UFreckles'))
        hmaster=hmaster(ih);
        break;
    end
end
handles=guidata(hmaster);
Umaster('reset_param',handles);
Umaster('loading_mat_file',hmaster,handles,param.result_file,[],1);
nim=1;
imdefs={param.deformed_image,param.deformed_image};
load(fullfile('TMP','params'),'param');
param.deformed_image=imdefs;
while run
    
    pause(0.1);
    listing=dir(sprintf('*%s',ext'));
    go=size(listing,1)>(nim);
    if go
        
        nim=nim+1;
        listing=listing(end);
        param.deformed_image{nim}=listing.name;
        save(fullfile('TMP','params'),'param');
        switch model.basis
            case 'fem'
        U=[U,zeros(2*prod(Nnodes),1)];
        [U]=SolvePreview(U,nscale,nmod,nim);
            case 'uni'
        U1=[U1,zeros(6,1)];
        [U1]=SolvePreview(U1,nscale,nmod,nim);
            Ux=phixo*U1;
            Uy=phiyo*U1;  
            exx=U1(4,:)/mean(sizeim);
                eyy=U1(5,:)/mean(sizeim);
                exy=U1(6,:)/mean(sizeim);

U=[Ux;Uy];
Up{1}=[U1;exx;eyy;exy];
        
        dlmwrite([strrep(filres,'.res',''),'-gage.csv'],[nim,exx(end),eyy(end),exy(end)],'-append');
        
        end
        handles=guidata(hmaster);
        handles.uvisu=U;
        handles.param.deformed_image=param.deformed_image;
        handles.animation.nbstep=size(U,2);
        handles.animation.frames=1:size(U,2);
        handles.animation.iim=handles.animation.nbstep;
        handles=Umaster('display_frame',handles);
        guidata(hmaster,handles);
        save(filres,'U','param','-append');
        switch model.basis
            case 'uni'
        save(filres,'Up','-append');  
        handles.uni_model.Up=Up;
        guidata(hmaster,handles);
        Umaster('plot_gage_data',handles,handles.uni_model.zone{2,1},1);
        end
    end
    rhandles=guidata(hrun);
    run=rhandles.run;
end
delete(fullfile('VTK',['camr*','-error.vtk']));
set(hrun,'String','Run');





% --- Executes on button press in run_button.
function run_button_Callback(hObject, eventdata, handles)
% hObject    handle to run_button (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
%for id=1:size(handles.jobs,1)
if handles.run==0
    handles.run=1;
    guidata(handles.figure1,handles);
    filname=handles.jobs{1};
    path=handles.jobs{2};
    cd(path);
    [pp,jobname,ext]=fileparts(filname);
    set(handles.runing_text,'String',sprintf('Running %s...',jobname))
    pause(0.1);
    if strcmp(ext,'.ufr')
        [param,model]=readINPFile(filname);
    else
        load(filname,'-mat','param','model');
    end
    if ~isfield(param,'da'), param.da=0;end
    if ~isfield(param,'ulive'), param.ulive=0;end
    if ~isfield(param,'thermo'), param.thermo=0;end
    if ~isfield(param,'psample'), param.psample=1;end
    if param.thermo==1
        if iscell(param.ir_calibration_data)
            pini=cd;
            
            T=UTcalib(param.ir_calibration_data{1},param.ir_calibration_data{2});
            cd(pini);
            param.ir_calibration_data=T;
            save(filname,'param','-append')
        end
        
    end
    switch model.basis
        case 'fem'
            if ~isfield(model,'phantom_nodes'), model.phantom_nodes=0;end
            switch param.analysis
                case 'correlation'
                    if param.ulive==1
                        run_live(handles.run_button,param,model);
                    else
                        if isfield(param,'detect')
                            if param.detect
                                run_crack_detect_job(param,model);
                            else
                                if isfield(model,'degree')
                                    run_nurbs_job(param,model);
                                else
                                    run_fem_job(param,model);
                                    
                                end
                            end
                        else
                            if isfield(model,'degree')
                                run_nurbs_job(param,model);
                            else
                                run_fem_job(param,model);
                            end
                        end
                        handles.run=0;
                        guidata(handles.figure1,handles);
                    end
                    
                case 'mechanics'
                    if param.da>0
                        run_crack_propa_job(param,model);
                    else
                        if model.phantom_nodes==0
                            run_fea_job(param,model);
                        else
                            clear run_fdfea_job
                            run_fdfea_job(param,model);
                        end
                    end
                    handles.run=0;
                    guidata(handles.figure1,handles);
            end
        case 'uni'
                    if param.ulive==1
                        run_live(handles.run_button,param,model);
                    else
            run_uni_job(param,model);
            handles.run=0;
            guidata(handles.figure1,handles);
                    end
        case 'vic'
            run_vic_job(param,model);
            handles.run=0;
            guidata(handles.figure1,handles);
        case {'beam','beam-nurbs'}
            run_beam_job(param,model);
            handles.run=0;
            guidata(handles.figure1,handles);
            
    end
    set(handles.runing_text,'String',sprintf('End running %s !',jobname))
else
    handles.run=0;
    guidata(handles.figure1,handles);
end
%set(handles.job_list,'String','');


%end