]> Creatis software - CreaPhase.git/blob - octave_packages/ga-0.10.0/rastriginsfcn.m
Add a useful package (from Source forge) for octave
[CreaPhase.git] / octave_packages / ga-0.10.0 / rastriginsfcn.m
1 ## Copyright (C) 2008, 2009, 2010, 2012 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 ## -*- texinfo -*-
17 ## @deftypefn{Function File} {@var{y} =} rastriginsfcn (@var{x})
18 ## Rastrigin's function.
19 ## @end deftypefn
20
21 ## Author: Luca Favatella <slackydeb@gmail.com>
22 ## Version: 2.0
23
24 function retval = rastriginsfcn (x)
25   if ((nargout != 1) ||
26       (nargin != 1) || (columns (x) != 2))
27     print_usage ();
28   else
29     x1 = x(:, 1);
30     x2 = x(:, 2);
31     retval = 20 + (x1 .** 2) + (x2 .** 2) - 10 .* (cos (2 .* pi .* x1) +
32                                                    cos (2 .* pi .* x2));
33   endif
34 endfunction
35
36
37 ## number of input arguments
38 %!error y = rastriginsfcn ()
39 %!error y = rastriginsfcn ([0, 0], "other argument")
40
41 ## number of output arguments
42 %!error [y1, y2] = rastriginsfcn ([0, 0])
43
44 ## type of arguments
45 %!error y = rastriginsfcn ([0; 0])
46 %!error y = rastriginsfcn (zeros (2, 3)) # TODO: document size of x
47
48 %!assert (rastriginsfcn ([0, 0]), 0)
49 %!assert (rastriginsfcn ([0, 0; 0, 0]), [0; 0])
50 %!assert (rastriginsfcn (zeros (3, 2)), [0; 0; 0])