Mercurial > hg > camir-aes2014
annotate toolboxes/FullBNT-1.0.7/bnt/potentials/@cgpot/set_domain_pot.m @ 0:e9a9cd732c1e tip
first hg version after svn
| author | wolffd |
|---|---|
| date | Tue, 10 Feb 2015 15:05:51 +0000 |
| parents | |
| children |
| rev | line source |
|---|---|
| wolffd@0 | 1 function pot = set_domain_pot(pot, domain) |
| wolffd@0 | 2 % SET_DOMAIN_POT Change the domain of a potential (cgpot) |
| wolffd@0 | 3 % pot = set_domain_pot(pot, domain) |
| wolffd@0 | 4 |
| wolffd@0 | 5 delta = domain(1) - pot.domain(1); |
| wolffd@0 | 6 assert(all(domain == pot.domain + delta)); |
| wolffd@0 | 7 pot.domain = pot.domain + delta; |
| wolffd@0 | 8 pot.ddom = pot.ddom + delta; |
| wolffd@0 | 9 pot.cdom = pot.cdom + delta; |
| wolffd@0 | 10 cdomain = pot.cdom; |
| wolffd@0 | 11 n = prod(pot.dsizes); |
| wolffd@0 | 12 if(pot.subtype == 'm') |
| wolffd@0 | 13 for i = 1: n |
| wolffd@0 | 14 pot.mom{i} = set_domain_pot(pot.mom{i}, cdomain); |
| wolffd@0 | 15 end |
| wolffd@0 | 16 end |
| wolffd@0 | 17 if(pot.subtype == 'c') |
| wolffd@0 | 18 for i = 1: n |
| wolffd@0 | 19 pot.can{i} = set_domain_pot(pot.can{i}, cdomain); |
| wolffd@0 | 20 end |
| wolffd@0 | 21 end |