octave_packages/ga-0.10.0/__ga_problem_update_state_at_each_generation__.m

   1 ## Copyright (C) 2008, 2010 Luca Favatella <slackydeb@gmail.com>
   2 ##
   3 ## This program is free software; you can redistribute it and/or modify
   4 ## it under the terms of the GNU General Public License as published by
   5 ## the Free Software Foundation; either version 2 of the License, or
   6 ## (at your option) any later version.
   7 ##
   8 ## This program is distributed in the hope that it will be useful,
   9 ## but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  11 ## GNU General Public License for more details.
  12 ##
  13 ## You should have received a copy of the GNU General Public License
  14 ## along with this program; If not, see <http://www.gnu.org/licenses/>.
  15
  16 ## Author: Luca Favatella <slackydeb@gmail.com>
  17 ## Version: 1.3.1
  18
  19 function state = \
  20       __ga_problem_update_state_at_each_generation__ (state, problem,
  21                                                       private_state)
  22   if ((state.Generation > 0) || isempty (problem.options.InitialScores))
  23     state.Score(1:problem.options.PopulationSize, 1) = \
  24         __ga_scores__ (problem, state.Population);
  25   else ## (Generation == 0) && (! isempty (InitialScores))
  26     nrIS = rows (problem.options.InitialScores);
  27     #assert (rows (problem.options.InitialPopulation) <= problem.options.PopulationSize); ## DEBUG
  28     if (nrIS <= rows (problem.options.InitialPopulation))
  29       missing_rows = (nrIS+1):problem.options.PopulationSize;
  30       state.Score(1:problem.options.PopulationSize, 1) = \
  31           [problem.options.InitialScores(:, 1);
  32            (__ga_scores__ (problem, state.Population(missing_rows, :)))
  33            ];
  34     else
  35       error ("rows (InitialScores) > rows (InitialPopulation)");
  36     endif
  37   endif
  38   state.Expectation(1, 1:problem.options.PopulationSize) = \
  39       problem.options.FitnessScalingFcn (state.Score, private_state.nParents);
  40   state.Best(state.Generation + 1, 1) = min (state.Score);
  41 endfunction
  42
  43
  44 %!error
  45 %! state.Generation = 0;
  46 %! problem = struct ("fitnessfcn", @rastriginsfcn, "nvars", 2, "options", gaoptimset ("Generations", 10, "InitialScores", [0; 0; 0]));
  47 %! unused = 0;
  48 %! __ga_problem_update_state_at_each_generation__ (state, problem, unused);