sourcecodes/bnt-master/BNT/potentials/@cgpot/Old/set_domain_pot.m


1
2
3
4
5
6
7
8
9

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;