function pot1 = multiply_by_gpot(pot1, pot2)
% MULTIPLY_BY_GPOT Compute pot1 *= pot2 for Gaussian potentials.
% pot1 = multiply_by_gpot(pot1, pot2)

if pot1.subtype == 'c' & pot2.subtype== 'c'
  [pot1.g, pot1.h, pot1.K] = multiply_by_canonical(pot1.g, pot1.h, pot1.K, pot2.g, pot2.h, pot2.K, ...
      pot1.domain, pot2.domain, pot1.sizes, pot2.sizes);
else
  error('bad type'); % multiplication is not defined for moment form
end


