comparison toolboxes/FullBNT-1.0.7/bnt/potentials/@cpot/multiply_by_pot.m @ 0:e9a9cd732c1e tip

first hg version after svn
author wolffd
date Tue, 10 Feb 2015 15:05:51 +0000
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:e9a9cd732c1e
1 function bigpot = multiply_by_pot(bigpot, smallpot, varargin)
2 % MULTIPLY_BY_POT bigpot *= smallpot for cpot
3 % bigpot = multiply_by_pot(bigpot, smallpot)
4 %
5 % smallpot's domain must be a subset of bigpot's domain.
6
7 bigpot.g = bigpot.g + smallpot.g;
8 if sum(smallpot.sizes) > 0
9 mask = find_equiv_posns(smallpot.domain, bigpot.domain);
10 u = block(mask, bigpot.sizes);
11 bigpot.h(u) = bigpot.h(u) + smallpot.h;
12 bigpot.K(u, u) = bigpot.K(u, u) + smallpot.K;
13 end