blob: d3c6e7f912c23f146f210208197cec4c47c13179 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
|
function bigpot = divide_by_pot(bigpot, smallpot)
% DIVIDE_BY_POT bigpot /= smallpot for cgpot
% bigpot = divide_by_pot(bigpot, smallpot)
%
% smallpot's domain must be a subset of bigpot's domain.
bigpot = cg_mom_to_can(bigpot);
smallpot = cg_mom_to_can(smallpot);
mask = find_equiv_posns(smallpot.ddom, bigpot.ddom);
for i=1:bigpot.dsize
if isempty(smallpot.ddom)
src = 1;
else
sub = ind2subv(bigpot.dsizes, i);
src = subv2ind(smallpot.dsizes, sub(mask));
end
bigpot.can{i} = divide_by_pot(bigpot.can{i}, smallpot.can{src});
end
|
