www.pudn.com > ConstrainedEM.zip > st_multiply_pots.m

function T = st_multiply_pots(T1, T2)
% same as multiply_pots of dpot of BNT
% MULTIPLY_POTS Multiply a pair of dpots together pointwise.
% T = multiply_pots(pots)

dom = myunion(T1.domain, T2.domain);
%ns = sparse(1, max(dom)); % causes problems in myreshape on NT
ns = zeros(1, max(dom));
ns(T1.domain) = T1.sizes;
ns(T2.domain) = T2.sizes;
T.domain=dom;
T.sizes=ns(dom);
T .T=ones(T.sizes); 
T = st_multiply_by_pot(T, T1);
T = st_multiply_by_pot(T, T2);