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