www.pudn.com > gmm_utilities.zip > kernel_convolve.m

function gr = kernel_convolve(g1, g2) 
 
D = size(g1.x, 1); 
M = size(g1.x, 2); 
N = size(g2.x, 2); 
R = M*N; 
 
gr.P = g1.P + g2.P; 
gr.x = zeros(D,R); 
gr.w = zeros(1,R); 
 
k=1; 
for i=1:M 
    xi = g1.x(:,i); 
    wi = g1.w(i); 
     
    for j=1:N 
        gr.x(:,k) = xi + g2.x(:,j); 
        gr.w(k) = wi * g2.w(j); 
         
        k = k+1; 
    end 
end