annotate toolboxes/FullBNT-1.0.7/bnt/potentials/@cgpot/set_domain_pot.m @ 0:cc4b1211e677 tip

initial commit to HG from Changeset: 646 (e263d8a21543) added further path and more save "camirversion.m"
author Daniel Wolff
date Fri, 19 Aug 2016 13:07:06 +0200
parents
children
rev   line source
Daniel@0 1 function pot = set_domain_pot(pot, domain)
Daniel@0 2 % SET_DOMAIN_POT Change the domain of a potential (cgpot)
Daniel@0 3 % pot = set_domain_pot(pot, domain)
Daniel@0 4
Daniel@0 5 delta = domain(1) - pot.domain(1);
Daniel@0 6 assert(all(domain == pot.domain + delta));
Daniel@0 7 pot.domain = pot.domain + delta;
Daniel@0 8 pot.ddom = pot.ddom + delta;
Daniel@0 9 pot.cdom = pot.cdom + delta;
Daniel@0 10 cdomain = pot.cdom;
Daniel@0 11 n = prod(pot.dsizes);
Daniel@0 12 if(pot.subtype == 'm')
Daniel@0 13 for i = 1: n
Daniel@0 14 pot.mom{i} = set_domain_pot(pot.mom{i}, cdomain);
Daniel@0 15 end
Daniel@0 16 end
Daniel@0 17 if(pot.subtype == 'c')
Daniel@0 18 for i = 1: n
Daniel@0 19 pot.can{i} = set_domain_pot(pot.can{i}, cdomain);
Daniel@0 20 end
Daniel@0 21 end