load_opt_data.m

Go to the documentation of this file.

function opt_data = load_opt_data(filestr, model)
%model = load_opt_data(filestr, model)
%
%Function loads optimization data from file named filestr.
%
%04.03.2010 Markus Dihlmann


Pfad=[getenv('RBMATLABTEMP'), '/',filestr];

fid = fopen(Pfad,'r');

str=fgetl(fid);
while(~strcmp(str(1:8),'Printing'))
    str=fgetl(fid);
    while(length(str)<8)
        str=fgetl(fid);
    end
end

str=fgetl(fid);
str=fgetl(fid);
index=1;
opt_data=[];
while(str~=-1)
    k = strfind(str, sprintf('\t'));
    opt_data.output(index)=str2double(str(1:k(1)));
    for i=1:length(k)-1
        opt_data.parameter_sets(index,i)=str2double(str(k(i):k(i+1)));
    end
    
    index=index+1;
    str=fgetl(fid);
end

opt_data.max_output = max(opt_data.output);

%where to reach this maximum value:
ind=1;
Pmax=[];
ind=search_pos(opt_data.max_output, opt_data.output,ind);
while(ind>0)
    Pmax=[Pmax;opt_data.parameter_sets(ind,:)];
    ind=search_pos(opt_data.max_output, opt_data.output,ind+1);
end

opt_data.max_output_paramsets=Pmax;


fclose(fid);


end



function ind=search_pos(value,vec,start)
    found=0;
    ind=start;
    while ((~found)&&(ind<=length(vec)))
        if(vec(ind)==value)
            found=1;
        else
            ind = ind+1;
        end
    end
    
    if(ind>length(vec))
        ind=-1;
    end
        
end