wolffd@0: function pot = set_domain_pot(pot, domain)
wolffd@0: % SET_DOMAIN_POT Change the domain of a potential (cgpot)
wolffd@0: % pot = set_domain_pot(pot, domain)
wolffd@0: 
wolffd@0: delta = domain(1) - pot.domain(1);
wolffd@0: assert(all(domain == pot.domain + delta));
wolffd@0: pot.domain = pot.domain + delta;
wolffd@0: pot.ddom = pot.ddom + delta;
wolffd@0: pot.cdom = pot.cdom + delta;