Daniel@0: function bigpot = divide_by_pot(bigpot, smallpot)
Daniel@0: % DIVIDE_BY_POT bigpot /= smallpot for cgpot
Daniel@0: % bigpot = divide_by_pot(bigpot, smallpot)
Daniel@0: %
Daniel@0: % smallpot's domain must be a subset of bigpot's domain.
Daniel@0: 
Daniel@0: bigpot = cg_mom_to_can(bigpot);
Daniel@0: smallpot = cg_mom_to_can(smallpot);
Daniel@0: 
Daniel@0: mask = find_equiv_posns(smallpot.ddom, bigpot.ddom);
Daniel@0: for i=1:bigpot.dsize
Daniel@0:   if isempty(smallpot.ddom)
Daniel@0:     src = 1;
Daniel@0:   else
Daniel@0:     sub = ind2subv(bigpot.dsizes, i);
Daniel@0:     src = subv2ind(smallpot.dsizes, sub(mask));
Daniel@0:   end
Daniel@0:   bigpot.can{i} = divide_by_pot(bigpot.can{i}, smallpot.can{src});
Daniel@0: end