Mercurial > hg > mauch-mirex-2010
view _FullBNT/BNT/potentials/Tables/divide_by_table.m @ 9:4ea6619cb3f5 tip
removed log files
author | matthiasm |
---|---|
date | Fri, 11 Apr 2014 15:55:11 +0100 |
parents | b5b38998ef3b |
children |
line wrap: on
line source
function bigT = divide_by_table(bigT, bigdom, bigsz, smallT, smalldom, smallsz) % DIVIDE_BY_TABLE % bigT = divide_by_table(bigT, bigdom, bigsz, smallT, smalldom, smallsz) % Ts = extend_domain_table(smallT, smalldom, smallsz, bigdom, bigsz); % Replace 0s by 1s before dividing. This is valid, Ts(i)=0 iff Tbig(i)=0. Ts = Ts + (Ts==0); %Tbig.T(:) = Tbig.T(:) ./ Ts(:); bigT(:) = bigT(:) ./ Ts(:);