1 ## Copyright (C) 2008 Luca Favatella <slackydeb@gmail.com>
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.
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.
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/>.
16 ## Author: Luca Favatella <slackydeb@gmail.com>
19 function parents = selectionstochunif (expectation, nParents, options)
20 nc_expectation = columns (expectation);
21 line(1, 1:nc_expectation) = cumsum (expectation(1, 1:nc_expectation));
22 max_step_size = line(1, nc_expectation);
23 step_size = max_step_size * rand ();
24 steps(1, 1:nParents) = rem (step_size * (1:nParents), max_step_size);
25 for index_steps = 1:nParents ## fix an entry of the steps (or parents) vector
26 #assert (steps(1, index_steps) < max_step_size); ## DEBUG
28 while (steps(1, index_steps) >= line(1, index_line))
29 #assert ((index_line >= 1) && (index_line < nc_expectation)); ## DEBUG
32 parents(1, index_steps) = index_line;
37 ## number of input arguments
40 ## number of output arguments