wolffd@0: function T = multiply_pots(T1, T2)
wolffd@0: % MULTIPLY_POTS Multiply a pair of dpots together pointwise (cgpot)
wolffd@0: % T = multiply_pots(pots)
wolffd@0: 
wolffd@0: ddom = myunion(T1.ddom, T2.ddom);
wolffd@0: cdom = myunion(T1.cdom, T2.cdom);
wolffd@0: dom = myunion(ddom, cdom);
wolffd@0: ns = zeros(1, max(dom));
wolffd@0: ns(T1.ddom) = T1.dsizes;
wolffd@0: ns(T2.ddom) = T2.dsizes;
wolffd@0: ns(T1.cdom) = T1.csizes;
wolffd@0: ns(T2.cdom) = T2.csizes;
wolffd@0: 
wolffd@0: T = cgpot(ddom, cdom, ns);
wolffd@0: T = multiply_by_pot(T, T1);
wolffd@0: T = multiply_by_pot(T, T2);