Daniel@0: function Tbig = multiply_by_pot(Tbig, Tsmall)
Daniel@0: % MULTIPLY_BY_POT Tbig *= Tsmall
Daniel@0: % Tbig = multiply_by_pot(Tbig, Tsmall)
Daniel@0: %
Daniel@0: % Tsmall's domain must be a subset of Tbig's domain.
Daniel@0: 
Daniel@0: %process sparse dpot, we do not consider only one of the two pots is sparse
Daniel@0: if issparse(Tbig.T) & issparse(Tsmall.T)
Daniel@0:    Tbig.T = mult_by_sparse_table(Tbig.T, Tbig.domain, Tbig.sizes, Tsmall.T, Tsmall.domain, Tsmall.sizes);
Daniel@0: else 
Daniel@0:    Tbig.T = mult_by_table(Tbig.T, Tbig.domain, Tbig.sizes, Tsmall.T, Tsmall.domain, Tsmall.sizes);
Daniel@0: end
Daniel@0: