blob: 1e2bb2b9d2f9be120839c8fd199a08b7a5096cfd (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
|
function bigpot = divide_by_pot(bigpot, smallpot)
% DIVIDE_BY_POT bigpot /= smallpot for cpot
% bigpot = divide_by_pot(bigpot, smallpot)
%
% smallpot's domain must be a subset of bigpot's domain.
bigpot.g = bigpot.g - smallpot.g;
if sum(smallpot.sizes) > 0
mask = find_equiv_posns(smallpot.domain, bigpot.domain);
u = block(mask, bigpot.sizes);
bigpot.h(u) = bigpot.h(u) - smallpot.h;
bigpot.K(u, u) = bigpot.K(u, u) - smallpot.K;
end
|
