annotate toolboxes/FullBNT-1.0.7/bnt/potentials/@upot/divide_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 Tbig = divide_by_pot(Tbig, Tsmall)
|
|
Daniel@0
|
2 % DIVIDE_BY_POT Tbig /= Tsmall
|
|
Daniel@0
|
3 % Tbig = divide_by_pot(Tbig, Tsmall)
|
|
Daniel@0
|
4 %
|
|
Daniel@0
|
5 % Tsmall's domain must be a subset of Tbig's domain.
|
|
Daniel@0
|
6
|
|
Daniel@0
|
7 smallp = extend_domain_table(Tsmall.p, Tsmall.domain, Tsmall.sizes, Tbig.domain, Tbig.sizes);
|
|
Daniel@0
|
8 smallp = smallp + (smallp==0);
|
|
Daniel@0
|
9 Tbig.p = Tbig.p ./ smallp;
|
|
Daniel@0
|
10
|
|
Daniel@0
|
11 smallu = extend_domain_table(Tsmall.u, Tsmall.domain, Tsmall.sizes, Tbig.domain, Tbig.sizes);
|
|
Daniel@0
|
12 Tbig.u = Tbig.u - smallu;
|
|
Daniel@0
|
13
|
