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;
Daniel@0: cdomain = pot.cdom;
Daniel@0: n = prod(pot.dsizes);
Daniel@0: if(pot.subtype == 'm')
Daniel@0:     for i = 1: n
Daniel@0:         pot.mom{i} = set_domain_pot(pot.mom{i}, cdomain);
Daniel@0:     end
Daniel@0: end
Daniel@0: if(pot.subtype == 'c')
Daniel@0:     for i = 1: n
Daniel@0:         pot.can{i} = set_domain_pot(pot.can{i}, cdomain);
Daniel@0:     end
Daniel@0: end