Daniel@0: function T = multiply_pots(T1, T2)
Daniel@0: % MULTIPLY_POTS Multiply a pair of dpots together pointwise.
Daniel@0: % T = multiply_pots(pots)
Daniel@0: 
Daniel@0: dom = myunion(T1.domain, T2.domain);
Daniel@0: %ns = sparse(1, max(dom)); % causes problems in myreshape on NT
Daniel@0: ns = zeros(1, max(dom));
Daniel@0: ns(T1.domain) = T1.sizes;
Daniel@0: ns(T2.domain) = T2.sizes;
Daniel@0: T = dpot(dom, ns(dom));
Daniel@0: T = multiply_by_pot(T, T1);
Daniel@0: T = multiply_by_pot(T, T2);