annotate toolboxes/FullBNT-1.0.7/bnt/potentials/@cpot/multiply_by_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 bigpot = multiply_by_pot(bigpot, smallpot, varargin)
|
|
Daniel@0
|
2 % MULTIPLY_BY_POT bigpot *= smallpot for cpot
|
|
Daniel@0
|
3 % bigpot = multiply_by_pot(bigpot, smallpot)
|
|
Daniel@0
|
4 %
|
|
Daniel@0
|
5 % smallpot's domain must be a subset of bigpot's domain.
|
|
Daniel@0
|
6
|
|
Daniel@0
|
7 bigpot.g = bigpot.g + smallpot.g;
|
|
Daniel@0
|
8 if sum(smallpot.sizes) > 0
|
|
Daniel@0
|
9 mask = find_equiv_posns(smallpot.domain, bigpot.domain);
|
|
Daniel@0
|
10 u = block(mask, bigpot.sizes);
|
|
Daniel@0
|
11 bigpot.h(u) = bigpot.h(u) + smallpot.h;
|
|
Daniel@0
|
12 bigpot.K(u, u) = bigpot.K(u, u) + smallpot.K;
|
|
Daniel@0
|
13 end
|
