blob: 073646482742e4c069a80b28caa8307ef43c0547 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
|
function pot = set_domain_pot(pot, domain)
% SET_DOMAIN_POT Change the domain of a potential (cgpot)
% pot = set_domain_pot(pot, domain)
delta = domain(1) - pot.domain(1);
assert(all(domain == pot.domain + delta));
pot.domain = pot.domain + delta;
pot.ddom = pot.ddom + delta;
pot.cdom = pot.cdom + delta;
cdomain = pot.cdom;
n = prod(pot.dsizes);
if(pot.subtype == 'm')
for i = 1: n
pot.mom{i} = set_domain_pot(pot.mom{i}, cdomain);
end
end
if(pot.subtype == 'c')
for i = 1: n
pot.can{i} = set_domain_pot(pot.can{i}, cdomain);
end
end
|
