annotate toolboxes/FullBNT-1.0.7/bnt/potentials/@cgpot/divide_by_pot.m @ 0:e9a9cd732c1e tip

first hg version after svn
author wolffd
date Tue, 10 Feb 2015 15:05:51 +0000
parents
children
rev   line source
wolffd@0 1 function bigpot = divide_by_pot(bigpot, smallpot)
wolffd@0 2 % DIVIDE_BY_POT bigpot /= smallpot for cgpot
wolffd@0 3 % bigpot = divide_by_pot(bigpot, smallpot)
wolffd@0 4 %
wolffd@0 5 % smallpot's domain must be a subset of bigpot's domain.
wolffd@0 6
wolffd@0 7 bigpot = cg_mom_to_can(bigpot);
wolffd@0 8 smallpot = cg_mom_to_can(smallpot);
wolffd@0 9
wolffd@0 10 mask = find_equiv_posns(smallpot.ddom, bigpot.ddom);
wolffd@0 11 for i=1:bigpot.dsize
wolffd@0 12 if isempty(smallpot.ddom)
wolffd@0 13 src = 1;
wolffd@0 14 else
wolffd@0 15 sub = ind2subv(bigpot.dsizes, i);
wolffd@0 16 src = subv2ind(smallpot.dsizes, sub(mask));
wolffd@0 17 end
wolffd@0 18 bigpot.can{i} = divide_by_pot(bigpot.can{i}, smallpot.can{src});
wolffd@0 19 end