annotate toolboxes/FullBNT-1.0.7/bnt/potentials/@upot/marginalize_pot.m @ 0:cc4b1211e677
tip
initial commit to HG from
Changeset:
646 (e263d8a21543) added further path and more save "camirversion.m"
| author |
Daniel Wolff |
| date |
Fri, 19 Aug 2016 13:07:06 +0200 |
| parents |
|
| children |
|
| rev |
line source |
|
Daniel@0
|
1 function smallpot = marginalize_pot(bigpot, onto, maximize)
|
|
Daniel@0
|
2 % MARGINALIZE_POT Marginalize a upot onto a smaller domain.
|
|
Daniel@0
|
3 % smallpot = marginalize_pot(bigpot, onto, maximize)
|
|
Daniel@0
|
4 %
|
|
Daniel@0
|
5 % The maximize argument is ignored
|
|
Daniel@0
|
6
|
|
Daniel@0
|
7 numer = marg_table(bigpot.p .* bigpot.u, bigpot.domain, bigpot.sizes, onto);
|
|
Daniel@0
|
8 denom = marg_table(bigpot.p, bigpot.domain, bigpot.sizes, onto);
|
|
Daniel@0
|
9
|
|
Daniel@0
|
10 p = denom;
|
|
Daniel@0
|
11 % replace 0s by 1s before dividing. This is valid since demon(i) = 0 => numer(i) = 0
|
|
Daniel@0
|
12 denom = denom + (denom == 0);
|
|
Daniel@0
|
13 u = numer ./ denom;
|
|
Daniel@0
|
14
|
|
Daniel@0
|
15 ns = zeros(1, max(bigpot.domain));
|
|
Daniel@0
|
16 ns(bigpot.domain) = bigpot.sizes;
|
|
Daniel@0
|
17
|
|
Daniel@0
|
18 smallpot = upot(onto, ns(onto), p, u);
|
