function [state,options,optchanged] = gaoutfun(options,state,flag)
persistent state_record
if isempty(state_record)
state_record = struct('Population', {}, 'Best', {}, 'Score', {});
end
if nargin == 0
state = state_record;
options = [];
optchanged = [];
else
state_record(end+1) = struct('Population', state.Population, 'Best', state.Best', 'Score', state.Score);
optchanged = false;
end
Before re-using this function, call
clear gaoutfun
Let the function be called as usual. Do not expect any changes outside the function when it is called.
After the ga run is finished, call
record = gaoutfun();
with no inputs. A structure of the state information will be returned.
However, this is likely to be a bit inefficient because the struct is expanded each iteration. I would suggest that when isempty() succeeds, that options be examined to find the maximum generations, and that state_record be initialized that large. Use state.Generation as the index for the output instead of end+1
