function ip = gaussian_inner_product(g1,g2)

sum  = g1.exp + g2.exp;
prod = g1.exp * g2.exp;
rat  = prod/sum;
deltaR = norm(g1.R-g2.R);

ip = 2^(3/2)*(rat/sum)^(3/4)*exp(-rat*deltaR^2);